#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 */
double dmean(double *, size_t);
double *dvector(size_t);
double sqrt(double);
double ess(double *values, size_t samples, size_t stepSize) {
  size_t maxLag = samples - 1;
  double *gammaStat = dvector(maxLag);
  double varStat = 0.0;
  double mean = dmean(values, samples);
  for (size_t lag = 0; lag < maxLag; lag++) {
    for (size_t j = 0; j < samples - lag; j++) {
      double del1 = values[j] - mean;
      double del2 = values[j + lag] - mean;
      gammaStat[lag] += (del1 * del2);
    }
    gammaStat[lag] /= ((double)(samples - lag));
    if (lag == 0) {
      varStat = gammaStat[0];
    } else if (lag % 2 == 0) {
      if (gammaStat[lag - 1] + gammaStat[lag] > 0) {
        varStat += 2.0 * (gammaStat[lag - 1] + gammaStat[lag]);
      } else
        maxLag = lag;
    }
  }
  double stdErrorOfMean = sqrt(varStat / samples);
  double ACT;
  double ESS;
  if (gammaStat[0] == 0)
    ACT = 0;
  else
    ACT = stepSize * varStat / gammaStat[0];
  if (ACT == 0)
    ESS = 1;
  else
    ESS = (stepSize * samples) / ACT;
  double stdErrOfACT =
      (2.0 * sqrt(2.0 * (2.0 * (double)(maxLag + 1)) / samples) *
       (varStat / gammaStat[0]) * stepSize);
  return ESS;
}


x86 assembly (clang, -O0):
<ess>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x20(%rbp)
mov    -0x20(%rbp),%rdi
call   1030 <dvector@plt>
mov    %rax,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
call   1050 <dmean@plt>
movsd  %xmm0,-0x38(%rbp)
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x20(%rbp),%rax
jae    1313 <ess+0x1e3>
movq   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%rax
mov    -0x10(%rbp),%rcx
sub    -0x40(%rbp),%rcx
cmp    %rcx,%rax
jae    1214 <ess+0xe4>
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
add    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    1198 <ess+0x68>
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
sub    %rcx,%rax
movq   %rax,%xmm1
movaps 0xdd5(%rip),%xmm0
punpckldq %xmm0,%xmm1
movapd 0xdd9(%rip),%xmm0
subpd  %xmm0,%xmm1
movaps %xmm1,%xmm0
unpckhpd %xmm1,%xmm1
addsd  %xmm0,%xmm1
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
cmpq   $0x0,-0x40(%rbp)
jne    1279 <ess+0x149>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    12fd <ess+0x1cd>
mov    -0x40(%rbp),%rax
and    $0x1,%rax
cmp    $0x0,%rax
jne    12f8 <ess+0x1c8>
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    12eb <ess+0x1bb>
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  -0x30(%rbp),%xmm1
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    12f3 <ess+0x1c3>
mov    -0x40(%rbp),%rax
mov    %rax,-0x20(%rbp)
jmp    12f8 <ess+0x1c8>
jmp    12fd <ess+0x1cd>
jmp    1302 <ess+0x1d2>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1182 <ess+0x52>
movsd  -0x30(%rbp),%xmm0
movaps 0xce1(%rip),%xmm2
movq   -0x10(%rbp),%xmm1
punpckldq %xmm2,%xmm1
movapd 0xce0(%rip),%xmm2
subpd  %xmm2,%xmm1
movaps %xmm1,%xmm2
unpckhpd %xmm1,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movsd  %xmm0,-0x60(%rbp)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1375 <ess+0x245>
jp     1375 <ess+0x245>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    13ae <ess+0x27e>
movaps 0xc84(%rip),%xmm1
movq   -0x18(%rbp),%xmm0
punpckldq %xmm1,%xmm0
movapd 0xc83(%rip),%xmm1
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
mulsd  -0x30(%rbp),%xmm0
mov    -0x28(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    13d8 <ess+0x2a8>
jp     13d8 <ess+0x2a8>
movsd  0xc52(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
jmp    1415 <ess+0x2e5>
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
imul   %rcx,%rax
movq   %rax,%xmm0
movaps 0xc10(%rip),%xmm1
punpckldq %xmm1,%xmm0
movapd 0xc14(%rip),%xmm1
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
movq   %rax,%xmm0
movaps 0xbd7(%rip),%xmm3
movaps %xmm3,-0xa0(%rbp)
punpckldq %xmm3,%xmm0
movapd 0xbd4(%rip),%xmm2
movaps %xmm2,-0x90(%rbp)
subpd  %xmm2,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
addsd  %xmm0,%xmm0
addsd  %xmm0,%xmm0
movq   -0x10(%rbp),%xmm1
punpckldq %xmm3,%xmm1
subpd  %xmm2,%xmm1
movaps %xmm1,%xmm2
unpckhpd %xmm1,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movaps -0xa0(%rbp),%xmm3
movaps -0x90(%rbp),%xmm2
addsd  %xmm0,%xmm0
movsd  -0x30(%rbp),%xmm1
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm4
divsd  %xmm4,%xmm1
mulsd  %xmm1,%xmm0
movq   -0x18(%rbp),%xmm1
punpckldq %xmm3,%xmm1
subpd  %xmm2,%xmm1
movaps %xmm1,%xmm2
unpckhpd %xmm1,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x70(%rbp),%xmm0
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ess>:
push   %r15
push   %r14
push   %r12
push   %rbx
push   %rax
mov    %rdx,%r14
mov    %rsi,%r15
mov    %rdi,%r12
lea    -0x1(%rsi),%rdi
call   1030 <dvector@plt>
mov    %rax,%rbx
mov    %r12,%rdi
mov    %r15,%rsi
call   1050 <dmean@plt>
mov    %r15,%r8
xorpd  %xmm8,%xmm8
xorpd  %xmm3,%xmm3
add    $0xffffffffffffffff,%r8
je     12b1 <ess+0x181>
lea    0x8(%r12),%rcx
xorpd  %xmm2,%xmm2
xor    %r9d,%r9d
movdqa 0xe7e(%rip),%xmm1
movapd 0xe86(%rip),%xmm4
mov    %r15,%r10
xorpd  %xmm3,%xmm3
xor    %edi,%edi
jmp    11b9 <ess+0x89>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rbx),%xmm3
add    $0x1,%rdi
add    $0x8,%r9
add    $0xffffffffffffffff,%r10
cmp    %r8,%rdi
jae    12b1 <ess+0x181>
mov    %r15,%r11
movsd  (%rbx,%rdi,8),%xmm5
sub    %rdi,%r11
je     1263 <ess+0x133>
cmp    %rdi,%r8
jne    11e0 <ess+0xb0>
xor    %edx,%edx
jmp    1239 <ess+0x109>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r10,%rsi
and    $0xfffffffffffffffe,%rsi
lea    (%rcx,%r9,1),%rax
xor    %edx,%edx
nopl   (%rax)
movsd  -0x8(%rcx,%rdx,8),%xmm6
subsd  %xmm0,%xmm6
movsd  -0x8(%rax,%rdx,8),%xmm7
subsd  %xmm0,%xmm7
mulsd  %xmm6,%xmm7
addsd  %xmm5,%xmm7
movsd  %xmm7,(%rbx,%rdi,8)
movsd  (%rcx,%rdx,8),%xmm6
subsd  %xmm0,%xmm6
movsd  (%rax,%rdx,8),%xmm5
subsd  %xmm0,%xmm5
mulsd  %xmm6,%xmm5
addsd  %xmm7,%xmm5
movsd  %xmm5,(%rbx,%rdi,8)
add    $0x2,%rdx
cmp    %rdx,%rsi
jne    11f0 <ess+0xc0>
test   $0x1,%r11b
je     1263 <ess+0x133>
movsd  (%r12,%rdx,8),%xmm6
subsd  %xmm0,%xmm6
add    %rdi,%rdx
movsd  (%r12,%rdx,8),%xmm7
subsd  %xmm0,%xmm7
mulsd  %xmm6,%xmm7
addsd  %xmm7,%xmm5
movsd  %xmm5,(%rbx,%rdi,8)
movq   %r11,%xmm6
punpckldq %xmm1,%xmm6
subpd  %xmm4,%xmm6
movapd %xmm6,%xmm7
unpckhpd %xmm6,%xmm7
addsd  %xmm6,%xmm7
divsd  %xmm7,%xmm5
movsd  %xmm5,(%rbx,%rdi,8)
test   %rdi,%rdi
je     11a0 <ess+0x70>
test   $0x1,%dil
jne    11a4 <ess+0x74>
addsd  -0x8(%rbx,%rdi,8),%xmm5
ucomisd %xmm2,%xmm5
jbe    12b1 <ess+0x181>
addsd  %xmm5,%xmm5
addsd  %xmm5,%xmm3
jmp    11a4 <ess+0x74>
movq   %r15,%xmm0
punpckldq 0xd42(%rip),%xmm0
subpd  0xd4a(%rip),%xmm0
movapd %xmm0,%xmm2
unpckhpd %xmm0,%xmm2
addsd  %xmm0,%xmm2
movapd %xmm3,%xmm0
divsd  %xmm2,%xmm0
ucomisd %xmm0,%xmm8
ja     135f <ess+0x22f>
movsd  (%rbx),%xmm1
movsd  0xd33(%rip),%xmm0
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm1
jne    12f9 <ess+0x1c9>
jnp    1353 <ess+0x223>
movq   %r14,%xmm5
punpckldq 0xcfa(%rip),%xmm5
subpd  0xd02(%rip),%xmm5
movapd %xmm5,%xmm4
unpckhpd %xmm5,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm4,%xmm3
divsd  %xmm1,%xmm3
ucomisd %xmm2,%xmm3
jne    132a <ess+0x1fa>
jnp    1353 <ess+0x223>
imul   %r15,%r14
movq   %r14,%xmm1
punpckldq 0xcc5(%rip),%xmm1
subpd  0xccd(%rip),%xmm1
movapd %xmm1,%xmm0
unpckhpd %xmm1,%xmm0
addsd  %xmm1,%xmm0
divsd  %xmm3,%xmm0
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret
pxor   %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jae    12e1 <ess+0x1b1>
movsd  %xmm3,(%rsp)
call   1040 <sqrt@plt>
movsd  (%rsp),%xmm3
jmp    12e1 <ess+0x1b1>

LLVM IR (-O0):
define double @ess(double* %0, i64 %1, i64 %2) {
%4 = alloca double*
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca double*
%9 = alloca double
%10 = alloca double
%11 = alloca i64
%12 = alloca i64
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
store double* %0, double** %4
store i64 %1, i64* %5
store i64 %2, i64* %6
%19 = load i64, i64* %5
%20 = sub i64 %19, 1
store i64 %20, i64* %7
%21 = load i64, i64* %7
%22 = call double* @dvector(i64 %21)
store double* %22, double** %8
store double 0.000000e+00, double* %9
%23 = load double*, double** %4
%24 = load i64, i64* %5
%25 = call double @dmean(double* %23, i64 %24)
store double %25, double* %10
store i64 0, i64* %11
br label %26
26: ; preds = %112, %3
%27 = load i64, i64* %11
%28 = load i64, i64* %7
%29 = icmp ult i64 %27, %28
br i1 %29, label %30, label %115
30: ; preds = %26
store i64 0, i64* %12
br label %31
31: ; preds = %59, %30
%32 = load i64, i64* %12
%33 = load i64, i64* %5
%34 = load i64, i64* %11
%35 = sub i64 %33, %34
%36 = icmp ult i64 %32, %35
br i1 %36, label %37, label %62
37: ; preds = %31
%38 = load double*, double** %4
%39 = load i64, i64* %12
%40 = getelementptr inbounds double, double* %38, i64 %39
%41 = load double, double* %40
%42 = load double, double* %10
%43 = fsub double %41, %42
store double %43, double* %13
%44 = load double*, double** %4
%45 = load i64, i64* %12
%46 = load i64, i64* %11
%47 = add i64 %45, %46
%48 = getelementptr inbounds double, double* %44, i64 %47
%49 = load double, double* %48
%50 = load double, double* %10
%51 = fsub double %49, %50
store double %51, double* %14
%52 = load double, double* %13
%53 = load double, double* %14
%54 = load double*, double** %8
%55 = load i64, i64* %11
%56 = getelementptr inbounds double, double* %54, i64 %55
%57 = load double, double* %56
%58 = call double @llvm.fmuladd.f64(double %52, double %53, double %57)
store double %58, double* %56
br label %59
59: ; preds = %37
%60 = load i64, i64* %12
%61 = add i64 %60, 1
store i64 %61, i64* %12
br label %31
62: ; preds = %31
%63 = load i64, i64* %5
%64 = load i64, i64* %11
%65 = sub i64 %63, %64
%66 = uitofp i64 %65 to double
%67 = load double*, double** %8
%68 = load i64, i64* %11
%69 = getelementptr inbounds double, double* %67, i64 %68
%70 = load double, double* %69
%71 = fdiv double %70, %66
store double %71, double* %69
%72 = load i64, i64* %11
%73 = icmp eq i64 %72, 0
br i1 %73, label %74, label %78
74: ; preds = %62
%75 = load double*, double** %8
%76 = getelementptr inbounds double, double* %75, i64 0
%77 = load double, double* %76
store double %77, double* %9
br label %111
78: ; preds = %62
%79 = load i64, i64* %11
%80 = urem i64 %79, 2
%81 = icmp eq i64 %80, 0
br i1 %81, label %82, label %110
82: ; preds = %78
%83 = load double*, double** %8
%84 = load i64, i64* %11
%85 = sub i64 %84, 1
%86 = getelementptr inbounds double, double* %83, i64 %85
%87 = load double, double* %86
%88 = load double*, double** %8
%89 = load i64, i64* %11
%90 = getelementptr inbounds double, double* %88, i64 %89
%91 = load double, double* %90
%92 = fadd double %87, %91
%93 = fcmp ogt double %92, 0.000000e+00
br i1 %93, label %94, label %107
94: ; preds = %82
%95 = load double*, double** %8
%96 = load i64, i64* %11
%97 = sub i64 %96, 1
%98 = getelementptr inbounds double, double* %95, i64 %97
%99 = load double, double* %98
%100 = load double*, double** %8
%101 = load i64, i64* %11
%102 = getelementptr inbounds double, double* %100, i64 %101
%103 = load double, double* %102
%104 = fadd double %99, %103
%105 = load double, double* %9
%106 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %104, double %105)
store double %106, double* %9
br label %109
107: ; preds = %82
%108 = load i64, i64* %11
store i64 %108, i64* %7
br label %109
109: ; preds = %107, %94
br label %110
110: ; preds = %109, %78
br label %111
111: ; preds = %110, %74
br label %112
112: ; preds = %111
%113 = load i64, i64* %11
%114 = add i64 %113, 1
store i64 %114, i64* %11
br label %26
115: ; preds = %26
%116 = load double, double* %9
%117 = load i64, i64* %5
%118 = uitofp i64 %117 to double
%119 = fdiv double %116, %118
%120 = call double @sqrt(double %119)
store double %120, double* %15
%121 = load double*, double** %8
%122 = getelementptr inbounds double, double* %121, i64 0
%123 = load double, double* %122
%124 = fcmp oeq double %123, 0.000000e+00
br i1 %124, label %125, label %126
125: ; preds = %115
store double 0.000000e+00, double* %16
br label %135
126: ; preds = %115
%127 = load i64, i64* %6
%128 = uitofp i64 %127 to double
%129 = load double, double* %9
%130 = fmul double %128, %129
%131 = load double*, double** %8
%132 = getelementptr inbounds double, double* %131, i64 0
%133 = load double, double* %132
%134 = fdiv double %130, %133
store double %134, double* %16
br label %135
135: ; preds = %126, %125
%136 = load double, double* %16
%137 = fcmp oeq double %136, 0.000000e+00
br i1 %137, label %138, label %139
138: ; preds = %135
store double 1.000000e+00, double* %17
br label %146
139: ; preds = %135
%140 = load i64, i64* %6
%141 = load i64, i64* %5
%142 = mul i64 %140, %141
%143 = uitofp i64 %142 to double
%144 = load double, double* %16
%145 = fdiv double %143, %144
store double %145, double* %17
br label %146
146: ; preds = %139, %138
%147 = load i64, i64* %7
%148 = add i64 %147, 1
%149 = uitofp i64 %148 to double
%150 = fmul double 2.000000e+00, %149
%151 = fmul double 2.000000e+00, %150
%152 = load i64, i64* %5
%153 = uitofp i64 %152 to double
%154 = fdiv double %151, %153
%155 = call double @sqrt(double %154)
%156 = fmul double 2.000000e+00, %155
%157 = load double, double* %9
%158 = load double*, double** %8
%159 = getelementptr inbounds double, double* %158, i64 0
%160 = load double, double* %159
%161 = fdiv double %157, %160
%162 = fmul double %156, %161
%163 = load i64, i64* %6
%164 = uitofp i64 %163 to double
%165 = fmul double %162, %164
store double %165, double* %18
%166 = load double, double* %17
ret double %166
}

LLVM IR (-O2):
define double @ess(double* %0, i64 %1, i64 %2) {
%4 = add i64 %1, -1
%5 = tail call double* @dvector(i64 %4)
%6 = tail call double @dmean(double* %0, i64 %1)
%7 = icmp eq i64 %4, 0
br i1 %7, label %21, label %8
8: ; preds = %3, %88
%9 = phi double [ %89, %88 ], [ 0.000000e+00, %3 ]
%10 = phi i64 [ %90, %88 ], [ 0, %3 ]
%11 = sub i64 %1, %10
%12 = sub i64 %1, %10
%13 = icmp eq i64 %12, 0
%14 = getelementptr inbounds double, double* %5, i64 %10
%15 = load double, double* %14
br i1 %13, label %45, label %16
16: ; preds = %8
%17 = and i64 %11, 1
%18 = icmp eq i64 %4, %10
br i1 %18, label %31, label %19
19: ; preds = %16
%20 = and i64 %11, -2
br label %51
21: ; preds = %80, %88, %3
%22 = phi double [ 0.000000e+00, %3 ], [ %9, %80 ], [ %89, %88 ]
%23 = uitofp i64 %1 to double
%24 = fdiv double %22, %23
%25 = fcmp olt double %24, 0.000000e+00
br i1 %25, label %26, label %28
26: ; preds = %21
%27 = tail call double @sqrt(double %24)
br label %28
28: ; preds = %21, %26
%29 = load double, double* %5
%30 = fcmp oeq double %29, 0.000000e+00
br i1 %30, label %101, label %92
31: ; preds = %51, %16
%32 = phi double [ undef, %16 ], [ %71, %51 ]
%33 = phi double [ %15, %16 ], [ %71, %51 ]
%34 = phi i64 [ 0, %16 ], [ %72, %51 ]
%35 = icmp eq i64 %17, 0
br i1 %35, label %45, label %36
36: ; preds = %31
%37 = getelementptr inbounds double, double* %0, i64 %34
%38 = load double, double* %37
%39 = fsub double %38, %6
%40 = add i64 %34, %10
%41 = getelementptr inbounds double, double* %0, i64 %40
%42 = load double, double* %41
%43 = fsub double %42, %6
%44 = tail call double @llvm.fmuladd.f64(double %39, double %43, double %33)
store double %44, double* %14
br label %45
45: ; preds = %36, %31, %8
%46 = phi double [ %15, %8 ], [ %32, %31 ], [ %44, %36 ]
%47 = uitofp i64 %12 to double
%48 = getelementptr inbounds double, double* %5, i64 %10
%49 = fdiv double %46, %47
store double %49, double* %48
%50 = icmp eq i64 %10, 0
br i1 %50, label %75, label %77
51: ; preds = %51, %19
%52 = phi double [ %15, %19 ], [ %71, %51 ]
%53 = phi i64 [ 0, %19 ], [ %72, %51 ]
%54 = phi i64 [ 0, %19 ], [ %73, %51 ]
%55 = getelementptr inbounds double, double* %0, i64 %53
%56 = load double, double* %55
%57 = fsub double %56, %6
%58 = add i64 %53, %10
%59 = getelementptr inbounds double, double* %0, i64 %58
%60 = load double, double* %59
%61 = fsub double %60, %6
%62 = tail call double @llvm.fmuladd.f64(double %57, double %61, double %52)
store double %62, double* %14
%63 = or i64 %53, 1
%64 = getelementptr inbounds double, double* %0, i64 %63
%65 = load double, double* %64
%66 = fsub double %65, %6
%67 = add i64 %63, %10
%68 = getelementptr inbounds double, double* %0, i64 %67
%69 = load double, double* %68
%70 = fsub double %69, %6
%71 = tail call double @llvm.fmuladd.f64(double %66, double %70, double %62)
store double %71, double* %14
%72 = add nuw i64 %53, 2
%73 = add i64 %54, 2
%74 = icmp eq i64 %73, %20
br i1 %74, label %31, label %51
75: ; preds = %45
%76 = load double, double* %5
br label %88
77: ; preds = %45
%78 = and i64 %10, 1
%79 = icmp eq i64 %78, 0
br i1 %79, label %80, label %88
80: ; preds = %77
%81 = add i64 %10, -1
%82 = getelementptr inbounds double, double* %5, i64 %81
%83 = load double, double* %82
%84 = fadd double %49, %83
%85 = fcmp ogt double %84, 0.000000e+00
br i1 %85, label %86, label %21
86: ; preds = %80
%87 = tail call double @llvm.fmuladd.f64(double %84, double 2.000000e+00, double %9)
br label %88
88: ; preds = %75, %86, %77
%89 = phi double [ %76, %75 ], [ %87, %86 ], [ %9, %77 ]
%90 = add nuw i64 %10, 1
%91 = icmp ult i64 %90, %4
br i1 %91, label %8, label %21
92: ; preds = %28
%93 = uitofp i64 %2 to double
%94 = fmul double %22, %93
%95 = fdiv double %94, %29
%96 = fcmp oeq double %95, 0.000000e+00
br i1 %96, label %101, label %97
97: ; preds = %92
%98 = mul i64 %2, %1
%99 = uitofp i64 %98 to double
%100 = fdiv double %99, %95
br label %101
101: ; preds = %97, %92, %28
%102 = phi double [ %100, %97 ], [ 1.000000e+00, %92 ], [ 1.000000e+00, %28 ]
ret double %102
}