int choldc_e(a, n, p)
float **a, *p;
int n;
{
  int i, j, k;
  float sum;
  for (i = 1; i <= n; i++) {
    for (j = i; j <= n; j++) {
      for (sum = a[i][j], k = i - 1; k >= 1; k--)
        sum -= a[i][k] * a[j][k];
      if (i == j) {
        if (sum <= 0.0)
          return (1);
        p[i] = sqrt(sum);
      } else
        a[j][i] = sum / p[i];
    }
  }
  return (0);
}


x86 assembly (clang, -O0):
<choldc_e>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     126b <choldc_e+0x15b>
mov    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1258 <choldc_e+0x148>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x1,-0x2c(%rbp)
jl     11cc <choldc_e+0xbc>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x30(%rbp),%xmm1
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    116b <choldc_e+0x5b>
mov    -0x24(%rbp),%eax
cmp    -0x28(%rbp),%eax
jne    121e <choldc_e+0x10e>
movss  -0x30(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jb     11fa <choldc_e+0xea>
movl   $0x1,-0x4(%rbp)
jmp    1272 <choldc_e+0x162>
movss  -0x30(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <sqrt@plt>
cvtsd2ss %xmm0,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1245 <choldc_e+0x135>
movss  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    124a <choldc_e+0x13a>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113c <choldc_e+0x2c>
jmp    125d <choldc_e+0x14d>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    112a <choldc_e+0x1a>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<choldc_e>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
test   %esi,%esi
jle    123a <choldc_e+0x12a>
mov    %rdx,%r14
mov    %esi,%r15d
mov    %rdi,%r12
add    $0x1,%r15d
mov    $0x1,%ebp
xor    %r13d,%r13d
xorps  %xmm2,%xmm2
jmp    1151 <choldc_e+0x41>
add    $0x1,%rbp
add    $0x1,%r13
cmp    %r15,%rbp
je     123a <choldc_e+0x12a>
lea    -0x1(%rbp),%rax
mov    %rax,0x10(%rsp)
mov    %eax,%eax
mov    %rax,0x8(%rsp)
mov    %rbp,%rbx
jmp    1188 <choldc_e+0x78>
cs nopw 0x0(%rax,%rax,1)
divss  (%r14,%rbp,4),%xmm0
mov    (%r12,%rbx,8),%rax
movss  %xmm0,(%rax,%rbp,4)
add    $0x1,%rbx
cmp    %r15,%rbx
je     1140 <choldc_e+0x30>
mov    (%r12,%rbp,8),%rax
movss  (%rax,%rbx,4),%xmm0
cmp    $0x2,%rbp
jb     11fe <choldc_e+0xee>
mov    (%r12,%rbx,8),%rcx
test   $0x1,%r13b
jne    11b0 <choldc_e+0xa0>
mov    %rbp,%rdx
cmp    $0x1,%r13
jne    11d0 <choldc_e+0xc0>
jmp    11fe <choldc_e+0xee>
nopl   0x0(%rax)
mov    0x8(%rsp),%rdx
movss  (%rax,%rdx,4),%xmm1
mulss  (%rcx,%rdx,4),%xmm1
subss  %xmm1,%xmm0
mov    0x10(%rsp),%rdx
cmp    $0x1,%r13
je     11fe <choldc_e+0xee>
xchg   %ax,%ax
lea    -0x1(%rdx),%esi
movss  (%rax,%rsi,4),%xmm1
mulss  (%rcx,%rsi,4),%xmm1
subss  %xmm1,%xmm0
lea    -0x2(%rdx),%rsi
mov    %esi,%edi
movss  (%rax,%rdi,4),%xmm1
mulss  (%rcx,%rdi,4),%xmm1
subss  %xmm1,%xmm0
cmp    $0x3,%rdx
mov    %rsi,%rdx
jg     11d0 <choldc_e+0xc0>
cmp    %rbx,%rbp
jne    1170 <choldc_e+0x60>
ucomiss %xmm0,%xmm2
jae    1233 <choldc_e+0x123>
ucomiss %xmm2,%xmm0
jb     1220 <choldc_e+0x110>
sqrtss %xmm0,%xmm0
movss  %xmm0,(%r14,%rbp,4)
jmp    117f <choldc_e+0x6f>
call   1030 <sqrtf@plt>
xorps  %xmm2,%xmm2
movss  %xmm0,(%r14,%rbp,4)
jmp    117f <choldc_e+0x6f>
mov    $0x1,%eax
jmp    123c <choldc_e+0x12c>
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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