#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 *allocate_vector(int /*<<< orphan*/, int);
int /*<<< orphan*/ free(double *);
double sqrt(double);
double f3_corr(double **P, int Ng) {
  int i, j;
  double sum_sqrx = 0, sum_sqry = 0, tmp, *px;
  double meanx = 0, meany = 0, stddevx, stddevy;
  px = allocate_vector(0, Ng);
  for (i = 0; i < Ng; ++i)
    px[i] = 0;
  for (i = 0; i < Ng; ++i)
    for (j = 0; j < Ng; ++j)
      px[i] += P[i][j];
  for (i = 0; i < Ng; ++i) {
    meanx += px[i] * i;
    sum_sqrx += px[i] * i * i;
  }
  meany = meanx;
  sum_sqry = sum_sqrx;
  stddevx = sqrt(sum_sqrx - (meanx * meanx));
  stddevy = stddevx;
  for (tmp = 0, i = 0; i < Ng; ++i)
    for (j = 0; j < Ng; ++j)
      tmp += i * j * P[i][j];
  free(px);
  return (tmp - meanx * meany) / (stddevx * stddevy);
}


x86 assembly (clang, -O0):
<f3_corr>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0xc(%rbp),%esi
xor    %edi,%edi
call   1040 <allocate_vector@plt>
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    119e <f3_corr+0x6e>
mov    -0x38(%rbp),%rax
movslq -0x10(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1174 <f3_corr+0x44>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    120c <f3_corr+0xdc>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11f9 <f3_corr+0xc9>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x10(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    11b8 <f3_corr+0x88>
jmp    11fe <f3_corr+0xce>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    11a5 <f3_corr+0x75>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    127e <f3_corr+0x14e>
mov    -0x38(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x10(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x10(%rbp),%xmm2
movsd  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1213 <f3_corr+0xe3>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <sqrt@plt>
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1355 <f3_corr+0x225>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1342 <f3_corr+0x212>
mov    -0x10(%rbp),%eax
imul   -0x14(%rbp),%eax
cvtsi2sd %eax,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    12f6 <f3_corr+0x1c6>
jmp    1347 <f3_corr+0x217>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    12e3 <f3_corr+0x1b3>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
movsd  -0x30(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x50(%rbp),%xmm1
mulsd  -0x58(%rbp),%xmm1
divsd  %xmm1,%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<f3_corr>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x30,%rsp
mov    %esi,%r15d
mov    %rdi,%r14
xor    %edi,%edi
call   1050 <allocate_vector@plt>
mov    %rax,%rbx
xorpd  %xmm0,%xmm0
mov    %r15d,%r12d
test   %r15d,%r15d
jle    12c7 <f3_corr+0x187>
lea    0x0(,%r12,8),%rdx
xor    %ebp,%ebp
mov    %rbx,%rdi
xor    %esi,%esi
call   1040 <memset@plt>
lea    -0x1(%r12),%rax
mov    %r12d,%ecx
and    $0x3,%ecx
mov    %r12d,%edx
and    $0xfffffffc,%edx
jmp    11a9 <f3_corr+0x69>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rbp
cmp    %r12,%rbp
je     1215 <f3_corr+0xd5>
mov    (%r14,%rbp,8),%rsi
movsd  (%rbx,%rbp,8),%xmm0
xor    %edi,%edi
cmp    $0x3,%rax
jb     11f4 <f3_corr+0xb4>
nopw   0x0(%rax,%rax,1)
addsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,(%rbx,%rbp,8)
addsd  0x8(%rsi,%rdi,8),%xmm0
movsd  %xmm0,(%rbx,%rbp,8)
addsd  0x10(%rsi,%rdi,8),%xmm0
movsd  %xmm0,(%rbx,%rbp,8)
addsd  0x18(%rsi,%rdi,8),%xmm0
movsd  %xmm0,(%rbx,%rbp,8)
add    $0x4,%rdi
cmp    %rdi,%rdx
jne    11c0 <f3_corr+0x80>
test   %rcx,%rcx
je     11a0 <f3_corr+0x60>
lea    (%rsi,%rdi,8),%rsi
xor    %edi,%edi
nop
addsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,(%rbx,%rbp,8)
add    $0x1,%rdi
cmp    %rdi,%rcx
jne    1200 <f3_corr+0xc0>
jmp    11a0 <f3_corr+0x60>
test   %r15d,%r15d
xorpd  %xmm0,%xmm0
jle    12c7 <f3_corr+0x187>
test   %rax,%rax
je     144a <f3_corr+0x30a>
mov    %r12d,%ecx
and    $0xfffffffe,%ecx
xorpd  %xmm5,%xmm5
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  (%rbx,%rax,8),%xmm2
movsd  0x8(%rbx,%rax,8),%xmm3
movapd %xmm2,%xmm4
mulsd  %xmm1,%xmm4
unpcklpd %xmm2,%xmm4
unpcklpd %xmm1,%xmm1
mulpd  %xmm4,%xmm1
lea    0x1(%rax),%edx
xorps  %xmm0,%xmm0
cvtsi2sd %edx,%xmm0
addpd  %xmm5,%xmm1
movapd %xmm3,%xmm2
mulsd  %xmm0,%xmm2
unpcklpd %xmm3,%xmm2
unpcklpd %xmm0,%xmm0
mulpd  %xmm2,%xmm0
movapd %xmm0,%xmm5
addpd  %xmm1,%xmm5
add    $0x2,%rax
cmp    %rcx,%rax
jne    1240 <f3_corr+0x100>
addpd  %xmm1,%xmm0
test   $0x1,%r12b
je     12c7 <f3_corr+0x187>
movsd  (%rbx,%rax,8),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movapd %xmm1,%xmm2
mulsd  %xmm0,%xmm2
unpcklpd %xmm1,%xmm2
unpcklpd %xmm0,%xmm0
mulpd  %xmm2,%xmm0
addpd  %xmm5,%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
movapd 0xd29(%rip),%xmm2
xorpd  %xmm1,%xmm2
movapd %xmm2,0x20(%rsp)
movapd %xmm1,0x10(%rsp)
mulsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     12ff <f3_corr+0x1bf>
sqrtsd %xmm0,%xmm0
jmp    1308 <f3_corr+0x1c8>
call   1060 <sqrt@plt>
xorpd  %xmm1,%xmm1
movsd  %xmm0,(%rsp)
test   %r15d,%r15d
jle    1410 <f3_corr+0x2d0>
lea    -0x1(%r12),%r8
mov    %r12d,%r10d
and    $0x3,%r10d
mov    %r12d,%r9d
and    $0xfffffffc,%r9d
xorpd  %xmm1,%xmm1
xor    %r11d,%r11d
xor    %r15d,%r15d
xor    %ebp,%ebp
xor    %eax,%eax
jmp    1358 <f3_corr+0x218>
nopl   0x0(%rax)
add    $0x1,%rax
add    $0x2,%ebp
add    $0x4,%r15d
add    $0x3,%r11d
cmp    %r12,%rax
je     1410 <f3_corr+0x2d0>
mov    (%r14,%rax,8),%rcx
cmp    $0x3,%r8
jae    1370 <f3_corr+0x230>
xor    %edx,%edx
jmp    13da <f3_corr+0x29a>
cs nopw 0x0(%rax,%rax,1)
xor    %esi,%esi
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
mulsd  (%rcx,%rdx,8),%xmm0
addsd  %xmm1,%xmm0
lea    (%rax,%rsi,1),%edi
xorps  %xmm1,%xmm1
cvtsi2sd %edi,%xmm1
mulsd  0x8(%rcx,%rdx,8),%xmm1
addsd  %xmm0,%xmm1
lea    0x0(%rbp,%rsi,1),%edi
xorps  %xmm0,%xmm0
cvtsi2sd %edi,%xmm0
mulsd  0x10(%rcx,%rdx,8),%xmm0
addsd  %xmm1,%xmm0
lea    (%r11,%rsi,1),%edi
xorps  %xmm1,%xmm1
cvtsi2sd %edi,%xmm1
mulsd  0x18(%rcx,%rdx,8),%xmm1
addsd  %xmm0,%xmm1
add    $0x4,%rdx
add    %r15d,%esi
cmp    %rdx,%r9
jne    1380 <f3_corr+0x240>
test   %r10,%r10
je     1340 <f3_corr+0x200>
lea    (%rcx,%rdx,8),%rcx
imul   %eax,%edx
xor    %esi,%esi
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
cvtsi2sd %edx,%xmm0
mulsd  (%rcx,%rsi,8),%xmm0
addsd  %xmm0,%xmm1
add    $0x1,%rsi
add    %eax,%edx
cmp    %rsi,%r10
jne    13f0 <f3_corr+0x2b0>
jmp    1340 <f3_corr+0x200>
movsd  %xmm1,0x8(%rsp)
mov    %rbx,%rdi
call   1030 <free@plt>
movapd 0x20(%rsp),%xmm0
mulsd  0x10(%rsp),%xmm0
addsd  0x8(%rsp),%xmm0
movsd  (%rsp),%xmm1
mulsd  %xmm1,%xmm1
divsd  %xmm1,%xmm0
add    $0x30,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
xorpd  %xmm5,%xmm5
xor    %eax,%eax
test   $0x1,%r12b
jne    12a3 <f3_corr+0x163>
jmp    12c7 <f3_corr+0x187>

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

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