void ukf_cholesky_solve(double *A, unsigned n, double *sigma, double *B,
                        unsigned m, double *X) {
  int i, j, k;
  double t;
  for (i = 0; i < m; i++) {
    for (j = 0; j < n; j++) {
      t = B[i * n + j];
      for (k = j - 1; k >= 0; k--)
        t -= A[j * n + k] * X[i * n + k];
      X[i * n + j] = t / sigma[j];
    }
    for (j = n - 1; j >= 0; j--) {
      t = X[i * n + j];
      for (k = j + 1; k < n; k++)
        t -= A[k * n + j] * X[i * n + k];
      X[i * n + j] = t / sigma[j];
    }
  }
}


x86 assembly (clang, -O0):
<ukf_cholesky_solve>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jae    12ef <ukf_cholesky_solve+0x1ef>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    1204 <ukf_cholesky_solve+0x104>
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x38(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x38(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x3c(%rbp)
cmpl   $0x0,-0x3c(%rbp)
jl     11cf <ukf_cholesky_solve+0xcf>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x3c(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x3c(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x3c(%rbp)
jmp    1164 <ukf_cholesky_solve+0x64>
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x38(%rbp),%ecx
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1135 <ukf_cholesky_solve+0x35>
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x0,-0x38(%rbp)
jl     12dc <ukf_cholesky_solve+0x1dc>
mov    -0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x38(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    12a7 <ukf_cholesky_solve+0x1a7>
mov    -0x8(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x38(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x3c(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    123a <ukf_cholesky_solve+0x13a>
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x38(%rbp),%ecx
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
jmp    120d <ukf_cholesky_solve+0x10d>
jmp    12e1 <ukf_cholesky_solve+0x1e1>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1122 <ukf_cholesky_solve+0x22>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ukf_cholesky_solve>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x8(%rsp)
mov    %r8d,-0x2c(%rsp)
test   %r8d,%r8d
je     139b <ukf_cholesky_solve+0x29b>
lea    -0x1(%rsi),%eax
mov    %eax,-0x30(%rsp)
cltq
mov    %rax,-0x18(%rsp)
mov    %esi,%eax
mov    %rax,-0x28(%rsp)
lea    (%rsi,%rsi,1),%r13d
mov    $0xffffffffffffffff,%r15
mov    $0xfffffffffffffffe,%r8
xor    %eax,%eax
mov    %rax,-0x40(%rsp)
movl   $0x0,-0x44(%rsp)
mov    %rdx,-0x10(%rsp)
mov    %rsi,-0x38(%rsp)
jmp    118c <ukf_cholesky_solve+0x8c>
nop
mov    -0x44(%rsp),%ecx
add    $0x1,%ecx
mov    -0x28(%rsp),%rax
add    %rax,%r15
add    %rax,%r8
mov    -0x40(%rsp),%rax
add    %esi,%eax
mov    %rax,-0x40(%rsp)
mov    %ecx,-0x44(%rsp)
cmp    -0x2c(%rsp),%ecx
je     139b <ukf_cholesky_solve+0x29b>
test   %esi,%esi
je     1290 <ukf_cholesky_solve+0x190>
mov    -0x44(%rsp),%eax
mov    %eax,%r10d
imul   %esi,%r10d
mov    $0xffffffffffffffff,%rbp
mov    $0xfffffffffffffffe,%rbx
xor    %r14d,%r14d
jmp    11e4 <ukf_cholesky_solve+0xe4>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
divsd  (%rdx,%r14,8),%xmm0
movsd  %xmm0,(%r9,%r11,8)
add    $0x1,%r14
mov    -0x28(%rsp),%rax
add    %rax,%rbp
add    %rax,%rbx
cmp    %rax,%r14
je     1290 <ukf_cholesky_solve+0x190>
lea    (%r10,%r14,1),%r11d
mov    -0x8(%rsp),%rax
movsd  (%rax,%r11,8),%xmm0
test   %r14,%r14
je     11c0 <ukf_cholesky_solve+0xc0>
test   $0x1,%r14b
jne    1210 <ukf_cholesky_solve+0x110>
mov    %r14,%rax
cmp    $0x1,%r14
je     11c0 <ukf_cholesky_solve+0xc0>
jmp    1235 <ukf_cholesky_solve+0x135>
nopl   0x0(%rax)
mov    %r14d,%ecx
imul   %esi,%ecx
lea    -0x1(%r14),%rax
add    %eax,%ecx
movsd  (%rdi,%rcx,8),%xmm1
lea    (%r10,%rax,1),%ecx
mulsd  (%r9,%rcx,8),%xmm1
subsd  %xmm1,%xmm0
cmp    $0x1,%r14
je     11c0 <ukf_cholesky_solve+0xc0>
add    $0x2,%rax
nopl   0x0(%rax)
lea    (%rax,%rbp,1),%ecx
add    $0xfffffffe,%ecx
movsd  (%rdi,%rcx,8),%xmm1
lea    -0x2(%r15,%rax,1),%ecx
mulsd  (%r9,%rcx,8),%xmm1
lea    -0x2(%rbx,%rax,1),%ecx
movsd  (%rdi,%rcx,8),%xmm2
lea    -0x2(%r8,%rax,1),%ecx
mulsd  (%r9,%rcx,8),%xmm2
subsd  %xmm1,%xmm0
subsd  %xmm2,%xmm0
add    $0xfffffffffffffffe,%rax
cmp    $0x2,%rax
jg     1240 <ukf_cholesky_solve+0x140>
jmp    11c0 <ukf_cholesky_solve+0xc0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
cmpl   $0x0,-0x30(%rsp)
js     1160 <ukf_cholesky_solve+0x60>
mov    -0x44(%rsp),%eax
imul   %esi,%eax
mov    %rax,-0x20(%rsp)
xor    %r10d,%r10d
mov    -0x18(%rsp),%r11
mov    %esi,%r14d
jmp    12ef <ukf_cholesky_solve+0x1ef>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    -0x10(%rsp),%rdx
divsd  (%rdx,%r11,8),%xmm0
movsd  %xmm0,(%r9,%rax,8)
add    $0xffffffffffffffff,%r11
add    $0xffffffff,%r14d
add    $0x1,%r10d
test   $0x80000000,%r11d
mov    -0x38(%rsp),%rsi
jne    1160 <ukf_cholesky_solve+0x60>
mov    -0x20(%rsp),%rax
add    %r11d,%eax
movsd  (%r9,%rax,8),%xmm0
cmp    %esi,%r14d
jae    12c0 <ukf_cholesky_solve+0x1c0>
mov    %r14d,%r12d
test   $0x1,%r10b
je     1331 <ukf_cholesky_solve+0x231>
mov    %r14d,%ecx
imul   -0x38(%rsp),%ecx
add    %r11d,%ecx
movsd  (%rdi,%rcx,8),%xmm1
mov    -0x20(%rsp),%rcx
add    %r14d,%ecx
mulsd  (%r9,%rcx,8),%xmm1
subsd  %xmm1,%xmm0
add    $0x1,%r12
cmp    $0x1,%r10d
je     12c0 <ukf_cholesky_solve+0x1c0>
lea    0x1(%r12),%ebx
mov    -0x38(%rsp),%rsi
imul   %esi,%ebx
add    %r11d,%ebx
mov    -0x40(%rsp),%rcx
add    %r12d,%ecx
mov    %esi,%edx
sub    %r12d,%edx
imul   %esi,%r12d
add    %r11d,%r12d
nopl   0x0(%rax,%rax,1)
mov    %r12d,%ebp
movsd  (%rdi,%rbp,8),%xmm1
mov    %ecx,%ebp
mulsd  (%r9,%rbp,8),%xmm1
subsd  %xmm1,%xmm0
mov    %ebx,%ebp
movsd  (%rdi,%rbp,8),%xmm1
lea    0x1(%rcx),%ebp
mulsd  (%r9,%rbp,8),%xmm1
subsd  %xmm1,%xmm0
add    %r13d,%ebx
add    $0x2,%ecx
add    %r13d,%r12d
add    $0xfffffffe,%edx
jne    1360 <ukf_cholesky_solve+0x260>
jmp    12c0 <ukf_cholesky_solve+0x1c0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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