#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ free(double *);
scalar_t__ mat(int, int);
int /*<<< orphan*/ memcpy(double *, double *, int);
double sqrt(double);
void cholesky_LL(double *P, double *L, int lda, int n, char lctl) {
  double *tmp = (double *)mat(n, n);
  int col, k, i;
  for (i = 0; i < n; i++)
    memcpy(tmp + i * n, P + lda * i, sizeof(double) * n);
  int iptx[n];
  for (i = 0; i < n; i++) {
    iptx[i] = (i + 1) * i / 2;
  }
  for (col = 0; col < n; col++) {
    if (lctl == 'H') {
      L[iptx[col] + col] = sqrt(tmp[col + col * n]);
    } else {
      L[col + col * n] = sqrt(tmp[col + col * n]);
    }
    for (k = col + 1; k < n; k++) {
      if (lctl == 'H') {
        L[iptx[k] + col] = tmp[k + col * n] / L[iptx[col] + col];
      } else {
        L[k + col * n] = tmp[k + col * n] / L[col + col * n];
      }
    }
    for (k = col + 1; k < n; k++) {
      for (i = k; i < n; i++) {
        if (lctl == 'H') {
          tmp[i + k * n] = tmp[i + k * n] - L[iptx[i] + col] * L[iptx[k] + col];
        } else {
          tmp[i + k * n] = tmp[i + k * n] - L[i + col * n] * L[k + col * n];
        }
      }
    }
  }
  free(tmp);
}


x86 assembly (clang, -O0):
<cholesky_LL>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %r8b,%al
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %al,-0x19(%rbp)
mov    -0x18(%rbp),%edi
mov    -0x18(%rbp),%esi
call   1040 <mat@plt>
mov    %rax,-0x28(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11c3 <cholesky_LL+0x83>
mov    -0x28(%rbp),%rdi
mov    -0x34(%rbp),%eax
imul   -0x18(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x8(%rbp),%rsi
mov    -0x14(%rbp),%eax
imul   -0x34(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
movslq -0x18(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edx
call   1050 <memcpy@plt>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1172 <cholesky_LL+0x32>
mov    -0x18(%rbp),%eax
mov    %rsp,%rcx
mov    %rcx,-0x40(%rbp)
lea    0xf(,%rax,4),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x50(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    122a <cholesky_LL+0xea>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
imul   -0x34(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%edx
mov    -0x50(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11f1 <cholesky_LL+0xb1>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    14f5 <cholesky_LL+0x3b5>
movsbl -0x19(%rbp),%eax
cmp    $0x48,%eax
jne    1286 <cholesky_LL+0x146>
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1060 <sqrt@plt>
mov    -0x50(%rbp),%rcx
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    12bb <cholesky_LL+0x17b>
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1060 <sqrt@plt>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1385 <cholesky_LL+0x245>
movsbl -0x19(%rbp),%eax
cmp    $0x48,%eax
jne    132a <cholesky_LL+0x1ea>
mov    -0x50(%rbp),%rcx
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%edx
mov    -0x2c(%rbp),%esi
imul   -0x18(%rbp),%esi
add    %esi,%edx
movslq %edx,%rdx
movsd  (%rax,%rdx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rdx
mov    (%rcx,%rdx,4),%edx
add    -0x2c(%rbp),%edx
movslq %edx,%rdx
divsd  (%rax,%rdx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1372 <cholesky_LL+0x232>
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1377 <cholesky_LL+0x237>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12c4 <cholesky_LL+0x184>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    14e2 <cholesky_LL+0x3a2>
mov    -0x30(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    14cf <cholesky_LL+0x38f>
movsbl -0x19(%rbp),%eax
cmp    $0x48,%eax
jne    143d <cholesky_LL+0x2fd>
mov    -0x50(%rbp),%rcx
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%edx
mov    -0x30(%rbp),%esi
imul   -0x18(%rbp),%esi
add    %esi,%edx
movslq %edx,%rdx
movsd  (%rax,%rdx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rdx
mov    (%rcx,%rdx,4),%edx
add    -0x2c(%rbp),%edx
movslq %edx,%rdx
movsd  (%rax,%rdx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    14bc <cholesky_LL+0x37c>
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    14c1 <cholesky_LL+0x381>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    13a0 <cholesky_LL+0x260>
jmp    14d4 <cholesky_LL+0x394>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    138e <cholesky_LL+0x24e>
jmp    14e7 <cholesky_LL+0x3a7>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1231 <cholesky_LL+0xf1>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
mov    -0x40(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cholesky_LL>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r8d,%r15d
mov    %ecx,%ebx
mov    %edx,%r12d
mov    %rsi,%r14
mov    %rdi,%r13
mov    %ecx,%edi
mov    %ecx,%esi
call   1040 <mat@plt>
mov    %rax,-0x50(%rbp)
mov    %rbx,-0x40(%rbp)
test   %ebx,%ebx
jle    151c <cholesky_LL+0x3dc>
mov    -0x40(%rbp),%rax
lea    0x0(,%rax,8),%ecx
mov    %ecx,-0x38(%rbp)
movslq %r12d,%rcx
mov    %eax,%r12d
shl    $0x3,%rcx
mov    %rcx,-0x58(%rbp)
lea    0x0(,%r12,8),%rax
mov    %rax,-0x48(%rbp)
mov    %r12,%rax
mov    -0x50(%rbp),%rbx
nopw   0x0(%rax,%rax,1)
mov    %rax,-0x30(%rbp)
mov    %rbx,%rdi
mov    %r13,%rsi
mov    -0x38(%rbp),%edx
call   1050 <memcpy@plt>
mov    -0x30(%rbp),%rax
add    -0x58(%rbp),%r13
add    -0x48(%rbp),%rbx
add    $0xffffffffffffffff,%rax
jne    11b0 <cholesky_LL+0x70>
mov    %rsp,%rcx
lea    0xf(,%r12,4),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%rcx
mov    %rcx,-0x30(%rbp)
mov    %rcx,%rsp
cmpl   $0x0,-0x40(%rbp)
jle    151c <cholesky_LL+0x3dc>
cmpl   $0x8,-0x40(%rbp)
jae    1208 <cholesky_LL+0xc8>
xor    %eax,%eax
mov    -0x30(%rbp),%rsi
jmp    1300 <cholesky_LL+0x1c0>
mov    %r12d,%eax
and    $0xfffffff8,%eax
lea    0x0(,%r12,4),%rcx
and    $0xffffffffffffffe0,%rcx
movaps 0xddf(%rip),%xmm0
movaps 0xde8(%rip),%xmm4
movdqa 0xdf0(%rip),%xmm3
xor    %edx,%edx
movdqa 0xdf5(%rip),%xmm8
pcmpeqd %xmm9,%xmm9
movdqa 0xdf7(%rip),%xmm10
movdqa 0xdfe(%rip),%xmm12
movdqa 0xe05(%rip),%xmm11
mov    -0x30(%rbp),%rsi
nop
movdqa %xmm3,%xmm2
paddd  %xmm8,%xmm2
movaps %xmm4,%xmm1
shufps $0x88,%xmm0,%xmm1
movaps %xmm1,%xmm5
psubd  %xmm9,%xmm5
paddd  %xmm10,%xmm1
pshufd $0xf5,%xmm5,%xmm7
pmuludq %xmm3,%xmm5
pshufd $0xe8,%xmm5,%xmm5
pshufd $0xf5,%xmm3,%xmm6
pmuludq %xmm7,%xmm6
pshufd $0xe8,%xmm6,%xmm6
punpckldq %xmm6,%xmm5
pshufd $0xf5,%xmm1,%xmm6
pmuludq %xmm2,%xmm1
pshufd $0xe8,%xmm1,%xmm1
pshufd $0xf5,%xmm2,%xmm2
pmuludq %xmm6,%xmm2
pshufd $0xe8,%xmm2,%xmm2
punpckldq %xmm2,%xmm1
psrld  $0x1,%xmm5
psrld  $0x1,%xmm1
movdqa %xmm5,(%rsi,%rdx,1)
movdqa %xmm1,0x10(%rsi,%rdx,1)
paddq  %xmm12,%xmm4
paddq  %xmm12,%xmm0
paddd  %xmm11,%xmm3
add    $0x20,%rdx
cmp    %rdx,%rcx
jne    1260 <cholesky_LL+0x120>
cmp    %r12,%rax
je     1316 <cholesky_LL+0x1d6>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
lea    0x1(%rax),%rcx
mov    %ecx,%edx
imul   %eax,%edx
shr    %edx
mov    %edx,(%rsi,%rax,4)
mov    %rcx,%rax
cmp    %rcx,%r12
jne    1300 <cholesky_LL+0x1c0>
cmpl   $0x0,-0x40(%rbp)
jle    151c <cholesky_LL+0x3dc>
lea    0x1(%r12),%rax
mov    %rax,-0x80(%rbp)
lea    0x8(,%r12,8),%rax
mov    %rax,-0x58(%rbp)
mov    -0x30(%rbp),%rax
lea    0x4(%rax),%rcx
lea    -0x1(%r12),%rdx
mov    -0x50(%rbp),%rax
lea    (%rax,%r12,8),%rsi
add    $0x8,%rsi
mov    %rsi,-0x68(%rbp)
mov    %rax,%rsi
add    $0x8,%rsi
mov    $0x1,%eax
mov    %rax,-0x38(%rbp)
xor    %r8d,%r8d
pxor   %xmm2,%xmm2
mov    $0x1,%eax
mov    %rax,-0x60(%rbp)
jmp    13b9 <cholesky_LL+0x279>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
addq   $0x1,-0x60(%rbp)
mov    -0x38(%rbp),%rax
add    -0x80(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x58(%rbp),%rax
mov    -0x70(%rbp),%rsi
add    %rax,%rsi
mov    -0x48(%rbp),%rcx
add    $0x4,%rcx
mov    -0x78(%rbp),%rdx
add    $0xffffffffffffffff,%rdx
add    %rax,-0x68(%rbp)
cmp    %r12,%r8
je     151c <cholesky_LL+0x3dc>
mov    %r8,%r13
mov    %r13d,%eax
imul   -0x40(%rbp),%eax
movslq %eax,%rbx
add    %r8,%rbx
mov    -0x50(%rbp),%rax
movsd  (%rax,%rbx,8),%xmm0
ucomisd %xmm2,%xmm0
mov    %rcx,-0x48(%rbp)
mov    %rdx,-0x78(%rbp)
mov    %rsi,-0x70(%rbp)
jb     13f0 <cholesky_LL+0x2b0>
sqrtsd %xmm0,%xmm0
cmp    $0x48,%r15b
je     13ff <cholesky_LL+0x2bf>
jmp    140a <cholesky_LL+0x2ca>
call   1060 <sqrt@plt>
xorpd  %xmm2,%xmm2
cmp    $0x48,%r15b
jne    140a <cholesky_LL+0x2ca>
mov    -0x30(%rbp),%rax
movslq (%rax,%r13,4),%rbx
add    %r13,%rbx
movsd  %xmm0,(%r14,%rbx,8)
lea    0x1(%r13),%r8
cmp    %r12,%r8
jae    1380 <cholesky_LL+0x240>
mov    %r13,%r9
imul   %r12,%r9
lea    (%r9,%r13,1),%r10
xor    %ecx,%ecx
mov    -0x30(%rbp),%r11
mov    -0x48(%rbp),%rbx
mov    -0x78(%rbp),%rax
mov    -0x70(%rbp),%rdi
jmp    1468 <cholesky_LL+0x328>
nopl   0x0(%rax)
movslq (%r11,%r13,4),%rsi
add    %r13,%rsi
movslq (%rbx,%rcx,4),%rdx
add    %r13,%rdx
movsd  (%rdi,%rcx,8),%xmm0
divsd  (%r14,%rsi,8),%xmm0
movsd  %xmm0,(%r14,%rdx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
je     1480 <cholesky_LL+0x340>
cmp    $0x48,%r15b
je     1440 <cholesky_LL+0x300>
mov    -0x38(%rbp),%rdx
add    %rcx,%rdx
mov    %r10,%rsi
jmp    144e <cholesky_LL+0x30e>
nopw   0x0(%rax,%rax,1)
cmp    %r12,%r8
jae    1380 <cholesky_LL+0x240>
mov    -0x48(%rbp),%r10
mov    -0x68(%rbp),%rsi
mov    -0x38(%rbp),%r11
mov    -0x60(%rbp),%rax
jmp    14b9 <cholesky_LL+0x379>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rax
add    $0x1,%r11
add    -0x58(%rbp),%rsi
add    $0x4,%r10
cmp    %r12,%rax
je     1380 <cholesky_LL+0x240>
lea    (%rax,%r9,1),%rdi
xor    %ebx,%ebx
jmp    1502 <cholesky_LL+0x3c2>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    (%r11,%rbx,1),%rdx
mov    %rdi,%rcx
movsd  (%rsi,%rbx,8),%xmm0
movsd  (%r14,%rdx,8),%xmm1
mulsd  (%r14,%rcx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rsi,%rbx,8)
lea    (%rax,%rbx,1),%rcx
add    $0x1,%rcx
add    $0x1,%rbx
cmp    %r12,%rcx
jae    14a0 <cholesky_LL+0x360>
cmp    $0x48,%r15b
jne    14d0 <cholesky_LL+0x390>
movslq (%r10,%rbx,4),%rdx
add    %r13,%rdx
mov    -0x30(%rbp),%rcx
movslq (%rcx,%rax,4),%rcx
add    %r13,%rcx
jmp    14d7 <cholesky_LL+0x397>
mov    -0x50(%rbp),%rdi
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @cholesky_LL(double* %0, double* %1, i32 %2, i32 %3, i8 signext %4) {
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i8
%11 = alloca double*
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i8*
%16 = alloca i64
store double* %0, double** %6
store double* %1, double** %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store i8 %4, i8* %10
%17 = load i32, i32* %9
%18 = load i32, i32* %9
%19 = call i64 @mat(i32 %17, i32 %18)
%20 = inttoptr i64 %19 to double*
store double* %20, double** %11
store i32 0, i32* %14
br label %21
21: ; preds = %43, %5
%22 = load i32, i32* %14
%23 = load i32, i32* %9
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %46
25: ; preds = %21
%26 = load double*, double** %11
%27 = load i32, i32* %14
%28 = load i32, i32* %9
%29 = mul nsw i32 %27, %28
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %26, i64 %30
%32 = load double*, double** %6
%33 = load i32, i32* %8
%34 = load i32, i32* %14
%35 = mul nsw i32 %33, %34
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %32, i64 %36
%38 = load i32, i32* %9
%39 = sext i32 %38 to i64
%40 = mul i64 8, %39
%41 = trunc i64 %40 to i32
%42 = call i32 @memcpy(double* %31, double* %37, i32 %41)
br label %43
43: ; preds = %25
%44 = load i32, i32* %14
%45 = add nsw i32 %44, 1
store i32 %45, i32* %14
br label %21
46: ; preds = %21
%47 = load i32, i32* %9
%48 = zext i32 %47 to i64
%49 = call i8* @llvm.stacksave()
store i8* %49, i8** %15
%50 = alloca i32, i64 %48
store i64 %48, i64* %16
store i32 0, i32* %14
br label %51
51: ; preds = %64, %46
%52 = load i32, i32* %14
%53 = load i32, i32* %9
%54 = icmp slt i32 %52, %53
br i1 %54, label %55, label %67
55: ; preds = %51
%56 = load i32, i32* %14
%57 = add nsw i32 %56, 1
%58 = load i32, i32* %14
%59 = mul nsw i32 %57, %58
%60 = sdiv i32 %59, 2
%61 = load i32, i32* %14
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds i32, i32* %50, i64 %62
store i32 %60, i32* %63
br label %64
64: ; preds = %55
%65 = load i32, i32* %14
%66 = add nsw i32 %65, 1
store i32 %66, i32* %14
br label %51
67: ; preds = %51
store i32 0, i32* %12
br label %68
68: ; preds = %292, %67
%69 = load i32, i32* %12
%70 = load i32, i32* %9
%71 = icmp slt i32 %69, %70
br i1 %71, label %72, label %295
72: ; preds = %68
%73 = load i8, i8* %10
%74 = sext i8 %73 to i32
%75 = icmp eq i32 %74, 72
br i1 %75, label %76, label %96
76: ; preds = %72
%77 = load double*, double** %11
%78 = load i32, i32* %12
%79 = load i32, i32* %12
%80 = load i32, i32* %9
%81 = mul nsw i32 %79, %80
%82 = add nsw i32 %78, %81
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds double, double* %77, i64 %83
%85 = load double, double* %84
%86 = call double @sqrt(double %85)
%87 = load double*, double** %7
%88 = load i32, i32* %12
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds i32, i32* %50, i64 %89
%91 = load i32, i32* %90
%92 = load i32, i32* %12
%93 = add nsw i32 %91, %92
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds double, double* %87, i64 %94
store double %86, double* %95
br label %115
96: ; preds = %72
%97 = load double*, double** %11
%98 = load i32, i32* %12
%99 = load i32, i32* %12
%100 = load i32, i32* %9
%101 = mul nsw i32 %99, %100
%102 = add nsw i32 %98, %101
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds double, double* %97, i64 %103
%105 = load double, double* %104
%106 = call double @sqrt(double %105)
%107 = load double*, double** %7
%108 = load i32, i32* %12
%109 = load i32, i32* %12
%110 = load i32, i32* %9
%111 = mul nsw i32 %109, %110
%112 = add nsw i32 %108, %111
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %107, i64 %113
store double %106, double* %114
br label %115
115: ; preds = %96, %76
%116 = load i32, i32* %12
%117 = add nsw i32 %116, 1
store i32 %117, i32* %13
br label %118
118: ; preds = %185, %115
%119 = load i32, i32* %13
%120 = load i32, i32* %9
%121 = icmp slt i32 %119, %120
br i1 %121, label %122, label %188
122: ; preds = %118
%123 = load i8, i8* %10
%124 = sext i8 %123 to i32
%125 = icmp eq i32 %124, 72
br i1 %125, label %126, label %156
126: ; preds = %122
%127 = load double*, double** %11
%128 = load i32, i32* %13
%129 = load i32, i32* %12
%130 = load i32, i32* %9
%131 = mul nsw i32 %129, %130
%132 = add nsw i32 %128, %131
%133 = sext i32 %132 to i64
%134 = getelementptr inbounds double, double* %127, i64 %133
%135 = load double, double* %134
%136 = load double*, double** %7
%137 = load i32, i32* %12
%138 = sext i32 %137 to i64
%139 = getelementptr inbounds i32, i32* %50, i64 %138
%140 = load i32, i32* %139
%141 = load i32, i32* %12
%142 = add nsw i32 %140, %141
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds double, double* %136, i64 %143
%145 = load double, double* %144
%146 = fdiv double %135, %145
%147 = load double*, double** %7
%148 = load i32, i32* %13
%149 = sext i32 %148 to i64
%150 = getelementptr inbounds i32, i32* %50, i64 %149
%151 = load i32, i32* %150
%152 = load i32, i32* %12
%153 = add nsw i32 %151, %152
%154 = sext i32 %153 to i64
%155 = getelementptr inbounds double, double* %147, i64 %154
store double %146, double* %155
br label %184
156: ; preds = %122
%157 = load double*, double** %11
%158 = load i32, i32* %13
%159 = load i32, i32* %12
%160 = load i32, i32* %9
%161 = mul nsw i32 %159, %160
%162 = add nsw i32 %158, %161
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds double, double* %157, i64 %163
%165 = load double, double* %164
%166 = load double*, double** %7
%167 = load i32, i32* %12
%168 = load i32, i32* %12
%169 = load i32, i32* %9
%170 = mul nsw i32 %168, %169
%171 = add nsw i32 %167, %170
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds double, double* %166, i64 %172
%174 = load double, double* %173
%175 = fdiv double %165, %174
%176 = load double*, double** %7
%177 = load i32, i32* %13
%178 = load i32, i32* %12
%179 = load i32, i32* %9
%180 = mul nsw i32 %178, %179
%181 = add nsw i32 %177, %180
%182 = sext i32 %181 to i64
%183 = getelementptr inbounds double, double* %176, i64 %182
store double %175, double* %183
br label %184
184: ; preds = %156, %126
br label %185
185: ; preds = %184
%186 = load i32, i32* %13
%187 = add nsw i32 %186, 1
store i32 %187, i32* %13
br label %118
188: ; preds = %118
%189 = load i32, i32* %12
%190 = add nsw i32 %189, 1
store i32 %190, i32* %13
br label %191
191: ; preds = %288, %188
%192 = load i32, i32* %13
%193 = load i32, i32* %9
%194 = icmp slt i32 %192, %193
br i1 %194, label %195, label %291
195: ; preds = %191
%196 = load i32, i32* %13
store i32 %196, i32* %14
br label %197
197: ; preds = %284, %195
%198 = load i32, i32* %14
%199 = load i32, i32* %9
%200 = icmp slt i32 %198, %199
br i1 %200, label %201, label %287
201: ; preds = %197
%202 = load i8, i8* %10
%203 = sext i8 %202 to i32
%204 = icmp eq i32 %203, 72
br i1 %204, label %205, label %245
205: ; preds = %201
%206 = load double*, double** %11
%207 = load i32, i32* %14
%208 = load i32, i32* %13
%209 = load i32, i32* %9
%210 = mul nsw i32 %208, %209
%211 = add nsw i32 %207, %210
%212 = sext i32 %211 to i64
%213 = getelementptr inbounds double, double* %206, i64 %212
%214 = load double, double* %213
%215 = load double*, double** %7
%216 = load i32, i32* %14
%217 = sext i32 %216 to i64
%218 = getelementptr inbounds i32, i32* %50, i64 %217
%219 = load i32, i32* %218
%220 = load i32, i32* %12
%221 = add nsw i32 %219, %220
%222 = sext i32 %221 to i64
%223 = getelementptr inbounds double, double* %215, i64 %222
%224 = load double, double* %223
%225 = load double*, double** %7
%226 = load i32, i32* %13
%227 = sext i32 %226 to i64
%228 = getelementptr inbounds i32, i32* %50, i64 %227
%229 = load i32, i32* %228
%230 = load i32, i32* %12
%231 = add nsw i32 %229, %230
%232 = sext i32 %231 to i64
%233 = getelementptr inbounds double, double* %225, i64 %232
%234 = load double, double* %233
%235 = fneg double %224
%236 = call double @llvm.fmuladd.f64(double %235, double %234, double %214)
%237 = load double*, double** %11
%238 = load i32, i32* %14
%239 = load i32, i32* %13
%240 = load i32, i32* %9
%241 = mul nsw i32 %239, %240
%242 = add nsw i32 %238, %241
%243 = sext i32 %242 to i64
%244 = getelementptr inbounds double, double* %237, i64 %243
store double %236, double* %244
br label %283
245: ; preds = %201
%246 = load double*, double** %11
%247 = load i32, i32* %14
%248 = load i32, i32* %13
%249 = load i32, i32* %9
%250 = mul nsw i32 %248, %249
%251 = add nsw i32 %247, %250
%252 = sext i32 %251 to i64
%253 = getelementptr inbounds double, double* %246, i64 %252
%254 = load double, double* %253
%255 = load double*, double** %7
%256 = load i32, i32* %14
%257 = load i32, i32* %12
%258 = load i32, i32* %9
%259 = mul nsw i32 %257, %258
%260 = add nsw i32 %256, %259
%261 = sext i32 %260 to i64
%262 = getelementptr inbounds double, double* %255, i64 %261
%263 = load double, double* %262
%264 = load double*, double** %7
%265 = load i32, i32* %13
%266 = load i32, i32* %12
%267 = load i32, i32* %9
%268 = mul nsw i32 %266, %267
%269 = add nsw i32 %265, %268
%270 = sext i32 %269 to i64
%271 = getelementptr inbounds double, double* %264, i64 %270
%272 = load double, double* %271
%273 = fneg double %263
%274 = call double @llvm.fmuladd.f64(double %273, double %272, double %254)
%275 = load double*, double** %11
%276 = load i32, i32* %14
%277 = load i32, i32* %13
%278 = load i32, i32* %9
%279 = mul nsw i32 %277, %278
%280 = add nsw i32 %276, %279
%281 = sext i32 %280 to i64
%282 = getelementptr inbounds double, double* %275, i64 %281
store double %274, double* %282
br label %283
283: ; preds = %245, %205
br label %284
284: ; preds = %283
%285 = load i32, i32* %14
%286 = add nsw i32 %285, 1
store i32 %286, i32* %14
br label %197
287: ; preds = %197
br label %288
288: ; preds = %287
%289 = load i32, i32* %13
%290 = add nsw i32 %289, 1
store i32 %290, i32* %13
br label %191
291: ; preds = %191
br label %292
292: ; preds = %291
%293 = load i32, i32* %12
%294 = add nsw i32 %293, 1
store i32 %294, i32* %12
br label %68
295: ; preds = %68
%296 = load double*, double** %11
%297 = call i32 @free(double* %296)
%298 = load i8*, i8** %15
ret void
}

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