void mod_chol_solve(const double *L_, const int *P, const double *b, double *x,
                    int n, double *y) {
  int i, k;
  for (k = 0; k < n; k++)
    x[k] = b[P[k]];
  for (k = 0; k < n; k++) {
    for (i = 0; i < k; i++)
      x[k] -= x[i] * L_[(((k) * ((k) + 1)) / 2 + (i))];
    x[k] /= L_[(((k) * ((k) + 1)) / 2 + (k))];
  }
  for (k = n - 1; k >= 0; k--) {
    for (i = k + 1; i < n; i++)
      x[k] -= x[i] * L_[(((i) * ((i) + 1)) / 2 + (k))];
    x[k] /= L_[(((k) * ((k) + 1)) / 2 + (k))];
  }
  for (k = 0; k < n; k++)
    y[P[k]] = x[k];
  for (k = 0; k < n; k++)
    x[k] = y[k];
}


x86 assembly (clang, -O0):
<mod_chol_solve>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    115f <mod_chol_solve+0x5f>
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x38(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1123 <mod_chol_solve+0x23>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    124f <mod_chol_solve+0x14f>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    11fe <mod_chol_solve+0xfe>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x40(%rbp),%rax
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1179 <mod_chol_solve+0x79>
mov    -0x8(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x48(%rbp),%rax
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1166 <mod_chol_solve+0x66>
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x0,-0x38(%rbp)
jl     1341 <mod_chol_solve+0x241>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    12f0 <mod_chol_solve+0x1f0>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x50(%rbp),%rax
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    126b <mod_chol_solve+0x16b>
mov    -0x8(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x58(%rbp),%rax
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
jmp    1258 <mod_chol_solve+0x158>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1384 <mod_chol_solve+0x284>
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x38(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1348 <mod_chol_solve+0x248>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    13bf <mod_chol_solve+0x2bf>
mov    -0x30(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    138b <mod_chol_solve+0x28b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<mod_chol_solve>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %r8d,%r8d
jle    1568 <mod_chol_solve+0x468>
mov    %r8d,%r14d
lea    -0x1(%r14),%r10
mov    %r14d,%r11d
and    $0x3,%r11d
cmp    $0x3,%r10
jae    112b <mod_chol_solve+0x2b>
xor    %ebx,%ebx
jmp    1187 <mod_chol_solve+0x87>
mov    %r14d,%r15d
and    $0xfffffffc,%r15d
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movslq (%rsi,%rbx,4),%rax
movsd  (%rdx,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%rbx,8)
movslq 0x4(%rsi,%rbx,4),%rax
movsd  (%rdx,%rax,8),%xmm0
movsd  %xmm0,0x8(%rcx,%rbx,8)
movslq 0x8(%rsi,%rbx,4),%rax
movsd  (%rdx,%rax,8),%xmm0
movsd  %xmm0,0x10(%rcx,%rbx,8)
movslq 0xc(%rsi,%rbx,4),%rax
movsd  (%rdx,%rax,8),%xmm0
movsd  %xmm0,0x18(%rcx,%rbx,8)
add    $0x4,%rbx
cmp    %rbx,%r15
jne    1140 <mod_chol_solve+0x40>
test   %r11,%r11
je     11b8 <mod_chol_solve+0xb8>
lea    (%rcx,%rbx,8),%r15
lea    (%rsi,%rbx,4),%r12
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
movslq (%r12,%rax,4),%rbx
movsd  (%rdx,%rbx,8),%xmm0
movsd  %xmm0,(%r15,%rax,8)
add    $0x1,%rax
cmp    %rax,%r11
jne    11a0 <mod_chol_solve+0xa0>
test   %r8d,%r8d
jle    1568 <mod_chol_solve+0x468>
lea    0x8(%rdi),%r11
mov    $0x2,%r12d
xor    %r15d,%r15d
xor    %eax,%eax
jmp    120f <mod_chol_solve+0x10f>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rcx),%xmm0
lea    0x1(%rax),%rdx
mov    %edx,%ebx
imul   %eax,%ebx
shr    %ebx
add    %rax,%rbx
divsd  (%rdi,%rbx,8),%xmm0
movsd  %xmm0,(%rcx,%rax,8)
add    %r12,%r15
add    $0x2,%r12
mov    %rdx,%rax
cmp    %r14,%rdx
je     12ac <mod_chol_solve+0x1ac>
test   %rax,%rax
je     11e0 <mod_chol_solve+0xe0>
movsd  (%rcx,%rax,8),%xmm0
cmp    $0x1,%rax
jne    1230 <mod_chol_solve+0x130>
xor    %ebx,%ebx
jmp    1281 <mod_chol_solve+0x181>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r15d,%edx
and    $0xfffffffe,%edx
lea    (%r11,%rdx,4),%rdx
mov    %rax,%r13
and    $0xfffffffffffffffe,%r13
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rcx,%rbx,8),%xmm1
mulsd  -0x8(%rdx,%rbx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%rax,8)
movsd  0x8(%rcx,%rbx,8),%xmm1
mulsd  (%rdx,%rbx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%rax,8)
add    $0x2,%rbx
cmp    %rbx,%r13
jne    1250 <mod_chol_solve+0x150>
test   $0x1,%al
je     11e4 <mod_chol_solve+0xe4>
lea    0x1(%rax),%edx
imul   %eax,%edx
shr    %edx
movsd  (%rcx,%rbx,8),%xmm1
add    %rbx,%rdx
mulsd  (%rdi,%rdx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%rax,8)
jmp    11e4 <mod_chol_solve+0xe4>
test   %r8d,%r8d
jle    1568 <mod_chol_solve+0x468>
xor    %r11d,%r11d
mov    %r14,%r15
jmp    12ee <mod_chol_solve+0x1ee>
nopl   (%rax)
mov    %r12d,%eax
imul   %r15d,%eax
mov    %eax,%edx
shr    $0x1f,%edx
add    %eax,%edx
sar    %edx
movslq %edx,%rax
add    %r12,%rax
divsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%r12,8)
add    $0x1,%r11d
cmp    $0x1,%r15
mov    %r12,%r15
jle    1347 <mod_chol_solve+0x247>
lea    -0x1(%r15),%r12
movsd  -0x8(%rcx,%r15,8),%xmm0
mov    %r11d,%eax
mov    %r15,%rbx
cmp    %r14,%r15
jge    12c0 <mod_chol_solve+0x1c0>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rcx,%rbx,8),%xmm1
lea    0x1(%rbx),%rdx
imul   %edx,%ebx
mov    %ebx,%ebp
shr    $0x1f,%ebp
add    %ebx,%ebp
sar    %ebp
movslq %ebp,%rbx
add    %r12,%rbx
mulsd  (%rdi,%rbx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%r12,8)
mov    %rdx,%rbx
add    $0xffffffff,%eax
jne    1310 <mod_chol_solve+0x210>
jmp    12c0 <mod_chol_solve+0x1c0>
test   %r8d,%r8d
jle    1568 <mod_chol_solve+0x468>
mov    %r14d,%eax
and    $0x3,%eax
cmp    $0x3,%r10
jae    1360 <mod_chol_solve+0x260>
xor    %edx,%edx
jmp    13bb <mod_chol_solve+0x2bb>
mov    %r14d,%edi
and    $0xfffffffc,%edi
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%rdx,8),%xmm0
movslq (%rsi,%rdx,4),%rbp
movsd  %xmm0,(%r9,%rbp,8)
movsd  0x8(%rcx,%rdx,8),%xmm0
movslq 0x4(%rsi,%rdx,4),%rbp
movsd  %xmm0,(%r9,%rbp,8)
movsd  0x10(%rcx,%rdx,8),%xmm0
movslq 0x8(%rsi,%rdx,4),%rbp
movsd  %xmm0,(%r9,%rbp,8)
movsd  0x18(%rcx,%rdx,8),%xmm0
movslq 0xc(%rsi,%rdx,4),%rbp
movsd  %xmm0,(%r9,%rbp,8)
add    $0x4,%rdx
cmp    %rdx,%rdi
jne    1370 <mod_chol_solve+0x270>
test   %rax,%rax
je     13e8 <mod_chol_solve+0x2e8>
lea    (%rsi,%rdx,4),%rsi
lea    (%rcx,%rdx,8),%rdx
xor    %edi,%edi
nopw   0x0(%rax,%rax,1)
movsd  (%rdx,%rdi,8),%xmm0
movslq (%rsi,%rdi,4),%rbp
movsd  %xmm0,(%r9,%rbp,8)
add    $0x1,%rdi
cmp    %rdi,%rax
jne    13d0 <mod_chol_solve+0x2d0>
test   %r8d,%r8d
jle    1568 <mod_chol_solve+0x468>
cmp    $0x4,%r8d
jb     1411 <mod_chol_solve+0x311>
lea    (%r9,%r14,8),%rax
cmp    %rcx,%rax
jbe    1490 <mod_chol_solve+0x390>
lea    (%rcx,%r14,8),%rax
cmp    %r9,%rax
jbe    1490 <mod_chol_solve+0x390>
xor    %edx,%edx
mov    %rdx,%rax
not    %rax
add    %r14,%rax
mov    %r14,%rsi
and    $0x3,%rsi
je     1445 <mod_chol_solve+0x345>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%r9,%rdx,8),%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
add    $0x1,%rdx
add    $0xffffffffffffffff,%rsi
jne    1430 <mod_chol_solve+0x330>
cmp    $0x3,%rax
jb     1568 <mod_chol_solve+0x468>
nop
movsd  (%r9,%rdx,8),%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
movsd  0x8(%r9,%rdx,8),%xmm0
movsd  %xmm0,0x8(%rcx,%rdx,8)
movsd  0x10(%r9,%rdx,8),%xmm0
movsd  %xmm0,0x10(%rcx,%rdx,8)
movsd  0x18(%r9,%rdx,8),%xmm0
movsd  %xmm0,0x18(%rcx,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%r14
jne    1450 <mod_chol_solve+0x350>
jmp    1568 <mod_chol_solve+0x468>
mov    %r14d,%edx
and    $0xfffffffc,%edx
lea    -0x4(%rdx),%rdi
mov    %rdi,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %esi,%eax
and    $0x3,%eax
cmp    $0xc,%rdi
jae    14b4 <mod_chol_solve+0x3b4>
xor    %edi,%edi
jmp    1524 <mod_chol_solve+0x424>
and    $0xfffffffffffffffc,%rsi
xor    %edi,%edi
nopw   0x0(%rax,%rax,1)
movups (%r9,%rdi,8),%xmm0
movups 0x10(%r9,%rdi,8),%xmm1
movups %xmm0,(%rcx,%rdi,8)
movups %xmm1,0x10(%rcx,%rdi,8)
movups 0x20(%r9,%rdi,8),%xmm0
movups 0x30(%r9,%rdi,8),%xmm1
movups %xmm0,0x20(%rcx,%rdi,8)
movups %xmm1,0x30(%rcx,%rdi,8)
movups 0x40(%r9,%rdi,8),%xmm0
movups 0x50(%r9,%rdi,8),%xmm1
movups %xmm0,0x40(%rcx,%rdi,8)
movups %xmm1,0x50(%rcx,%rdi,8)
movupd 0x60(%r9,%rdi,8),%xmm0
movupd 0x70(%r9,%rdi,8),%xmm1
movupd %xmm0,0x60(%rcx,%rdi,8)
movupd %xmm1,0x70(%rcx,%rdi,8)
add    $0x10,%rdi
add    $0xfffffffffffffffc,%rsi
jne    14c0 <mod_chol_solve+0x3c0>
test   %rax,%rax
je     155f <mod_chol_solve+0x45f>
lea    (%rcx,%rdi,8),%rsi
add    $0x10,%rsi
lea    (%r9,%rdi,8),%rdi
add    $0x10,%rdi
shl    $0x5,%rax
xor    %ebx,%ebx
nop
movupd -0x10(%rdi,%rbx,1),%xmm0
movupd (%rdi,%rbx,1),%xmm1
movupd %xmm0,-0x10(%rsi,%rbx,1)
movupd %xmm1,(%rsi,%rbx,1)
add    $0x20,%rbx
cmp    %rbx,%rax
jne    1540 <mod_chol_solve+0x440>
cmp    %r14,%rdx
jne    1413 <mod_chol_solve+0x313>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @mod_chol_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
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* %14
br label %15
15: ; preds = %33, %6
%16 = load i32, i32* %14
%17 = load i32, i32* %11
%18 = icmp slt i32 %16, %17
br i1 %18, label %19, label %36
19: ; preds = %15
%20 = load double*, double** %9
%21 = load i32*, i32** %8
%22 = load i32, i32* %14
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds i32, i32* %21, i64 %23
%25 = load i32, i32* %24
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double, double* %20, i64 %26
%28 = load double, double* %27
%29 = load double*, double** %10
%30 = load i32, i32* %14
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds double, double* %29, i64 %31
store double %28, double* %32
br label %33
33: ; preds = %19
%34 = load i32, i32* %14
%35 = add nsw i32 %34, 1
store i32 %35, i32* %14
br label %15
36: ; preds = %15
store i32 0, i32* %14
br label %37
37: ; preds = %91, %36
%38 = load i32, i32* %14
%39 = load i32, i32* %11
%40 = icmp slt i32 %38, %39
br i1 %40, label %41, label %94
41: ; preds = %37
store i32 0, i32* %13
br label %42
42: ; preds = %70, %41
%43 = load i32, i32* %13
%44 = load i32, i32* %14
%45 = icmp slt i32 %43, %44
br i1 %45, label %46, label %73
46: ; preds = %42
%47 = load double*, double** %10
%48 = load i32, i32* %13
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %47, i64 %49
%51 = load double, double* %50
%52 = load double*, double** %7
%53 = load i32, i32* %14
%54 = load i32, i32* %14
%55 = add nsw i32 %54, 1
%56 = mul nsw i32 %53, %55
%57 = sdiv i32 %56, 2
%58 = load i32, i32* %13
%59 = add nsw i32 %57, %58
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %52, i64 %60
%62 = load double, double* %61
%63 = load double*, double** %10
%64 = load i32, i32* %14
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
%67 = load double, double* %66
%68 = fneg double %51
%69 = call double @llvm.fmuladd.f64(double %68, double %62, double %67)
store double %69, double* %66
br label %70
70: ; preds = %46
%71 = load i32, i32* %13
%72 = add nsw i32 %71, 1
store i32 %72, i32* %13
br label %42
73: ; preds = %42
%74 = load double*, double** %7
%75 = load i32, i32* %14
%76 = load i32, i32* %14
%77 = add nsw i32 %76, 1
%78 = mul nsw i32 %75, %77
%79 = sdiv i32 %78, 2
%80 = load i32, i32* %14
%81 = add nsw i32 %79, %80
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds double, double* %74, i64 %82
%84 = load double, double* %83
%85 = load double*, double** %10
%86 = load i32, i32* %14
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %85, i64 %87
%89 = load double, double* %88
%90 = fdiv double %89, %84
store double %90, double* %88
br label %91
91: ; preds = %73
%92 = load i32, i32* %14
%93 = add nsw i32 %92, 1
store i32 %93, i32* %14
br label %37
94: ; preds = %37
%95 = load i32, i32* %11
%96 = sub nsw i32 %95, 1
store i32 %96, i32* %14
br label %97
97: ; preds = %152, %94
%98 = load i32, i32* %14
%99 = icmp sge i32 %98, 0
br i1 %99, label %100, label %155
100: ; preds = %97
%101 = load i32, i32* %14
%102 = add nsw i32 %101, 1
store i32 %102, i32* %13
br label %103
103: ; preds = %131, %100
%104 = load i32, i32* %13
%105 = load i32, i32* %11
%106 = icmp slt i32 %104, %105
br i1 %106, label %107, label %134
107: ; preds = %103
%108 = load double*, double** %10
%109 = load i32, i32* %13
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds double, double* %108, i64 %110
%112 = load double, double* %111
%113 = load double*, double** %7
%114 = load i32, i32* %13
%115 = load i32, i32* %13
%116 = add nsw i32 %115, 1
%117 = mul nsw i32 %114, %116
%118 = sdiv i32 %117, 2
%119 = load i32, i32* %14
%120 = add nsw i32 %118, %119
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds double, double* %113, i64 %121
%123 = load double, double* %122
%124 = load double*, double** %10
%125 = load i32, i32* %14
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds double, double* %124, i64 %126
%128 = load double, double* %127
%129 = fneg double %112
%130 = call double @llvm.fmuladd.f64(double %129, double %123, double %128)
store double %130, double* %127
br label %131
131: ; preds = %107
%132 = load i32, i32* %13
%133 = add nsw i32 %132, 1
store i32 %133, i32* %13
br label %103
134: ; preds = %103
%135 = load double*, double** %7
%136 = load i32, i32* %14
%137 = load i32, i32* %14
%138 = add nsw i32 %137, 1
%139 = mul nsw i32 %136, %138
%140 = sdiv i32 %139, 2
%141 = load i32, i32* %14
%142 = add nsw i32 %140, %141
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds double, double* %135, i64 %143
%145 = load double, double* %144
%146 = load double*, double** %10
%147 = load i32, i32* %14
%148 = sext i32 %147 to i64
%149 = getelementptr inbounds double, double* %146, i64 %148
%150 = load double, double* %149
%151 = fdiv double %150, %145
store double %151, double* %149
br label %152
152: ; preds = %134
%153 = load i32, i32* %14
%154 = add nsw i32 %153, -1
store i32 %154, i32* %14
br label %97
155: ; preds = %97
store i32 0, i32* %14
br label %156
156: ; preds = %174, %155
%157 = load i32, i32* %14
%158 = load i32, i32* %11
%159 = icmp slt i32 %157, %158
br i1 %159, label %160, label %177
160: ; preds = %156
%161 = load double*, double** %10
%162 = load i32, i32* %14
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds double, double* %161, i64 %163
%165 = load double, double* %164
%166 = load double*, double** %12
%167 = load i32*, i32** %8
%168 = load i32, i32* %14
%169 = sext i32 %168 to i64
%170 = getelementptr inbounds i32, i32* %167, i64 %169
%171 = load i32, i32* %170
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds double, double* %166, i64 %172
store double %165, double* %173
br label %174
174: ; preds = %160
%175 = load i32, i32* %14
%176 = add nsw i32 %175, 1
store i32 %176, i32* %14
br label %156
177: ; preds = %156
store i32 0, i32* %14
br label %178
178: ; preds = %192, %177
%179 = load i32, i32* %14
%180 = load i32, i32* %11
%181 = icmp slt i32 %179, %180
br i1 %181, label %182, label %195
182: ; preds = %178
%183 = load double*, double** %12
%184 = load i32, i32* %14
%185 = sext i32 %184 to i64
%186 = getelementptr inbounds double, double* %183, i64 %185
%187 = load double, double* %186
%188 = load double*, double** %10
%189 = load i32, i32* %14
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds double, double* %188, i64 %190
store double %187, double* %191
br label %192
192: ; preds = %182
%193 = load i32, i32* %14
%194 = add nsw i32 %193, 1
store i32 %194, i32* %14
br label %178
195: ; preds = %178
ret void
}

LLVM IR (-O2):
define void @mod_chol_solve(double* nocapture %0, i32* nocapture %1, double* nocapture %2, double* nocapture %3, i32 %4, double* nocapture %5) {
%7 = icmp sgt i32 %4, 0
br i1 %7, label %8, label %349
8: ; preds = %6
%9 = zext i32 %4 to i64
%10 = add nsw i64 %9, -1
%11 = and i64 %9, 3
%12 = icmp ult i64 %10, 3
br i1 %12, label %15, label %13
13: ; preds = %8
%14 = and i64 %9, 4294967292
br label %33
15: ; preds = %33, %8
%16 = phi i64 [ 0, %8 ], [ %63, %33 ]
%17 = icmp eq i64 %11, 0
br i1 %17, label %30, label %18
18: ; preds = %15, %18
%19 = phi i64 [ %27, %18 ], [ %16, %15 ]
%20 = phi i64 [ %28, %18 ], [ 0, %15 ]
%21 = getelementptr inbounds i32, i32* %1, i64 %19
%22 = load i32, i32* %21
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %2, i64 %23
%25 = load double, double* %24
%26 = getelementptr inbounds double, double* %3, i64 %19
store double %25, double* %26
%27 = add nuw nsw i64 %19, 1
%28 = add i64 %20, 1
%29 = icmp eq i64 %28, %11
br i1 %29, label %30, label %18
30: ; preds = %18, %15
br i1 %7, label %31, label %349
31: ; preds = %30
%32 = zext i32 %4 to i64
br label %66
33: ; preds = %33, %13
%34 = phi i64 [ 0, %13 ], [ %63, %33 ]
%35 = phi i64 [ 0, %13 ], [ %64, %33 ]
%36 = getelementptr inbounds i32, i32* %1, i64 %34
%37 = load i32, i32* %36
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double, double* %2, i64 %38
%40 = load double, double* %39
%41 = getelementptr inbounds double, double* %3, i64 %34
store double %40, double* %41
%42 = or i64 %34, 1
%43 = getelementptr inbounds i32, i32* %1, i64 %42
%44 = load i32, i32* %43
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double, double* %2, i64 %45
%47 = load double, double* %46
%48 = getelementptr inbounds double, double* %3, i64 %42
store double %47, double* %48
%49 = or i64 %34, 2
%50 = getelementptr inbounds i32, i32* %1, i64 %49
%51 = load i32, i32* %50
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds double, double* %2, i64 %52
%54 = load double, double* %53
%55 = getelementptr inbounds double, double* %3, i64 %49
store double %54, double* %55
%56 = or i64 %34, 3
%57 = getelementptr inbounds i32, i32* %1, i64 %56
%58 = load i32, i32* %57
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %2, i64 %59
%61 = load double, double* %60
%62 = getelementptr inbounds double, double* %3, i64 %56
store double %61, double* %62
%63 = add nuw nsw i64 %34, 4
%64 = add i64 %35, 4
%65 = icmp eq i64 %64, %14
br i1 %65, label %15, label %33
66: ; preds = %31, %120
%67 = phi i64 [ 0, %31 ], [ %122, %120 ]
%68 = icmp eq i64 %67, 0
br i1 %68, label %69, label %71
69: ; preds = %66
%70 = load double, double* %3
br label %120
71: ; preds = %66
%72 = add nuw i64 %67, 1
%73 = mul i64 %72, %67
%74 = lshr i64 %73, 1
%75 = getelementptr inbounds double, double* %3, i64 %67
%76 = and i64 %74, 2147483647
%77 = load double, double* %75
%78 = and i64 %67, 1
%79 = icmp eq i64 %67, 1
br i1 %79, label %107, label %80
80: ; preds = %71
%81 = and i64 %67, 9223372036854775806
br label %85
82: ; preds = %120
br i1 %7, label %83, label %349
83: ; preds = %82
%84 = zext i32 %4 to i64
br label %132
85: ; preds = %85, %80
%86 = phi double [ %77, %80 ], [ %103, %85 ]
%87 = phi i64 [ 0, %80 ], [ %104, %85 ]
%88 = phi i64 [ 0, %80 ], [ %105, %85 ]
%89 = getelementptr inbounds double, double* %3, i64 %87
%90 = load double, double* %89
%91 = add nuw nsw i64 %87, %76
%92 = getelementptr inbounds double, double* %0, i64 %91
%93 = load double, double* %92
%94 = fneg double %90
%95 = tail call double @llvm.fmuladd.f64(double %94, double %93, double %86)
store double %95, double* %75
%96 = or i64 %87, 1
%97 = getelementptr inbounds double, double* %3, i64 %96
%98 = load double, double* %97
%99 = add nuw nsw i64 %96, %76
%100 = getelementptr inbounds double, double* %0, i64 %99
%101 = load double, double* %100
%102 = fneg double %98
%103 = tail call double @llvm.fmuladd.f64(double %102, double %101, double %95)
store double %103, double* %75
%104 = add nuw nsw i64 %87, 2
%105 = add i64 %88, 2
%106 = icmp eq i64 %105, %81
br i1 %106, label %107, label %85
107: ; preds = %85, %71
%108 = phi double [ undef, %71 ], [ %103, %85 ]
%109 = phi double [ %77, %71 ], [ %103, %85 ]
%110 = phi i64 [ 0, %71 ], [ %104, %85 ]
%111 = icmp eq i64 %78, 0
br i1 %111, label %120, label %112
112: ; preds = %107
%113 = getelementptr inbounds double, double* %3, i64 %110
%114 = load double, double* %113
%115 = add nuw nsw i64 %110, %76
%116 = getelementptr inbounds double, double* %0, i64 %115
%117 = load double, double* %116
%118 = fneg double %114
%119 = tail call double @llvm.fmuladd.f64(double %118, double %117, double %109)
store double %119, double* %75
br label %120
120: ; preds = %112, %107, %69
%121 = phi double [ %70, %69 ], [ %108, %107 ], [ %119, %112 ]
%122 = add nuw nsw i64 %67, 1
%123 = mul i64 %122, %67
%124 = lshr i64 %123, 1
%125 = and i64 %124, 2147483647
%126 = add nuw nsw i64 %125, %67
%127 = getelementptr inbounds double, double* %0, i64 %126
%128 = load double, double* %127
%129 = getelementptr inbounds double, double* %3, i64 %67
%130 = fdiv double %121, %128
store double %130, double* %129
%131 = icmp eq i64 %122, %32
br i1 %131, label %82, label %66
132: ; preds = %83, %162
%133 = phi i64 [ %84, %83 ], [ %134, %162 ]
%134 = add nsw i64 %133, -1
%135 = icmp slt i64 %133, %84
%136 = getelementptr inbounds double, double* %3, i64 %134
%137 = load double, double* %136
br i1 %135, label %144, label %162
138: ; preds = %162
br i1 %7, label %139, label %349
139: ; preds = %138
%140 = and i64 %9, 3
%141 = icmp ult i64 %10, 3
br i1 %141, label %175, label %142
142: ; preds = %139
%143 = and i64 %9, 4294967292
br label %297
144: ; preds = %132, %144
%145 = phi double [ %160, %144 ], [ %137, %132 ]
%146 = phi i64 [ %149, %144 ], [ %133, %132 ]
%147 = getelementptr inbounds double, double* %3, i64 %146
%148 = load double, double* %147
%149 = add nsw i64 %146, 1
%150 = trunc i64 %149 to i32
%151 = and i64 %149, 4294967295
%152 = mul i64 %151, %146
%153 = trunc i64 %152 to i32
%154 = sdiv i32 %153, 2
%155 = sext i32 %154 to i64
%156 = add nsw i64 %134, %155
%157 = getelementptr inbounds double, double* %0, i64 %156
%158 = load double, double* %157
%159 = fneg double %148
%160 = tail call double @llvm.fmuladd.f64(double %159, double %158, double %145)
store double %160, double* %136
%161 = icmp eq i32 %150, %4
br i1 %161, label %162, label %144
162: ; preds = %144, %132
%163 = phi double [ %137, %132 ], [ %160, %144 ]
%164 = and i64 %134, 4294967295
%165 = mul i64 %164, %133
%166 = trunc i64 %165 to i32
%167 = sdiv i32 %166, 2
%168 = sext i32 %167 to i64
%169 = add nsw i64 %134, %168
%170 = getelementptr inbounds double, double* %0, i64 %169
%171 = load double, double* %170
%172 = getelementptr inbounds double, double* %3, i64 %134
%173 = fdiv double %163, %171
store double %173, double* %172
%174 = icmp sgt i64 %133, 1
br i1 %174, label %132, label %138
175: ; preds = %297, %139
%176 = phi i64 [ 0, %139 ], [ %327, %297 ]
%177 = icmp eq i64 %140, 0
br i1 %177, label %190, label %178
178: ; preds = %175, %178
%179 = phi i64 [ %187, %178 ], [ %176, %175 ]
%180 = phi i64 [ %188, %178 ], [ 0, %175 ]
%181 = getelementptr inbounds double, double* %3, i64 %179
%182 = load double, double* %181
%183 = getelementptr inbounds i32, i32* %1, i64 %179
%184 = load i32, i32* %183
%185 = sext i32 %184 to i64
%186 = getelementptr inbounds double, double* %5, i64 %185
store double %182, double* %186
%187 = add nuw nsw i64 %179, 1
%188 = add i64 %180, 1
%189 = icmp eq i64 %188, %140
br i1 %189, label %190, label %178
190: ; preds = %178, %175
br i1 %7, label %191, label %349
191: ; preds = %190
%192 = zext i32 %4 to i64
%193 = icmp ult i32 %4, 4
br i1 %193, label %279, label %194
194: ; preds = %191
%195 = getelementptr double, double* %3, i64 %9
%196 = getelementptr double, double* %5, i64 %9
%197 = icmp ugt double* %196, %3
%198 = icmp ugt double* %195, %5
%199 = and i1 %197, %198
br i1 %199, label %279, label %200
200: ; preds = %194
%201 = and i64 %9, 4294967292
%202 = add nsw i64 %201, -4
%203 = lshr exact i64 %202, 2
%204 = add nuw nsw i64 %203, 1
%205 = and i64 %204, 3
%206 = icmp ult i64 %202, 12
br i1 %206, label %258, label %207
207: ; preds = %200
%208 = and i64 %204, 9223372036854775804
br label %209
209: ; preds = %209, %207
%210 = phi i64 [ 0, %207 ], [ %255, %209 ]
%211 = phi i64 [ 0, %207 ], [ %256, %209 ]
%212 = getelementptr inbounds double, double* %5, i64 %210
%213 = bitcast double* %212 to <2 x double>*
%214 = load <2 x double>, <2 x double>* %213
%215 = getelementptr inbounds double, double* %212, i64 2
%216 = bitcast double* %215 to <2 x double>*
%217 = load <2 x double>, <2 x double>* %216
%218 = getelementptr inbounds double, double* %3, i64 %210
%219 = bitcast double* %218 to <2 x double>*
store <2 x double> %214, <2 x double>* %219
%220 = getelementptr inbounds double, double* %218, i64 2
%221 = bitcast double* %220 to <2 x double>*
store <2 x double> %217, <2 x double>* %221
%222 = or i64 %210, 4
%223 = getelementptr inbounds double, double* %5, i64 %222
%224 = bitcast double* %223 to <2 x double>*
%225 = load <2 x double>, <2 x double>* %224
%226 = getelementptr inbounds double, double* %223, i64 2
%227 = bitcast double* %226 to <2 x double>*
%228 = load <2 x double>, <2 x double>* %227
%229 = getelementptr inbounds double, double* %3, i64 %222
%230 = bitcast double* %229 to <2 x double>*
store <2 x double> %225, <2 x double>* %230
%231 = getelementptr inbounds double, double* %229, i64 2
%232 = bitcast double* %231 to <2 x double>*
store <2 x double> %228, <2 x double>* %232
%233 = or i64 %210, 8
%234 = getelementptr inbounds double, double* %5, i64 %233
%235 = bitcast double* %234 to <2 x double>*
%236 = load <2 x double>, <2 x double>* %235
%237 = getelementptr inbounds double, double* %234, i64 2
%238 = bitcast double* %237 to <2 x double>*
%239 = load <2 x double>, <2 x double>* %238
%240 = getelementptr inbounds double, double* %3, i64 %233
%241 = bitcast double* %240 to <2 x double>*
store <2 x double> %236, <2 x double>* %241
%242 = getelementptr inbounds double, double* %240, i64 2
%243 = bitcast double* %242 to <2 x double>*
store <2 x double> %239, <2 x double>* %243
%244 = or i64 %210, 12
%245 = getelementptr inbounds double, double* %5, i64 %244
%246 = bitcast double* %245 to <2 x double>*
%247 = load <2 x double>, <2 x double>* %246
%248 = getelementptr inbounds double, double* %245, i64 2
%249 = bitcast double* %248 to <2 x double>*
%250 = load <2 x double>, <2 x double>* %249
%251 = getelementptr inbounds double, double* %3, i64 %244
%252 = bitcast double* %251 to <2 x double>*
store <2 x double> %247, <2 x double>* %252
%253 = getelementptr inbounds double, double* %251, i64 2
%254 = bitcast double* %253 to <2 x double>*
store <2 x double> %250, <2 x double>* %254
%255 = add nuw i64 %210, 16
%256 = add nuw i64 %211, 4
%257 = icmp eq i64 %256, %208
br i1 %257, label %258, label %209
258: ; preds = %209, %200
%259 = phi i64 [ 0, %200 ], [ %255, %209 ]
%260 = icmp eq i64 %205, 0
br i1 %260, label %277, label %261
261: ; preds = %258, %261
%262 = phi i64 [ %274, %261 ], [ %259, %258 ]
%263 = phi i64 [ %275, %261 ], [ 0, %258 ]
%264 = getelementptr inbounds double, double* %5, i64 %262
%265 = bitcast double* %264 to <2 x double>*
%266 = load <2 x double>, <2 x double>* %265
%267 = getelementptr inbounds double, double* %264, i64 2
%268 = bitcast double* %267 to <2 x double>*
%269 = load <2 x double>, <2 x double>* %268
%270 = getelementptr inbounds double, double* %3, i64 %262
%271 = bitcast double* %270 to <2 x double>*
store <2 x double> %266, <2 x double>* %271
%272 = getelementptr inbounds double, double* %270, i64 2
%273 = bitcast double* %272 to <2 x double>*
store <2 x double> %269, <2 x double>* %273
%274 = add nuw i64 %262, 4
%275 = add i64 %263, 1
%276 = icmp eq i64 %275, %205
br i1 %276, label %277, label %261
277: ; preds = %261, %258
%278 = icmp eq i64 %201, %9
br i1 %278, label %349, label %279
279: ; preds = %194, %191, %277
%280 = phi i64 [ 0, %194 ], [ 0, %191 ], [ %201, %277 ]
%281 = xor i64 %280, -1
%282 = add nsw i64 %281, %9
%283 = and i64 %9, 3
%284 = icmp eq i64 %283, 0
br i1 %284, label %294, label %285
285: ; preds = %279, %285
%286 = phi i64 [ %291, %285 ], [ %280, %279 ]
%287 = phi i64 [ %292, %285 ], [ 0, %279 ]
%288 = getelementptr inbounds double, double* %5, i64 %286
%289 = load double, double* %288
%290 = getelementptr inbounds double, double* %3, i64 %286
store double %289, double* %290
%291 = add nuw nsw i64 %286, 1
%292 = add i64 %287, 1
%293 = icmp eq i64 %292, %283
br i1 %293, label %294, label %285
294: ; preds = %285, %279
%295 = phi i64 [ %280, %279 ], [ %291, %285 ]
%296 = icmp ult i64 %282, 3
br i1 %296, label %349, label %330
297: ; preds = %297, %142
%298 = phi i64 [ 0, %142 ], [ %327, %297 ]
%299 = phi i64 [ 0, %142 ], [ %328, %297 ]
%300 = getelementptr inbounds double, double* %3, i64 %298
%301 = load double, double* %300
%302 = getelementptr inbounds i32, i32* %1, i64 %298
%303 = load i32, i32* %302
%304 = sext i32 %303 to i64
%305 = getelementptr inbounds double, double* %5, i64 %304
store double %301, double* %305
%306 = or i64 %298, 1
%307 = getelementptr inbounds double, double* %3, i64 %306
%308 = load double, double* %307
%309 = getelementptr inbounds i32, i32* %1, i64 %306
%310 = load i32, i32* %309
%311 = sext i32 %310 to i64
%312 = getelementptr inbounds double, double* %5, i64 %311
store double %308, double* %312
%313 = or i64 %298, 2
%314 = getelementptr inbounds double, double* %3, i64 %313
%315 = load double, double* %314
%316 = getelementptr inbounds i32, i32* %1, i64 %313
%317 = load i32, i32* %316
%318 = sext i32 %317 to i64
%319 = getelementptr inbounds double, double* %5, i64 %318
store double %315, double* %319
%320 = or i64 %298, 3
%321 = getelementptr inbounds double, double* %3, i64 %320
%322 = load double, double* %321
%323 = getelementptr inbounds i32, i32* %1, i64 %320
%324 = load i32, i32* %323
%325 = sext i32 %324 to i64
%326 = getelementptr inbounds double, double* %5, i64 %325
store double %322, double* %326
%327 = add nuw nsw i64 %298, 4
%328 = add i64 %299, 4
%329 = icmp eq i64 %328, %143
br i1 %329, label %175, label %297
330: ; preds = %294, %330
%331 = phi i64 [ %347, %330 ], [ %295, %294 ]
%332 = getelementptr inbounds double, double* %5, i64 %331
%333 = load double, double* %332
%334 = getelementptr inbounds double, double* %3, i64 %331
store double %333, double* %334
%335 = add nuw nsw i64 %331, 1
%336 = getelementptr inbounds double, double* %5, i64 %335
%337 = load double, double* %336
%338 = getelementptr inbounds double, double* %3, i64 %335
store double %337, double* %338
%339 = add nuw nsw i64 %331, 2
%340 = getelementptr inbounds double, double* %5, i64 %339
%341 = load double, double* %340
%342 = getelementptr inbounds double, double* %3, i64 %339
store double %341, double* %342
%343 = add nuw nsw i64 %331, 3
%344 = getelementptr inbounds double, double* %5, i64 %343
%345 = load double, double* %344
%346 = getelementptr inbounds double, double* %3, i64 %343
store double %345, double* %346
%347 = add nuw nsw i64 %331, 4
%348 = icmp eq i64 %347, %192
br i1 %348, label %349, label %330
349: ; preds = %294, %330, %277, %6, %30, %82, %138, %190
ret void
}