#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 */
double **dmatrix(int /*<<< orphan*/, int, int /*<<< orphan*/, int);
double *dvector(int /*<<< orphan*/, int);
int /*<<< orphan*/ free_dmatrix(double **, int /*<<< orphan*/, int,
                                int /*<<< orphan*/, int);
int /*<<< orphan*/ free_dvector(double *, int /*<<< orphan*/, int);
int /*<<< orphan*/ free_ivector(int *, int /*<<< orphan*/, int);
int *ivector(int /*<<< orphan*/, int);
int /*<<< orphan*/ lu_backsub(double **, int, int *, double *);
scalar_t__ lu_decomp(double **, int, int *, double *);
int /*<<< orphan*/ lu_polish(double **, double **, int, double *, double *,
                             int *);
int polished_solve_se(double **a, double *b, int n) {
  double rip;
  int *pivx, PIVX[10];
  double **sa;
  double *sb;
  int i, j;
  if (n <= 10)
    pivx = PIVX;
  else
    pivx = ivector(0, n - 1);
  sa = dmatrix(0, n - 1, 0, n - 1);
  sb = dvector(0, n - 1);
  for (i = 0; i < n; i++) {
    sb[i] = b[i];
    for (j = 0; j < n; j++)
      sa[i][j] = a[i][j];
  }
  if (lu_decomp(a, n, pivx, &rip)) {
    free_dvector(sb, 0, n - 1);
    free_dmatrix(sa, 0, n - 1, 0, n - 1);
    if (pivx != PIVX)
      free_ivector(pivx, 0, n - 1);
    return 1;
  }
  lu_backsub(a, n, pivx, b);
  lu_polish(sa, a, n, sb, b, pivx);
  free_dvector(sb, 0, n - 1);
  free_dmatrix(sa, 0, n - 1, 0, n - 1);
  if (pivx != PIVX)
    free_ivector(pivx, 0, n - 1);
  return 0;
}


x86 assembly (clang, -O0):
<polished_solve_se>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
cmpl   $0xa,-0x1c(%rbp)
jg     11bd <polished_solve_se+0x2d>
lea    -0x60(%rbp),%rax
mov    %rax,-0x30(%rbp)
jmp    11ce <polished_solve_se+0x3e>
mov    -0x1c(%rbp),%esi
sub    $0x1,%esi
xor    %edi,%edi
call   1090 <ivector@plt>
mov    %rax,-0x30(%rbp)
mov    -0x1c(%rbp),%esi
sub    $0x1,%esi
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
xor    %edx,%edx
mov    %edx,%edi
call   1030 <dmatrix@plt>
mov    %rax,-0x68(%rbp)
mov    -0x1c(%rbp),%esi
sub    $0x1,%esi
xor    %edi,%edi
call   1050 <dvector@plt>
mov    %rax,-0x70(%rbp)
movl   $0x0,-0x74(%rbp)
mov    -0x74(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1283 <polished_solve_se+0xf3>
mov    -0x18(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x70(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x78(%rbp)
mov    -0x78(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1270 <polished_solve_se+0xe0>
mov    -0x10(%rbp),%rax
movslq -0x74(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x78(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x68(%rbp),%rax
movslq -0x74(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x78(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x78(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x78(%rbp)
jmp    122c <polished_solve_se+0x9c>
jmp    1275 <polished_solve_se+0xe5>
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    11ff <polished_solve_se+0x6f>
mov    -0x10(%rbp),%rdi
mov    -0x1c(%rbp),%esi
mov    -0x30(%rbp),%rdx
lea    -0x28(%rbp),%rcx
call   1070 <lu_decomp@plt>
cmp    $0x0,%rax
je     12fb <polished_solve_se+0x16b>
mov    -0x70(%rbp),%rdi
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
xor    %esi,%esi
call   1080 <free_dvector@plt>
mov    -0x68(%rbp),%rdi
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
mov    -0x1c(%rbp),%r8d
sub    $0x1,%r8d
xor    %ecx,%ecx
mov    %ecx,%esi
call   10a0 <free_dmatrix@plt>
mov    -0x30(%rbp),%rax
lea    -0x60(%rbp),%rcx
cmp    %rcx,%rax
je     12ef <polished_solve_se+0x15f>
mov    -0x30(%rbp),%rdi
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
xor    %esi,%esi
call   1060 <free_ivector@plt>
movl   $0x1,-0x4(%rbp)
jmp    1380 <polished_solve_se+0x1f0>
mov    -0x10(%rbp),%rdi
mov    -0x1c(%rbp),%esi
mov    -0x30(%rbp),%rdx
mov    -0x18(%rbp),%rcx
call   1040 <lu_backsub@plt>
mov    -0x68(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x1c(%rbp),%edx
mov    -0x70(%rbp),%rcx
mov    -0x18(%rbp),%r8
mov    -0x30(%rbp),%r9
call   10b0 <lu_polish@plt>
mov    -0x70(%rbp),%rdi
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
xor    %esi,%esi
call   1080 <free_dvector@plt>
mov    -0x68(%rbp),%rdi
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
mov    -0x1c(%rbp),%r8d
sub    $0x1,%r8d
xor    %ecx,%ecx
mov    %ecx,%esi
call   10a0 <free_dmatrix@plt>
mov    -0x30(%rbp),%rax
lea    -0x60(%rbp),%rcx
cmp    %rcx,%rax
je     1379 <polished_solve_se+0x1e9>
mov    -0x30(%rbp),%rdi
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
xor    %esi,%esi
call   1060 <free_ivector@plt>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<polished_solve_se>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %edx,%r13d
mov    %rsi,0x18(%rsp)
mov    %rdi,0x8(%rsp)
cmp    $0xa,%edx
jg     11c0 <polished_solve_se+0x30>
lea    0x40(%rsp),%rax
mov    %rax,0x10(%rsp)
lea    -0x1(%r13),%ebx
jmp    11d2 <polished_solve_se+0x42>
lea    -0x1(%r13),%ebx
xor    %edi,%edi
mov    %ebx,%esi
call   1090 <ivector@plt>
mov    %rax,0x10(%rsp)
xor    %edi,%edi
mov    %ebx,%esi
xor    %edx,%edx
mov    %ebx,%ecx
call   1030 <dmatrix@plt>
mov    %rax,%r12
xor    %edi,%edi
mov    %ebx,0x4(%rsp)
mov    %ebx,%esi
call   1050 <dvector@plt>
mov    %rax,%r14
test   %r13d,%r13d
mov    0x8(%rsp),%rdi
jle    13f4 <polished_solve_se+0x264>
mov    %r13d,%eax
mov    %eax,%r8d
and    $0xfffffffc,%r8d
lea    -0x4(%r8),%rcx
mov    %rcx,0x30(%rsp)
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %ecx,%r10d
and    $0x3,%r10d
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0x20(%rsp)
mov    %eax,%r9d
and    $0x3,%r9d
mov    %r10,0x28(%rsp)
shl    $0x5,%r10
xor    %esi,%esi
jmp    124d <polished_solve_se+0xbd>
nop
add    $0x1,%rsi
cmp    %rax,%rsi
je     13f4 <polished_solve_se+0x264>
mov    0x18(%rsp),%rcx
movsd  (%rcx,%rsi,8),%xmm0
movsd  %xmm0,(%r14,%rsi,8)
mov    (%rdi,%rsi,8),%rcx
mov    (%r12,%rsi,8),%r15
cmp    $0x4,%r13d
jae    1280 <polished_solve_se+0xf0>
xor    %r11d,%r11d
jmp    1370 <polished_solve_se+0x1e0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
lea    (%rcx,%rax,8),%rdx
cmp    %rdx,%r15
jae    12a0 <polished_solve_se+0x110>
lea    (%r15,%rax,8),%rdx
cmp    %rdx,%rcx
jae    12a0 <polished_solve_se+0x110>
xor    %r11d,%r11d
mov    0x8(%rsp),%rdi
jmp    1370 <polished_solve_se+0x1e0>
nop
cmpq   $0xc,0x30(%rsp)
jae    12ac <polished_solve_se+0x11c>
xor    %ebx,%ebx
jmp    1320 <polished_solve_se+0x190>
mov    0x20(%rsp),%rbp
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movups (%rcx,%rbx,8),%xmm0
movups 0x10(%rcx,%rbx,8),%xmm1
movups %xmm0,(%r15,%rbx,8)
movups %xmm1,0x10(%r15,%rbx,8)
movups 0x20(%rcx,%rbx,8),%xmm0
movups 0x30(%rcx,%rbx,8),%xmm1
movups %xmm0,0x20(%r15,%rbx,8)
movups %xmm1,0x30(%r15,%rbx,8)
movups 0x40(%rcx,%rbx,8),%xmm0
movups 0x50(%rcx,%rbx,8),%xmm1
movups %xmm0,0x40(%r15,%rbx,8)
movups %xmm1,0x50(%r15,%rbx,8)
movups 0x60(%rcx,%rbx,8),%xmm0
movups 0x70(%rcx,%rbx,8),%xmm1
movups %xmm0,0x60(%r15,%rbx,8)
movups %xmm1,0x70(%r15,%rbx,8)
add    $0x10,%rbx
add    $0xfffffffffffffffc,%rbp
jne    12c0 <polished_solve_se+0x130>
cmpq   $0x0,0x28(%rsp)
je     135c <polished_solve_se+0x1cc>
lea    (%r15,%rbx,8),%rbp
add    $0x10,%rbp
lea    (%rcx,%rbx,8),%rbx
add    $0x10,%rbx
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
movups -0x10(%rbx,%rdx,1),%xmm0
movups (%rbx,%rdx,1),%xmm1
movups %xmm0,-0x10(%rbp,%rdx,1)
movups %xmm1,0x0(%rbp,%rdx,1)
add    $0x20,%rdx
cmp    %rdx,%r10
jne    1340 <polished_solve_se+0x1b0>
mov    %r8,%r11
cmp    %rax,%r8
mov    0x8(%rsp),%rdi
je     1240 <polished_solve_se+0xb0>
nopl   (%rax)
mov    %r11,%rbx
not    %rbx
add    %rax,%rbx
test   %r9,%r9
je     13a6 <polished_solve_se+0x216>
mov    %r9,%rbp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%r11,8),%xmm0
movsd  %xmm0,(%r15,%r11,8)
add    $0x1,%r11
add    $0xffffffffffffffff,%rbp
jne    1390 <polished_solve_se+0x200>
cmp    $0x3,%rbx
jb     1240 <polished_solve_se+0xb0>
movsd  (%rcx,%r11,8),%xmm0
movsd  %xmm0,(%r15,%r11,8)
movsd  0x8(%rcx,%r11,8),%xmm0
movsd  %xmm0,0x8(%r15,%r11,8)
movsd  0x10(%rcx,%r11,8),%xmm0
movsd  %xmm0,0x10(%r15,%r11,8)
movsd  0x18(%rcx,%r11,8),%xmm0
movsd  %xmm0,0x18(%r15,%r11,8)
add    $0x4,%r11
cmp    %r11,%rax
jne    13b0 <polished_solve_se+0x220>
jmp    1240 <polished_solve_se+0xb0>
lea    0x38(%rsp),%rcx
mov    %r13d,%esi
mov    0x10(%rsp),%rbp
mov    %rdi,%rbx
mov    %rbp,%rdx
call   1070 <lu_decomp@plt>
test   %rax,%rax
je     143a <polished_solve_se+0x2aa>
mov    %r14,%rdi
xor    %esi,%esi
mov    0x4(%rsp),%ebx
mov    %ebx,%edx
call   1080 <free_dvector@plt>
mov    %r12,%rdi
xor    %esi,%esi
mov    %ebx,%edx
xor    %ecx,%ecx
mov    %ebx,%r8d
call   10a0 <free_dmatrix@plt>
mov    $0x1,%r15d
jmp    148b <polished_solve_se+0x2fb>
mov    %rbx,%rdi
mov    %r13d,%esi
mov    %rbp,%rdx
mov    0x18(%rsp),%r15
mov    %r15,%rcx
call   1040 <lu_backsub@plt>
mov    %r12,%rdi
mov    %rbx,%rsi
mov    %r13d,%edx
mov    %r14,%rcx
mov    %r15,%r8
mov    %rbp,%r9
call   10b0 <lu_polish@plt>
xor    %r15d,%r15d
mov    %r14,%rdi
xor    %esi,%esi
mov    0x4(%rsp),%ebx
mov    %ebx,%edx
call   1080 <free_dvector@plt>
mov    %r12,%rdi
xor    %esi,%esi
mov    %ebx,%edx
xor    %ecx,%ecx
mov    %ebx,%r8d
call   10a0 <free_dmatrix@plt>
lea    0x40(%rsp),%rax
cmp    %rax,%rbp
je     14a1 <polished_solve_se+0x311>
mov    %rbp,%rdi
xor    %esi,%esi
mov    %ebx,%edx
call   1060 <free_ivector@plt>
mov    %r15d,%eax
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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