#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 sqrt(double);
int cholesky_decomp(int n, double l[n][n], double a[n][n]) {
  int i, j, k;
  double sum;
  for (i = 0; i < n; i++) {
    for (j = 0; j < i; j++) {
      sum = 0;
      for (k = 0; k < j; k++)
        sum += l[i][k] * l[j][k];
      l[i][j] = (a[i][j] - sum) / l[j][j];
    }
    sum = 0;
    for (k = 0; k < i; k++)
      sum += l[i][k] * l[i][k];
    sum = a[i][i] - sum;
    if (sum < 1e-30)
      return 1;
    l[i][i] = sqrt(sum);
  }
  return 0;
}


x86 assembly (clang, -O0):
<cholesky_decomp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x8(%rbp),%eax
mov    %rax,-0x40(%rbp)
mov    -0x8(%rbp),%eax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1335 <cholesky_decomp+0x225>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    123d <cholesky_decomp+0x12d>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11ce <cholesky_decomp+0xbe>
mov    -0x40(%rbp),%rdx
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1166 <cholesky_decomp+0x56>
mov    -0x40(%rbp),%rdx
mov    -0x38(%rbp),%rsi
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   %rsi,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    114b <cholesky_decomp+0x3b>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12b4 <cholesky_decomp+0x1a4>
mov    -0x40(%rbp),%rdx
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    124c <cholesky_decomp+0x13c>
mov    -0x38(%rbp),%rdx
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xd1a(%rip),%xmm0
ucomisd -0x30(%rbp),%xmm0
jbe    12fd <cholesky_decomp+0x1ed>
movl   $0x1,-0x4(%rbp)
jmp    133c <cholesky_decomp+0x22c>
movsd  -0x30(%rbp),%xmm0
call   1030 <sqrt@plt>
mov    -0x40(%rbp),%rdx
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1138 <cholesky_decomp+0x28>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cholesky_decomp>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rdx,0x18(%rsp)
test   %edi,%edi
jle    13dc <cholesky_decomp+0x2cc>
mov    %rsi,%r13
mov    %edi,%eax
lea    0x18(%rsi),%r14
mov    %rax,(%rsp)
lea    0x0(,%rax,8),%rbx
xor    %r15d,%r15d
movsd  0xeb5(%rip),%xmm4
xorpd  %xmm5,%xmm5
mov    %rsi,%r12
mov    %r14,0x20(%rsp)
jmp    119c <cholesky_decomp+0x8c>
nopl   0x0(%rax)
mov    %rdx,%rbp
call   1030 <sqrt@plt>
mov    %rbp,%rdx
xorpd  %xmm5,%xmm5
movsd  0xe89(%rip),%xmm4
lea    0x0(,%rdx,8),%rax
add    %r13,%rax
movsd  %xmm0,(%rax,%r15,8)
add    $0x1,%r15
add    %rbx,%r14
add    %rbx,%r12
cmp    (%rsp),%r15
je     13dc <cholesky_decomp+0x2cc>
mov    %r15,%rax
sub    $0x1,%rax
jae    11b0 <cholesky_decomp+0xa0>
xor    %edx,%edx
xorpd  %xmm1,%xmm1
jmp    13b0 <cholesky_decomp+0x2a0>
mov    %rax,0x30(%rsp)
mov    %r15d,%eax
and    $0x3,%eax
mov    %rax,0x8(%rsp)
mov    %r15,%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x28(%rsp)
mov    %r15,%rcx
imul   (%rsp),%rcx
lea    0x0(,%rcx,8),%r11
add    %r13,%r11
mov    0x18(%rsp),%rax
mov    %rcx,0x10(%rsp)
lea    (%rax,%rcx,8),%rdi
mov    %r13,%rdx
mov    0x20(%rsp),%rcx
xor    %r8d,%r8d
jmp    123a <cholesky_decomp+0x12a>
nopw   0x0(%rax,%rax,1)
xorpd  %xmm0,%xmm0
xor    %esi,%esi
movsd  (%rdi,%r8,8),%xmm1
subsd  %xmm0,%xmm1
lea    0x0(,%rsi,8),%rax
add    %r13,%rax
divsd  (%rax,%r8,8),%xmm1
movsd  %xmm1,(%r11,%r8,8)
add    $0x1,%r8
add    %rbx,%rcx
add    %rbx,%rdx
cmp    %r15,%r8
je     1300 <cholesky_decomp+0x1f0>
test   %r8,%r8
je     1200 <cholesky_decomp+0xf0>
lea    -0x1(%r8),%rax
cmp    $0x3,%rax
jae    1260 <cholesky_decomp+0x150>
xorpd  %xmm0,%xmm0
xor    %ebp,%ebp
jmp    12bb <cholesky_decomp+0x1ab>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r8,%rsi
and    $0xfffffffffffffffc,%rsi
xorpd  %xmm0,%xmm0
xor    %ebp,%ebp
nopl   (%rax)
movsd  -0x18(%r14,%rbp,8),%xmm1
movsd  -0x10(%r14,%rbp,8),%xmm2
mulsd  -0x18(%rcx,%rbp,8),%xmm1
mulsd  -0x10(%rcx,%rbp,8),%xmm2
addsd  %xmm0,%xmm1
movsd  -0x8(%r14,%rbp,8),%xmm3
mulsd  -0x8(%rcx,%rbp,8),%xmm3
addsd  %xmm1,%xmm2
movsd  (%r14,%rbp,8),%xmm0
mulsd  (%rcx,%rbp,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rbp
cmp    %rbp,%rsi
jne    1270 <cholesky_decomp+0x160>
mov    %r8,%rsi
imul   (%rsp),%rsi
test   $0x3,%r8b
je     1206 <cholesky_decomp+0xf6>
mov    %r8d,%eax
and    $0x3,%eax
lea    (%rdx,%rbp,8),%r9
lea    (%r12,%rbp,8),%rbp
xor    %r10d,%r10d
xchg   %ax,%ax
movsd  0x0(%rbp,%r10,8),%xmm1
mulsd  (%r9,%r10,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%r10
cmp    %r10,%rax
jne    12e0 <cholesky_decomp+0x1d0>
jmp    1206 <cholesky_decomp+0xf6>
nop
test   %r15,%r15
je     131f <cholesky_decomp+0x20f>
cmpq   $0x3,0x30(%rsp)
jae    132a <cholesky_decomp+0x21a>
xorpd  %xmm1,%xmm1
xor    %eax,%eax
mov    0x10(%rsp),%rdx
mov    0x8(%rsp),%rsi
jmp    1384 <cholesky_decomp+0x274>
xorpd  %xmm1,%xmm1
xor    %edx,%edx
jmp    13b0 <cholesky_decomp+0x2a0>
xorpd  %xmm1,%xmm1
xor    %eax,%eax
mov    0x10(%rsp),%rdx
mov    0x8(%rsp),%rsi
mov    0x28(%rsp),%rcx
nop
movsd  -0x18(%r14,%rax,8),%xmm0
movsd  -0x10(%r14,%rax,8),%xmm2
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
mulsd  %xmm2,%xmm2
addsd  %xmm0,%xmm2
movsd  -0x8(%r14,%rax,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
movsd  (%r14,%rax,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
add    $0x4,%rax
cmp    %rax,%rcx
jne    1340 <cholesky_decomp+0x230>
test   $0x3,%r15b
je     13b0 <cholesky_decomp+0x2a0>
lea    (%r12,%rax,8),%rax
xor    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm0,%xmm1
add    $0x1,%rcx
cmp    %rcx,%rsi
jne    1390 <cholesky_decomp+0x280>
cs nopw 0x0(%rax,%rax,1)
mov    0x18(%rsp),%rax
lea    (%rax,%rdx,8),%rax
movsd  (%rax,%r15,8),%xmm0
subsd  %xmm1,%xmm0
ucomisd %xmm0,%xmm4
ja     13e0 <cholesky_decomp+0x2d0>
ucomisd %xmm5,%xmm0
jb     1160 <cholesky_decomp+0x50>
sqrtsd %xmm0,%xmm0
jmp    1177 <cholesky_decomp+0x67>
xor    %eax,%eax
jmp    13e5 <cholesky_decomp+0x2d5>
mov    $0x1,%eax
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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