#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 const mat_Dabs(double);
int /*<<< orphan*/ mat_copy_matrix_d3(double **, double **);
double mat_get_determinant_d3(double **);
int /*<<< orphan*/ warning_print(char *, double);
int mat_inverse_matrix_d3(double m[3][3], double a[3][3],
                          const double precision) {
  double det;
  double c[3][3];
  det = mat_get_determinant_d3(a);
  if (mat_Dabs(det) < precision) {
    warning_print("spglib: No inverse matrix (det=%f)\n", det);
    return 0;
  }
  c[0][0] = (a[1][1] * a[2][2] - a[1][2] * a[2][1]) / det;
  c[1][0] = (a[1][2] * a[2][0] - a[1][0] * a[2][2]) / det;
  c[2][0] = (a[1][0] * a[2][1] - a[1][1] * a[2][0]) / det;
  c[0][1] = (a[2][1] * a[0][2] - a[2][2] * a[0][1]) / det;
  c[1][1] = (a[2][2] * a[0][0] - a[2][0] * a[0][2]) / det;
  c[2][1] = (a[2][0] * a[0][1] - a[2][1] * a[0][0]) / det;
  c[0][2] = (a[0][1] * a[1][2] - a[0][2] * a[1][1]) / det;
  c[1][2] = (a[0][2] * a[1][0] - a[0][0] * a[1][2]) / det;
  c[2][2] = (a[0][0] * a[1][1] - a[0][1] * a[1][0]) / det;
  mat_copy_matrix_d3(m, c);
  return 1;
}


x86 assembly (clang, -O0):
<mat_inverse_matrix_d3>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%rdi
call   1030 <mat_get_determinant_d3@plt>
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1050 <mat_Dabs@plt>
movaps %xmm0,%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    119c <mat_inverse_matrix_d3+0x5c>
movsd  -0x28(%rbp),%xmm0
lea    0xe75(%rip),%rdi
call   1040 <warning_print@plt>
movl   $0x0,-0x4(%rbp)
jmp    1461 <mat_inverse_matrix_d3+0x321>
mov    -0x18(%rbp),%rax
movsd  0x20(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x40(%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x38(%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x30(%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x40(%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x38(%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x20(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x30(%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x38(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x40(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x8(%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x40(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x30(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x10(%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x30(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x38(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  (%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x20(%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x28(%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x20(%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x18(%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rdi
lea    -0x70(%rbp),%rsi
call   1060 <mat_copy_matrix_d3@plt>
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mat_inverse_matrix_d3>:
push   %r14
push   %rbx
sub    $0x68,%rsp
movsd  %xmm0,0x8(%rsp)
mov    %rsi,%rbx
mov    %rdi,%r14
mov    %rsi,%rdi
call   1030 <mat_get_determinant_d3@plt>
movapd %xmm0,0x10(%rsp)
call   1050 <mat_Dabs@plt>
movsd  0x8(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jbe    118b <mat_inverse_matrix_d3+0x4b>
lea    0xe87(%rip),%rdi
movapd 0x10(%rsp),%xmm0
call   1040 <warning_print@plt>
xor    %eax,%eax
jmp    129d <mat_inverse_matrix_d3+0x15d>
movsd  0x20(%rbx),%xmm8
movupd 0x28(%rbx),%xmm3
movupd 0x10(%rbx),%xmm4
movapd %xmm4,%xmm10
unpckhpd %xmm4,%xmm10
movsd  0x8(%rbx),%xmm7
movsd  0x40(%rbx),%xmm9
movsd  0x38(%rbx),%xmm6
movapd %xmm8,%xmm12
unpcklpd %xmm6,%xmm12
unpcklpd %xmm7,%xmm6
movapd %xmm3,%xmm11
movapd %xmm4,%xmm1
unpckhpd %xmm3,%xmm1
mulsd  %xmm7,%xmm10
unpcklpd %xmm3,%xmm7
mulpd  %xmm3,%xmm7
movapd %xmm4,%xmm2
unpcklpd %xmm3,%xmm2
unpcklpd %xmm9,%xmm3
mulpd  %xmm6,%xmm3
movapd %xmm8,%xmm0
unpcklpd %xmm9,%xmm0
unpcklpd %xmm4,%xmm9
movapd %xmm12,%xmm5
mulpd  %xmm9,%xmm5
subpd  %xmm3,%xmm5
movapd 0x10(%rsp),%xmm13
movapd %xmm13,%xmm3
unpcklpd %xmm13,%xmm3
divpd  %xmm3,%xmm5
movapd %xmm5,0x20(%rsp)
movsd  (%rbx),%xmm5
shufpd $0x1,%xmm5,%xmm11
mulpd  %xmm11,%xmm12
mulpd  %xmm6,%xmm1
subpd  %xmm12,%xmm1
divpd  %xmm3,%xmm1
movapd %xmm1,0x50(%rsp)
mulpd  %xmm4,%xmm0
subpd  %xmm0,%xmm7
divpd  %xmm3,%xmm7
movapd %xmm7,0x30(%rsp)
mulpd  %xmm11,%xmm2
movsd  %xmm5,%xmm4
mulpd  %xmm9,%xmm4
subpd  %xmm2,%xmm4
divpd  %xmm3,%xmm4
movapd %xmm4,0x40(%rsp)
mulsd  %xmm8,%xmm5
subsd  %xmm10,%xmm5
divsd  %xmm13,%xmm5
movsd  %xmm5,0x60(%rsp)
lea    0x20(%rsp),%rsi
mov    %r14,%rdi
call   1060 <mat_copy_matrix_d3@plt>
mov    $0x1,%eax
add    $0x68,%rsp
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
@.str = [36 x i8] c"spglib: No inverse matrix (det=%f)\0A\00"
define i32 @mat_inverse_matrix_d3([3 x double]* %0, [3 x double]* %1, double %2) {
%4 = alloca [3 x [3 x double]]
%5 = bitcast [3 x [3 x double]]* %4 to i8*
%6 = bitcast [3 x double]* %1 to double**
%7 = tail call double @mat_get_determinant_d3(double** %6)
%8 = tail call double @mat_Dabs(double %7)
%9 = fcmp olt double %8, %2
br i1 %9, label %10, label %12
10: ; preds = %3
%11 = tail call i32 @warning_print(i8* getelementptr inbounds ([36 x i8], [36 x i8]* @.str, i64 0, i64 0), double %7)
br label %79
12: ; preds = %3
%13 = getelementptr inbounds [3 x double], [3 x double]* %1, i64 1, i64 1
%14 = load double, double* %13
%15 = getelementptr inbounds [3 x double], [3 x double]* %1, i64 2, i64 2
%16 = load double, double* %15
%17 = getelementptr inbounds [3 x double], [3 x double]* %1, i64 1, i64 2
%18 = getelementptr inbounds [3 x double], [3 x double]* %1, i64 2, i64 1
%19 = load double, double* %18
%20 = bitcast double* %17 to <2 x double>*
%21 = load <2 x double>, <2 x double>* %20
%22 = extractelement <2 x double> %21, i64 1
%23 = getelementptr inbounds [3 x double], [3 x double]* %1, i64 0, i64 2
%24 = bitcast double* %23 to <2 x double>*
%25 = load <2 x double>, <2 x double>* %24
%26 = extractelement <2 x double> %25, i64 1
%27 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %4, i64 0, i64 2, i64 0
%28 = getelementptr inbounds [3 x double], [3 x double]* %1, i64 0, i64 1
%29 = load double, double* %28
%30 = insertelement <2 x double> %21, double %16, i64 1
%31 = fneg <2 x double> %30
%32 = insertelement <2 x double> poison, double %19, i64 0
%33 = insertelement <2 x double> %32, double %29, i64 1
%34 = fmul <2 x double> %33, %31
%35 = insertelement <2 x double> poison, double %14, i64 0
%36 = insertelement <2 x double> %35, double %19, i64 1
%37 = insertelement <2 x double> poison, double %16, i64 0
%38 = shufflevector <2 x double> %37, <2 x double> %25, <2 x i32> <i32 0, i32 2>
%39 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %36, <2 x double> %38, <2 x double> %34)
%40 = insertelement <2 x double> poison, double %7, i64 0
%41 = shufflevector <2 x double> %40, <2 x double> poison, <2 x i32> zeroinitializer
%42 = fdiv <2 x double> %39, %41
%43 = bitcast [3 x [3 x double]]* %4 to <2 x double>*
store <2 x double> %42, <2 x double>* %43
%44 = getelementptr inbounds [3 x double], [3 x double]* %1, i64 0, i64 0
%45 = load double, double* %44
%46 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %4, i64 0, i64 1, i64 1
%47 = fneg <2 x double> %36
%48 = insertelement <2 x double> poison, double %22, i64 0
%49 = insertelement <2 x double> %48, double %45, i64 1
%50 = fmul <2 x double> %49, %47
%51 = shufflevector <2 x double> %25, <2 x double> %21, <2 x i32> <i32 1, i32 3>
%52 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %51, <2 x double> %33, <2 x double> %50)
%53 = fdiv <2 x double> %52, %41
%54 = bitcast double* %27 to <2 x double>*
store <2 x double> %53, <2 x double>* %54
%55 = fneg <2 x double> %25
%56 = insertelement <2 x double> %35, double %16, i64 1
%57 = fmul <2 x double> %56, %55
%58 = insertelement <2 x double> poison, double %29, i64 0
%59 = shufflevector <2 x double> %58, <2 x double> %21, <2 x i32> <i32 0, i32 2>
%60 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %59, <2 x double> %21, <2 x double> %57)
%61 = fdiv <2 x double> %60, %41
%62 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %4, i64 0, i64 0, i64 2
%63 = bitcast double* %62 to <2 x double>*
store <2 x double> %61, <2 x double>* %63
%64 = fneg <2 x double> %49
%65 = shufflevector <2 x double> %25, <2 x double> %21, <2 x i32> <i32 0, i32 2>
%66 = fmul <2 x double> %65, %64
%67 = insertelement <2 x double> %25, double %45, i64 0
%68 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %38, <2 x double> %67, <2 x double> %66)
%69 = fdiv <2 x double> %68, %41
%70 = bitcast double* %46 to <2 x double>*
store <2 x double> %69, <2 x double>* %70
%71 = fneg double %29
%72 = fmul double %26, %71
%73 = tail call double @llvm.fmuladd.f64(double %45, double %14, double %72)
%74 = fdiv double %73, %7
%75 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %4, i64 0, i64 2, i64 2
store double %74, double* %75
%76 = bitcast [3 x double]* %0 to double**
%77 = bitcast [3 x [3 x double]]* %4 to double**
%78 = call i32 @mat_copy_matrix_d3(double** %76, double** %77)
br label %79
79: ; preds = %12, %10
%80 = phi i32 [ 0, %10 ], [ 1, %12 ]
ret i32 %80
}