#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*/ fprintf(int /*<<< orphan*/, char *, ...);
int /*<<< orphan*/ stderr;
void lu_dcmp(double **A, int n, double *b, int reduce, int solve, int *pd) {
  double pivot;
  int i, j, k;
  *pd = 1;
  if (reduce) {
    for (k = 1; k <= n; k++) {
      if (0.0 == (pivot = A[k][k])) {
        fprintf(stderr, " lu_dcmp: zero found on the diagonal\n");
        fprintf(stderr, " A[%d][%d] = %11.4e\n", k, k, A[k][k]);
        *pd = 0;
        return;
      }
      for (i = k + 1; i <= n; i++) {
        A[i][k] /= pivot;
        for (j = k + 1; j <= n; j++)
          A[i][j] -= A[i][k] * A[k][j];
      }
    }
  }
  if (solve) {
    for (k = 1; k <= n; k++)
      for (i = k + 1; i <= n; i++)
        b[i] -= A[i][k] * b[k];
    for (j = n; j >= 2; j--)
      for (i = 1; i <= j - 1; i++)
        b[i] -= b[j] * A[i][j] / A[j][j];
    for (i = 1; i <= n; i++)
      b[i] /= A[i][i];
  }
  return;
}


x86 assembly (clang, -O0):
<lu_dcmp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0x28(%rbp),%rax
movl   $0x1,(%rax)
cmpl   $0x0,-0x1c(%rbp)
je     12c3 <lu_dcmp+0x1b3>
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     12be <lu_dcmp+0x1ae>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  %xmm1,-0x30(%rbp)
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jne    11da <lu_dcmp+0xca>
jp     11da <lu_dcmp+0xca>
mov    0x2e57(%rip),%rax
mov    (%rax),%edi
lea    0xe6e(%rip),%rsi
mov    $0x0,%al
call   1030 <fprintf@plt>
mov    0x2e40(%rip),%rax
mov    (%rax),%edi
mov    -0x3c(%rbp),%edx
mov    -0x3c(%rbp),%ecx
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rsi
mov    (%rax,%rsi,8),%rax
movslq -0x3c(%rbp),%rsi
movsd  (%rax,%rsi,8),%xmm0
lea    0xe62(%rip),%rsi
mov    $0x1,%al
call   1030 <fprintf@plt>
mov    -0x28(%rbp),%rax
movl   $0x0,(%rax)
jmp    1455 <lu_dcmp+0x345>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     12ab <lu_dcmp+0x19b>
movsd  -0x30(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1298 <lu_dcmp+0x188>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%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    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    121b <lu_dcmp+0x10b>
jmp    129d <lu_dcmp+0x18d>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11e3 <lu_dcmp+0xd3>
jmp    12b0 <lu_dcmp+0x1a0>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1149 <lu_dcmp+0x39>
jmp    12c3 <lu_dcmp+0x1b3>
cmpl   $0x0,-0x20(%rbp)
je     1450 <lu_dcmp+0x340>
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1369 <lu_dcmp+0x259>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1356 <lu_dcmp+0x246>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x34(%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    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    12e9 <lu_dcmp+0x1d9>
jmp    135b <lu_dcmp+0x24b>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    12d4 <lu_dcmp+0x1c4>
mov    -0xc(%rbp),%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x2,-0x38(%rbp)
jl     13ff <lu_dcmp+0x2ef>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     13ec <lu_dcmp+0x2dc>
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1380 <lu_dcmp+0x270>
jmp    13f1 <lu_dcmp+0x2e1>
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
jmp    136f <lu_dcmp+0x25f>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     144b <lu_dcmp+0x33b>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1406 <lu_dcmp+0x2f6>
jmp    1450 <lu_dcmp+0x340>
jmp    1455 <lu_dcmp+0x345>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<lu_dcmp>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r8d,0x14(%rsp)
mov    %rsi,0x8(%rsp)
mov    %rdi,%r15
mov    %r9,0x18(%rsp)
movl   $0x1,(%r9)
test   %ecx,%ecx
je     1286 <lu_dcmp+0x176>
cmpl   $0x0,0x8(%rsp)
jle    1286 <lu_dcmp+0x176>
mov    0x8(%rsp),%rcx
lea    0x1(%rcx),%eax
mov    %ecx,%esi
mov    %rsi,0x20(%rsp)
lea    -0x3(%rcx),%r10d
mov    $0x2,%r9d
mov    $0x1,%r11d
xor    %r13d,%r13d
xorpd  %xmm0,%xmm0
jmp    1191 <lu_dcmp+0x81>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r9
add    $0x1,%r13d
cmp    %rax,%r11
je     1286 <lu_dcmp+0x176>
mov    %r11,%rbx
mov    (%r15,%r11,8),%rdi
movsd  (%rdi,%r11,8),%xmm1
ucomisd %xmm0,%xmm1
jne    11aa <lu_dcmp+0x9a>
jnp    1427 <lu_dcmp+0x317>
lea    0x1(%rbx),%r11
cmp    0x20(%rsp),%rbx
jae    1180 <lu_dcmp+0x70>
mov    0x8(%rsp),%rcx
mov    %ecx,%r12d
sub    %r13d,%r12d
lea    0x1(%r9),%rsi
lea    0x2(%r9),%r14
mov    %r9,%rcx
jmp    11d8 <lu_dcmp+0xc8>
nopl   (%rax)
add    $0x1,%rcx
cmp    %ecx,%eax
je     1180 <lu_dcmp+0x70>
mov    (%r15,%rcx,8),%r8
movsd  (%r8,%rbx,8),%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,(%r8,%rbx,8)
mulsd  (%rdi,%r9,8),%xmm2
movsd  (%r8,%r9,8),%xmm3
subsd  %xmm2,%xmm3
movsd  %xmm3,(%r8,%r9,8)
cmp    %esi,%eax
je     11d0 <lu_dcmp+0xc0>
mov    %rsi,%rbp
test   $0x1,%r12b
je     1231 <lu_dcmp+0x121>
movsd  (%r8,%rbx,8),%xmm2
movsd  0x8(%r8,%r9,8),%xmm3
mulsd  0x8(%rdi,%r9,8),%xmm2
subsd  %xmm2,%xmm3
movsd  %xmm3,0x8(%r8,%r9,8)
mov    %r14,%rbp
cmp    %r13d,%r10d
je     11d0 <lu_dcmp+0xc0>
cs nopw 0x0(%rax,%rax,1)
movsd  (%r8,%rbx,8),%xmm2
mulsd  (%rdi,%rbp,8),%xmm2
movsd  (%r8,%rbp,8),%xmm3
subsd  %xmm2,%xmm3
movsd  0x8(%r8,%rbp,8),%xmm2
movsd  %xmm3,(%r8,%rbp,8)
movsd  (%r8,%rbx,8),%xmm3
mulsd  0x8(%rdi,%rbp,8),%xmm3
subsd  %xmm3,%xmm2
movsd  %xmm2,0x8(%r8,%rbp,8)
add    $0x2,%rbp
cmp    %ebp,%eax
jne    1240 <lu_dcmp+0x130>
jmp    11d0 <lu_dcmp+0xc0>
cmpl   $0x0,0x14(%rsp)
mov    0x8(%rsp),%r11
je     146e <lu_dcmp+0x35e>
test   %r11d,%r11d
jle    146e <lu_dcmp+0x35e>
lea    0x1(%r11),%r10d
mov    %r11d,%ecx
lea    -0x2(%r11),%r8d
mov    $0x2,%ebp
mov    $0x1,%r9d
xor    %eax,%eax
jmp    12d0 <lu_dcmp+0x1c0>
nopl   0x0(%rax)
add    $0x1,%rbp
add    $0x1,%eax
cmp    %r10,%r9
je     1359 <lu_dcmp+0x249>
mov    %r9,%rdi
add    $0x1,%r9
cmp    %rcx,%rdi
jae    12c0 <lu_dcmp+0x1b0>
mov    %eax,%ebx
not    %ebx
add    %r11d,%ebx
mov    %rbp,%rsi
test   $0x1,%bl
je     130b <lu_dcmp+0x1fb>
mov    (%r15,%rbp,8),%rsi
movsd  (%rsi,%rdi,8),%xmm0
mulsd  (%rdx,%rdi,8),%xmm0
movsd  (%rdx,%rbp,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rbp,8)
lea    0x1(%rbp),%rsi
cmp    %eax,%r8d
je     12c0 <lu_dcmp+0x1b0>
mov    (%r15,%rsi,8),%rbx
movsd  (%rbx,%rdi,8),%xmm0
mulsd  (%rdx,%rdi,8),%xmm0
movsd  (%rdx,%rsi,8),%xmm1
subsd  %xmm0,%xmm1
movsd  0x8(%rdx,%rsi,8),%xmm0
movsd  %xmm1,(%rdx,%rsi,8)
mov    0x8(%r15,%rsi,8),%rbx
movsd  (%rbx,%rdi,8),%xmm1
mulsd  (%rdx,%rdi,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rdx,%rsi,8)
add    $0x2,%rsi
cmp    %esi,%r10d
jne    1310 <lu_dcmp+0x200>
jmp    12c0 <lu_dcmp+0x1b0>
cmp    $0x2,%r11d
jl     13a7 <lu_dcmp+0x297>
nop
mov    (%r15,%rcx,8),%rsi
mov    $0x1,%edi
nopl   0x0(%rax)
movsd  (%rdx,%rcx,8),%xmm0
mov    (%r15,%rdi,8),%rax
mulsd  (%rax,%rcx,8),%xmm0
divsd  (%rsi,%rcx,8),%xmm0
movsd  (%rdx,%rdi,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rdi,8)
add    $0x1,%rdi
cmp    %rdi,%rcx
jne    1370 <lu_dcmp+0x260>
lea    -0x1(%rcx),%rax
cmp    $0x2,%rcx
mov    %rax,%rcx
jg     1360 <lu_dcmp+0x250>
test   %r11d,%r11d
jle    146e <lu_dcmp+0x35e>
mov    %r10d,%eax
add    $0xffffffffffffffff,%r10
mov    $0x1,%ecx
cmp    $0x2,%eax
je     140c <lu_dcmp+0x2fc>
mov    %r10,%rsi
and    $0xfffffffffffffffe,%rsi
neg    %rsi
mov    $0x1,%ecx
mov    (%r15,%rcx,8),%rax
movsd  (%rdx,%rcx,8),%xmm0
movsd  0x8(%rdx,%rcx,8),%xmm1
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rdx,%rcx,8)
mov    0x8(%r15,%rcx,8),%rax
divsd  0x8(%rax,%rcx,8),%xmm1
movsd  %xmm1,0x8(%rdx,%rcx,8)
lea    (%rsi,%rcx,1),%rax
add    $0x2,%rax
add    $0x2,%rcx
cmp    $0x1,%rax
jne    13d0 <lu_dcmp+0x2c0>
test   $0x1,%r10b
je     146e <lu_dcmp+0x35e>
mov    (%r15,%rcx,8),%rax
movsd  (%rdx,%rcx,8),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rdx,%rcx,8)
jmp    146e <lu_dcmp+0x35e>
mov    %ebx,%ebp
shl    $0x3,%rbp
mov    0x2bac(%rip),%r14
mov    (%r14),%edi
lea    0xbc2(%rip),%rsi
xor    %eax,%eax
call   1030 <fprintf@plt>
mov    (%r14),%edi
mov    (%r15,%rbp,1),%rax
movsd  (%rax,%rbp,1),%xmm0
lea    0xbce(%rip),%rsi
mov    %ebx,%edx
mov    %ebx,%ecx
mov    $0x1,%al
call   1030 <fprintf@plt>
mov    0x18(%rsp),%rax
movl   $0x0,(%rax)
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@stderr = global i32 0
@.str = [38 x i8] c" lu_dcmp: zero found on the diagonal\0A\00"
@.str.1 = [21 x i8] c" A[%d][%d] = %11.4e\0A\00"
define void @lu_dcmp(double** %0, i32 %1, double* %2, i32 %3, i32 %4, i32* %5) {
%7 = alloca double**
%8 = alloca i32
%9 = alloca double*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32*
%13 = alloca double
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
store double** %0, double*** %7
store i32 %1, i32* %8
store double* %2, double** %9
store i32 %3, i32* %10
store i32 %4, i32* %11
store i32* %5, i32** %12
%17 = load i32*, i32** %12
store i32 1, i32* %17
%18 = load i32, i32* %10
%19 = icmp ne i32 %18, 0
br i1 %19, label %20, label %120
20: ; preds = %6
store i32 1, i32* %16
br label %21
21: ; preds = %116, %20
%22 = load i32, i32* %16
%23 = load i32, i32* %8
%24 = icmp sle i32 %22, %23
br i1 %24, label %25, label %119
25: ; preds = %21
%26 = load double**, double*** %7
%27 = load i32, i32* %16
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds double*, double** %26, i64 %28
%30 = load double*, double** %29
%31 = load i32, i32* %16
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %30, i64 %32
%34 = load double, double* %33
store double %34, double* %13
%35 = fcmp oeq double 0.000000e+00, %34
br i1 %35, label %36, label %53
36: ; preds = %25
%37 = load i32, i32* @stderr
%38 = call i32 (i32, i8*, ...) @fprintf(i32 %37, i8* getelementptr inbounds ([38 x i8], [38 x i8]* @.str, i64 0, i64 0))
%39 = load i32, i32* @stderr
%40 = load i32, i32* %16
%41 = load i32, i32* %16
%42 = load double**, double*** %7
%43 = load i32, i32* %16
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double*, double** %42, i64 %44
%46 = load double*, double** %45
%47 = load i32, i32* %16
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %46, i64 %48
%50 = load double, double* %49
%51 = call i32 (i32, i8*, ...) @fprintf(i32 %39, i8* getelementptr inbounds ([21 x i8], [21 x i8]* @.str.1, i64 0, i64 0), i32 %40, i32 %41, double %50)
%52 = load i32*, i32** %12
store i32 0, i32* %52
br label %240
53: ; preds = %25
%54 = load i32, i32* %16
%55 = add nsw i32 %54, 1
store i32 %55, i32* %14
br label %56
56: ; preds = %112, %53
%57 = load i32, i32* %14
%58 = load i32, i32* %8
%59 = icmp sle i32 %57, %58
br i1 %59, label %60, label %115
60: ; preds = %56
%61 = load double, double* %13
%62 = load double**, double*** %7
%63 = load i32, i32* %14
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds double*, double** %62, i64 %64
%66 = load double*, double** %65
%67 = load i32, i32* %16
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %66, i64 %68
%70 = load double, double* %69
%71 = fdiv double %70, %61
store double %71, double* %69
%72 = load i32, i32* %16
%73 = add nsw i32 %72, 1
store i32 %73, i32* %15
br label %74
74: ; preds = %108, %60
%75 = load i32, i32* %15
%76 = load i32, i32* %8
%77 = icmp sle i32 %75, %76
br i1 %77, label %78, label %111
78: ; preds = %74
%79 = load double**, double*** %7
%80 = load i32, i32* %14
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double*, double** %79, i64 %81
%83 = load double*, double** %82
%84 = load i32, i32* %16
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %83, i64 %85
%87 = load double, double* %86
%88 = load double**, double*** %7
%89 = load i32, i32* %16
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds double*, double** %88, i64 %90
%92 = load double*, double** %91
%93 = load i32, i32* %15
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds double, double* %92, i64 %94
%96 = load double, double* %95
%97 = load double**, double*** %7
%98 = load i32, i32* %14
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double*, double** %97, i64 %99
%101 = load double*, double** %100
%102 = load i32, i32* %15
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds double, double* %101, i64 %103
%105 = load double, double* %104
%106 = fneg double %87
%107 = call double @llvm.fmuladd.f64(double %106, double %96, double %105)
store double %107, double* %104
br label %108
108: ; preds = %78
%109 = load i32, i32* %15
%110 = add nsw i32 %109, 1
store i32 %110, i32* %15
br label %74
111: ; preds = %74
br label %112
112: ; preds = %111
%113 = load i32, i32* %14
%114 = add nsw i32 %113, 1
store i32 %114, i32* %14
br label %56
115: ; preds = %56
br label %116
116: ; preds = %115
%117 = load i32, i32* %16
%118 = add nsw i32 %117, 1
store i32 %118, i32* %16
br label %21
119: ; preds = %21
br label %120
120: ; preds = %119, %6
%121 = load i32, i32* %11
%122 = icmp ne i32 %121, 0
br i1 %122, label %123, label %239
123: ; preds = %120
store i32 1, i32* %16
br label %124
124: ; preds = %161, %123
%125 = load i32, i32* %16
%126 = load i32, i32* %8
%127 = icmp sle i32 %125, %126
br i1 %127, label %128, label %164
128: ; preds = %124
%129 = load i32, i32* %16
%130 = add nsw i32 %129, 1
store i32 %130, i32* %14
br label %131
131: ; preds = %157, %128
%132 = load i32, i32* %14
%133 = load i32, i32* %8
%134 = icmp sle i32 %132, %133
br i1 %134, label %135, label %160
135: ; preds = %131
%136 = load double**, double*** %7
%137 = load i32, i32* %14
%138 = sext i32 %137 to i64
%139 = getelementptr inbounds double*, double** %136, i64 %138
%140 = load double*, double** %139
%141 = load i32, i32* %16
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds double, double* %140, i64 %142
%144 = load double, double* %143
%145 = load double*, double** %9
%146 = load i32, i32* %16
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds double, double* %145, i64 %147
%149 = load double, double* %148
%150 = load double*, double** %9
%151 = load i32, i32* %14
%152 = sext i32 %151 to i64
%153 = getelementptr inbounds double, double* %150, i64 %152
%154 = load double, double* %153
%155 = fneg double %144
%156 = call double @llvm.fmuladd.f64(double %155, double %149, double %154)
store double %156, double* %153
br label %157
157: ; preds = %135
%158 = load i32, i32* %14
%159 = add nsw i32 %158, 1
store i32 %159, i32* %14
br label %131
160: ; preds = %131
br label %161
161: ; preds = %160
%162 = load i32, i32* %16
%163 = add nsw i32 %162, 1
store i32 %163, i32* %16
br label %124
164: ; preds = %124
%165 = load i32, i32* %8
store i32 %165, i32* %15
br label %166
166: ; preds = %211, %164
%167 = load i32, i32* %15
%168 = icmp sge i32 %167, 2
br i1 %168, label %169, label %214
169: ; preds = %166
store i32 1, i32* %14
br label %170
170: ; preds = %207, %169
%171 = load i32, i32* %14
%172 = load i32, i32* %15
%173 = sub nsw i32 %172, 1
%174 = icmp sle i32 %171, %173
br i1 %174, label %175, label %210
175: ; preds = %170
%176 = load double*, double** %9
%177 = load i32, i32* %15
%178 = sext i32 %177 to i64
%179 = getelementptr inbounds double, double* %176, i64 %178
%180 = load double, double* %179
%181 = load double**, double*** %7
%182 = load i32, i32* %14
%183 = sext i32 %182 to i64
%184 = getelementptr inbounds double*, double** %181, i64 %183
%185 = load double*, double** %184
%186 = load i32, i32* %15
%187 = sext i32 %186 to i64
%188 = getelementptr inbounds double, double* %185, i64 %187
%189 = load double, double* %188
%190 = fmul double %180, %189
%191 = load double**, double*** %7
%192 = load i32, i32* %15
%193 = sext i32 %192 to i64
%194 = getelementptr inbounds double*, double** %191, i64 %193
%195 = load double*, double** %194
%196 = load i32, i32* %15
%197 = sext i32 %196 to i64
%198 = getelementptr inbounds double, double* %195, i64 %197
%199 = load double, double* %198
%200 = fdiv double %190, %199
%201 = load double*, double** %9
%202 = load i32, i32* %14
%203 = sext i32 %202 to i64
%204 = getelementptr inbounds double, double* %201, i64 %203
%205 = load double, double* %204
%206 = fsub double %205, %200
store double %206, double* %204
br label %207
207: ; preds = %175
%208 = load i32, i32* %14
%209 = add nsw i32 %208, 1
store i32 %209, i32* %14
br label %170
210: ; preds = %170
br label %211
211: ; preds = %210
%212 = load i32, i32* %15
%213 = add nsw i32 %212, -1
store i32 %213, i32* %15
br label %166
214: ; preds = %166
store i32 1, i32* %14
br label %215
215: ; preds = %235, %214
%216 = load i32, i32* %14
%217 = load i32, i32* %8
%218 = icmp sle i32 %216, %217
br i1 %218, label %219, label %238
219: ; preds = %215
%220 = load double**, double*** %7
%221 = load i32, i32* %14
%222 = sext i32 %221 to i64
%223 = getelementptr inbounds double*, double** %220, i64 %222
%224 = load double*, double** %223
%225 = load i32, i32* %14
%226 = sext i32 %225 to i64
%227 = getelementptr inbounds double, double* %224, i64 %226
%228 = load double, double* %227
%229 = load double*, double** %9
%230 = load i32, i32* %14
%231 = sext i32 %230 to i64
%232 = getelementptr inbounds double, double* %229, i64 %231
%233 = load double, double* %232
%234 = fdiv double %233, %228
store double %234, double* %232
br label %235
235: ; preds = %219
%236 = load i32, i32* %14
%237 = add nsw i32 %236, 1
store i32 %237, i32* %14
br label %215
238: ; preds = %215
br label %239
239: ; preds = %238, %120
br label %240
240: ; preds = %239, %36
ret void
}

LLVM IR (-O2):
@stderr = global i32 0
@.str = [38 x i8] c" lu_dcmp: zero found on the diagonal\0A\00"
@.str.1 = [21 x i8] c" A[%d][%d] = %11.4e\0A\00"
define void @lu_dcmp(double** nocapture %0, i32 %1, double* nocapture %2, i32 %3, i32 %4, i32* nocapture writeonly %5) {
store i32 1, i32* %5
%7 = icmp eq i32 %3, 0
%8 = icmp slt i32 %1, 1
%9 = or i1 %7, %8
br i1 %9, label %99, label %10
10: ; preds = %6
%11 = add nuw i32 %1, 1
%12 = zext i32 %1 to i64
%13 = zext i32 %11 to i64
%14 = add i32 %1, -3
br label %19
15: ; preds = %95, %40
%16 = add nuw nsw i64 %22, 1
%17 = icmp eq i64 %41, %13
%18 = add i32 %20, 1
br i1 %17, label %99, label %19
19: ; preds = %10, %15
%20 = phi i32 [ 0, %10 ], [ %18, %15 ]
%21 = phi i64 [ 1, %10 ], [ %41, %15 ]
%22 = phi i64 [ 2, %10 ], [ %16, %15 ]
%23 = sub i32 %1, %20
%24 = getelementptr inbounds double*, double** %0, i64 %21
%25 = load double*, double** %24
%26 = getelementptr inbounds double, double* %25, i64 %21
%27 = load double, double* %26
%28 = fcmp oeq double %27, 0.000000e+00
br i1 %28, label %29, label %40
29: ; preds = %19
%30 = trunc i64 %21 to i32
%31 = and i64 %21, 4294967295
%32 = getelementptr inbounds double*, double** %0, i64 %31
%33 = load i32, i32* @stderr
%34 = tail call i32 (i32, i8*, ...) @fprintf(i32 %33, i8* getelementptr inbounds ([38 x i8], [38 x i8]* @.str, i64 0, i64 0))
%35 = load i32, i32* @stderr
%36 = load double*, double** %32
%37 = getelementptr inbounds double, double* %36, i64 %31
%38 = load double, double* %37
%39 = tail call i32 (i32, i8*, ...) @fprintf(i32 %35, i8* getelementptr inbounds ([21 x i8], [21 x i8]* @.str.1, i64 0, i64 0), i32 %30, i32 %30, double %38)
store i32 0, i32* %5
br label %228
40: ; preds = %19
%41 = add nuw nsw i64 %21, 1
%42 = icmp ult i64 %21, %12
br i1 %42, label %43, label %15
43: ; preds = %40
%44 = getelementptr inbounds double, double* %25, i64 %22
%45 = add nuw nsw i64 %22, 1
%46 = trunc i64 %45 to i32
%47 = icmp eq i32 %11, %46
%48 = and i32 %23, 1
%49 = icmp eq i32 %48, 0
%50 = getelementptr inbounds double, double* %25, i64 %45
%51 = add nuw nsw i64 %22, 2
%52 = icmp eq i32 %14, %20
br label %53
53: ; preds = %43, %95
%54 = phi i64 [ %96, %95 ], [ %22, %43 ]
%55 = getelementptr inbounds double*, double** %0, i64 %54
%56 = load double*, double** %55
%57 = getelementptr inbounds double, double* %56, i64 %21
%58 = load double, double* %57
%59 = fdiv double %58, %27
store double %59, double* %57
%60 = load double, double* %44
%61 = getelementptr inbounds double, double* %56, i64 %22
%62 = load double, double* %61
%63 = fneg double %59
%64 = tail call double @llvm.fmuladd.f64(double %63, double %60, double %62)
store double %64, double* %61
br i1 %47, label %95, label %65
65: ; preds = %53
br i1 %49, label %73, label %66
66: ; preds = %65
%67 = load double, double* %57
%68 = load double, double* %50
%69 = getelementptr inbounds double, double* %56, i64 %45
%70 = load double, double* %69
%71 = fneg double %67
%72 = tail call double @llvm.fmuladd.f64(double %71, double %68, double %70)
store double %72, double* %69
br label %73
73: ; preds = %66, %65
%74 = phi i64 [ %45, %65 ], [ %51, %66 ]
br i1 %52, label %95, label %75
75: ; preds = %73, %75
%76 = phi i64 [ %92, %75 ], [ %74, %73 ]
%77 = load double, double* %57
%78 = getelementptr inbounds double, double* %25, i64 %76
%79 = load double, double* %78
%80 = getelementptr inbounds double, double* %56, i64 %76
%81 = load double, double* %80
%82 = fneg double %77
%83 = tail call double @llvm.fmuladd.f64(double %82, double %79, double %81)
store double %83, double* %80
%84 = add nuw nsw i64 %76, 1
%85 = load double, double* %57
%86 = getelementptr inbounds double, double* %25, i64 %84
%87 = load double, double* %86
%88 = getelementptr inbounds double, double* %56, i64 %84
%89 = load double, double* %88
%90 = fneg double %85
%91 = tail call double @llvm.fmuladd.f64(double %90, double %87, double %89)
store double %91, double* %88
%92 = add nuw nsw i64 %76, 2
%93 = trunc i64 %92 to i32
%94 = icmp eq i32 %11, %93
br i1 %94, label %95, label %75
95: ; preds = %73, %75, %53
%96 = add nuw nsw i64 %54, 1
%97 = trunc i64 %96 to i32
%98 = icmp eq i32 %11, %97
br i1 %98, label %15, label %53
99: ; preds = %15, %6
%100 = icmp eq i32 %4, 0
%101 = or i1 %100, %8
br i1 %101, label %228, label %102
102: ; preds = %99
%103 = add nuw i32 %1, 1
%104 = zext i32 %1 to i64
%105 = zext i32 %103 to i64
%106 = add i32 %1, -2
br label %115
107: ; preds = %138, %141, %115
%108 = add nuw nsw i64 %118, 1
%109 = icmp eq i64 %119, %105
%110 = add i32 %116, 1
br i1 %109, label %111, label %115
111: ; preds = %107
%112 = icmp sgt i32 %1, 1
br i1 %112, label %113, label %171
113: ; preds = %111
%114 = zext i32 %1 to i64
br label %165
115: ; preds = %102, %107
%116 = phi i32 [ 0, %102 ], [ %110, %107 ]
%117 = phi i64 [ 1, %102 ], [ %119, %107 ]
%118 = phi i64 [ 2, %102 ], [ %108, %107 ]
%119 = add nuw nsw i64 %117, 1
%120 = icmp ult i64 %117, %104
br i1 %120, label %121, label %107
121: ; preds = %115
%122 = xor i32 %116, -1
%123 = add i32 %122, %1
%124 = getelementptr inbounds double, double* %2, i64 %117
%125 = and i32 %123, 1
%126 = icmp eq i32 %125, 0
br i1 %126, label %138, label %127
127: ; preds = %121
%128 = getelementptr inbounds double*, double** %0, i64 %118
%129 = load double*, double** %128
%130 = getelementptr inbounds double, double* %129, i64 %117
%131 = load double, double* %130
%132 = load double, double* %124
%133 = getelementptr inbounds double, double* %2, i64 %118
%134 = load double, double* %133
%135 = fneg double %131
%136 = tail call double @llvm.fmuladd.f64(double %135, double %132, double %134)
store double %136, double* %133
%137 = add nuw nsw i64 %118, 1
br label %138
138: ; preds = %127, %121
%139 = phi i64 [ %118, %121 ], [ %137, %127 ]
%140 = icmp eq i32 %106, %116
br i1 %140, label %107, label %141
141: ; preds = %138, %141
%142 = phi i64 [ %162, %141 ], [ %139, %138 ]
%143 = getelementptr inbounds double*, double** %0, i64 %142
%144 = load double*, double** %143
%145 = getelementptr inbounds double, double* %144, i64 %117
%146 = load double, double* %145
%147 = load double, double* %124
%148 = getelementptr inbounds double, double* %2, i64 %142
%149 = load double, double* %148
%150 = fneg double %146
%151 = tail call double @llvm.fmuladd.f64(double %150, double %147, double %149)
store double %151, double* %148
%152 = add nuw nsw i64 %142, 1
%153 = getelementptr inbounds double*, double** %0, i64 %152
%154 = load double*, double** %153
%155 = getelementptr inbounds double, double* %154, i64 %117
%156 = load double, double* %155
%157 = load double, double* %124
%158 = getelementptr inbounds double, double* %2, i64 %152
%159 = load double, double* %158
%160 = fneg double %156
%161 = tail call double @llvm.fmuladd.f64(double %160, double %157, double %159)
store double %161, double* %158
%162 = add nuw nsw i64 %142, 2
%163 = trunc i64 %162 to i32
%164 = icmp eq i32 %103, %163
br i1 %164, label %107, label %141
165: ; preds = %193, %113
%166 = phi i64 [ %114, %113 ], [ %194, %193 ]
%167 = getelementptr inbounds double, double* %2, i64 %166
%168 = getelementptr inbounds double*, double** %0, i64 %166
%169 = load double*, double** %168
%170 = getelementptr inbounds double, double* %169, i64 %166
br label %178
171: ; preds = %193, %111
br i1 %8, label %228, label %172
172: ; preds = %171
%173 = add nsw i64 %105, -1
%174 = and i64 %173, 1
%175 = icmp eq i32 %103, 2
br i1 %175, label %217, label %176
176: ; preds = %172
%177 = and i64 %173, -2
br label %196
178: ; preds = %165, %178
%179 = phi i64 [ 1, %165 ], [ %191, %178 ]
%180 = load double, double* %167
%181 = getelementptr inbounds double*, double** %0, i64 %179
%182 = load double*, double** %181
%183 = getelementptr inbounds double, double* %182, i64 %166
%184 = load double, double* %183
%185 = fmul double %180, %184
%186 = load double, double* %170
%187 = fdiv double %185, %186
%188 = getelementptr inbounds double, double* %2, i64 %179
%189 = load double, double* %188
%190 = fsub double %189, %187
store double %190, double* %188
%191 = add nuw nsw i64 %179, 1
%192 = icmp eq i64 %191, %166
br i1 %192, label %193, label %178
193: ; preds = %178
%194 = add nsw i64 %166, -1
%195 = icmp sgt i64 %166, 2
br i1 %195, label %165, label %171
196: ; preds = %196, %176
%197 = phi i64 [ 1, %176 ], [ %214, %196 ]
%198 = phi i64 [ 0, %176 ], [ %215, %196 ]
%199 = getelementptr inbounds double*, double** %0, i64 %197
%200 = load double*, double** %199
%201 = getelementptr inbounds double, double* %200, i64 %197
%202 = load double, double* %201
%203 = getelementptr inbounds double, double* %2, i64 %197
%204 = load double, double* %203
%205 = fdiv double %204, %202
store double %205, double* %203
%206 = add nuw nsw i64 %197, 1
%207 = getelementptr inbounds double*, double** %0, i64 %206
%208 = load double*, double** %207
%209 = getelementptr inbounds double, double* %208, i64 %206
%210 = load double, double* %209
%211 = getelementptr inbounds double, double* %2, i64 %206
%212 = load double, double* %211
%213 = fdiv double %212, %210
store double %213, double* %211
%214 = add nuw nsw i64 %197, 2
%215 = add i64 %198, 2
%216 = icmp eq i64 %215, %177
br i1 %216, label %217, label %196
217: ; preds = %196, %172
%218 = phi i64 [ 1, %172 ], [ %214, %196 ]
%219 = icmp eq i64 %174, 0
br i1 %219, label %228, label %220
220: ; preds = %217
%221 = getelementptr inbounds double*, double** %0, i64 %218
%222 = load double*, double** %221
%223 = getelementptr inbounds double, double* %222, i64 %218
%224 = load double, double* %223
%225 = getelementptr inbounds double, double* %2, i64 %218
%226 = load double, double* %225
%227 = fdiv double %226, %224
store double %227, double* %225
br label %228
228: ; preds = %220, %217, %99, %171, %29
ret void
}