#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
scalar_t__ calloc(int, int);
scalar_t__ ceil(double);
scalar_t__ floor(double);
int /*<<< orphan*/ free(double *);
double log(double);
scalar_t__ pow(double, double);
void DyadUpdate(double GAMMA_p[], double Y_p[], double C_p[], const int n) {
  int L, ii, ell, k, pos, scale, s_length;
  double *s_p;
  L = (int)ceil(log((double)n) / log((double)2));
  s_length = (int)pow(((double)2), ((double)L + 1));
  s_p = (double *)calloc(s_length, sizeof(double));
  for (ii = 2; ii <= n; ii++) {
    for (ell = 0; ell <= L - 1; ell++) {
      k = (int)ceil(Y_p[ii - 1 - 1] / ((int)pow(((double)2), ((double)ell))));
      pos = k;
      if (ell > 0) {
        for (scale = ell - 1; scale >= 0; scale--) {
          pos = pos + (int)pow(((double)2), ((double)L - scale));
        }
      }
      s_p[pos - 1] = s_p[pos - 1] + C_p[ii - 1 - 1];
    }
    for (ell = 0; ell <= L - 1; ell++) {
      k = (int)floor((double)(Y_p[ii - 1] - 1) /
                     ((int)pow(((double)2), ((double)ell))));
      if ((double)k / 2 > (int)floor(((double)k) / 2)) {
        pos = k;
        if (ell > 0) {
          for (scale = ell - 1; scale >= 0; scale--) {
            pos = pos + (int)pow(((double)2), ((double)L - scale));
          }
        }
        GAMMA_p[ii - 1] = GAMMA_p[ii - 1] + s_p[pos - 1];
      }
    }
  }
  free(s_p);
}


x86 assembly (clang, -O0):
<DyadUpdate>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
cvtsi2sdl -0x1c(%rbp),%xmm0
call   1060 <log@plt>
movsd  %xmm0,-0x48(%rbp)
movsd  0xe72(%rip),%xmm0
call   1060 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
divsd  %xmm1,%xmm0
call   1040 <ceil@plt>
mov    %eax,-0x20(%rbp)
cvtsi2sdl -0x20(%rbp),%xmm1
movsd  0xe54(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  0xe40(%rip),%xmm0
call   1070 <pow@plt>
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x40(%rbp)
movl   $0x2,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     144d <DyadUpdate+0x2ed>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     12ee <DyadUpdate+0x18e>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
cvtsi2sdl -0x28(%rbp),%xmm1
movsd  0xdd5(%rip),%xmm0
call   1070 <pow@plt>
movsd  -0x50(%rbp),%xmm0
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
call   1040 <ceil@plt>
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x28(%rbp)
jle    12a7 <DyadUpdate+0x147>
mov    -0x28(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jl     12a2 <DyadUpdate+0x142>
mov    -0x30(%rbp),%eax
mov    %eax,-0x54(%rbp)
cvtsi2sdl -0x20(%rbp),%xmm1
cvtsi2sdl -0x34(%rbp),%xmm0
subsd  %xmm0,%xmm1
movsd  0xd7c(%rip),%xmm0
call   1070 <pow@plt>
mov    %rax,%rcx
mov    -0x54(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x30(%rbp)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    125e <DyadUpdate+0xfe>
jmp    12a7 <DyadUpdate+0x147>
mov    -0x40(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11f3 <DyadUpdate+0x93>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     143a <DyadUpdate+0x2da>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xce8(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
cvtsi2sdl -0x28(%rbp),%xmm1
movsd  0xcca(%rip),%xmm0
call   1070 <pow@plt>
movsd  -0x68(%rbp),%xmm0
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
call   1080 <floor@plt>
mov    %eax,-0x2c(%rbp)
cvtsi2sdl -0x2c(%rbp),%xmm0
movsd  0xca3(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
cvtsi2sdl -0x2c(%rbp),%xmm0
movsd  0xc8d(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1080 <floor@plt>
movsd  -0x60(%rbp),%xmm0
cvtsi2sd %eax,%xmm1
ucomisd %xmm1,%xmm0
jbe    1427 <DyadUpdate+0x2c7>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x28(%rbp)
jle    13f1 <DyadUpdate+0x291>
mov    -0x28(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jl     13ec <DyadUpdate+0x28c>
mov    -0x30(%rbp),%eax
mov    %eax,-0x6c(%rbp)
cvtsi2sdl -0x20(%rbp),%xmm1
cvtsi2sdl -0x34(%rbp),%xmm0
subsd  %xmm0,%xmm1
movsd  0xc32(%rip),%xmm0
call   1070 <pow@plt>
mov    %rax,%rcx
mov    -0x6c(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x30(%rbp)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    13a8 <DyadUpdate+0x248>
jmp    13f1 <DyadUpdate+0x291>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    142c <DyadUpdate+0x2cc>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    12f5 <DyadUpdate+0x195>
jmp    143f <DyadUpdate+0x2df>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11e0 <DyadUpdate+0x80>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<DyadUpdate>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %ecx,%ebx
mov    %rdx,0x38(%rsp)
mov    %rsi,0x48(%rsp)
mov    %rdi,0x50(%rsp)
cvtsi2sd %ecx,%xmm0
call   1060 <log@plt>
divsd  0xe70(%rip),%xmm0
call   1040 <ceil@plt>
mov    %rax,0x18(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0xe5f(%rip),%xmm1
movsd  %xmm0,0x8(%rsp)
addsd  %xmm0,%xmm1
movsd  0xe55(%rip),%xmm0
call   1070 <pow@plt>
mov    %eax,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,0x30(%rsp)
mov    %rbx,0x28(%rsp)
cmp    $0x2,%ebx
jge    11f3 <DyadUpdate+0x93>
mov    0x30(%rsp),%rdi
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>
mov    0x28(%rsp),%rax
add    $0x1,%eax
mov    %rax,0x28(%rsp)
mov    $0x2,%ebx
movsd  0xe03(%rip),%xmm0
jmp    121b <DyadUpdate+0xbb>
nop
add    $0x1,%rbx
cmp    0x28(%rsp),%rbx
je     11db <DyadUpdate+0x7b>
cmpl   $0x0,0x18(%rsp)
jle    1210 <DyadUpdate+0xb0>
xorpd  %xmm1,%xmm1
movl   $0x0,0x4(%rsp)
mov    $0xffffffff,%r15d
xor    %r14d,%r14d
mov    0x38(%rsp),%r12
mov    %rbx,0x40(%rsp)
jmp    129e <DyadUpdate+0x13e>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  0xdb8(%rip),%xmm0
movslq %ebx,%rax
mov    0x30(%rsp),%rcx
movsd  -0x8(%rcx,%rax,8),%xmm1
mov    0x40(%rsp),%rbx
addsd  -0x10(%r12,%rbx,8),%xmm1
movsd  %xmm1,-0x8(%rcx,%rax,8)
movsd  0x20(%rsp),%xmm1
addsd  0xd82(%rip),%xmm1
add    $0x1,%r14d
addl   $0xffffffff,0x4(%rsp)
add    $0x1,%r15d
cmp    0x18(%rsp),%r14d
je     1400 <DyadUpdate+0x2a0>
mov    0x48(%rsp),%rax
movsd  -0x10(%rax,%rbx,8),%xmm2
movsd  %xmm2,0x10(%rsp)
movsd  %xmm1,0x20(%rsp)
call   1070 <pow@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0x10(%rsp),%xmm0
divsd  %xmm1,%xmm0
call   1040 <ceil@plt>
mov    %rax,%rbx
mov    %r14d,%r13d
sub    $0x1,%r13d
jb     1250 <DyadUpdate+0xf0>
test   $0x3,%r14b
je     1337 <DyadUpdate+0x1d7>
mov    %r14d,%r12d
and    $0x3,%r12d
neg    %r12d
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    (%r15,%rbp,1),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  0xcf3(%rip),%xmm0
call   1070 <pow@plt>
add    %eax,%ebx
add    $0xffffffff,%ebp
cmp    %ebp,%r12d
jne    1300 <DyadUpdate+0x1a0>
sub    0x4(%rsp),%ebp
mov    0x38(%rsp),%r12
jmp    133a <DyadUpdate+0x1da>
mov    %r14d,%ebp
cmp    $0x3,%r13d
movsd  0xcca(%rip),%xmm0
jb     1258 <DyadUpdate+0xf8>
add    $0x4,%ebp
nop
lea    -0x5(%rbp),%eax
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movsd  0x8(%rsp),%xmm1
subsd  %xmm2,%xmm1
call   1070 <pow@plt>
mov    %rax,%r13
add    %ebx,%r13d
lea    -0x6(%rbp),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  0xc85(%rip),%xmm0
call   1070 <pow@plt>
mov    %rax,%rbx
add    %r13d,%ebx
lea    -0x7(%rbp),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  0xc5e(%rip),%xmm0
call   1070 <pow@plt>
mov    %rax,%r13
add    %ebx,%r13d
lea    -0x8(%rbp),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  0xc37(%rip),%xmm0
call   1070 <pow@plt>
movsd  0xc2a(%rip),%xmm0
mov    %rax,%rbx
add    %r13d,%ebx
add    $0xfffffffc,%ebp
cmp    $0x4,%ebp
jg     1350 <DyadUpdate+0x1f0>
jmp    1258 <DyadUpdate+0xf8>
nopl   (%rax)
cmpl   $0x0,0x18(%rsp)
jle    1210 <DyadUpdate+0xb0>
xorpd  %xmm1,%xmm1
xor    %ebp,%ebp
mov    $0xffffffff,%r15d
xor    %r12d,%r12d
jmp    1474 <DyadUpdate+0x314>
nopl   0x0(%rax)
mov    0x50(%rsp),%rcx
mov    0x40(%rsp),%rbx
movsd  -0x8(%rcx,%rbx,8),%xmm0
movslq %r14d,%rax
mov    0x30(%rsp),%rdx
addsd  -0x8(%rdx,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rcx,%rbx,8)
movsd  0x20(%rsp),%xmm1
mov    0x4(%rsp),%ebp
addsd  0xbb2(%rip),%xmm1
add    $0x1,%r12d
add    $0xffffffff,%ebp
add    $0x1,%r15d
cmp    0x18(%rsp),%r12d
movsd  0xba2(%rip),%xmm0
je     1210 <DyadUpdate+0xb0>
mov    0x48(%rsp),%rax
movsd  -0x8(%rax,%rbx,8),%xmm2
addsd  0xb91(%rip),%xmm2
movsd  %xmm2,0x10(%rsp)
movsd  %xmm1,0x20(%rsp)
call   1070 <pow@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0x10(%rsp),%xmm0
divsd  %xmm1,%xmm0
call   1080 <floor@plt>
mov    %rax,%r14
xorps  %xmm0,%xmm0
cvtsi2sd %r14d,%xmm0
mulsd  0xb5f(%rip),%xmm0
movsd  %xmm0,0x10(%rsp)
call   1080 <floor@plt>
movsd  0x20(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x10(%rsp),%xmm2
ucomisd %xmm0,%xmm2
jbe    144e <DyadUpdate+0x2ee>
test   %r12d,%r12d
mov    %ebp,0x4(%rsp)
je     1420 <DyadUpdate+0x2c0>
lea    -0x1(%r12),%ebx
test   $0x3,%r12b
je     1545 <DyadUpdate+0x3e5>
mov    %r12d,%ebp
and    $0x3,%ebp
neg    %ebp
xor    %r13d,%r13d
nopl   0x0(%rax)
lea    (%r15,%r13,1),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  0xae3(%rip),%xmm0
call   1070 <pow@plt>
add    %eax,%r14d
add    $0xffffffff,%r13d
cmp    %r13d,%ebp
jne    1510 <DyadUpdate+0x3b0>
sub    0x4(%rsp),%r13d
jmp    1548 <DyadUpdate+0x3e8>
mov    %r12d,%r13d
cmp    $0x3,%ebx
movsd  0xabd(%rip),%xmm0
jb     1420 <DyadUpdate+0x2c0>
add    $0x4,%r13d
nopl   (%rax)
lea    -0x5(%r13),%eax
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movsd  0x8(%rsp),%xmm1
subsd  %xmm2,%xmm1
call   1070 <pow@plt>
mov    %rax,%rbx
add    %r14d,%ebx
lea    -0x6(%r13),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  0xa73(%rip),%xmm0
call   1070 <pow@plt>
mov    %rax,%rbp
add    %ebx,%ebp
lea    -0x7(%r13),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  0xa4c(%rip),%xmm0
call   1070 <pow@plt>
mov    %rax,%rbx
add    %ebp,%ebx
lea    -0x8(%r13),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  0xa25(%rip),%xmm0
call   1070 <pow@plt>
movsd  0xa18(%rip),%xmm0
mov    %rax,%r14
add    %ebx,%r14d
add    $0xfffffffc,%r13d
cmp    $0x4,%r13d
jg     1560 <DyadUpdate+0x400>
jmp    1420 <DyadUpdate+0x2c0>

LLVM IR (-O0):
define void @DyadUpdate(double* %0, double* %1, double* %2, i32 %3) {
%5 = alloca double*
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca double*
store double* %0, double** %5
store double* %1, double** %6
store double* %2, double** %7
store i32 %3, i32* %8
%17 = load i32, i32* %8
%18 = sitofp i32 %17 to double
%19 = call double @log(double %18)
%20 = call double @log(double 2.000000e+00)
%21 = fdiv double %19, %20
%22 = call i64 @ceil(double %21)
%23 = trunc i64 %22 to i32
store i32 %23, i32* %9
%24 = load i32, i32* %9
%25 = sitofp i32 %24 to double
%26 = fadd double %25, 1.000000e+00
%27 = call i64 @pow(double 2.000000e+00, double %26)
%28 = trunc i64 %27 to i32
store i32 %28, i32* %15
%29 = load i32, i32* %15
%30 = call i64 @calloc(i32 %29, i32 8)
%31 = inttoptr i64 %30 to double*
store double* %31, double** %16
store i32 2, i32* %10
br label %32
32: ; preds = %184, %4
%33 = load i32, i32* %10
%34 = load i32, i32* %8
%35 = icmp sle i32 %33, %34
br i1 %35, label %36, label %187
36: ; preds = %32
store i32 0, i32* %11
br label %37
37: ; preds = %101, %36
%38 = load i32, i32* %11
%39 = load i32, i32* %9
%40 = sub nsw i32 %39, 1
%41 = icmp sle i32 %38, %40
br i1 %41, label %42, label %104
42: ; preds = %37
%43 = load double*, double** %6
%44 = load i32, i32* %10
%45 = sub nsw i32 %44, 1
%46 = sub nsw i32 %45, 1
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %43, i64 %47
%49 = load double, double* %48
%50 = load i32, i32* %11
%51 = sitofp i32 %50 to double
%52 = call i64 @pow(double 2.000000e+00, double %51)
%53 = trunc i64 %52 to i32
%54 = sitofp i32 %53 to double
%55 = fdiv double %49, %54
%56 = call i64 @ceil(double %55)
%57 = trunc i64 %56 to i32
store i32 %57, i32* %12
%58 = load i32, i32* %12
store i32 %58, i32* %13
%59 = load i32, i32* %11
%60 = icmp sgt i32 %59, 0
br i1 %60, label %61, label %81
61: ; preds = %42
%62 = load i32, i32* %11
%63 = sub nsw i32 %62, 1
store i32 %63, i32* %14
br label %64
64: ; preds = %77, %61
%65 = load i32, i32* %14
%66 = icmp sge i32 %65, 0
br i1 %66, label %67, label %80
67: ; preds = %64
%68 = load i32, i32* %13
%69 = load i32, i32* %9
%70 = sitofp i32 %69 to double
%71 = load i32, i32* %14
%72 = sitofp i32 %71 to double
%73 = fsub double %70, %72
%74 = call i64 @pow(double 2.000000e+00, double %73)
%75 = trunc i64 %74 to i32
%76 = add nsw i32 %68, %75
store i32 %76, i32* %13
br label %77
77: ; preds = %67
%78 = load i32, i32* %14
%79 = add nsw i32 %78, -1
store i32 %79, i32* %14
br label %64
80: ; preds = %64
br label %81
81: ; preds = %80, %42
%82 = load double*, double** %16
%83 = load i32, i32* %13
%84 = sub nsw i32 %83, 1
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %82, i64 %85
%87 = load double, double* %86
%88 = load double*, double** %7
%89 = load i32, i32* %10
%90 = sub nsw i32 %89, 1
%91 = sub nsw i32 %90, 1
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds double, double* %88, i64 %92
%94 = load double, double* %93
%95 = fadd double %87, %94
%96 = load double*, double** %16
%97 = load i32, i32* %13
%98 = sub nsw i32 %97, 1
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double, double* %96, i64 %99
store double %95, double* %100
br label %101
101: ; preds = %81
%102 = load i32, i32* %11
%103 = add nsw i32 %102, 1
store i32 %103, i32* %11
br label %37
104: ; preds = %37
store i32 0, i32* %11
br label %105
105: ; preds = %180, %104
%106 = load i32, i32* %11
%107 = load i32, i32* %9
%108 = sub nsw i32 %107, 1
%109 = icmp sle i32 %106, %108
br i1 %109, label %110, label %183
110: ; preds = %105
%111 = load double*, double** %6
%112 = load i32, i32* %10
%113 = sub nsw i32 %112, 1
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds double, double* %111, i64 %114
%116 = load double, double* %115
%117 = fsub double %116, 1.000000e+00
%118 = load i32, i32* %11
%119 = sitofp i32 %118 to double
%120 = call i64 @pow(double 2.000000e+00, double %119)
%121 = trunc i64 %120 to i32
%122 = sitofp i32 %121 to double
%123 = fdiv double %117, %122
%124 = call i64 @floor(double %123)
%125 = trunc i64 %124 to i32
store i32 %125, i32* %12
%126 = load i32, i32* %12
%127 = sitofp i32 %126 to double
%128 = fdiv double %127, 2.000000e+00
%129 = load i32, i32* %12
%130 = sitofp i32 %129 to double
%131 = fdiv double %130, 2.000000e+00
%132 = call i64 @floor(double %131)
%133 = trunc i64 %132 to i32
%134 = sitofp i32 %133 to double
%135 = fcmp ogt double %128, %134
br i1 %135, label %136, label %179
136: ; preds = %110
%137 = load i32, i32* %12
store i32 %137, i32* %13
%138 = load i32, i32* %11
%139 = icmp sgt i32 %138, 0
br i1 %139, label %140, label %160
140: ; preds = %136
%141 = load i32, i32* %11
%142 = sub nsw i32 %141, 1
store i32 %142, i32* %14
br label %143
143: ; preds = %156, %140
%144 = load i32, i32* %14
%145 = icmp sge i32 %144, 0
br i1 %145, label %146, label %159
146: ; preds = %143
%147 = load i32, i32* %13
%148 = load i32, i32* %9
%149 = sitofp i32 %148 to double
%150 = load i32, i32* %14
%151 = sitofp i32 %150 to double
%152 = fsub double %149, %151
%153 = call i64 @pow(double 2.000000e+00, double %152)
%154 = trunc i64 %153 to i32
%155 = add nsw i32 %147, %154
store i32 %155, i32* %13
br label %156
156: ; preds = %146
%157 = load i32, i32* %14
%158 = add nsw i32 %157, -1
store i32 %158, i32* %14
br label %143
159: ; preds = %143
br label %160
160: ; preds = %159, %136
%161 = load double*, double** %5
%162 = load i32, i32* %10
%163 = sub nsw i32 %162, 1
%164 = sext i32 %163 to i64
%165 = getelementptr inbounds double, double* %161, i64 %164
%166 = load double, double* %165
%167 = load double*, double** %16
%168 = load i32, i32* %13
%169 = sub nsw i32 %168, 1
%170 = sext i32 %169 to i64
%171 = getelementptr inbounds double, double* %167, i64 %170
%172 = load double, double* %171
%173 = fadd double %166, %172
%174 = load double*, double** %5
%175 = load i32, i32* %10
%176 = sub nsw i32 %175, 1
%177 = sext i32 %176 to i64
%178 = getelementptr inbounds double, double* %174, i64 %177
store double %173, double* %178
br label %179
179: ; preds = %160, %110
br label %180
180: ; preds = %179
%181 = load i32, i32* %11
%182 = add nsw i32 %181, 1
store i32 %182, i32* %11
br label %105
183: ; preds = %105
br label %184
184: ; preds = %183
%185 = load i32, i32* %10
%186 = add nsw i32 %185, 1
store i32 %186, i32* %10
br label %32
187: ; preds = %32
%188 = load double*, double** %16
%189 = call i32 @free(double* %188)
ret void
}

LLVM IR (-O2):
define void @DyadUpdate(double* nocapture %0, double* nocapture %1, double* nocapture %2, i32 %3) {
%5 = sitofp i32 %3 to double
%6 = tail call double @log(double %5)
%7 = fdiv double %6, 0x3FE62E42FEFA39EF
%8 = tail call i64 @ceil(double %7)
%9 = trunc i64 %8 to i32
%10 = sitofp i32 %9 to double
%11 = fadd double %10, 1.000000e+00
%12 = tail call i64 @pow(double 2.000000e+00, double %11)
%13 = trunc i64 %12 to i32
%14 = tail call i64 @calloc(i32 %13, i32 8)
%15 = inttoptr i64 %14 to double*
%16 = icmp slt i32 %3, 2
br i1 %16, label %184, label %17
17: ; preds = %4
%18 = icmp sgt i32 %9, 0
%19 = add nuw i32 %3, 1
%20 = zext i32 %19 to i64
br label %21
21: ; preds = %17, %181
%22 = phi i64 [ 2, %17 ], [ %182, %181 ]
br i1 %18, label %23, label %181
23: ; preds = %21
%24 = add nsw i64 %22, -2
%25 = getelementptr inbounds double, double* %1, i64 %24
%26 = getelementptr inbounds double, double* %2, i64 %24
br label %32
27: ; preds = %92
br i1 %18, label %28, label %181
28: ; preds = %27
%29 = add nsw i64 %22, -1
%30 = getelementptr inbounds double, double* %1, i64 %29
%31 = getelementptr inbounds double, double* %0, i64 %29
br label %102
32: ; preds = %23, %92
%33 = phi i32 [ 0, %23 ], [ %100, %92 ]
%34 = add nsw i32 %33, -1
%35 = load double, double* %25
%36 = sitofp i32 %33 to double
%37 = tail call i64 @pow(double 2.000000e+00, double %36)
%38 = trunc i64 %37 to i32
%39 = sitofp i32 %38 to double
%40 = fdiv double %35, %39
%41 = tail call i64 @ceil(double %40)
%42 = trunc i64 %41 to i32
%43 = icmp eq i32 %33, 0
br i1 %43, label %92, label %44
44: ; preds = %32
%45 = and i32 %33, 3
%46 = icmp eq i32 %45, 0
br i1 %46, label %59, label %47
47: ; preds = %44, %47
%48 = phi i32 [ %51, %47 ], [ %33, %44 ]
%49 = phi i32 [ %56, %47 ], [ %42, %44 ]
%50 = phi i32 [ %57, %47 ], [ 0, %44 ]
%51 = add nsw i32 %48, -1
%52 = sitofp i32 %51 to double
%53 = fsub double %10, %52
%54 = tail call i64 @pow(double 2.000000e+00, double %53)
%55 = trunc i64 %54 to i32
%56 = add nsw i32 %49, %55
%57 = add i32 %50, 1
%58 = icmp eq i32 %57, %45
br i1 %58, label %59, label %47
59: ; preds = %47, %44
%60 = phi i32 [ %33, %44 ], [ %51, %47 ]
%61 = phi i32 [ %42, %44 ], [ %56, %47 ]
%62 = phi i32 [ undef, %44 ], [ %56, %47 ]
%63 = icmp ult i32 %34, 3
br i1 %63, label %92, label %64
64: ; preds = %59, %64
%65 = phi i32 [ %85, %64 ], [ %60, %59 ]
%66 = phi i32 [ %90, %64 ], [ %61, %59 ]
%67 = add nsw i32 %65, -1
%68 = sitofp i32 %67 to double
%69 = fsub double %10, %68
%70 = tail call i64 @pow(double 2.000000e+00, double %69)
%71 = trunc i64 %70 to i32
%72 = add nsw i32 %66, %71
%73 = add nsw i32 %65, -2
%74 = sitofp i32 %73 to double
%75 = fsub double %10, %74
%76 = tail call i64 @pow(double 2.000000e+00, double %75)
%77 = trunc i64 %76 to i32
%78 = add nsw i32 %72, %77
%79 = add nsw i32 %65, -3
%80 = sitofp i32 %79 to double
%81 = fsub double %10, %80
%82 = tail call i64 @pow(double 2.000000e+00, double %81)
%83 = trunc i64 %82 to i32
%84 = add nsw i32 %78, %83
%85 = add nsw i32 %65, -4
%86 = sitofp i32 %85 to double
%87 = fsub double %10, %86
%88 = tail call i64 @pow(double 2.000000e+00, double %87)
%89 = trunc i64 %88 to i32
%90 = add nsw i32 %84, %89
%91 = icmp sgt i32 %65, 4
br i1 %91, label %64, label %92
92: ; preds = %59, %64, %32
%93 = phi i32 [ %42, %32 ], [ %62, %59 ], [ %90, %64 ]
%94 = add nsw i32 %93, -1
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds double, double* %15, i64 %95
%97 = load double, double* %96
%98 = load double, double* %26
%99 = fadd double %97, %98
store double %99, double* %96
%100 = add nuw nsw i32 %33, 1
%101 = icmp eq i32 %100, %9
br i1 %101, label %27, label %32
102: ; preds = %28, %178
%103 = phi i32 [ 0, %28 ], [ %179, %178 ]
%104 = add nsw i32 %103, -1
%105 = load double, double* %30
%106 = fadd double %105, -1.000000e+00
%107 = sitofp i32 %103 to double
%108 = tail call i64 @pow(double 2.000000e+00, double %107)
%109 = trunc i64 %108 to i32
%110 = sitofp i32 %109 to double
%111 = fdiv double %106, %110
%112 = tail call i64 @floor(double %111)
%113 = trunc i64 %112 to i32
%114 = sitofp i32 %113 to double
%115 = fmul double %114, 5.000000e-01
%116 = tail call i64 @floor(double %115)
%117 = trunc i64 %116 to i32
%118 = sitofp i32 %117 to double
%119 = fcmp ogt double %115, %118
br i1 %119, label %120, label %178
120: ; preds = %102
%121 = icmp eq i32 %103, 0
br i1 %121, label %170, label %122
122: ; preds = %120
%123 = and i32 %103, 3
%124 = icmp eq i32 %123, 0
br i1 %124, label %137, label %125
125: ; preds = %122, %125
%126 = phi i32 [ %129, %125 ], [ %103, %122 ]
%127 = phi i32 [ %134, %125 ], [ %113, %122 ]
%128 = phi i32 [ %135, %125 ], [ 0, %122 ]
%129 = add nsw i32 %126, -1
%130 = sitofp i32 %129 to double
%131 = fsub double %10, %130
%132 = tail call i64 @pow(double 2.000000e+00, double %131)
%133 = trunc i64 %132 to i32
%134 = add nsw i32 %127, %133
%135 = add i32 %128, 1
%136 = icmp eq i32 %135, %123
br i1 %136, label %137, label %125
137: ; preds = %125, %122
%138 = phi i32 [ %103, %122 ], [ %129, %125 ]
%139 = phi i32 [ %113, %122 ], [ %134, %125 ]
%140 = phi i32 [ undef, %122 ], [ %134, %125 ]
%141 = icmp ult i32 %104, 3
br i1 %141, label %170, label %142
142: ; preds = %137, %142
%143 = phi i32 [ %163, %142 ], [ %138, %137 ]
%144 = phi i32 [ %168, %142 ], [ %139, %137 ]
%145 = add nsw i32 %143, -1
%146 = sitofp i32 %145 to double
%147 = fsub double %10, %146
%148 = tail call i64 @pow(double 2.000000e+00, double %147)
%149 = trunc i64 %148 to i32
%150 = add nsw i32 %144, %149
%151 = add nsw i32 %143, -2
%152 = sitofp i32 %151 to double
%153 = fsub double %10, %152
%154 = tail call i64 @pow(double 2.000000e+00, double %153)
%155 = trunc i64 %154 to i32
%156 = add nsw i32 %150, %155
%157 = add nsw i32 %143, -3
%158 = sitofp i32 %157 to double
%159 = fsub double %10, %158
%160 = tail call i64 @pow(double 2.000000e+00, double %159)
%161 = trunc i64 %160 to i32
%162 = add nsw i32 %156, %161
%163 = add nsw i32 %143, -4
%164 = sitofp i32 %163 to double
%165 = fsub double %10, %164
%166 = tail call i64 @pow(double 2.000000e+00, double %165)
%167 = trunc i64 %166 to i32
%168 = add nsw i32 %162, %167
%169 = icmp sgt i32 %143, 4
br i1 %169, label %142, label %170
170: ; preds = %137, %142, %120
%171 = phi i32 [ %113, %120 ], [ %140, %137 ], [ %168, %142 ]
%172 = load double, double* %31
%173 = add nsw i32 %171, -1
%174 = sext i32 %173 to i64
%175 = getelementptr inbounds double, double* %15, i64 %174
%176 = load double, double* %175
%177 = fadd double %172, %176
store double %177, double* %31
br label %178
178: ; preds = %102, %170
%179 = add nuw nsw i32 %103, 1
%180 = icmp eq i32 %179, %9
br i1 %180, label %181, label %102
181: ; preds = %178, %21, %27
%182 = add nuw nsw i64 %22, 1
%183 = icmp eq i64 %182, %20
br i1 %183, label %184, label %21
184: ; preds = %181, %4
%185 = tail call i32 @free(double* %15)
ret void
}