#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 pow(double, double);
double sqrt(int);
double twtest(double);
void tw(double *values, double *pvalues, double *twstat, double *effectn,
        int M) {
  int Mp, i;
  double nhat, lambda, mu, sigma, s, s2;
  s = 0;
  s2 = 0;
  for (i = 0; i < M; i++) {
    s += values[i];
    s2 += values[i] * values[i];
  }
  for (i = 0; i < M; i++) {
    Mp = M - i;
    nhat = (Mp + 2) * s * s / (Mp * s2 - s * s);
    lambda = values[i] * Mp / s;
    mu = (sqrt(nhat - 1) + sqrt(Mp)) * (sqrt(nhat - 1) + sqrt(Mp)) / nhat;
    sigma = (sqrt(nhat - 1) + sqrt(Mp)) / nhat *
            pow(1.0 / sqrt(nhat - 1) + 1.0 / sqrt(Mp), (double)(1.0 / 3.0));
    lambda = (lambda - mu) / sigma;
    pvalues[i] = twtest(lambda);
    twstat[i] = lambda;
    effectn[i] = nhat;
    s -= values[i];
    s2 -= values[i] * values[i];
  }
}


x86 assembly (clang, -O0):
<tw>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11c3 <tw+0x93>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1166 <tw+0x36>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    143d <tw+0x30d>
mov    -0x24(%rbp),%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
add    $0x2,%eax
cvtsi2sd %eax,%xmm0
mulsd  -0x58(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
cvtsi2sdl -0x28(%rbp),%xmm1
movsd  -0x60(%rbp),%xmm3
movsd  -0x58(%rbp),%xmm2
mulsd  -0x58(%rbp),%xmm2
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x28(%rbp),%xmm1
mulsd  %xmm1,%xmm0
divsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  0xdac(%rip),%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sqrt@plt>
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x28(%rbp),%edi
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x90(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  0xd6e(%rip),%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sqrt@plt>
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x28(%rbp),%edi
call   1050 <sqrt@plt>
movsd  -0x88(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x80(%rbp),%xmm0
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  0xd22(%rip),%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sqrt@plt>
movsd  %xmm0,-0x78(%rbp)
mov    -0x28(%rbp),%edi
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
addsd  %xmm1,%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  0xce5(%rip),%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xccd(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x28(%rbp),%edi
call   1050 <sqrt@plt>
movaps %xmm0,%xmm2
movsd  -0x70(%rbp),%xmm0
movsd  0xcac(%rip),%xmm1
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  0xc94(%rip),%xmm1
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
divsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
call   1030 <twtest@plt>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x58(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x60(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11ca <tw+0x9a>
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<tw>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
test   %r8d,%r8d
jle    1393 <tw+0x263>
mov    %rdx,%r12
mov    %rsi,%r13
mov    %rdi,%r15
mov    %r8d,%ebx
lea    -0x1(%rbx),%rdi
mov    %ebx,%eax
and    $0x3,%eax
cmp    $0x3,%rdi
jae    116e <tw+0x3e>
xorpd  %xmm3,%xmm3
xor    %ebp,%ebp
xorpd  %xmm2,%xmm2
jmp    11d4 <tw+0xa4>
mov    %ebx,%edi
and    $0xfffffffc,%edi
xorpd  %xmm3,%xmm3
xor    %ebp,%ebp
xorpd  %xmm2,%xmm2
nopl   (%rax)
movsd  (%r15,%rbp,8),%xmm0
movsd  0x8(%r15,%rbp,8),%xmm1
addsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm0
addsd  %xmm3,%xmm0
addsd  %xmm1,%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
movsd  0x10(%r15,%rbp,8),%xmm0
addsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  0x18(%r15,%rbp,8),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm3,%xmm3
addsd  %xmm0,%xmm3
add    $0x4,%rbp
cmp    %rbp,%rdi
jne    1180 <tw+0x50>
test   %rax,%rax
je     11fb <tw+0xcb>
lea    (%r15,%rbp,8),%rbp
xor    %edi,%edi
nop
movsd  0x0(%rbp,%rdi,8),%xmm0
addsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm0
addsd  %xmm0,%xmm3
add    $0x1,%rdi
cmp    %rdi,%rax
jne    11e0 <tw+0xb0>
test   %r8d,%r8d
jle    1393 <tw+0x263>
xor    %r14d,%r14d
mov    %rcx,0x20(%rsp)
nopl   0x0(%rax)
movsd  %xmm3,0x28(%rsp)
movsd  %xmm2,0x30(%rsp)
lea    0x2(%rbx),%eax
cvtsi2sd %eax,%xmm4
mulsd  %xmm2,%xmm4
mulsd  %xmm2,%xmm4
cvtsi2sd %ebx,%xmm5
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movapd %xmm5,%xmm1
mulsd  %xmm3,%xmm1
subsd  %xmm0,%xmm1
divsd  %xmm1,%xmm4
movsd  %xmm4,0x18(%rsp)
mulsd  (%r15,%r14,1),%xmm5
divsd  %xmm2,%xmm5
movsd  %xmm5,0x8(%rsp)
movapd %xmm4,%xmm0
addsd  0xd97(%rip),%xmm0
cvttsd2si %xmm0,%ebp
mov    %ebp,%edi
call   1050 <sqrt@plt>
movsd  %xmm0,(%rsp)
mov    %ebx,%edi
call   1050 <sqrt@plt>
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
mov    %ebp,%edi
call   1050 <sqrt@plt>
movsd  %xmm0,0x10(%rsp)
mov    %ebx,%edi
call   1050 <sqrt@plt>
addsd  0x10(%rsp),%xmm0
mulsd  (%rsp),%xmm0
divsd  0x18(%rsp),%xmm0
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
mov    %ebp,%edi
call   1050 <sqrt@plt>
movsd  %xmm0,(%rsp)
mov    %ebx,%edi
call   1050 <sqrt@plt>
addsd  (%rsp),%xmm0
divsd  0x18(%rsp),%xmm0
movsd  %xmm0,(%rsp)
mov    %ebp,%edi
call   1050 <sqrt@plt>
movsd  0xd17(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
mov    %ebx,%edi
call   1050 <sqrt@plt>
movsd  0xcfe(%rip),%xmm1
divsd  %xmm0,%xmm1
addsd  0x10(%rsp),%xmm1
movapd %xmm1,%xmm0
movsd  0xcf0(%rip),%xmm1
call   1040 <pow@plt>
mulsd  (%rsp),%xmm0
movsd  0x8(%rsp),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
movapd %xmm1,%xmm0
call   1030 <twtest@plt>
movsd  0x28(%rsp),%xmm3
movsd  0x30(%rsp),%xmm2
mov    0x20(%rsp),%rcx
movsd  %xmm0,0x0(%r13,%r14,1)
movsd  0x8(%rsp),%xmm0
movsd  %xmm0,(%r12,%r14,1)
movsd  0x18(%rsp),%xmm0
movsd  %xmm0,(%rcx,%r14,1)
movsd  (%r15,%r14,1),%xmm0
subsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm0
subsd  %xmm0,%xmm3
add    $0x8,%r14
add    $0xffffffffffffffff,%rbx
jne    1210 <tw+0xe0>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @tw(double* nocapture %0, double* nocapture writeonly %1, double* nocapture writeonly %2, double* nocapture writeonly %3, i32 %4) {
%6 = icmp sgt i32 %4, 0
br i1 %6, label %7, label %117
7: ; preds = %5
%8 = zext i32 %4 to i64
%9 = add nsw i64 %8, -1
%10 = and i64 %8, 3
%11 = icmp ult i64 %9, 3
br i1 %11, label %14, label %12
12: ; preds = %7
%13 = and i64 %8, 4294967292
br label %38
14: ; preds = %38, %7
%15 = phi double [ undef, %7 ], [ %60, %38 ]
%16 = phi double [ undef, %7 ], [ %61, %38 ]
%17 = phi i64 [ 0, %7 ], [ %62, %38 ]
%18 = phi double [ 0.000000e+00, %7 ], [ %61, %38 ]
%19 = phi double [ 0.000000e+00, %7 ], [ %60, %38 ]
%20 = icmp eq i64 %10, 0
br i1 %20, label %33, label %21
21: ; preds = %14, %21
%22 = phi i64 [ %30, %21 ], [ %17, %14 ]
%23 = phi double [ %29, %21 ], [ %18, %14 ]
%24 = phi double [ %28, %21 ], [ %19, %14 ]
%25 = phi i64 [ %31, %21 ], [ 0, %14 ]
%26 = getelementptr inbounds double, double* %0, i64 %22
%27 = load double, double* %26
%28 = fadd double %24, %27
%29 = tail call double @llvm.fmuladd.f64(double %27, double %27, double %23)
%30 = add nuw nsw i64 %22, 1
%31 = add i64 %25, 1
%32 = icmp eq i64 %31, %10
br i1 %32, label %33, label %21
33: ; preds = %21, %14
%34 = phi double [ %15, %14 ], [ %28, %21 ]
%35 = phi double [ %16, %14 ], [ %29, %21 ]
br i1 %6, label %36, label %117
36: ; preds = %33
%37 = zext i32 %4 to i64
br label %65
38: ; preds = %38, %12
%39 = phi i64 [ 0, %12 ], [ %62, %38 ]
%40 = phi double [ 0.000000e+00, %12 ], [ %61, %38 ]
%41 = phi double [ 0.000000e+00, %12 ], [ %60, %38 ]
%42 = phi i64 [ 0, %12 ], [ %63, %38 ]
%43 = getelementptr inbounds double, double* %0, i64 %39
%44 = load double, double* %43
%45 = fadd double %41, %44
%46 = tail call double @llvm.fmuladd.f64(double %44, double %44, double %40)
%47 = or i64 %39, 1
%48 = getelementptr inbounds double, double* %0, i64 %47
%49 = load double, double* %48
%50 = fadd double %45, %49
%51 = tail call double @llvm.fmuladd.f64(double %49, double %49, double %46)
%52 = or i64 %39, 2
%53 = getelementptr inbounds double, double* %0, i64 %52
%54 = load double, double* %53
%55 = fadd double %50, %54
%56 = tail call double @llvm.fmuladd.f64(double %54, double %54, double %51)
%57 = or i64 %39, 3
%58 = getelementptr inbounds double, double* %0, i64 %57
%59 = load double, double* %58
%60 = fadd double %55, %59
%61 = tail call double @llvm.fmuladd.f64(double %59, double %59, double %56)
%62 = add nuw nsw i64 %39, 4
%63 = add i64 %42, 4
%64 = icmp eq i64 %63, %13
br i1 %64, label %14, label %38
65: ; preds = %36, %65
%66 = phi i64 [ 0, %36 ], [ %115, %65 ]
%67 = phi double [ %35, %36 ], [ %114, %65 ]
%68 = phi double [ %34, %36 ], [ %112, %65 ]
%69 = trunc i64 %66 to i32
%70 = sub nsw i32 %4, %69
%71 = add nsw i32 %70, 2
%72 = sitofp i32 %71 to double
%73 = fmul double %68, %72
%74 = fmul double %68, %73
%75 = sitofp i32 %70 to double
%76 = fneg double %68
%77 = fmul double %68, %76
%78 = tail call double @llvm.fmuladd.f64(double %75, double %67, double %77)
%79 = fdiv double %74, %78
%80 = getelementptr inbounds double, double* %0, i64 %66
%81 = load double, double* %80
%82 = fmul double %81, %75
%83 = fdiv double %82, %68
%84 = fadd double %79, -1.000000e+00
%85 = fptosi double %84 to i32
%86 = tail call double @sqrt(i32 %85)
%87 = tail call double @sqrt(i32 %70)
%88 = fadd double %86, %87
%89 = tail call double @sqrt(i32 %85)
%90 = tail call double @sqrt(i32 %70)
%91 = fadd double %89, %90
%92 = fmul double %88, %91
%93 = fdiv double %92, %79
%94 = tail call double @sqrt(i32 %85)
%95 = tail call double @sqrt(i32 %70)
%96 = fadd double %94, %95
%97 = fdiv double %96, %79
%98 = tail call double @sqrt(i32 %85)
%99 = fdiv double 1.000000e+00, %98
%100 = tail call double @sqrt(i32 %70)
%101 = fdiv double 1.000000e+00, %100
%102 = fadd double %99, %101
%103 = tail call double @pow(double %102, double 0x3FD5555555555555)
%104 = fmul double %97, %103
%105 = fsub double %83, %93
%106 = fdiv double %105, %104
%107 = tail call double @twtest(double %106)
%108 = getelementptr inbounds double, double* %1, i64 %66
store double %107, double* %108
%109 = getelementptr inbounds double, double* %2, i64 %66
store double %106, double* %109
%110 = getelementptr inbounds double, double* %3, i64 %66
store double %79, double* %110
%111 = load double, double* %80
%112 = fsub double %68, %111
%113 = fneg double %111
%114 = tail call double @llvm.fmuladd.f64(double %113, double %111, double %67)
%115 = add nuw nsw i64 %66, 1
%116 = icmp eq i64 %115, %37
br i1 %116, label %117, label %65
117: ; preds = %65, %5, %33
ret void
}