double *cholesky(double *A, int n) {
  register int i, j, k;
  double *L = (double *)calloc(n * n, sizeof(double));
  for (i = 0; i < n; ++i) {
    for (j = 0; j < i; ++j) {
      L[((((i) * (i + 1)) / 2) + (j))] = A[((i * n) + j)];
      for (k = 0; k < j; ++k) {
        L[((((i) * (i + 1)) / 2) + (j))] -=
            L[((((i) * (i + 1)) / 2) + (k))] * L[((((j) * (j + 1)) / 2) + (k))];
      }
      L[((((i) * (i + 1)) / 2) + (j))] /= L[((((j) * (j + 1)) / 2) + (j))];
    }
    L[((((j) * (j + 1)) / 2) + (j))] = A[((j * n) + j)];
    for (k = 0; k < i; ++k) {
      L[((((i) * (i + 1)) / 2) + (i))] -=
          L[((((i) * (i + 1)) / 2) + (k))] * L[((((i) * (i + 1)) / 2) + (k))];
    }
    L[((((i) * (i + 1)) / 2) + (i))] = sqrt(L[((((i) * (i + 1)) / 2) + (i))]);
  }
  return L;
}


x86 assembly (clang, -O0):
<cholesky>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0xc(%rbp),%eax
imul   -0xc(%rbp),%eax
movslq %eax,%rdi
mov    $0x8,%esi
call   1030 <calloc@plt>
mov    %rax,-0x20(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    146c <cholesky+0x34c>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    12f0 <cholesky+0x1d0>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x28(%rbp),%rax
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    127f <cholesky+0x15f>
mov    -0x20(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x40(%rbp),%rax
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x38(%rbp),%rax
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x30(%rbp),%rax
add    -0x14(%rbp),%ecx
movslq %ecx,%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    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    11ba <cholesky+0x9a>
mov    -0x20(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x50(%rbp),%rax
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x48(%rbp),%rax
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1164 <cholesky+0x44>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x58(%rbp),%rax
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    13ff <cholesky+0x2df>
mov    -0x20(%rbp),%rax
mov    %rax,-0x70(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x70(%rbp),%rax
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    %rax,-0x68(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x68(%rbp),%rax
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x60(%rbp),%rax
add    -0x10(%rbp),%ecx
movslq %ecx,%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    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    133a <cholesky+0x21a>
mov    -0x20(%rbp),%rax
mov    %rax,-0x80(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x80(%rbp),%rax
add    -0x10(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1040 <sqrt@plt>
mov    -0x20(%rbp),%rax
mov    %rax,-0x78(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x78(%rbp),%rax
add    -0x10(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1151 <cholesky+0x31>
mov    -0x20(%rbp),%rax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cholesky>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %esi,%ebx
mov    %rdi,%r15
mov    %esi,%edi
imul   %edi,%edi
mov    $0x8,%esi
call   1030 <calloc@plt>
mov    %rax,%rbp
mov    %ebx,0x8(%rsp)
test   %ebx,%ebx
jle    13b6 <cholesky+0x296>
mov    0x8(%rsp),%eax
mov    %rax,0x10(%rsp)
mov    %rbp,%r12
add    $0x8,%r12
mov    $0x2,%eax
xor    %r14d,%r14d
xorpd  %xmm2,%xmm2
xor    %r8d,%r8d
jmp    11a7 <cholesky+0x87>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
call   1040 <sqrt@plt>
xorpd  %xmm2,%xmm2
mov    0xc(%rsp),%eax
movsd  %xmm0,0x0(%rbp,%rbx,8)
add    %eax,%r14d
add    $0x2,%eax
mov    %r13,%r8
cmp    0x10(%rsp),%r13
je     13b6 <cholesky+0x296>
mov    %eax,0xc(%rsp)
mov    %r14d,%ecx
and    $0xfffffffe,%ecx
shl    $0x2,%rcx
add    %r12,%rcx
test   %r8,%r8
je     12e0 <cholesky+0x1c0>
mov    %r8,%r9
imul   0x10(%rsp),%r9
lea    0x1(%r8),%r10d
imul   %r8d,%r10d
shr    %r10d
mov    $0x2,%edi
xor    %r11d,%r11d
xor    %r13d,%r13d
jmp    1221 <cholesky+0x101>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %r13,%rax
add    $0x1,%rax
mov    %eax,%esi
imul   %r13d,%esi
shr    %esi
add    %r13d,%esi
divsd  0x0(%rbp,%rsi,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
add    %rdi,%r11
add    $0x2,%rdi
mov    %rax,%r13
cmp    %r8,%rax
je     12e0 <cholesky+0x1c0>
lea    (%r9,%r13,1),%rax
movsd  (%r15,%rax,8),%xmm0
lea    (%r10,%r13,1),%rdx
movsd  %xmm0,0x0(%rbp,%rdx,8)
test   %r13,%r13
je     11f0 <cholesky+0xd0>
cmp    $0x1,%r13
jne    1250 <cholesky+0x130>
xor    %esi,%esi
jmp    12a3 <cholesky+0x183>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %r11d,%eax
and    $0xfffffffe,%eax
lea    (%r12,%rax,4),%rbx
mov    %r13,%rax
and    $0xfffffffffffffffe,%rax
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  -0x8(%rcx,%rsi,8),%xmm1
mulsd  -0x8(%rbx,%rsi,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
movsd  (%rcx,%rsi,8),%xmm1
mulsd  (%rbx,%rsi,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
add    $0x2,%rsi
cmp    %rsi,%rax
jne    1270 <cholesky+0x150>
test   $0x1,%r13b
je     11f0 <cholesky+0xd0>
lea    0x1(%r13),%eax
imul   %r13d,%eax
shr    %eax
lea    (%rsi,%r10,1),%rbx
movsd  0x0(%rbp,%rbx,8),%xmm1
add    %rsi,%rax
mulsd  0x0(%rbp,%rax,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
jmp    11f0 <cholesky+0xd0>
nopl   0x0(%rax)
mov    %r8d,%eax
imul   0x8(%rsp),%eax
add    %r8d,%eax
lea    0x1(%r8),%edx
imul   %r8d,%edx
shr    %edx
add    %r8d,%edx
cltq
movsd  (%r15,%rax,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
test   %r8,%r8
je     138d <cholesky+0x26d>
lea    0x1(%r8),%eax
imul   %r8d,%eax
shr    %eax
lea    (%rax,%r8,1),%edx
movsd  0x0(%rbp,%rdx,8),%xmm0
cmp    $0x1,%r8
jne    1330 <cholesky+0x210>
xor    %esi,%esi
jmp    1370 <cholesky+0x250>
nopl   (%rax)
mov    %r8,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %esi,%esi
nopl   0x0(%rax)
movsd  -0x8(%rcx,%rsi,8),%xmm1
mulsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
movsd  (%rcx,%rsi,8),%xmm1
mulsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
add    $0x2,%rsi
cmp    %rsi,%rdi
jne    1340 <cholesky+0x220>
test   $0x1,%r8b
je     138d <cholesky+0x26d>
add    %rax,%rsi
movsd  0x0(%rbp,%rsi,8),%xmm1
mulsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
lea    0x1(%r8),%r13
mov    %r8d,%ebx
imul   %r13d,%ebx
shr    %ebx
add    %r8d,%ebx
movsd  0x0(%rbp,%rbx,8),%xmm0
ucomisd %xmm2,%xmm0
jb     1180 <cholesky+0x60>
sqrtsd %xmm0,%xmm0
jmp    1189 <cholesky+0x69>
mov    %rbp,%rax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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