void solve_axb_cholesky(int n, double l[n][n], double x[n], double b[n]) {
  int i, j;
  double sum;
  for (i = 0; i < n; i++) {
    sum = 0;
    for (j = 0; j < i; j++)
      sum += l[i][j] * x[j];
    x[i] = (b[i] - sum) / l[i][i];
  }
  for (i = n - 1; i >= 0; i--) {
    sum = 0;
    for (j = i + 1; j < n; j++)
      sum += l[j][i] * x[j];
    x[i] = (x[i] - sum) / l[i][i];
  }
}


x86 assembly (clang, -O0):
<solve_axb_cholesky>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x4(%rbp),%eax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11e2 <solve_axb_cholesky+0xe2>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1195 <solve_axb_cholesky+0x95>
mov    -0x38(%rbp),%rdx
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113c <solve_axb_cholesky+0x3c>
mov    -0x38(%rbp),%rdx
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1121 <solve_axb_cholesky+0x21>
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     12ac <solve_axb_cholesky+0x1ac>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    125f <solve_axb_cholesky+0x15f>
mov    -0x38(%rbp),%rdx
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1206 <solve_axb_cholesky+0x106>
mov    -0x38(%rbp),%rdx
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    11eb <solve_axb_cholesky+0xeb>
pop    %rbp
ret

x86 assembly (clang, -O2):
<solve_axb_cholesky>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rsi,-0x20(%rsp)
test   %edi,%edi
jle    13c0 <solve_axb_cholesky+0x2c0>
mov    %edi,%r14d
mov    -0x20(%rsp),%r12
lea    0x18(%r12),%r15
lea    0x0(,%r14,8),%r8
xor    %r10d,%r10d
mov    %r12,%r9
jmp    1174 <solve_axb_cholesky+0x74>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
xorpd  %xmm0,%xmm0
xor    %r11d,%r11d
movsd  (%rcx,%r10,8),%xmm1
subsd  %xmm0,%xmm1
lea    (%r12,%r11,8),%rax
divsd  (%rax,%r10,8),%xmm1
movsd  %xmm1,(%rdx,%r10,8)
add    $0x1,%r10
add    %r8,%r15
add    %r8,%r9
cmp    %r14,%r10
je     122c <solve_axb_cholesky+0x12c>
test   %r10,%r10
je     1140 <solve_axb_cholesky+0x40>
lea    -0x1(%r10),%rax
cmp    $0x3,%rax
jae    1190 <solve_axb_cholesky+0x90>
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
jmp    11eb <solve_axb_cholesky+0xeb>
nopl   0x0(%rax,%rax,1)
mov    %r10,%rsi
and    $0xfffffffffffffffc,%rsi
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
nopl   (%rax)
movsd  -0x18(%r15,%rbx,8),%xmm1
movsd  -0x10(%r15,%rbx,8),%xmm2
mulsd  (%rdx,%rbx,8),%xmm1
mulsd  0x8(%rdx,%rbx,8),%xmm2
addsd  %xmm0,%xmm1
movsd  -0x8(%r15,%rbx,8),%xmm3
mulsd  0x10(%rdx,%rbx,8),%xmm3
addsd  %xmm1,%xmm2
movsd  (%r15,%rbx,8),%xmm0
mulsd  0x18(%rdx,%rbx,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rbx
cmp    %rbx,%rsi
jne    11a0 <solve_axb_cholesky+0xa0>
mov    %r10,%r11
imul   %r14,%r11
test   $0x3,%r10b
je     1147 <solve_axb_cholesky+0x47>
mov    %r10d,%ebp
and    $0x3,%ebp
lea    (%rdx,%rbx,8),%rsi
lea    (%r9,%rbx,8),%rbx
xor    %eax,%eax
nopl   0x0(%rax)
movsd  (%rbx,%rax,8),%xmm1
mulsd  (%rsi,%rax,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rax
cmp    %rax,%rbp
jne    1210 <solve_axb_cholesky+0x110>
jmp    1147 <solve_axb_cholesky+0x47>
test   %edi,%edi
jle    13c0 <solve_axb_cholesky+0x2c0>
mov    %r14,%rax
neg    %rax
mov    %rax,-0x18(%rsp)
lea    (%rdx,%r14,8),%rdi
lea    0x0(,%r14,8),%r12
lea    0x8(,%r14,8),%rax
imul   %r14,%rax
mov    -0x20(%rsp),%rcx
lea    (%rcx,%rax,1),%r9
add    $0xfffffffffffffff8,%r9
mov    %r12,%rax
xor    $0xfffffffffffffff8,%rax
mov    %rax,-0x8(%rsp)
lea    (%rcx,%r14,8),%r8
add    $0xfffffffffffffff8,%r8
mov    %r14,%rbx
shl    $0x5,%rbx
xor    %r13d,%r13d
mov    %r14,%r15
jmp    12d9 <solve_axb_cholesky+0x1d9>
nopw   0x0(%rax,%rax,1)
movsd  (%rdx,%r15,8),%xmm1
subsd  %xmm0,%xmm1
mov    %r15,%rax
imul   %r14,%rax
mov    -0x20(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
divsd  (%rax,%r15,8),%xmm1
movsd  %xmm1,(%rdx,%r15,8)
add    $0x1,%r13d
addq   $0x1,-0x18(%rsp)
add    $0xfffffffffffffff8,%rdi
add    -0x8(%rsp),%r9
add    $0xfffffffffffffff8,%r8
cmpq   $0x1,-0x10(%rsp)
jle    13c0 <solve_axb_cholesky+0x2c0>
mov    %r15,%rax
add    $0xffffffffffffffff,%r15
xorpd  %xmm0,%xmm0
mov    %rax,-0x10(%rsp)
cmp    %r14,%rax
jge    1290 <solve_axb_cholesky+0x190>
lea    -0x1(%r13),%ebp
test   $0x3,%r13b
je     1338 <solve_axb_cholesky+0x238>
mov    %r13d,%ecx
and    $0x3,%ecx
xorpd  %xmm0,%xmm0
mov    %r9,%rsi
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movsd  (%rsi),%xmm1
mulsd  (%rdi,%rax,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rax
add    %r12,%rsi
cmp    %eax,%ecx
jne    1310 <solve_axb_cholesky+0x210>
sub    -0x18(%rsp),%rax
cmp    $0x3,%ebp
jb     1290 <solve_axb_cholesky+0x190>
jmp    134a <solve_axb_cholesky+0x24a>
xorpd  %xmm0,%xmm0
mov    -0x10(%rsp),%rax
cmp    $0x3,%ebp
jb     1290 <solve_axb_cholesky+0x190>
lea    0x3(%rax),%rbp
imul   %r12,%rbp
lea    0x2(%rax),%r10
imul   %r12,%r10
mov    %r12,%rcx
imul   %rax,%rcx
lea    0x1(%rax),%rsi
imul   %r12,%rsi
mov    %r8,%r11
nopl   0x0(%rax)
movsd  (%r11,%rcx,1),%xmm1
mulsd  (%rdx,%rax,8),%xmm1
addsd  %xmm0,%xmm1
movsd  (%r11,%rsi,1),%xmm0
mulsd  0x8(%rdx,%rax,8),%xmm0
addsd  %xmm1,%xmm0
movsd  (%r11,%r10,1),%xmm1
mulsd  0x10(%rdx,%rax,8),%xmm1
addsd  %xmm0,%xmm1
movsd  (%r11,%rbp,1),%xmm0
mulsd  0x18(%rdx,%rax,8),%xmm0
addsd  %xmm1,%xmm0
add    $0x4,%rax
add    %rbx,%r11
cmp    %eax,%r14d
jne    1370 <solve_axb_cholesky+0x270>
jmp    1290 <solve_axb_cholesky+0x190>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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