#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 sqrt(float);
void MatToQuat(float srcMatrix[3][3], float destQuat[4]) {
  double trace, s;
  int i, j, k;
  static int next[3] = {1, 2, 0};
  trace = srcMatrix[0][0] + srcMatrix[1][1] + srcMatrix[2][2];
  if (trace > 0.0) {
    s = sqrt(trace + 1.0);
    destQuat[3] = s * 0.5;
    s = 0.5 / s;
    destQuat[0] = (srcMatrix[2][1] - srcMatrix[1][2]) * s;
    destQuat[1] = (srcMatrix[0][2] - srcMatrix[2][0]) * s;
    destQuat[2] = (srcMatrix[1][0] - srcMatrix[0][1]) * s;
  } else {
    i = 0;
    if (srcMatrix[1][1] > srcMatrix[0][0])
      i = 1;
    if (srcMatrix[2][2] > srcMatrix[i][i])
      i = 2;
    j = next[i];
    k = next[j];
    s = sqrt((srcMatrix[i][i] - (srcMatrix[j][j] + srcMatrix[k][k])) + 1.0);
    destQuat[i] = s * 0.5;
    s = 0.5 / s;
    destQuat[3] = (srcMatrix[k][j] - srcMatrix[j][k]) * s;
    destQuat[j] = (srcMatrix[j][i] + srcMatrix[i][j]) * s;
    destQuat[k] = (srcMatrix[k][i] + srcMatrix[i][k]) * s;
  }
}


x86 assembly (clang, -O0):
<MatToQuat>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x8(%rbp),%rax
addss  0x10(%rax),%xmm0
mov    -0x8(%rbp),%rax
addss  0x20(%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1218 <MatToQuat+0x108>
movsd  0xeab(%rip),%xmm0
addsd  -0x18(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x20(%rbp)
movsd  0xe88(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0xc(%rax)
movsd  0xe6e(%rip),%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movss  0x1c(%rax),%xmm0
mov    -0x8(%rbp),%rax
subss  0x14(%rax),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  -0x20(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
movss  0x8(%rax),%xmm0
mov    -0x8(%rbp),%rax
subss  0x18(%rax),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  -0x20(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x4(%rax)
mov    -0x8(%rbp),%rax
movss  0xc(%rax),%xmm0
mov    -0x8(%rbp),%rax
subss  0x4(%rax),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  -0x20(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x8(%rax)
jmp    13ff <MatToQuat+0x2ef>
movl   $0x0,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movss  0x10(%rax),%xmm0
mov    -0x8(%rbp),%rax
ucomiss (%rax),%xmm0
jbe    123c <MatToQuat+0x12c>
movl   $0x1,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movss  0x20(%rax),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
ucomiss (%rax,%rcx,4),%xmm0
jbe    1269 <MatToQuat+0x159>
movl   $0x2,-0x24(%rbp)
movslq -0x24(%rbp),%rcx
lea    0x2db4(%rip),%rax
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x28(%rbp)
movslq -0x28(%rbp),%rcx
lea    0x2da3(%rip),%rax
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm1
subss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xd25(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x20(%rbp)
movsd  0xd03(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movsd  0xce5(%rip),%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  -0x20(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0xc(%rax)
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  -0x20(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  -0x20(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<MatToQuat>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %rsi,%r14
mov    %rdi,%rbx
movss  (%rdi),%xmm2
movss  0x10(%rdi),%xmm3
movaps %xmm2,%xmm0
addss  %xmm3,%xmm0
movss  0x20(%rdi),%xmm1
addss  %xmm1,%xmm0
xorps  %xmm4,%xmm4
ucomiss %xmm4,%xmm0
jbe    11b8 <MatToQuat+0xa8>
addss  0xeb6(%rip),%xmm0
call   1030 <sqrt@plt>
movsd  0xeb1(%rip),%xmm1
movaps %xmm0,%xmm2
mulsd  %xmm1,%xmm2
cvtsd2ss %xmm2,%xmm2
movss  %xmm2,0xc(%r14)
divsd  %xmm0,%xmm1
movss  0x1c(%rbx),%xmm0
subss  0x14(%rbx),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%r14)
movss  0x8(%rbx),%xmm0
subss  0x18(%rbx),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x4(%r14)
movss  0xc(%rbx),%xmm0
subss  0x4(%rbx),%xmm0
mov    $0x2,%r15d
jmp    1291 <MatToQuat+0x181>
xor    %eax,%eax
ucomiss %xmm2,%xmm3
seta   %al
lea    (%rax,%rax,2),%rcx
lea    (%rbx,%rcx,4),%rcx
ucomiss (%rcx,%rax,4),%xmm1
mov    $0x2,%r12d
cmovbe %rax,%r12
lea    0xe33(%rip),%rax
movslq (%rax,%r12,4),%r13
movslq (%rax,%r13,4),%r15
lea    (%r12,%r12,2),%rax
lea    (%rbx,%rax,4),%rcx
mov    %rcx,(%rsp)
lea    0x0(%r13,%r13,2),%rax
lea    (%rbx,%rax,4),%rbp
movss  0x0(%rbp,%r13,4),%xmm1
lea    (%r15,%r15,2),%rax
lea    (%rbx,%rax,4),%rbx
addss  (%rbx,%r15,4),%xmm1
movss  (%rcx,%r12,4),%xmm0
subss  %xmm1,%xmm0
addss  0xddf(%rip),%xmm0
call   1030 <sqrt@plt>
movsd  0xdda(%rip),%xmm1
movaps %xmm0,%xmm2
mulsd  %xmm1,%xmm2
cvtsd2ss %xmm2,%xmm2
movss  %xmm2,(%r14,%r12,4)
divsd  %xmm0,%xmm1
movss  (%rbx,%r13,4),%xmm0
subss  0x0(%rbp,%r15,4),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0xc(%r14)
movss  0x0(%rbp,%r12,4),%xmm0
mov    (%rsp),%rax
addss  (%rax,%r13,4),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%r14,%r13,4)
movss  (%rbx,%r12,4),%xmm0
addss  (%rax,%r15,4),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%r14,%r15,4)
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@MatToQuat.next = internal unnamed_addr constant [3 x i32] [i32 1, i32 2, i32 0]
define void @MatToQuat([3 x float]* nocapture %0, float* nocapture writeonly %1) {
%3 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 0, i64 0
%4 = load float, float* %3
%5 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 1, i64 1
%6 = load float, float* %5
%7 = fadd float %4, %6
%8 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 2, i64 2
%9 = load float, float* %8
%10 = fadd float %7, %9
%11 = fcmp ogt float %10, 0.000000e+00
br i1 %11, label %12, label %41
12: ; preds = %2
%13 = fadd float %10, 1.000000e+00
%14 = tail call double @sqrt(float %13)
%15 = fmul double %14, 5.000000e-01
%16 = fptrunc double %15 to float
%17 = getelementptr inbounds float, float* %1, i64 3
store float %16, float* %17
%18 = fdiv double 5.000000e-01, %14
%19 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 2, i64 1
%20 = load float, float* %19
%21 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 1, i64 2
%22 = load float, float* %21
%23 = fsub float %20, %22
%24 = fpext float %23 to double
%25 = fmul double %18, %24
%26 = fptrunc double %25 to float
store float %26, float* %1
%27 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 0, i64 2
%28 = load float, float* %27
%29 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 2, i64 0
%30 = load float, float* %29
%31 = fsub float %28, %30
%32 = fpext float %31 to double
%33 = fmul double %18, %32
%34 = fptrunc double %33 to float
%35 = getelementptr inbounds float, float* %1, i64 1
store float %34, float* %35
%36 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 1, i64 0
%37 = load float, float* %36
%38 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 0, i64 1
%39 = load float, float* %38
%40 = fsub float %37, %39
br label %91
41: ; preds = %2
%42 = fcmp ogt float %6, %4
%43 = zext i1 %42 to i64
%44 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %43, i64 %43
%45 = load float, float* %44
%46 = fcmp ogt float %9, %45
%47 = select i1 %46, i64 2, i64 %43
%48 = getelementptr inbounds [3 x i32], [3 x i32]* @MatToQuat.next, i64 0, i64 %47
%49 = load i32, i32* %48
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds [3 x i32], [3 x i32]* @MatToQuat.next, i64 0, i64 %50
%52 = load i32, i32* %51
%53 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %47, i64 %47
%54 = load float, float* %53
%55 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %50, i64 %50
%56 = load float, float* %55
%57 = sext i32 %52 to i64
%58 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %57, i64 %57
%59 = load float, float* %58
%60 = fadd float %56, %59
%61 = fsub float %54, %60
%62 = fadd float %61, 1.000000e+00
%63 = tail call double @sqrt(float %62)
%64 = fmul double %63, 5.000000e-01
%65 = fptrunc double %64 to float
%66 = getelementptr inbounds float, float* %1, i64 %47
store float %65, float* %66
%67 = fdiv double 5.000000e-01, %63
%68 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %57, i64 %50
%69 = load float, float* %68
%70 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %50, i64 %57
%71 = load float, float* %70
%72 = fsub float %69, %71
%73 = fpext float %72 to double
%74 = fmul double %67, %73
%75 = fptrunc double %74 to float
%76 = getelementptr inbounds float, float* %1, i64 3
store float %75, float* %76
%77 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %50, i64 %47
%78 = load float, float* %77
%79 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %47, i64 %50
%80 = load float, float* %79
%81 = fadd float %78, %80
%82 = fpext float %81 to double
%83 = fmul double %67, %82
%84 = fptrunc double %83 to float
%85 = getelementptr inbounds float, float* %1, i64 %50
store float %84, float* %85
%86 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %57, i64 %47
%87 = load float, float* %86
%88 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %47, i64 %57
%89 = load float, float* %88
%90 = fadd float %87, %89
br label %91
91: ; preds = %41, %12
%92 = phi float [ %90, %41 ], [ %40, %12 ]
%93 = phi double [ %67, %41 ], [ %18, %12 ]
%94 = phi i64 [ %57, %41 ], [ 2, %12 ]
%95 = fpext float %92 to double
%96 = fmul double %93, %95
%97 = fptrunc double %96 to float
%98 = getelementptr inbounds float, float* %1, i64 %94
store float %97, float* %98
ret void
}