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])) {
        *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
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     1266 <lu_dcmp+0x166>
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1261 <lu_dcmp+0x161>
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    117d <lu_dcmp+0x7d>
jp     117d <lu_dcmp+0x7d>
mov    -0x28(%rbp),%rax
movl   $0x0,(%rax)
jmp    13f8 <lu_dcmp+0x2f8>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     124e <lu_dcmp+0x14e>
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     123b <lu_dcmp+0x13b>
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    11be <lu_dcmp+0xbe>
jmp    1240 <lu_dcmp+0x140>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1186 <lu_dcmp+0x86>
jmp    1253 <lu_dcmp+0x153>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1135 <lu_dcmp+0x35>
jmp    1266 <lu_dcmp+0x166>
cmpl   $0x0,-0x20(%rbp)
je     13f3 <lu_dcmp+0x2f3>
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     130c <lu_dcmp+0x20c>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     12f9 <lu_dcmp+0x1f9>
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    128c <lu_dcmp+0x18c>
jmp    12fe <lu_dcmp+0x1fe>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1277 <lu_dcmp+0x177>
mov    -0xc(%rbp),%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x2,-0x38(%rbp)
jl     13a2 <lu_dcmp+0x2a2>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     138f <lu_dcmp+0x28f>
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    1323 <lu_dcmp+0x223>
jmp    1394 <lu_dcmp+0x294>
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
jmp    1312 <lu_dcmp+0x212>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     13ee <lu_dcmp+0x2ee>
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    13a9 <lu_dcmp+0x2a9>
jmp    13f3 <lu_dcmp+0x2f3>
jmp    13f8 <lu_dcmp+0x2f8>
pop    %rbp
ret

x86 assembly (clang, -O2):
<lu_dcmp>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r8d,-0x14(%rsp)
mov    %rsi,-0x20(%rsp)
mov    %r9,-0x10(%rsp)
movl   $0x1,(%r9)
test   %ecx,%ecx
je     126a <lu_dcmp+0x16a>
cmpl   $0x0,-0x20(%rsp)
jle    126a <lu_dcmp+0x16a>
mov    -0x20(%rsp),%rcx
lea    0x1(%rcx),%eax
mov    %ecx,%esi
mov    %rsi,-0x8(%rsp)
lea    -0x3(%rcx),%r11d
mov    $0x2,%r8d
mov    $0x1,%r14d
xor    %r15d,%r15d
xorpd  %xmm0,%xmm0
jmp    1171 <lu_dcmp+0x71>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r8
add    $0x1,%r15d
cmp    %rax,%r14
je     126a <lu_dcmp+0x16a>
mov    (%rdi,%r14,8),%rbp
movsd  0x0(%rbp,%r14,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1188 <lu_dcmp+0x88>
jnp    1437 <lu_dcmp+0x337>
mov    %r14,%rbx
add    $0x1,%r14
cmp    -0x8(%rsp),%rbx
jae    1160 <lu_dcmp+0x60>
mov    -0x20(%rsp),%rcx
mov    %ecx,%r12d
sub    %r15d,%r12d
lea    0x1(%r8),%rsi
lea    0x2(%r8),%r10
mov    %r8,%rcx
jmp    11b8 <lu_dcmp+0xb8>
xchg   %ax,%ax
add    $0x1,%rcx
cmp    %ecx,%eax
je     1160 <lu_dcmp+0x60>
mov    (%rdi,%rcx,8),%r9
movsd  (%r9,%rbx,8),%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,(%r9,%rbx,8)
mulsd  0x0(%rbp,%r8,8),%xmm2
movsd  (%r9,%r8,8),%xmm3
subsd  %xmm2,%xmm3
movsd  %xmm3,(%r9,%r8,8)
cmp    %esi,%eax
je     11b0 <lu_dcmp+0xb0>
mov    %rsi,%r13
test   $0x1,%r12b
je     1212 <lu_dcmp+0x112>
movsd  (%r9,%rbx,8),%xmm2
movsd  0x8(%r9,%r8,8),%xmm3
mulsd  0x8(%rbp,%r8,8),%xmm2
subsd  %xmm2,%xmm3
movsd  %xmm3,0x8(%r9,%r8,8)
mov    %r10,%r13
cmp    %r15d,%r11d
je     11b0 <lu_dcmp+0xb0>
nopw   0x0(%rax,%rax,1)
movsd  (%r9,%rbx,8),%xmm2
mulsd  0x0(%rbp,%r13,8),%xmm2
movsd  (%r9,%r13,8),%xmm3
subsd  %xmm2,%xmm3
movsd  0x8(%r9,%r13,8),%xmm2
movsd  %xmm3,(%r9,%r13,8)
movsd  (%r9,%rbx,8),%xmm3
mulsd  0x8(%rbp,%r13,8),%xmm3
subsd  %xmm3,%xmm2
movsd  %xmm2,0x8(%r9,%r13,8)
add    $0x2,%r13
cmp    %r13d,%eax
jne    1220 <lu_dcmp+0x120>
jmp    11b0 <lu_dcmp+0xb0>
cmpl   $0x0,-0x14(%rsp)
mov    -0x20(%rsp),%r14
je     1442 <lu_dcmp+0x342>
test   %r14d,%r14d
jle    1442 <lu_dcmp+0x342>
lea    0x1(%r14),%r10d
mov    %r14d,%ecx
lea    -0x2(%r14),%r8d
mov    $0x2,%ebx
mov    $0x1,%r9d
xor    %r11d,%r11d
jmp    12b1 <lu_dcmp+0x1b1>
xchg   %ax,%ax
add    $0x1,%rbx
add    $0x1,%r11d
cmp    %r10,%r9
je     134b <lu_dcmp+0x24b>
mov    %r9,%rax
add    $0x1,%r9
cmp    %rcx,%rax
jae    12a0 <lu_dcmp+0x1a0>
mov    %r11d,%ebp
not    %ebp
add    %r14d,%ebp
mov    %rbx,%rsi
test   $0x1,%bpl
je     12ee <lu_dcmp+0x1ee>
mov    (%rdi,%rbx,8),%rsi
movsd  (%rsi,%rax,8),%xmm0
mulsd  (%rdx,%rax,8),%xmm0
movsd  (%rdx,%rbx,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rbx,8)
lea    0x1(%rbx),%rsi
cmp    %r11d,%r8d
je     12a0 <lu_dcmp+0x1a0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    (%rdi,%rsi,8),%rbp
movsd  0x0(%rbp,%rax,8),%xmm0
mulsd  (%rdx,%rax,8),%xmm0
movsd  (%rdx,%rsi,8),%xmm1
subsd  %xmm0,%xmm1
movsd  0x8(%rdx,%rsi,8),%xmm0
movsd  %xmm1,(%rdx,%rsi,8)
mov    0x8(%rdi,%rsi,8),%rbp
movsd  0x0(%rbp,%rax,8),%xmm1
mulsd  (%rdx,%rax,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rdx,%rsi,8)
add    $0x2,%rsi
cmp    %esi,%r10d
jne    1300 <lu_dcmp+0x200>
jmp    12a0 <lu_dcmp+0x1a0>
cmp    $0x2,%r14d
jl     13a8 <lu_dcmp+0x2a8>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    (%rdi,%rcx,8),%rax
mov    $0x1,%esi
nopl   0x0(%rax)
movsd  (%rdx,%rcx,8),%xmm0
mov    (%rdi,%rsi,8),%rbp
mulsd  0x0(%rbp,%rcx,8),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  (%rdx,%rsi,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    1370 <lu_dcmp+0x270>
lea    -0x1(%rcx),%rax
cmp    $0x2,%rcx
mov    %rax,%rcx
jg     1360 <lu_dcmp+0x260>
test   %r14d,%r14d
jle    1442 <lu_dcmp+0x342>
mov    %r10d,%ecx
add    $0xffffffffffffffff,%r10
mov    $0x1,%eax
cmp    $0x2,%ecx
je     141c <lu_dcmp+0x31c>
mov    %r10,%rcx
and    $0xfffffffffffffffe,%rcx
neg    %rcx
mov    $0x1,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    (%rdi,%rax,8),%rsi
movsd  (%rdx,%rax,8),%xmm0
movsd  0x8(%rdx,%rax,8),%xmm1
divsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
mov    0x8(%rdi,%rax,8),%rsi
divsd  0x8(%rsi,%rax,8),%xmm1
movsd  %xmm1,0x8(%rdx,%rax,8)
lea    (%rcx,%rax,1),%rsi
add    $0x2,%rsi
add    $0x2,%rax
cmp    $0x1,%rsi
jne    13e0 <lu_dcmp+0x2e0>
test   $0x1,%r10b
je     1442 <lu_dcmp+0x342>
mov    (%rdi,%rax,8),%rcx
movsd  (%rdx,%rax,8),%xmm0
divsd  (%rcx,%rax,8),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
jmp    1442 <lu_dcmp+0x342>
mov    -0x10(%rsp),%rax
movl   $0x0,(%rax)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
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 %105
20: ; preds = %6
store i32 1, i32* %16
br label %21
21: ; preds = %101, %20
%22 = load i32, i32* %16
%23 = load i32, i32* %8
%24 = icmp sle i32 %22, %23
br i1 %24, label %25, label %104
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 %38
36: ; preds = %25
%37 = load i32*, i32** %12
store i32 0, i32* %37
br label %225
38: ; preds = %25
%39 = load i32, i32* %16
%40 = add nsw i32 %39, 1
store i32 %40, i32* %14
br label %41
41: ; preds = %97, %38
%42 = load i32, i32* %14
%43 = load i32, i32* %8
%44 = icmp sle i32 %42, %43
br i1 %44, label %45, label %100
45: ; preds = %41
%46 = load double, double* %13
%47 = load double**, double*** %7
%48 = load i32, i32* %14
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double*, double** %47, i64 %49
%51 = load double*, double** %50
%52 = load i32, i32* %16
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %51, i64 %53
%55 = load double, double* %54
%56 = fdiv double %55, %46
store double %56, double* %54
%57 = load i32, i32* %16
%58 = add nsw i32 %57, 1
store i32 %58, i32* %15
br label %59
59: ; preds = %93, %45
%60 = load i32, i32* %15
%61 = load i32, i32* %8
%62 = icmp sle i32 %60, %61
br i1 %62, label %63, label %96
63: ; preds = %59
%64 = load double**, double*** %7
%65 = load i32, i32* %14
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double*, double** %64, i64 %66
%68 = load double*, double** %67
%69 = load i32, i32* %16
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds double, double* %68, i64 %70
%72 = load double, double* %71
%73 = load double**, double*** %7
%74 = load i32, i32* %16
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds double*, double** %73, i64 %75
%77 = load double*, double** %76
%78 = load i32, i32* %15
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds double, double* %77, i64 %79
%81 = load double, double* %80
%82 = load double**, double*** %7
%83 = load i32, i32* %14
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double*, double** %82, i64 %84
%86 = load double*, double** %85
%87 = load i32, i32* %15
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds double, double* %86, i64 %88
%90 = load double, double* %89
%91 = fneg double %72
%92 = call double @llvm.fmuladd.f64(double %91, double %81, double %90)
store double %92, double* %89
br label %93
93: ; preds = %63
%94 = load i32, i32* %15
%95 = add nsw i32 %94, 1
store i32 %95, i32* %15
br label %59
96: ; preds = %59
br label %97
97: ; preds = %96
%98 = load i32, i32* %14
%99 = add nsw i32 %98, 1
store i32 %99, i32* %14
br label %41
100: ; preds = %41
br label %101
101: ; preds = %100
%102 = load i32, i32* %16
%103 = add nsw i32 %102, 1
store i32 %103, i32* %16
br label %21
104: ; preds = %21
br label %105
105: ; preds = %104, %6
%106 = load i32, i32* %11
%107 = icmp ne i32 %106, 0
br i1 %107, label %108, label %224
108: ; preds = %105
store i32 1, i32* %16
br label %109
109: ; preds = %146, %108
%110 = load i32, i32* %16
%111 = load i32, i32* %8
%112 = icmp sle i32 %110, %111
br i1 %112, label %113, label %149
113: ; preds = %109
%114 = load i32, i32* %16
%115 = add nsw i32 %114, 1
store i32 %115, i32* %14
br label %116
116: ; preds = %142, %113
%117 = load i32, i32* %14
%118 = load i32, i32* %8
%119 = icmp sle i32 %117, %118
br i1 %119, label %120, label %145
120: ; preds = %116
%121 = load double**, double*** %7
%122 = load i32, i32* %14
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds double*, double** %121, i64 %123
%125 = load double*, double** %124
%126 = load i32, i32* %16
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds double, double* %125, i64 %127
%129 = load double, double* %128
%130 = load double*, double** %9
%131 = load i32, i32* %16
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double, double* %130, i64 %132
%134 = load double, double* %133
%135 = load double*, double** %9
%136 = load i32, i32* %14
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double, double* %135, i64 %137
%139 = load double, double* %138
%140 = fneg double %129
%141 = call double @llvm.fmuladd.f64(double %140, double %134, double %139)
store double %141, double* %138
br label %142
142: ; preds = %120
%143 = load i32, i32* %14
%144 = add nsw i32 %143, 1
store i32 %144, i32* %14
br label %116
145: ; preds = %116
br label %146
146: ; preds = %145
%147 = load i32, i32* %16
%148 = add nsw i32 %147, 1
store i32 %148, i32* %16
br label %109
149: ; preds = %109
%150 = load i32, i32* %8
store i32 %150, i32* %15
br label %151
151: ; preds = %196, %149
%152 = load i32, i32* %15
%153 = icmp sge i32 %152, 2
br i1 %153, label %154, label %199
154: ; preds = %151
store i32 1, i32* %14
br label %155
155: ; preds = %192, %154
%156 = load i32, i32* %14
%157 = load i32, i32* %15
%158 = sub nsw i32 %157, 1
%159 = icmp sle i32 %156, %158
br i1 %159, label %160, label %195
160: ; preds = %155
%161 = load double*, double** %9
%162 = load i32, i32* %15
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds double, double* %161, i64 %163
%165 = load double, double* %164
%166 = load double**, double*** %7
%167 = load i32, i32* %14
%168 = sext i32 %167 to i64
%169 = getelementptr inbounds double*, double** %166, i64 %168
%170 = load double*, double** %169
%171 = load i32, i32* %15
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds double, double* %170, i64 %172
%174 = load double, double* %173
%175 = fmul double %165, %174
%176 = load double**, double*** %7
%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 i32, i32* %15
%182 = sext i32 %181 to i64
%183 = getelementptr inbounds double, double* %180, i64 %182
%184 = load double, double* %183
%185 = fdiv double %175, %184
%186 = load double*, double** %9
%187 = load i32, i32* %14
%188 = sext i32 %187 to i64
%189 = getelementptr inbounds double, double* %186, i64 %188
%190 = load double, double* %189
%191 = fsub double %190, %185
store double %191, double* %189
br label %192
192: ; preds = %160
%193 = load i32, i32* %14
%194 = add nsw i32 %193, 1
store i32 %194, i32* %14
br label %155
195: ; preds = %155
br label %196
196: ; preds = %195
%197 = load i32, i32* %15
%198 = add nsw i32 %197, -1
store i32 %198, i32* %15
br label %151
199: ; preds = %151
store i32 1, i32* %14
br label %200
200: ; preds = %220, %199
%201 = load i32, i32* %14
%202 = load i32, i32* %8
%203 = icmp sle i32 %201, %202
br i1 %203, label %204, label %223
204: ; preds = %200
%205 = load double**, double*** %7
%206 = load i32, i32* %14
%207 = sext i32 %206 to i64
%208 = getelementptr inbounds double*, double** %205, i64 %207
%209 = load double*, double** %208
%210 = load i32, i32* %14
%211 = sext i32 %210 to i64
%212 = getelementptr inbounds double, double* %209, i64 %211
%213 = load double, double* %212
%214 = load double*, double** %9
%215 = load i32, i32* %14
%216 = sext i32 %215 to i64
%217 = getelementptr inbounds double, double* %214, i64 %216
%218 = load double, double* %217
%219 = fdiv double %218, %213
store double %219, double* %217
br label %220
220: ; preds = %204
%221 = load i32, i32* %14
%222 = add nsw i32 %221, 1
store i32 %222, i32* %14
br label %200
223: ; preds = %200
br label %224
224: ; preds = %223, %105
br label %225
225: ; preds = %224, %36
ret void
}

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