#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 */
int sqrt(int);
void ComputeHs(int *s, int M, int K, double *h0, double *h1) {
  int i;
  double iC, iL, iR;
  i = 0;
  if (K == 1) {
    iL = 0;
    iC = (double)(s[i] + 1);
    iR = M;
    h0[i] = ((M - iC) * iC / M) * (iR - iL) / ((iR - iC) * ((iC - iL)));
  } else {
    iL = 0;
    iC = (double)(s[i] + 1);
    iR = (double)(s[i + 1] + 1);
    h0[i] = ((M - iC) * iC / M) * (iR - iL) / ((iR - iC) * ((iC - iL)));
    h1[i] = -sqrt((M - iC) * iC * (M - iR) * iR) / (M * (iR - iC));
    for (i = 1; i < K - 1; i++) {
      iL = iC;
      iC = iR;
      iR = (double)(s[i + 1] + 1);
      h0[i] = ((M - iC) * iC / M) * (iR - iL) / ((iR - iC) * ((iC - iL)));
      h1[i] = -sqrt((M - iC) * iC * (M - iR) * iR) / (M * (iR - iC));
    }
    iL = iC;
    iC = iR;
    iR = M;
    h0[i] = ((M - iC) * iC / M) * (iR - iL) / ((iR - iC) * ((iC - iL)));
  }
}


x86 assembly (clang, -O0):
<ComputeHs>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
cmpl   $0x1,-0x10(%rbp)
jne    11b8 <ComputeHs+0xa8>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x30(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x30(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm2
subsd  -0x38(%rbp),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1406 <ComputeHs+0x2f6>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x40(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x30(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm2
subsd  -0x38(%rbp),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0xc(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
subsd  -0x40(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mulsd  -0x40(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sqrt@plt>
mov    %eax,%ecx
xor    %eax,%eax
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm2
subsd  -0x30(%rbp),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1399 <ComputeHs+0x289>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x40(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x30(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm2
subsd  -0x38(%rbp),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0xc(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
subsd  -0x40(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mulsd  -0x40(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sqrt@plt>
mov    %eax,%ecx
xor    %eax,%eax
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm2
subsd  -0x30(%rbp),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    12a2 <ComputeHs+0x192>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x30(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm2
subsd  -0x38(%rbp),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ComputeHs>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rcx,%r14
mov    (%rdi),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm4
cmp    $0x1,%edx
jne    1159 <ComputeHs+0x49>
cvtsi2sd %esi,%xmm0
movapd %xmm0,%xmm1
subsd  %xmm4,%xmm1
mulsd  %xmm4,%xmm1
movapd %xmm1,%xmm2
divsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,(%r14)
jmp    1326 <ComputeHs+0x216>
mov    %r8,%r15
mov    %edx,%ebp
mov    %rdi,%r13
mov    0x4(%rdi),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm3
cvtsi2sd %esi,%xmm2
movapd %xmm2,%xmm0
subsd  %xmm4,%xmm0
movsd  %xmm0,0x30(%rsp)
mulsd  %xmm4,%xmm0
movapd %xmm0,%xmm1
divsd  %xmm2,%xmm1
movapd %xmm2,%xmm5
movsd  %xmm2,0x20(%rsp)
mulsd  %xmm3,%xmm1
movapd %xmm3,%xmm2
subsd  %xmm4,%xmm2
movsd  %xmm2,0x18(%rsp)
movsd  %xmm4,0x10(%rsp)
mulsd  %xmm4,%xmm2
divsd  %xmm2,%xmm1
movsd  %xmm1,(%r14)
movapd %xmm5,%xmm1
subsd  %xmm3,%xmm1
movsd  %xmm1,0x28(%rsp)
mulsd  %xmm1,%xmm0
movsd  %xmm3,0x8(%rsp)
mulsd  %xmm3,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sqrt@plt>
movsd  0x18(%rsp),%xmm3
movsd  0x20(%rsp),%xmm2
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movapd %xmm3,%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%r15)
mov    $0x1,%r12d
cmp    $0x3,%ebp
jl     12fe <ComputeHs+0x1ee>
add    $0xffffffff,%ebp
mov    %ebp,%r12d
lea    -0x1(%r12),%rbp
xor    %ebx,%ebx
movsd  0x8(%rsp),%xmm0
movsd  0x10(%rsp),%xmm4
nopl   0x0(%rax,%rax,1)
movapd %xmm0,%xmm3
movsd  %xmm0,0x10(%rsp)
mov    0x8(%r13,%rbx,4),%eax
add    $0x1,%eax
xorps  %xmm5,%xmm5
cvtsi2sd %eax,%xmm5
movapd %xmm2,%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x30(%rsp)
mulsd  %xmm1,%xmm0
movapd %xmm0,%xmm1
divsd  %xmm2,%xmm1
movapd %xmm5,%xmm2
subsd  %xmm4,%xmm2
mulsd  %xmm1,%xmm2
movapd %xmm5,%xmm6
subsd  %xmm3,%xmm6
movsd  %xmm6,0x18(%rsp)
movapd %xmm3,%xmm1
subsd  %xmm4,%xmm1
mulsd  %xmm6,%xmm1
divsd  %xmm1,%xmm2
movsd  %xmm2,0x8(%r14,%rbx,8)
movsd  0x20(%rsp),%xmm1
movsd  %xmm5,0x8(%rsp)
subsd  %xmm5,%xmm1
movsd  %xmm1,0x28(%rsp)
mulsd  %xmm1,%xmm0
mulsd  %xmm5,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sqrt@plt>
movsd  0x18(%rsp),%xmm3
movsd  0x20(%rsp),%xmm2
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movapd %xmm3,%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%r15,%rbx,8)
movsd  0x8(%rsp),%xmm0
add    $0x1,%rbx
movsd  0x10(%rsp),%xmm4
cmp    %rbx,%rbp
jne    1230 <ComputeHs+0x120>
jmp    1304 <ComputeHs+0x1f4>
movsd  0x8(%rsp),%xmm0
movsd  0x28(%rsp),%xmm1
mulsd  %xmm1,%xmm0
divsd  %xmm2,%xmm0
mulsd  0x30(%rsp),%xmm0
mulsd  %xmm1,%xmm3
divsd  %xmm3,%xmm0
movsd  %xmm0,(%r14,%r12,8)
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @ComputeHs(i32* nocapture %0, i32 %1, i32 %2, double* nocapture writeonly %3, double* nocapture writeonly %4) {
%6 = icmp eq i32 %2, 1
%7 = load i32, i32* %0
%8 = add nsw i32 %7, 1
%9 = sitofp i32 %8 to double
br i1 %6, label %10, label %17
10: ; preds = %5
%11 = sitofp i32 %1 to double
%12 = fsub double %11, %9
%13 = fmul double %12, %9
%14 = fdiv double %13, %11
%15 = fmul double %14, %11
%16 = fdiv double %15, %13
store double %16, double* %3
br label %87
17: ; preds = %5
%18 = getelementptr inbounds i32, i32* %0, i64 1
%19 = load i32, i32* %18
%20 = add nsw i32 %19, 1
%21 = sitofp i32 %20 to double
%22 = sitofp i32 %1 to double
%23 = fsub double %22, %9
%24 = fmul double %23, %9
%25 = fdiv double %24, %22
%26 = fmul double %25, %21
%27 = fsub double %21, %9
%28 = fmul double %27, %9
%29 = fdiv double %26, %28
store double %29, double* %3
%30 = fsub double %22, %21
%31 = fmul double %24, %30
%32 = fmul double %31, %21
%33 = fptosi double %32 to i32
%34 = tail call i32 @sqrt(i32 %33)
%35 = sub nsw i32 0, %34
%36 = sitofp i32 %35 to double
%37 = fmul double %27, %22
%38 = fdiv double %36, %37
store double %38, double* %4
%39 = add i32 %2, -1
%40 = icmp sgt i32 %2, 2
br i1 %40, label %41, label %75
41: ; preds = %17
%42 = zext i32 %39 to i64
br label %43
43: ; preds = %41, %43
%44 = phi i64 [ 1, %41 ], [ %47, %43 ]
%45 = phi double [ %21, %41 ], [ %51, %43 ]
%46 = phi double [ %9, %41 ], [ %45, %43 ]
%47 = add nuw nsw i64 %44, 1
%48 = getelementptr inbounds i32, i32* %0, i64 %47
%49 = load i32, i32* %48
%50 = add nsw i32 %49, 1
%51 = sitofp i32 %50 to double
%52 = fsub double %22, %45
%53 = fmul double %45, %52
%54 = fdiv double %53, %22
%55 = fsub double %51, %46
%56 = fmul double %54, %55
%57 = fsub double %51, %45
%58 = fsub double %45, %46
%59 = fmul double %58, %57
%60 = fdiv double %56, %59
%61 = getelementptr inbounds double, double* %3, i64 %44
store double %60, double* %61
%62 = fsub double %22, %51
%63 = fmul double %53, %62
%64 = fmul double %63, %51
%65 = fptosi double %64 to i32
%66 = tail call i32 @sqrt(i32 %65)
%67 = sub nsw i32 0, %66
%68 = sitofp i32 %67 to double
%69 = fmul double %57, %22
%70 = fdiv double %68, %69
%71 = getelementptr inbounds double, double* %4, i64 %44
store double %70, double* %71
%72 = icmp eq i64 %47, %42
br i1 %72, label %73, label %43
73: ; preds = %43
%74 = zext i32 %39 to i64
br label %75
75: ; preds = %17, %73
%76 = phi double [ %57, %73 ], [ %27, %17 ]
%77 = phi double [ %52, %73 ], [ %23, %17 ]
%78 = phi double [ %62, %73 ], [ %30, %17 ]
%79 = phi i64 [ %74, %73 ], [ 1, %17 ]
%80 = phi double [ %51, %73 ], [ %21, %17 ]
%81 = fmul double %80, %78
%82 = fdiv double %81, %22
%83 = fmul double %77, %82
%84 = fmul double %78, %76
%85 = fdiv double %83, %84
%86 = getelementptr inbounds double, double* %3, i64 %79
store double %85, double* %86
br label %87
87: ; preds = %75, %10
ret void
}