#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*/ lubksb(double **, int, int *, double *);
int /*<<< orphan*/ ludcmp(double **, int, int *, double *);
double calc_gyakugyouretu(double **ans, double **a, int n) {
  double y[n][n];
  double col[n];
  double d;
  int indx[n];
  int i, j;
  ludcmp(a, n, indx, &d);
  for (j = 0; j < n; j++)
    d *= (*((double *)a + j * n + j));
  for (j = 0; j < n; j++) {
    for (i = 0; i < n; i++)
      col[i] = 0.0;
    col[j] = 1.0;
    lubksb(a, n, indx, col);
    for (i = 0; i < n; i++)
      y[i][j] = col[i];
  }
  for (i = 0; i < n; i++)
    for (j = 0; j < n; j++)
      (*((double *)ans + i * n + j)) = y[i][j];
  return d;
}


x86 assembly (clang, -O0):
<calc_gyakugyouretu>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    %rax,-0x70(%rbp)
mov    %rsp,%rcx
mov    %rcx,-0x20(%rbp)
mov    %rax,%rcx
imul   %rcx,%rcx
lea    0xf(,%rcx,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x68(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x28(%rbp)
mov    %rax,-0x30(%rbp)
mov    -0x14(%rbp),%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x60(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x38(%rbp)
mov    -0x14(%rbp),%eax
lea    0xf(,%rax,4),%rcx
and    $0xfffffffffffffff0,%rcx
mov    %rsp,%rdx
sub    %rcx,%rdx
mov    %rdx,-0x58(%rbp)
mov    %rdx,%rsp
mov    %rax,-0x48(%rbp)
mov    -0x10(%rbp),%rdi
mov    -0x14(%rbp),%esi
lea    -0x40(%rbp),%rcx
call   1030 <ludcmp@plt>
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1202 <calc_gyakugyouretu+0xe2>
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    11c0 <calc_gyakugyouretu+0xa0>
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12cc <calc_gyakugyouretu+0x1ac>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1246 <calc_gyakugyouretu+0x126>
mov    -0x60(%rbp),%rax
movslq -0x4c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    121c <calc_gyakugyouretu+0xfc>
mov    -0x60(%rbp),%rcx
mov    -0x58(%rbp),%rdx
movslq -0x50(%rbp),%rax
movsd  0xda6(%rip),%xmm0
movsd  %xmm0,(%rcx,%rax,8)
mov    -0x10(%rbp),%rdi
mov    -0x14(%rbp),%esi
call   1040 <lubksb@plt>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12b9 <calc_gyakugyouretu+0x199>
mov    -0x68(%rbp),%rax
mov    -0x70(%rbp),%rdx
mov    -0x60(%rbp),%rcx
movslq -0x4c(%rbp),%rsi
movsd  (%rcx,%rsi,8),%xmm0
movslq -0x4c(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1272 <calc_gyakugyouretu+0x152>
jmp    12be <calc_gyakugyouretu+0x19e>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1209 <calc_gyakugyouretu+0xe9>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1351 <calc_gyakugyouretu+0x231>
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    133e <calc_gyakugyouretu+0x21e>
mov    -0x68(%rbp),%rax
mov    -0x70(%rbp),%rdx
movslq -0x4c(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x4c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    12e6 <calc_gyakugyouretu+0x1c6>
jmp    1343 <calc_gyakugyouretu+0x223>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    12d3 <calc_gyakugyouretu+0x1b3>
movsd  -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calc_gyakugyouretu>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %edx,%r15d
mov    %rdi,%r14
mov    %edx,%r12d
mov    %r12,%rax
imul   %r12,%rax
mov    %rsp,%rcx
lea    0xf(,%rax,8),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%rcx
mov    %rcx,-0x48(%rbp)
mov    %rcx,%rsp
mov    %rsp,%rbx
lea    0xf(,%r12,8),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%rbx
mov    %rbx,%rsp
mov    %rsp,%rdx
lea    0xf(,%r12,4),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%rdx
mov    %rdx,%rsp
lea    -0x60(%rbp),%rcx
mov    %rsi,-0x58(%rbp)
mov    %rsi,%rdi
mov    %r15d,%esi
mov    %rdx,-0x88(%rbp)
call   1030 <ludcmp@plt>
test   %r15d,%r15d
jle    14a1 <calc_gyakugyouretu+0x361>
mov    %r12,%r8
movsd  -0x60(%rbp),%xmm0
lea    -0x1(%r12),%rcx
mov    %r8d,%eax
and    $0x3,%eax
mov    %rcx,-0x80(%rbp)
cmp    $0x3,%rcx
jae    11f0 <calc_gyakugyouretu+0xb0>
xor    %ecx,%ecx
test   %rax,%rax
jne    1266 <calc_gyakugyouretu+0x126>
jmp    128d <calc_gyakugyouretu+0x14d>
mov    %r8,%rsi
shl    $0x4,%rsi
mov    %r8,%r10
shl    $0x5,%r10
add    $0x20,%r10
lea    0x0(,%r8,8),%rcx
lea    (%rcx,%rcx,2),%r9
mov    %r8d,%r11d
and    $0xfffffffc,%r11d
neg    %r11
mov    $0x1,%ecx
mov    -0x58(%rbp),%rdi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mulsd  (%rdi),%xmm0
mulsd  0x8(%rdi,%r8,8),%xmm0
mulsd  0x10(%rdi,%rsi,1),%xmm0
mulsd  0x18(%rdi,%r9,1),%xmm0
add    %r10,%rdi
lea    (%r11,%rcx,1),%rdx
add    $0x4,%rdx
add    $0x4,%rcx
cmp    $0x1,%rdx
jne    1230 <calc_gyakugyouretu+0xf0>
add    $0xffffffffffffffff,%rcx
test   %rax,%rax
je     128d <calc_gyakugyouretu+0x14d>
lea    0x8(,%r8,8),%rdx
imul   %rdx,%rcx
add    -0x58(%rbp),%rcx
cs nopw 0x0(%rax,%rax,1)
mulsd  (%rcx),%xmm0
add    %rdx,%rcx
add    $0xffffffffffffffff,%rax
jne    1280 <calc_gyakugyouretu+0x140>
movsd  %xmm0,-0x60(%rbp)
test   %r15d,%r15d
jle    14a1 <calc_gyakugyouretu+0x361>
mov    %r15,-0x78(%rbp)
mov    %r14,-0x68(%rbp)
lea    0x0(,%r8,8),%rdx
mov    %r8d,%r12d
and    $0x3,%r12d
mov    %r8d,%r13d
and    $0xfffffffc,%r13d
mov    %r8,%rsi
shl    $0x4,%rsi
mov    %r8,%r15
shl    $0x5,%r15
lea    (%rdx,%rdx,2),%r14
xor    %eax,%eax
mov    -0x48(%rbp),%r9
mov    %r8,-0x50(%rbp)
mov    %rdx,-0x38(%rbp)
mov    %rsi,-0x70(%rbp)
jmp    12f5 <calc_gyakugyouretu+0x1b5>
nop
mov    -0x40(%rbp),%rax
add    $0x1,%rax
add    $0x8,%r9
cmp    %r8,%rax
je     13ea <calc_gyakugyouretu+0x2aa>
mov    %rax,-0x40(%rbp)
mov    %r9,-0x30(%rbp)
mov    %rbx,%rdi
xor    %esi,%esi
call   1040 <memset@plt>
movabs $0x3ff0000000000000,%rax
mov    -0x40(%rbp),%rcx
mov    %rax,(%rbx,%rcx,8)
mov    -0x58(%rbp),%rdi
mov    -0x78(%rbp),%rsi
mov    -0x88(%rbp),%rdx
mov    %rbx,%rcx
call   1060 <lubksb@plt>
cmpq   $0x3,-0x80(%rbp)
jae    1350 <calc_gyakugyouretu+0x210>
xor    %eax,%eax
mov    -0x50(%rbp),%r8
mov    -0x30(%rbp),%r9
jmp    13a8 <calc_gyakugyouretu+0x268>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    -0x30(%rbp),%r9
mov    %r9,%rcx
xor    %eax,%eax
mov    -0x50(%rbp),%r8
mov    -0x70(%rbp),%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,(%rcx)
movsd  0x8(%rbx,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%r8,8)
movsd  0x10(%rbx,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%rsi,1)
movsd  0x18(%rbx,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%r14,1)
add    $0x4,%rax
add    %r15,%rcx
cmp    %rax,%r13
jne    1370 <calc_gyakugyouretu+0x230>
test   %r12,%r12
mov    -0x38(%rbp),%rdx
je     12e0 <calc_gyakugyouretu+0x1a0>
mov    %r8,%rcx
imul   %rax,%rcx
lea    (%r9,%rcx,8),%rcx
lea    (%rbx,%rax,8),%rax
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
movsd  (%rax,%rdi,8),%xmm0
movsd  %xmm0,(%rcx)
add    $0x1,%rdi
add    %rdx,%rcx
cmp    %rdi,%r12
jne    13d0 <calc_gyakugyouretu+0x290>
jmp    12e0 <calc_gyakugyouretu+0x1a0>
mov    -0x78(%rbp),%r12
test   %r12d,%r12d
jle    14a1 <calc_gyakugyouretu+0x361>
xor    %ebx,%ebx
cmpq   $0x0,-0x80(%rbp)
je     146f <calc_gyakugyouretu+0x32f>
and    $0xfffffffe,%r8d
mov    %r8,-0x40(%rbp)
lea    (%r12,%r12,1),%eax
mov    %eax,-0x30(%rbp)
mov    -0x48(%rbp),%r15
xor    %r14d,%r14d
mov    -0x68(%rbp),%r13
nopw   0x0(%rax,%rax,1)
movslq %ebx,%rbx
lea    0x0(,%rbx,8),%rdi
add    %r13,%rdi
mov    %r15,%rsi
call   1050 <memcpy@plt>
mov    -0x38(%rbp),%rax
lea    (%r15,%rax,1),%rsi
lea    (%r12,%rbx,1),%eax
cltq
lea    0x0(,%rax,8),%rdi
add    %r13,%rdi
mov    -0x38(%rbp),%rdx
call   1050 <memcpy@plt>
mov    -0x38(%rbp),%rdx
add    $0x2,%r14
add    -0x70(%rbp),%r15
add    -0x30(%rbp),%ebx
cmp    %r14,-0x40(%rbp)
jne    1420 <calc_gyakugyouretu+0x2e0>
jmp    1476 <calc_gyakugyouretu+0x336>
xor    %r14d,%r14d
mov    -0x68(%rbp),%r13
mov    -0x50(%rbp),%rax
test   $0x1,%al
je     14a1 <calc_gyakugyouretu+0x361>
imul   %r14,%rax
mov    -0x48(%rbp),%rcx
lea    (%rcx,%rax,8),%rsi
imul   %r12d,%r14d
movslq %r14d,%rax
lea    0x0(,%rax,8),%rdi
add    %r13,%rdi
call   1050 <memcpy@plt>
movsd  -0x60(%rbp),%xmm0
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @calc_gyakugyouretu(double** nocapture writeonly %0, double** %1, i32 %2) {
%4 = alloca double
%5 = zext i32 %2 to i64
%6 = mul nuw i64 %5, %5
%7 = alloca double, i64 %6
%8 = alloca double, i64 %5
%9 = bitcast double* %8 to i8*
%10 = bitcast double* %4 to i8*
%11 = alloca i32, i64 %5
%12 = call i32 @ludcmp(double** %1, i32 %2, i32* %11, double* %4)
%13 = icmp sgt i32 %2, 0
br i1 %13, label %14, label %169
14: ; preds = %3
%15 = load double, double* %4
%16 = zext i32 %2 to i64
%17 = add nsw i64 %5, -1
%18 = and i64 %5, 3
%19 = icmp ult i64 %17, 3
br i1 %19, label %22, label %20
20: ; preds = %14
%21 = and i64 %5, 4294967292
br label %48
22: ; preds = %48, %14
%23 = phi double [ undef, %14 ], [ %78, %48 ]
%24 = phi i64 [ 0, %14 ], [ %79, %48 ]
%25 = phi double [ %15, %14 ], [ %78, %48 ]
%26 = icmp eq i64 %18, 0
br i1 %26, label %40, label %27
27: ; preds = %22, %27
%28 = phi i64 [ %37, %27 ], [ %24, %22 ]
%29 = phi double [ %36, %27 ], [ %25, %22 ]
%30 = phi i64 [ %38, %27 ], [ 0, %22 ]
%31 = mul nsw i64 %28, %16
%32 = getelementptr inbounds double*, double** %1, i64 %31
%33 = getelementptr inbounds double*, double** %32, i64 %28
%34 = bitcast double** %33 to double*
%35 = load double, double* %34
%36 = fmul double %35, %29
%37 = add nuw nsw i64 %28, 1
%38 = add i64 %30, 1
%39 = icmp eq i64 %38, %18
br i1 %39, label %40, label %27
40: ; preds = %27, %22
%41 = phi double [ %23, %22 ], [ %36, %27 ]
store double %41, double* %4
br i1 %13, label %42, label %169
42: ; preds = %40
%43 = shl nuw nsw i64 %5, 3
%44 = and i64 %5, 3
%45 = icmp ult i64 %17, 3
%46 = and i64 %5, 4294967292
%47 = icmp eq i64 %44, 0
br label %89
48: ; preds = %48, %20
%49 = phi i64 [ 0, %20 ], [ %79, %48 ]
%50 = phi double [ %15, %20 ], [ %78, %48 ]
%51 = phi i64 [ 0, %20 ], [ %80, %48 ]
%52 = mul nsw i64 %49, %16
%53 = getelementptr inbounds double*, double** %1, i64 %52
%54 = getelementptr inbounds double*, double** %53, i64 %49
%55 = bitcast double** %54 to double*
%56 = load double, double* %55
%57 = fmul double %56, %50
%58 = or i64 %49, 1
%59 = mul nsw i64 %58, %16
%60 = getelementptr inbounds double*, double** %1, i64 %59
%61 = getelementptr inbounds double*, double** %60, i64 %58
%62 = bitcast double** %61 to double*
%63 = load double, double* %62
%64 = fmul double %63, %57
%65 = or i64 %49, 2
%66 = mul nsw i64 %65, %16
%67 = getelementptr inbounds double*, double** %1, i64 %66
%68 = getelementptr inbounds double*, double** %67, i64 %65
%69 = bitcast double** %68 to double*
%70 = load double, double* %69
%71 = fmul double %70, %64
%72 = or i64 %49, 3
%73 = mul nsw i64 %72, %16
%74 = getelementptr inbounds double*, double** %1, i64 %73
%75 = getelementptr inbounds double*, double** %74, i64 %72
%76 = bitcast double** %75 to double*
%77 = load double, double* %76
%78 = fmul double %77, %71
%79 = add nuw nsw i64 %49, 4
%80 = add i64 %51, 4
%81 = icmp eq i64 %80, %21
br i1 %81, label %22, label %48
82: ; preds = %131
br i1 %13, label %83, label %169
83: ; preds = %82
%84 = shl nuw nsw i64 %5, 3
%85 = and i64 %5, 1
%86 = icmp eq i64 %17, 0
br i1 %86, label %157, label %87
87: ; preds = %83
%88 = and i64 %5, 4294967294
br label %134
89: ; preds = %131, %42
%90 = phi i64 [ 0, %42 ], [ %132, %131 ]
%91 = getelementptr inbounds double, double* %8, i64 %90
store double 1.000000e+00, double* %91
%92 = call i32 @lubksb(double** %1, i32 %2, i32* %11, double* %8)
%93 = getelementptr inbounds double, double* %7, i64 %90
br i1 %45, label %119, label %94
94: ; preds = %89, %94
%95 = phi i64 [ %116, %94 ], [ 0, %89 ]
%96 = phi i64 [ %117, %94 ], [ 0, %89 ]
%97 = getelementptr inbounds double, double* %8, i64 %95
%98 = load double, double* %97
%99 = mul nuw nsw i64 %95, %5
%100 = getelementptr inbounds double, double* %93, i64 %99
store double %98, double* %100
%101 = or i64 %95, 1
%102 = getelementptr inbounds double, double* %8, i64 %101
%103 = load double, double* %102
%104 = mul nuw nsw i64 %101, %5
%105 = getelementptr inbounds double, double* %93, i64 %104
store double %103, double* %105
%106 = or i64 %95, 2
%107 = getelementptr inbounds double, double* %8, i64 %106
%108 = load double, double* %107
%109 = mul nuw nsw i64 %106, %5
%110 = getelementptr inbounds double, double* %93, i64 %109
store double %108, double* %110
%111 = or i64 %95, 3
%112 = getelementptr inbounds double, double* %8, i64 %111
%113 = load double, double* %112
%114 = mul nuw nsw i64 %111, %5
%115 = getelementptr inbounds double, double* %93, i64 %114
store double %113, double* %115
%116 = add nuw nsw i64 %95, 4
%117 = add i64 %96, 4
%118 = icmp eq i64 %117, %46
br i1 %118, label %119, label %94
119: ; preds = %94, %89
%120 = phi i64 [ 0, %89 ], [ %116, %94 ]
br i1 %47, label %131, label %121
121: ; preds = %119, %121
%122 = phi i64 [ %128, %121 ], [ %120, %119 ]
%123 = phi i64 [ %129, %121 ], [ 0, %119 ]
%124 = getelementptr inbounds double, double* %8, i64 %122
%125 = load double, double* %124
%126 = mul nuw nsw i64 %122, %5
%127 = getelementptr inbounds double, double* %93, i64 %126
store double %125, double* %127
%128 = add nuw nsw i64 %122, 1
%129 = add i64 %123, 1
%130 = icmp eq i64 %129, %44
br i1 %130, label %131, label %121
131: ; preds = %121, %119
%132 = add nuw nsw i64 %90, 1
%133 = icmp eq i64 %132, %5
br i1 %133, label %82, label %89
134: ; preds = %134, %87
%135 = phi i64 [ 0, %87 ], [ %154, %134 ]
%136 = phi i64 [ 0, %87 ], [ %155, %134 ]
%137 = mul nuw nsw i64 %135, %5
%138 = getelementptr double, double* %7, i64 %137
%139 = bitcast double* %138 to i8*
%140 = trunc i64 %135 to i32
%141 = mul i32 %140, %2
%142 = sext i32 %141 to i64
%143 = getelementptr double*, double** %0, i64 %142
%144 = bitcast double** %143 to i8*
%145 = or i64 %135, 1
%146 = mul nuw nsw i64 %145, %5
%147 = getelementptr double, double* %7, i64 %146
%148 = bitcast double* %147 to i8*
%149 = trunc i64 %145 to i32
%150 = mul i32 %149, %2
%151 = sext i32 %150 to i64
%152 = getelementptr double*, double** %0, i64 %151
%153 = bitcast double** %152 to i8*
%154 = add nuw nsw i64 %135, 2
%155 = add i64 %136, 2
%156 = icmp eq i64 %155, %88
br i1 %156, label %157, label %134
157: ; preds = %134, %83
%158 = phi i64 [ 0, %83 ], [ %154, %134 ]
%159 = icmp eq i64 %85, 0
br i1 %159, label %169, label %160
160: ; preds = %157
%161 = mul nuw nsw i64 %158, %5
%162 = getelementptr double, double* %7, i64 %161
%163 = bitcast double* %162 to i8*
%164 = trunc i64 %158 to i32
%165 = mul i32 %164, %2
%166 = sext i32 %165 to i64
%167 = getelementptr double*, double** %0, i64 %166
%168 = bitcast double** %167 to i8*
br label %169
169: ; preds = %160, %157, %3, %40, %82
%170 = load double, double* %4
ret double %170
}