void factoriza_lu(double **A, long n) {
  long i, j, k;
  for (k = 0; k < n; k++) {
    for (j = k + 1; j < n; j++) {
      A[k][j] = A[k][j] / A[k][k];
    }
    for (i = k + 1; i < n; i++) {
      for (j = k + 1; j < n; j++) {
        A[i][j] = A[i][j] - A[i][k] * A[k][j];
      }
    }
  }
}


x86 assembly (clang, -O0):
<factoriza_lu>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    1270 <factoriza_lu+0x170>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    118c <factoriza_lu+0x8c>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x28(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    112e <factoriza_lu+0x2e>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x18(%rbp)
mov    -0x18(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    125a <factoriza_lu+0x15a>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    1244 <factoriza_lu+0x144>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%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    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    11b2 <factoriza_lu+0xb2>
jmp    1249 <factoriza_lu+0x149>
mov    -0x18(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x18(%rbp)
jmp    1198 <factoriza_lu+0x98>
jmp    125f <factoriza_lu+0x15f>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    1114 <factoriza_lu+0x14>
pop    %rbp
ret

x86 assembly (clang, -O2):
<factoriza_lu>:
test   %rsi,%rsi
jle    11a0 <factoriza_lu+0xa0>
xor    %r8d,%r8d
jmp    1119 <factoriza_lu+0x19>
xchg   %ax,%ax
cmp    %rsi,%r8
je     11a0 <factoriza_lu+0xa0>
mov    %r8,%r10
add    $0x1,%r8
cmp    %rsi,%r8
jge    1110 <factoriza_lu+0x10>
mov    (%rdi,%r10,8),%rax
mov    %r8,%rcx
nopl   0x0(%rax)
movsd  (%rax,%rcx,8),%xmm0
divsd  (%rax,%r10,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
add    $0x1,%rcx
cmp    %rsi,%rcx
jl     1130 <factoriza_lu+0x30>
cmp    %rsi,%r8
jge    1110 <factoriza_lu+0x10>
mov    (%rdi,%r10,8),%rdx
mov    %r8,%r9
cs nopw 0x0(%rax,%rax,1)
nop
mov    (%rdi,%r9,8),%rax
mov    %r8,%rcx
nopw   0x0(%rax,%rax,1)
movsd  (%rax,%rcx,8),%xmm0
movsd  (%rax,%r10,8),%xmm1
mulsd  (%rdx,%rcx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
add    $0x1,%rcx
cmp    %rsi,%rcx
jl     1170 <factoriza_lu+0x70>
add    $0x1,%r9
cmp    %rsi,%r9
jl     1160 <factoriza_lu+0x60>
jmp    1110 <factoriza_lu+0x10>
ret

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

LLVM IR (-O2):
define void @factoriza_lu(double** nocapture %0, i64 %1) {
%3 = icmp sgt i64 %1, 0
br i1 %3, label %6, label %45
4: ; preds = %42, %6, %14
%5 = icmp eq i64 %8, %1
br i1 %5, label %45, label %6
6: ; preds = %2, %4
%7 = phi i64 [ %8, %4 ], [ 0, %2 ]
%8 = add nuw nsw i64 %7, 1
%9 = icmp slt i64 %8, %1
br i1 %9, label %10, label %4
10: ; preds = %6
%11 = getelementptr inbounds double*, double** %0, i64 %7
%12 = load double*, double** %11
%13 = getelementptr inbounds double, double* %12, i64 %7
br label %18
14: ; preds = %18
br i1 %9, label %15, label %4
15: ; preds = %14
%16 = getelementptr inbounds double*, double** %0, i64 %7
%17 = load double*, double** %16
br label %26
18: ; preds = %10, %18
%19 = phi i64 [ %8, %10 ], [ %24, %18 ]
%20 = getelementptr inbounds double, double* %12, i64 %19
%21 = load double, double* %20
%22 = load double, double* %13
%23 = fdiv double %21, %22
store double %23, double* %20
%24 = add nuw nsw i64 %19, 1
%25 = icmp slt i64 %24, %1
br i1 %25, label %18, label %14
26: ; preds = %42, %15
%27 = phi i64 [ %8, %15 ], [ %43, %42 ]
%28 = getelementptr inbounds double*, double** %0, i64 %27
%29 = load double*, double** %28
%30 = getelementptr inbounds double, double* %29, i64 %7
br label %31
31: ; preds = %26, %31
%32 = phi i64 [ %8, %26 ], [ %40, %31 ]
%33 = getelementptr inbounds double, double* %29, i64 %32
%34 = load double, double* %33
%35 = load double, double* %30
%36 = getelementptr inbounds double, double* %17, i64 %32
%37 = load double, double* %36
%38 = fneg double %35
%39 = tail call double @llvm.fmuladd.f64(double %38, double %37, double %34)
store double %39, double* %33
%40 = add nuw nsw i64 %32, 1
%41 = icmp slt i64 %40, %1
br i1 %41, label %31, label %42
42: ; preds = %31
%43 = add nuw nsw i64 %27, 1
%44 = icmp slt i64 %43, %1
br i1 %44, label %26, label %4
45: ; preds = %4, %2
ret void
}