#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 */
double sqrt(double);
int choleskyDecomposition(double **A, int n) {
  int i, j, k;
  for (k = 0; k < n; k++) {
    if (A[k][k] <= 0.0) {
      return 0;
    }
    A[k][k] = sqrt(A[k][k]);
    for (i = k + 1; i < n; i++) {
      A[i][k] /= A[k][k];
    }
    for (j = k + 1; j < n; j++) {
      for (i = j; i < n; i++) {
        A[i][j] -= A[i][k] * A[j][k];
      }
    }
  }
  return 1;
}


x86 assembly (clang, -O0):
<choleskyDecomposition>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    129f <choleskyDecomposition+0x18f>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jb     115c <choleskyDecomposition+0x4c>
movl   $0x0,-0x4(%rbp)
jmp    12a6 <choleskyDecomposition+0x196>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <sqrt@plt>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11e1 <choleskyDecomposition+0xd1>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1194 <choleskyDecomposition+0x84>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    128c <choleskyDecomposition+0x17c>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1279 <choleskyDecomposition+0x169>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1c(%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    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    11fc <choleskyDecomposition+0xec>
jmp    127e <choleskyDecomposition+0x16e>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11ea <choleskyDecomposition+0xda>
jmp    1291 <choleskyDecomposition+0x181>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1126 <choleskyDecomposition+0x16>
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<choleskyDecomposition>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    $0x1,%eax
test   %esi,%esi
jle    1281 <choleskyDecomposition+0x171>
mov    %esi,%ebp
mov    %rdi,%rbx
mov    %esi,%r14d
mov    $0x1,%r8d
xor    %eax,%eax
xorpd  %xmm2,%xmm2
mov    %rdi,%r13
mov    %r14,%r15
jmp    1165 <choleskyDecomposition+0x55>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r8
add    $0xffffffffffffffff,%r15
add    $0x8,%r13
cmp    %r14,%rax
je     1278 <choleskyDecomposition+0x168>
mov    %rax,%r12
mov    (%rbx,%rax,8),%rax
movsd  (%rax,%r12,8),%xmm0
ucomisd %xmm0,%xmm2
jae    127f <choleskyDecomposition+0x16f>
ucomisd %xmm2,%xmm0
jb     1190 <choleskyDecomposition+0x80>
sqrtsd %xmm0,%xmm0
jmp    11a1 <choleskyDecomposition+0x91>
nopl   0x0(%rax,%rax,1)
mov    %r8,(%rsp)
call   1030 <sqrt@plt>
xorpd  %xmm2,%xmm2
mov    (%rsp),%r8
mov    (%rbx,%r12,8),%rcx
movsd  %xmm0,(%rcx,%r12,8)
lea    0x1(%r12),%rax
cmp    %r14,%rax
jae    1150 <choleskyDecomposition+0x40>
mov    (%rbx,%r8,8),%rdx
movsd  (%rdx,%r12,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%r12,8)
lea    0x1(%r8),%rdx
cmp    %r14,%rdx
je     1200 <choleskyDecomposition+0xf0>
mov    $0x2,%edx
nopw   0x0(%rax,%rax,1)
mov    0x0(%r13,%rdx,8),%rsi
movsd  (%rsi,%r12,8),%xmm0
divsd  (%rcx,%r12,8),%xmm0
movsd  %xmm0,(%rsi,%r12,8)
add    $0x1,%rdx
cmp    %rdx,%r15
jne    11e0 <choleskyDecomposition+0xd0>
mov    %r8,%rcx
cmp    %r14,%rax
jb     121d <choleskyDecomposition+0x10d>
jmp    1150 <choleskyDecomposition+0x40>
nopl   (%rax)
add    $0x1,%rcx
cmp    %r14,%rcx
je     1150 <choleskyDecomposition+0x40>
mov    (%rbx,%rcx,8),%rdx
movsd  (%rdx,%r12,8),%xmm0
movsd  (%rdx,%rcx,8),%xmm1
mulsd  %xmm0,%xmm0
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rcx,8)
lea    0x1(%rcx),%rsi
cmp    %ebp,%esi
jge    1210 <choleskyDecomposition+0x100>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    (%rbx,%rsi,8),%rdi
movsd  (%rdi,%r12,8),%xmm0
movsd  (%rdi,%rcx,8),%xmm1
mulsd  (%rdx,%r12,8),%xmm0
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%rcx,8)
add    $0x1,%rsi
cmp    %ebp,%esi
jl     1250 <choleskyDecomposition+0x140>
jmp    1210 <choleskyDecomposition+0x100>
mov    $0x1,%eax
jmp    1281 <choleskyDecomposition+0x171>
xor    %eax,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @choleskyDecomposition(double** nocapture %0, i32 %1) {
%3 = icmp sgt i32 %1, 0
br i1 %3, label %4, label %74
4: ; preds = %2
%5 = zext i32 %1 to i64
%6 = zext i32 %1 to i64
br label %10
7: ; preds = %71, %18, %32
%8 = add nuw nsw i64 %12, 1
%9 = icmp eq i64 %22, %6
br i1 %9, label %74, label %10
10: ; preds = %4, %7
%11 = phi i64 [ 0, %4 ], [ %22, %7 ]
%12 = phi i64 [ 1, %4 ], [ %8, %7 ]
%13 = getelementptr inbounds double*, double** %0, i64 %11
%14 = load double*, double** %13
%15 = getelementptr inbounds double, double* %14, i64 %11
%16 = load double, double* %15
%17 = fcmp ugt double %16, 0.000000e+00
br i1 %17, label %18, label %74
18: ; preds = %10
%19 = tail call double @sqrt(double %16)
%20 = load double*, double** %13
%21 = getelementptr inbounds double, double* %20, i64 %11
store double %19, double* %21
%22 = add nuw nsw i64 %11, 1
%23 = icmp ult i64 %22, %5
br i1 %23, label %24, label %7
24: ; preds = %18
%25 = getelementptr inbounds double*, double** %0, i64 %12
%26 = load double*, double** %25
%27 = getelementptr inbounds double, double* %26, i64 %11
%28 = load double, double* %27
%29 = fdiv double %28, %19
store double %29, double* %27
%30 = add nuw nsw i64 %12, 1
%31 = icmp eq i64 %30, %6
br i1 %31, label %32, label %33
32: ; preds = %33, %24
br i1 %23, label %43, label %7
33: ; preds = %24, %33
%34 = phi i64 [ %41, %33 ], [ %30, %24 ]
%35 = load double, double* %21
%36 = getelementptr inbounds double*, double** %0, i64 %34
%37 = load double*, double** %36
%38 = getelementptr inbounds double, double* %37, i64 %11
%39 = load double, double* %38
%40 = fdiv double %39, %35
store double %40, double* %38
%41 = add nuw nsw i64 %34, 1
%42 = icmp eq i64 %41, %6
br i1 %42, label %32, label %33
43: ; preds = %32, %71
%44 = phi i64 [ %72, %71 ], [ %12, %32 ]
%45 = getelementptr inbounds double*, double** %0, i64 %44
%46 = load double*, double** %45
%47 = getelementptr inbounds double, double* %46, i64 %11
%48 = getelementptr inbounds double, double* %46, i64 %11
%49 = load double, double* %48
%50 = getelementptr inbounds double, double* %46, i64 %44
%51 = load double, double* %50
%52 = fneg double %49
%53 = tail call double @llvm.fmuladd.f64(double %52, double %49, double %51)
store double %53, double* %50
%54 = add nuw nsw i64 %44, 1
%55 = trunc i64 %54 to i32
%56 = icmp slt i32 %55, %1
br i1 %56, label %57, label %71
57: ; preds = %43, %57
%58 = phi i64 [ %68, %57 ], [ %54, %43 ]
%59 = getelementptr inbounds double*, double** %0, i64 %58
%60 = load double*, double** %59
%61 = getelementptr inbounds double, double* %60, i64 %11
%62 = load double, double* %61
%63 = load double, double* %47
%64 = getelementptr inbounds double, double* %60, i64 %44
%65 = load double, double* %64
%66 = fneg double %62
%67 = tail call double @llvm.fmuladd.f64(double %66, double %63, double %65)
store double %67, double* %64
%68 = add nuw nsw i64 %58, 1
%69 = trunc i64 %68 to i32
%70 = icmp slt i32 %69, %1
br i1 %70, label %57, label %71
71: ; preds = %57, %43
%72 = add nuw nsw i64 %44, 1
%73 = icmp eq i64 %72, %6
br i1 %73, label %7, label %43
74: ; preds = %10, %7, %2
%75 = phi i32 [ 1, %2 ], [ 1, %7 ], [ 0, %10 ]
ret i32 %75
}