#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 abs(int /*<<< orphan*/);
int delaunay_reduce_basis(double **, double const);
int /*<<< orphan*/ get_delaunay_shortest_vectors(double **, double const);
int /*<<< orphan*/ get_exteneded_basis(double **, double **);
double const mat_Dabs(double);
int /*<<< orphan*/ mat_cast_matrix_3d_to_3i(int **, double **);
int /*<<< orphan*/ mat_copy_matrix_d3(double **, double **);
double mat_get_determinant_d3(double **);
int /*<<< orphan*/ mat_get_determinant_i3(int **);
int /*<<< orphan*/ mat_inverse_matrix_d3(double **, double **, double const);
int /*<<< orphan*/ mat_multiply_matrix_d3(double **, double **, double **);
int delaunay_reduce(double red_lattice[3][3], double lattice[3][3],
                    const double symprec) {
  int i, j, attempt, succeeded;
  int tmp_mat_int[3][3];
  double volume;
  double orig_lattice[3][3], tmp_mat[3][3], basis[4][3];
  mat_copy_matrix_d3(orig_lattice, lattice);
  get_exteneded_basis(basis, lattice);
  for (attempt = 0; attempt < 100; attempt++) {
    succeeded = delaunay_reduce_basis(basis, symprec);
    if (succeeded) {
      break;
    }
  }
  if (!succeeded) {
    goto err;
  }
  get_delaunay_shortest_vectors(basis, symprec);
  for (i = 0; i < 3; i++) {
    for (j = 0; j < 3; j++) {
      red_lattice[i][j] = basis[j][i];
    }
  }
  volume = mat_get_determinant_d3(red_lattice);
  if (mat_Dabs(volume) < symprec) {

    ;
    goto err;
  }
  if (volume < 0) {
    for (i = 0; i < 3; i++) {
      for (j = 0; j < 3; j++) {
        red_lattice[i][j] = -red_lattice[i][j];
      }
    }
  }
  mat_inverse_matrix_d3(tmp_mat, red_lattice, symprec);
  mat_multiply_matrix_d3(tmp_mat, tmp_mat, orig_lattice);
  mat_cast_matrix_3d_to_3i(tmp_mat_int, tmp_mat);
  if (abs(mat_get_determinant_i3(tmp_mat_int)) != 1) {

    ;
    goto err;
  }
  return 1;
err:
  return 0;
}


x86 assembly (clang, -O0):
<delaunay_reduce>:
push   %rbp
mov    %rsp,%rbp
sub    $0x160,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
lea    -0xb0(%rbp),%rdi
mov    -0x18(%rbp),%rsi
call   10c0 <mat_copy_matrix_d3@plt>
lea    -0x160(%rbp),%rdi
mov    -0x18(%rbp),%rsi
call   1070 <get_exteneded_basis@plt>
movl   $0x0,-0x2c(%rbp)
cmpl   $0x64,-0x2c(%rbp)
jge    122f <delaunay_reduce+0x7f>
lea    -0x160(%rbp),%rdi
movsd  -0x20(%rbp),%xmm0
call   1060 <delaunay_reduce_basis@plt>
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x30(%rbp)
je     121c <delaunay_reduce+0x6c>
jmp    122f <delaunay_reduce+0x7f>
jmp    1221 <delaunay_reduce+0x71>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11ef <delaunay_reduce+0x3f>
cmpl   $0x0,-0x30(%rbp)
jne    123e <delaunay_reduce+0x8e>
jmp    13fa <delaunay_reduce+0x24a>
lea    -0x160(%rbp),%rdi
movsd  -0x20(%rbp),%xmm0
call   1040 <get_delaunay_shortest_vectors@plt>
movl   $0x0,-0x24(%rbp)
cmpl   $0x3,-0x24(%rbp)
jge    12c5 <delaunay_reduce+0x115>
movl   $0x0,-0x28(%rbp)
cmpl   $0x3,-0x28(%rbp)
jge    12b2 <delaunay_reduce+0x102>
movslq -0x28(%rbp),%rcx
lea    -0x160(%rbp),%rax
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1267 <delaunay_reduce+0xb7>
jmp    12b7 <delaunay_reduce+0x107>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1256 <delaunay_reduce+0xa6>
mov    -0x10(%rbp),%rdi
call   1050 <mat_get_determinant_d3@plt>
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
call   10b0 <mat_Dabs@plt>
movaps %xmm0,%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    12f4 <delaunay_reduce+0x144>
jmp    13fa <delaunay_reduce+0x24a>
xorps  %xmm0,%xmm0
ucomisd -0x68(%rbp),%xmm0
jbe    1391 <delaunay_reduce+0x1e1>
movl   $0x0,-0x24(%rbp)
cmpl   $0x3,-0x24(%rbp)
jge    138c <delaunay_reduce+0x1dc>
movl   $0x0,-0x28(%rbp)
cmpl   $0x3,-0x28(%rbp)
jge    1379 <delaunay_reduce+0x1c9>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    131a <delaunay_reduce+0x16a>
jmp    137e <delaunay_reduce+0x1ce>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1309 <delaunay_reduce+0x159>
jmp    1391 <delaunay_reduce+0x1e1>
lea    -0x100(%rbp),%rdi
mov    -0x10(%rbp),%rsi
movsd  -0x20(%rbp),%xmm0
call   1080 <mat_inverse_matrix_d3@plt>
lea    -0x100(%rbp),%rdi
lea    -0x100(%rbp),%rsi
lea    -0xb0(%rbp),%rdx
call   10d0 <mat_multiply_matrix_d3@plt>
lea    -0x60(%rbp),%rdi
lea    -0x100(%rbp),%rsi
call   10a0 <mat_cast_matrix_3d_to_3i@plt>
lea    -0x60(%rbp),%rdi
call   1090 <mat_get_determinant_i3@plt>
mov    %eax,%edi
call   1030 <abs@plt>
cmp    $0x1,%eax
je     13ee <delaunay_reduce+0x23e>
jmp    13fa <delaunay_reduce+0x24a>
movl   $0x1,-0x4(%rbp)
jmp    1401 <delaunay_reduce+0x251>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x160,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<delaunay_reduce>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x138,%rsp
movsd  %xmm0,(%rsp)
mov    %rsi,%r15
mov    %rdi,%r14
lea    0x70(%rsp),%rdi
call   10b0 <mat_copy_matrix_d3@plt>
lea    0x10(%rsp),%rbx
mov    %rbx,%rdi
mov    %r15,%rsi
call   1060 <get_exteneded_basis@plt>
mov    $0x64,%ebp
nopw   0x0(%rax,%rax,1)
mov    %rbx,%rdi
movsd  (%rsp),%xmm0
call   1050 <delaunay_reduce_basis@plt>
test   %eax,%eax
jne    11fb <delaunay_reduce+0x5b>
add    $0xffffffff,%ebp
jne    11e0 <delaunay_reduce+0x40>
jmp    134d <delaunay_reduce+0x1ad>
lea    0x10(%rsp),%rdi
movsd  (%rsp),%xmm0
call   1030 <get_delaunay_shortest_vectors@plt>
movsd  0x10(%rsp),%xmm0
movsd  %xmm0,(%r14)
movsd  0x28(%rsp),%xmm0
movsd  %xmm0,0x8(%r14)
movsd  0x40(%rsp),%xmm0
movsd  %xmm0,0x10(%r14)
movsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x18(%r14)
movsd  0x30(%rsp),%xmm0
movsd  %xmm0,0x20(%r14)
movsd  0x48(%rsp),%xmm0
movsd  %xmm0,0x28(%r14)
movsd  0x20(%rsp),%xmm0
movsd  %xmm0,0x30(%r14)
movsd  0x38(%rsp),%xmm0
movsd  %xmm0,0x38(%r14)
movsd  0x50(%rsp),%xmm0
movsd  %xmm0,0x40(%r14)
mov    %r14,%rdi
call   1040 <mat_get_determinant_d3@plt>
movsd  %xmm0,0x8(%rsp)
call   10a0 <mat_Dabs@plt>
movsd  (%rsp),%xmm1
ucomisd %xmm0,%xmm1
ja     134d <delaunay_reduce+0x1ad>
xorpd  %xmm0,%xmm0
ucomisd 0x8(%rsp),%xmm0
jbe    12f9 <delaunay_reduce+0x159>
movapd 0xd55(%rip),%xmm0
movupd (%r14),%xmm1
movupd 0x10(%r14),%xmm2
movupd 0x20(%r14),%xmm3
movupd 0x30(%r14),%xmm4
xorpd  %xmm0,%xmm1
movupd %xmm1,(%r14)
xorpd  %xmm0,%xmm2
movupd %xmm2,0x10(%r14)
xorpd  %xmm0,%xmm3
movupd %xmm3,0x20(%r14)
xorpd  %xmm0,%xmm4
movupd %xmm4,0x30(%r14)
movsd  0x40(%r14),%xmm1
xorpd  %xmm0,%xmm1
movlpd %xmm1,0x40(%r14)
lea    0xf0(%rsp),%rbx
mov    %rbx,%rdi
mov    %r14,%rsi
movsd  (%rsp),%xmm0
call   1070 <mat_inverse_matrix_d3@plt>
lea    0x70(%rsp),%rdx
mov    %rbx,%rdi
mov    %rbx,%rsi
call   10c0 <mat_multiply_matrix_d3@plt>
lea    0xc0(%rsp),%r14
mov    %r14,%rdi
mov    %rbx,%rsi
call   1090 <mat_cast_matrix_3d_to_3i@plt>
mov    %r14,%rdi
call   1080 <mat_get_determinant_i3@plt>
mov    %eax,%ecx
neg    %ecx
cmovs  %eax,%ecx
mov    $0x1,%eax
cmp    $0x1,%ecx
je     134f <delaunay_reduce+0x1af>
xor    %eax,%eax
add    $0x138,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @delaunay_reduce([3 x double]* %0, [3 x double]* %1, double %2) {
%4 = alloca i32
%5 = alloca [3 x double]*
%6 = alloca [3 x double]*
%7 = alloca double
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca [3 x [3 x i32]]
%13 = alloca double
%14 = alloca [3 x [3 x double]]
%15 = alloca [3 x [3 x double]]
%16 = alloca [4 x [3 x double]]
store [3 x double]* %0, [3 x double]** %5
store [3 x double]* %1, [3 x double]** %6
store double %2, double* %7
%17 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %14, i64 0, i64 0
%18 = bitcast [3 x double]* %17 to double**
%19 = load [3 x double]*, [3 x double]** %6
%20 = bitcast [3 x double]* %19 to double**
%21 = call i32 @mat_copy_matrix_d3(double** %18, double** %20)
%22 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %16, i64 0, i64 0
%23 = bitcast [3 x double]* %22 to double**
%24 = load [3 x double]*, [3 x double]** %6
%25 = bitcast [3 x double]* %24 to double**
%26 = call i32 @get_exteneded_basis(double** %23, double** %25)
store i32 0, i32* %10
br label %27
27: ; preds = %39, %3
%28 = load i32, i32* %10
%29 = icmp slt i32 %28, 100
br i1 %29, label %30, label %42
30: ; preds = %27
%31 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %16, i64 0, i64 0
%32 = bitcast [3 x double]* %31 to double**
%33 = load double, double* %7
%34 = call i32 @delaunay_reduce_basis(double** %32, double %33)
store i32 %34, i32* %11
%35 = load i32, i32* %11
%36 = icmp ne i32 %35, 0
br i1 %36, label %37, label %38
37: ; preds = %30
br label %42
38: ; preds = %30
br label %39
39: ; preds = %38
%40 = load i32, i32* %10
%41 = add nsw i32 %40, 1
store i32 %41, i32* %10
br label %27
42: ; preds = %37, %27
%43 = load i32, i32* %11
%44 = icmp ne i32 %43, 0
br i1 %44, label %46, label %45
45: ; preds = %42
br label %151
46: ; preds = %42
%47 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %16, i64 0, i64 0
%48 = bitcast [3 x double]* %47 to double**
%49 = load double, double* %7
%50 = call i32 @get_delaunay_shortest_vectors(double** %48, double %49)
store i32 0, i32* %8
br label %51
51: ; preds = %77, %46
%52 = load i32, i32* %8
%53 = icmp slt i32 %52, 3
br i1 %53, label %54, label %80
54: ; preds = %51
store i32 0, i32* %9
br label %55
55: ; preds = %73, %54
%56 = load i32, i32* %9
%57 = icmp slt i32 %56, 3
br i1 %57, label %58, label %76
58: ; preds = %55
%59 = load i32, i32* %9
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %16, i64 0, i64 %60
%62 = load i32, i32* %8
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds [3 x double], [3 x double]* %61, i64 0, i64 %63
%65 = load double, double* %64
%66 = load [3 x double]*, [3 x double]** %5
%67 = load i32, i32* %8
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds [3 x double], [3 x double]* %66, i64 %68
%70 = load i32, i32* %9
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds [3 x double], [3 x double]* %69, i64 0, i64 %71
store double %65, double* %72
br label %73
73: ; preds = %58
%74 = load i32, i32* %9
%75 = add nsw i32 %74, 1
store i32 %75, i32* %9
br label %55
76: ; preds = %55
br label %77
77: ; preds = %76
%78 = load i32, i32* %8
%79 = add nsw i32 %78, 1
store i32 %79, i32* %8
br label %51
80: ; preds = %51
%81 = load [3 x double]*, [3 x double]** %5
%82 = bitcast [3 x double]* %81 to double**
%83 = call double @mat_get_determinant_d3(double** %82)
store double %83, double* %13
%84 = load double, double* %13
%85 = call double @mat_Dabs(double %84)
%86 = load double, double* %7
%87 = fcmp olt double %85, %86
br i1 %87, label %88, label %89
88: ; preds = %80
br label %151
89: ; preds = %80
%90 = load double, double* %13
%91 = fcmp olt double %90, 0.000000e+00
br i1 %91, label %92, label %125
92: ; preds = %89
store i32 0, i32* %8
br label %93
93: ; preds = %121, %92
%94 = load i32, i32* %8
%95 = icmp slt i32 %94, 3
br i1 %95, label %96, label %124
96: ; preds = %93
store i32 0, i32* %9
br label %97
97: ; preds = %117, %96
%98 = load i32, i32* %9
%99 = icmp slt i32 %98, 3
br i1 %99, label %100, label %120
100: ; preds = %97
%101 = load [3 x double]*, [3 x double]** %5
%102 = load i32, i32* %8
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds [3 x double], [3 x double]* %101, i64 %103
%105 = load i32, i32* %9
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds [3 x double], [3 x double]* %104, i64 0, i64 %106
%108 = load double, double* %107
%109 = fneg double %108
%110 = load [3 x double]*, [3 x double]** %5
%111 = load i32, i32* %8
%112 = sext i32 %111 to i64
%113 = getelementptr inbounds [3 x double], [3 x double]* %110, i64 %112
%114 = load i32, i32* %9
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds [3 x double], [3 x double]* %113, i64 0, i64 %115
store double %109, double* %116
br label %117
117: ; preds = %100
%118 = load i32, i32* %9
%119 = add nsw i32 %118, 1
store i32 %119, i32* %9
br label %97
120: ; preds = %97
br label %121
121: ; preds = %120
%122 = load i32, i32* %8
%123 = add nsw i32 %122, 1
store i32 %123, i32* %8
br label %93
124: ; preds = %93
br label %125
125: ; preds = %124, %89
%126 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %15, i64 0, i64 0
%127 = bitcast [3 x double]* %126 to double**
%128 = load [3 x double]*, [3 x double]** %5
%129 = bitcast [3 x double]* %128 to double**
%130 = load double, double* %7
%131 = call i32 @mat_inverse_matrix_d3(double** %127, double** %129, double %130)
%132 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %15, i64 0, i64 0
%133 = bitcast [3 x double]* %132 to double**
%134 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %15, i64 0, i64 0
%135 = bitcast [3 x double]* %134 to double**
%136 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %14, i64 0, i64 0
%137 = bitcast [3 x double]* %136 to double**
%138 = call i32 @mat_multiply_matrix_d3(double** %133, double** %135, double** %137)
%139 = getelementptr inbounds [3 x [3 x i32]], [3 x [3 x i32]]* %12, i64 0, i64 0
%140 = bitcast [3 x i32]* %139 to i32**
%141 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %15, i64 0, i64 0
%142 = bitcast [3 x double]* %141 to double**
%143 = call i32 @mat_cast_matrix_3d_to_3i(i32** %140, double** %142)
%144 = getelementptr inbounds [3 x [3 x i32]], [3 x [3 x i32]]* %12, i64 0, i64 0
%145 = bitcast [3 x i32]* %144 to i32**
%146 = call i32 @mat_get_determinant_i3(i32** %145)
%147 = call i32 @abs(i32 %146)
%148 = icmp ne i32 %147, 1
br i1 %148, label %149, label %150
149: ; preds = %125
br label %151
150: ; preds = %125
store i32 1, i32* %4
br label %152
151: ; preds = %149, %88, %45
store i32 0, i32* %4
br label %152
152: ; preds = %151, %150
%153 = load i32, i32* %4
ret i32 %153
}

LLVM IR (-O2):
define i32 @delaunay_reduce([3 x double]* %0, [3 x double]* %1, double %2) {
%4 = alloca [3 x [3 x i32]]
%5 = alloca [3 x [3 x double]]
%6 = alloca [3 x [3 x double]]
%7 = alloca [4 x [3 x double]]
%8 = bitcast [3 x [3 x i32]]* %4 to i8*
%9 = bitcast [3 x [3 x double]]* %5 to i8*
%10 = bitcast [3 x [3 x double]]* %6 to i8*
%11 = bitcast [4 x [3 x double]]* %7 to i8*
%12 = bitcast [3 x [3 x double]]* %5 to double**
%13 = bitcast [3 x double]* %1 to double**
%14 = call i32 @mat_copy_matrix_d3(double** %12, double** %13)
%15 = bitcast [4 x [3 x double]]* %7 to double**
%16 = call i32 @get_exteneded_basis(double** %15, double** %13)
br label %20
17: ; preds = %20
%18 = add nuw nsw i32 %21, 1
%19 = icmp eq i32 %18, 100
br i1 %19, label %87, label %20
20: ; preds = %3, %17
%21 = phi i32 [ 0, %3 ], [ %18, %17 ]
%22 = call i32 @delaunay_reduce_basis(double** %15, double %2)
%23 = icmp eq i32 %22, 0
br i1 %23, label %17, label %24
24: ; preds = %20
%25 = call i32 @get_delaunay_shortest_vectors(double** %15, double %2)
%26 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %7, i64 0, i64 0, i64 0
%27 = load double, double* %26
%28 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 0, i64 0
store double %27, double* %28
%29 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %7, i64 0, i64 1, i64 0
%30 = load double, double* %29
%31 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 0, i64 1
store double %30, double* %31
%32 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %7, i64 0, i64 2, i64 0
%33 = load double, double* %32
%34 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 0, i64 2
store double %33, double* %34
%35 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %7, i64 0, i64 0, i64 1
%36 = load double, double* %35
%37 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 1, i64 0
store double %36, double* %37
%38 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %7, i64 0, i64 1, i64 1
%39 = load double, double* %38
%40 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 1, i64 1
store double %39, double* %40
%41 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %7, i64 0, i64 2, i64 1
%42 = load double, double* %41
%43 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 1, i64 2
store double %42, double* %43
%44 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %7, i64 0, i64 0, i64 2
%45 = load double, double* %44
%46 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 2, i64 0
store double %45, double* %46
%47 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %7, i64 0, i64 1, i64 2
%48 = load double, double* %47
%49 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 2, i64 1
store double %48, double* %49
%50 = getelementptr inbounds [4 x [3 x double]], [4 x [3 x double]]* %7, i64 0, i64 2, i64 2
%51 = load double, double* %50
%52 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 2, i64 2
store double %51, double* %52
%53 = bitcast [3 x double]* %0 to double**
%54 = call double @mat_get_determinant_d3(double** %53)
%55 = call double @mat_Dabs(double %54)
%56 = fcmp olt double %55, %2
br i1 %56, label %87, label %57
57: ; preds = %24
%58 = fcmp olt double %54, 0.000000e+00
br i1 %58, label %59, label %78
59: ; preds = %57
%60 = bitcast [3 x double]* %0 to <2 x double>*
%61 = load <2 x double>, <2 x double>* %60
%62 = fneg <2 x double> %61
%63 = bitcast [3 x double]* %0 to <2 x double>*
store <2 x double> %62, <2 x double>* %63
%64 = bitcast double* %34 to <2 x double>*
%65 = load <2 x double>, <2 x double>* %64
%66 = fneg <2 x double> %65
%67 = bitcast double* %34 to <2 x double>*
store <2 x double> %66, <2 x double>* %67
%68 = bitcast double* %40 to <2 x double>*
%69 = load <2 x double>, <2 x double>* %68
%70 = fneg <2 x double> %69
%71 = bitcast double* %40 to <2 x double>*
store <2 x double> %70, <2 x double>* %71
%72 = bitcast double* %46 to <2 x double>*
%73 = load <2 x double>, <2 x double>* %72
%74 = fneg <2 x double> %73
%75 = bitcast double* %46 to <2 x double>*
store <2 x double> %74, <2 x double>* %75
%76 = load double, double* %52
%77 = fneg double %76
store double %77, double* %52
br label %78
78: ; preds = %59, %57
%79 = bitcast [3 x [3 x double]]* %6 to double**
%80 = call i32 @mat_inverse_matrix_d3(double** %79, double** %53, double %2)
%81 = call i32 @mat_multiply_matrix_d3(double** %79, double** %79, double** %12)
%82 = bitcast [3 x [3 x i32]]* %4 to i32**
%83 = call i32 @mat_cast_matrix_3d_to_3i(i32** %82, double** %79)
%84 = call i32 @mat_get_determinant_i3(i32** %82)
%85 = call i32 @llvm.abs.i32(i32 %84, i1 true)
%86 = icmp eq i32 %85, 1
br i1 %86, label %88, label %87
87: ; preds = %17, %78, %24
br label %88
88: ; preds = %78, %87
%89 = phi i32 [ 0, %87 ], [ 1, %78 ]
ret i32 %89
}