#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 M_PI;
double cos(double);
double sin(double);
double sqrt(double);
int matMakeO2T(double a, double b, double c, double alpha, double beta,
               double gamma, double *rm) {
  double ax, bx, cx;
  double sa, sb, sc, ca, cb, cc;
  double sax, sbx, scx, cax, cbx, ccx;
  double v, v_1;
  double p[3], q[3], r[3];
  alpha *= M_PI / 180.0;
  beta *= M_PI / 180.0;
  gamma *= M_PI / 180.0;
  sa = sin(alpha);
  sb = sin(beta);
  sc = sin(gamma);
  ca = cos(alpha);
  cb = cos(beta);
  cc = cos(gamma);
  v = a * b * c * sqrt(1 - ca * ca - cb * cb - cc * cc + 2 * ca * cb * cc);
  v_1 = 1.0 / v;
  ax = b * c * sa * v_1;
  bx = a * c * sb * v_1;
  cx = a * b * sc * v_1;
  sax = v / (a * b * c * sb * sc);
  sbx = v / (a * b * c * sa * sc);
  scx = v / (a * b * c * sa * sb);
  cax = (cb * cc - ca) / (sb * sc);
  cbx = (ca * cc - cb) / (sa * sc);
  ccx = (ca * cb - cc) / (sa * sb);
  rm[0] = a;
  rm[1] = 0.0;
  rm[2] = 0.0;
  rm[3] = 0.0;
  rm[4] = b * cc;
  rm[5] = b * sc;
  rm[6] = 0.0;
  rm[7] = 0.0;
  rm[8] = c * cb;
  rm[9] = -c * sb * cax;
  rm[10] = 1.0 / cx;
  rm[11] = 0.0;
  rm[12] = 0.0;
  rm[13] = 0.0;
  rm[14] = 0.0;
  rm[15] = 1.0;
  p[0] = a;
  p[1] = 0.0;
  p[2] = 0.0;
  q[0] = b * cc;
  q[1] = b * sc;
  q[2] = 0.0;
  r[0] = c * cb;
  r[1] = (b * c * ca - b * c * cc * cb) / (b * sc);
  r[2] = sqrt(c * c - r[0] * r[0] - r[1] * r[1]);
  return 0;
}


x86 assembly (clang, -O0):
<matMakeO2T>:
push   %rbp
mov    %rsp,%rbp
sub    $0x130,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movsd  %xmm3,-0x20(%rbp)
movsd  %xmm4,-0x28(%rbp)
movsd  %xmm5,-0x30(%rbp)
mov    %rdi,-0x38(%rbp)
mov    0x2e84(%rip),%rax
movsd  (%rax),%xmm0
movsd  0xea0(%rip),%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    0x2e63(%rip),%rax
movsd  (%rax),%xmm0
movsd  0xe7f(%rip),%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    0x2e42(%rip),%rax
movsd  (%rax),%xmm0
movsd  0xe5e(%rip),%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x20(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x60(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x68(%rbp)
movsd  -0x20(%rbp),%xmm0
call   1040 <cos@plt>
movsd  %xmm0,-0x70(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1040 <cos@plt>
movsd  %xmm0,-0x78(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1040 <cos@plt>
movsd  %xmm0,-0x80(%rbp)
movsd  -0x8(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x70(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm0
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm0,%xmm1
movsd  0xda2(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x80(%rbp),%xmm1
movsd  -0x80(%rbp),%xmm2
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  0xd4c(%rip),%xmm0
mulsd  -0x70(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x128(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  0xd09(%rip),%xmm0
divsd  -0xb8(%rbp),%xmm0
ff
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
mulsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x40(%rbp)
movsd  -0x8(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x48(%rbp)
movsd  -0x8(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
mulsd  -0x68(%rbp),%xmm0
mulsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x50(%rbp)
movsd  -0xb8(%rbp),%xmm0
ff
movsd  -0x8(%rbp),%xmm1
mulsd  -0x10(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
movsd  -0x8(%rbp),%xmm1
mulsd  -0x10(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
mulsd  -0x58(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
movsd  -0x8(%rbp),%xmm1
mulsd  -0x10(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
mulsd  -0x58(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
movsd  -0x70(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x60(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
movsd  -0x78(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x58(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm2
movsd  -0x80(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x58(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x8(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x38(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x8(%rax)
mov    -0x38(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x10(%rax)
mov    -0x38(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x18(%rax)
movsd  -0x10(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,0x20(%rax)
movsd  -0x10(%rbp),%xmm0
mulsd  -0x68(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,0x28(%rax)
mov    -0x38(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x30(%rax)
mov    -0x38(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x38(%rax)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,0x40(%rax)
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
mov    -0x38(%rbp),%rax
movsd  %xmm0,0x48(%rax)
movsd  0xa96(%rip),%xmm0
divsd  -0x50(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,0x50(%rax)
mov    -0x38(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x58(%rax)
mov    -0x38(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x60(%rax)
mov    -0x38(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x68(%rax)
mov    -0x38(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x70(%rax)
mov    -0x38(%rbp),%rax
movsd  0xa4c(%rip),%xmm0
movsd  %xmm0,0x78(%rax)
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
mulsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm2
movsd  -0x10(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
mulsd  -0x80(%rbp),%xmm1
mulsd  -0x78(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x10(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
movsd  -0x18(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm2
movsd  -0x120(%rbp),%xmm0
ff
mulsd  -0x120(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x118(%rbp),%xmm0
ff
movsd  -0x118(%rbp),%xmm2
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <sqrt@plt>
movsd  %xmm0,-0x110(%rbp)
ff
xor    %eax,%eax
add    $0x130,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<matMakeO2T>:
push   %rbx
sub    $0x80,%rsp
mov    %rdi,%rbx
movaps %xmm2,0x60(%rsp)
movsd  %xmm1,0x28(%rsp)
movsd  %xmm0,0x20(%rsp)
mov    0x2e95(%rip),%rax
movsd  (%rax),%xmm0
divsd  0xea1(%rip),%xmm0
mulsd  %xmm0,%xmm3
mulsd  %xmm0,%xmm4
mulsd  %xmm5,%xmm0
movsd  %xmm0,0x8(%rsp)
ucomisd 0xe8f(%rip),%xmm3
movsd  %xmm4,0x10(%rsp)
movsd  %xmm3,0x30(%rsp)
jae    13f4 <matMakeO2T+0x2c4>
movsd  0xe7d(%rip),%xmm0
ucomisd %xmm3,%xmm0
jae    13f4 <matMakeO2T+0x2c4>
movapd %xmm4,%xmm0
call   1030 <sin@plt>
movsd  %xmm0,0x50(%rsp)
movsd  0x8(%rsp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,0x58(%rsp)
movsd  0x30(%rsp),%xmm0
call   1040 <cos@plt>
movaps %xmm0,0x30(%rsp)
movsd  0x10(%rsp),%xmm0
call   1040 <cos@plt>
movaps %xmm0,0x10(%rsp)
movsd  0x8(%rsp),%xmm0
call   1040 <cos@plt>
movapd %xmm0,%xmm3
movsd  0x20(%rsp),%xmm7
mulsd  0x28(%rsp),%xmm7
movapd %xmm7,%xmm4
movapd 0x60(%rsp),%xmm5
mulsd  %xmm5,%xmm4
movapd 0xe11(%rip),%xmm8
movapd 0x30(%rsp),%xmm1
xorpd  %xmm1,%xmm8
movapd %xmm1,%xmm0
movapd %xmm1,%xmm6
mulsd  %xmm1,%xmm0
movsd  0xde9(%rip),%xmm12
movapd %xmm12,%xmm1
subsd  %xmm0,%xmm1
movapd 0x10(%rsp),%xmm2
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
subsd  %xmm0,%xmm1
movapd %xmm3,%xmm0
mulsd  %xmm3,%xmm0
subsd  %xmm0,%xmm1
movapd %xmm6,%xmm0
addsd  %xmm6,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm1,%xmm0
xorpd  %xmm11,%xmm11
ucomisd %xmm11,%xmm0
jb     127c <matMakeO2T+0x14c>
sqrtsd %xmm0,%xmm0
jmp    12c7 <matMakeO2T+0x197>
movsd  %xmm3,0x8(%rsp)
movsd  %xmm7,0x48(%rsp)
movapd %xmm8,0x70(%rsp)
movsd  %xmm4,0x40(%rsp)
call   1050 <sqrt@plt>
movsd  0x40(%rsp),%xmm4
movapd 0x70(%rsp),%xmm8
movsd  0xd68(%rip),%xmm12
xorpd  %xmm11,%xmm11
movsd  0x48(%rsp),%xmm7
movsd  0x8(%rsp),%xmm3
movapd 0x60(%rsp),%xmm5
mulsd  %xmm0,%xmm4
movapd %xmm12,%xmm1
divsd  %xmm4,%xmm1
movsd  0x28(%rsp),%xmm4
movapd %xmm4,%xmm0
mulsd  %xmm5,%xmm0
movsd  0x58(%rsp),%xmm6
mulsd  %xmm6,%xmm7
mulsd  %xmm1,%xmm7
movapd 0x10(%rsp),%xmm9
movapd %xmm9,%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm8,%xmm1
movsd  0x50(%rsp),%xmm10
movapd %xmm10,%xmm2
mulsd  %xmm6,%xmm2
divsd  %xmm2,%xmm1
movsd  0x20(%rsp),%xmm2
movsd  %xmm2,(%rbx)
movapd %xmm3,%xmm2
mulsd  %xmm4,%xmm2
xorpd  %xmm8,%xmm8
movupd %xmm8,0x8(%rbx)
movq   $0x0,0x18(%rbx)
movsd  %xmm2,0x20(%rbx)
mulsd  %xmm4,%xmm6
movsd  %xmm6,0x28(%rbx)
movupd %xmm8,0x30(%rbx)
movapd 0xcc6(%rip),%xmm2
xorpd  %xmm5,%xmm2
mulsd  %xmm10,%xmm2
movapd %xmm9,%xmm4
unpcklpd %xmm2,%xmm4
movapd %xmm5,%xmm2
unpcklpd %xmm1,%xmm2
mulpd  %xmm4,%xmm2
movupd %xmm2,0x40(%rbx)
divsd  %xmm7,%xmm12
movsd  %xmm12,0x50(%rbx)
movupd %xmm8,0x58(%rbx)
movupd %xmm8,0x68(%rbx)
movabs $0x3ff0000000000000,%rax
mov    %rax,0x78(%rbx)
mulsd  %xmm0,%xmm3
mulsd  %xmm9,%xmm3
mulsd  0x30(%rsp),%xmm0
subsd  %xmm3,%xmm0
divsd  %xmm6,%xmm0
mulsd  %xmm2,%xmm2
mulsd  %xmm5,%xmm5
subsd  %xmm2,%xmm5
mulsd  %xmm0,%xmm0
subsd  %xmm0,%xmm5
ucomisd %xmm5,%xmm11
ja     13df <matMakeO2T+0x2af>
xor    %eax,%eax
add    $0x80,%rsp
pop    %rbx
ret
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm5
jae    13d4 <matMakeO2T+0x2a4>
movapd %xmm5,%xmm0
call   1050 <sqrt@plt>
jmp    13d4 <matMakeO2T+0x2a4>
movapd %xmm3,%xmm0
call   1030 <sin@plt>
movsd  0x10(%rsp),%xmm4
jmp    119d <matMakeO2T+0x6d>

LLVM IR (-O0):
@M_PI = global double 0.000000e+00
define i32 @matMakeO2T(double %0, double %1, double %2, double %3, double %4, double %5, double* %6) {
%8 = alloca double
%9 = alloca double
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double*
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca [3 x double]
%33 = alloca [3 x double]
%34 = alloca [3 x double]
store double %0, double* %8
store double %1, double* %9
store double %2, double* %10
store double %3, double* %11
store double %4, double* %12
store double %5, double* %13
store double* %6, double** %14
%35 = load double, double* @M_PI
%36 = fdiv double %35, 1.800000e+02
%37 = load double, double* %11
%38 = fmul double %37, %36
store double %38, double* %11
%39 = load double, double* @M_PI
%40 = fdiv double %39, 1.800000e+02
%41 = load double, double* %12
%42 = fmul double %41, %40
store double %42, double* %12
%43 = load double, double* @M_PI
%44 = fdiv double %43, 1.800000e+02
%45 = load double, double* %13
%46 = fmul double %45, %44
store double %46, double* %13
%47 = load double, double* %11
%48 = call double @sin(double %47)
store double %48, double* %18
%49 = load double, double* %12
%50 = call double @sin(double %49)
store double %50, double* %19
%51 = load double, double* %13
%52 = call double @sin(double %51)
store double %52, double* %20
%53 = load double, double* %11
%54 = call double @cos(double %53)
store double %54, double* %21
%55 = load double, double* %12
%56 = call double @cos(double %55)
store double %56, double* %22
%57 = load double, double* %13
%58 = call double @cos(double %57)
store double %58, double* %23
%59 = load double, double* %8
%60 = load double, double* %9
%61 = fmul double %59, %60
%62 = load double, double* %10
%63 = fmul double %61, %62
%64 = load double, double* %21
%65 = load double, double* %21
%66 = fneg double %64
%67 = call double @llvm.fmuladd.f64(double %66, double %65, double 1.000000e+00)
%68 = load double, double* %22
%69 = load double, double* %22
%70 = fneg double %68
%71 = call double @llvm.fmuladd.f64(double %70, double %69, double %67)
%72 = load double, double* %23
%73 = load double, double* %23
%74 = fneg double %72
%75 = call double @llvm.fmuladd.f64(double %74, double %73, double %71)
%76 = load double, double* %21
%77 = fmul double 2.000000e+00, %76
%78 = load double, double* %22
%79 = fmul double %77, %78
%80 = load double, double* %23
%81 = call double @llvm.fmuladd.f64(double %79, double %80, double %75)
%82 = call double @sqrt(double %81)
%83 = fmul double %63, %82
store double %83, double* %30
%84 = load double, double* %30
%85 = fdiv double 1.000000e+00, %84
store double %85, double* %31
%86 = load double, double* %9
%87 = load double, double* %10
%88 = fmul double %86, %87
%89 = load double, double* %18
%90 = fmul double %88, %89
%91 = load double, double* %31
%92 = fmul double %90, %91
store double %92, double* %15
%93 = load double, double* %8
%94 = load double, double* %10
%95 = fmul double %93, %94
%96 = load double, double* %19
%97 = fmul double %95, %96
%98 = load double, double* %31
%99 = fmul double %97, %98
store double %99, double* %16
%100 = load double, double* %8
%101 = load double, double* %9
%102 = fmul double %100, %101
%103 = load double, double* %20
%104 = fmul double %102, %103
%105 = load double, double* %31
%106 = fmul double %104, %105
store double %106, double* %17
%107 = load double, double* %30
%108 = load double, double* %8
%109 = load double, double* %9
%110 = fmul double %108, %109
%111 = load double, double* %10
%112 = fmul double %110, %111
%113 = load double, double* %19
%114 = fmul double %112, %113
%115 = load double, double* %20
%116 = fmul double %114, %115
%117 = fdiv double %107, %116
store double %117, double* %24
%118 = load double, double* %30
%119 = load double, double* %8
%120 = load double, double* %9
%121 = fmul double %119, %120
%122 = load double, double* %10
%123 = fmul double %121, %122
%124 = load double, double* %18
%125 = fmul double %123, %124
%126 = load double, double* %20
%127 = fmul double %125, %126
%128 = fdiv double %118, %127
store double %128, double* %25
%129 = load double, double* %30
%130 = load double, double* %8
%131 = load double, double* %9
%132 = fmul double %130, %131
%133 = load double, double* %10
%134 = fmul double %132, %133
%135 = load double, double* %18
%136 = fmul double %134, %135
%137 = load double, double* %19
%138 = fmul double %136, %137
%139 = fdiv double %129, %138
store double %139, double* %26
%140 = load double, double* %22
%141 = load double, double* %23
%142 = load double, double* %21
%143 = fneg double %142
%144 = call double @llvm.fmuladd.f64(double %140, double %141, double %143)
%145 = load double, double* %19
%146 = load double, double* %20
%147 = fmul double %145, %146
%148 = fdiv double %144, %147
store double %148, double* %27
%149 = load double, double* %21
%150 = load double, double* %23
%151 = load double, double* %22
%152 = fneg double %151
%153 = call double @llvm.fmuladd.f64(double %149, double %150, double %152)
%154 = load double, double* %18
%155 = load double, double* %20
%156 = fmul double %154, %155
%157 = fdiv double %153, %156
store double %157, double* %28
%158 = load double, double* %21
%159 = load double, double* %22
%160 = load double, double* %23
%161 = fneg double %160
%162 = call double @llvm.fmuladd.f64(double %158, double %159, double %161)
%163 = load double, double* %18
%164 = load double, double* %19
%165 = fmul double %163, %164
%166 = fdiv double %162, %165
store double %166, double* %29
%167 = load double, double* %8
%168 = load double*, double** %14
%169 = getelementptr inbounds double, double* %168, i64 0
store double %167, double* %169
%170 = load double*, double** %14
%171 = getelementptr inbounds double, double* %170, i64 1
store double 0.000000e+00, double* %171
%172 = load double*, double** %14
%173 = getelementptr inbounds double, double* %172, i64 2
store double 0.000000e+00, double* %173
%174 = load double*, double** %14
%175 = getelementptr inbounds double, double* %174, i64 3
store double 0.000000e+00, double* %175
%176 = load double, double* %9
%177 = load double, double* %23
%178 = fmul double %176, %177
%179 = load double*, double** %14
%180 = getelementptr inbounds double, double* %179, i64 4
store double %178, double* %180
%181 = load double, double* %9
%182 = load double, double* %20
%183 = fmul double %181, %182
%184 = load double*, double** %14
%185 = getelementptr inbounds double, double* %184, i64 5
store double %183, double* %185
%186 = load double*, double** %14
%187 = getelementptr inbounds double, double* %186, i64 6
store double 0.000000e+00, double* %187
%188 = load double*, double** %14
%189 = getelementptr inbounds double, double* %188, i64 7
store double 0.000000e+00, double* %189
%190 = load double, double* %10
%191 = load double, double* %22
%192 = fmul double %190, %191
%193 = load double*, double** %14
%194 = getelementptr inbounds double, double* %193, i64 8
store double %192, double* %194
%195 = load double, double* %10
%196 = fneg double %195
%197 = load double, double* %19
%198 = fmul double %196, %197
%199 = load double, double* %27
%200 = fmul double %198, %199
%201 = load double*, double** %14
%202 = getelementptr inbounds double, double* %201, i64 9
store double %200, double* %202
%203 = load double, double* %17
%204 = fdiv double 1.000000e+00, %203
%205 = load double*, double** %14
%206 = getelementptr inbounds double, double* %205, i64 10
store double %204, double* %206
%207 = load double*, double** %14
%208 = getelementptr inbounds double, double* %207, i64 11
store double 0.000000e+00, double* %208
%209 = load double*, double** %14
%210 = getelementptr inbounds double, double* %209, i64 12
store double 0.000000e+00, double* %210
%211 = load double*, double** %14
%212 = getelementptr inbounds double, double* %211, i64 13
store double 0.000000e+00, double* %212
%213 = load double*, double** %14
%214 = getelementptr inbounds double, double* %213, i64 14
store double 0.000000e+00, double* %214
%215 = load double*, double** %14
%216 = getelementptr inbounds double, double* %215, i64 15
store double 1.000000e+00, double* %216
%217 = load double, double* %8
%218 = getelementptr inbounds [3 x double], [3 x double]* %32, i64 0, i64 0
store double %217, double* %218
%219 = getelementptr inbounds [3 x double], [3 x double]* %32, i64 0, i64 1
store double 0.000000e+00, double* %219
%220 = getelementptr inbounds [3 x double], [3 x double]* %32, i64 0, i64 2
store double 0.000000e+00, double* %220
%221 = load double, double* %9
%222 = load double, double* %23
%223 = fmul double %221, %222
%224 = getelementptr inbounds [3 x double], [3 x double]* %33, i64 0, i64 0
store double %223, double* %224
%225 = load double, double* %9
%226 = load double, double* %20
%227 = fmul double %225, %226
%228 = getelementptr inbounds [3 x double], [3 x double]* %33, i64 0, i64 1
store double %227, double* %228
%229 = getelementptr inbounds [3 x double], [3 x double]* %33, i64 0, i64 2
store double 0.000000e+00, double* %229
%230 = load double, double* %10
%231 = load double, double* %22
%232 = fmul double %230, %231
%233 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 0
store double %232, double* %233
%234 = load double, double* %9
%235 = load double, double* %10
%236 = fmul double %234, %235
%237 = load double, double* %21
%238 = load double, double* %9
%239 = load double, double* %10
%240 = fmul double %238, %239
%241 = load double, double* %23
%242 = fmul double %240, %241
%243 = load double, double* %22
%244 = fmul double %242, %243
%245 = fneg double %244
%246 = call double @llvm.fmuladd.f64(double %236, double %237, double %245)
%247 = load double, double* %9
%248 = load double, double* %20
%249 = fmul double %247, %248
%250 = fdiv double %246, %249
%251 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 1
store double %250, double* %251
%252 = load double, double* %10
%253 = load double, double* %10
%254 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 0
%255 = load double, double* %254
%256 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 0
%257 = load double, double* %256
%258 = fmul double %255, %257
%259 = fneg double %258
%260 = call double @llvm.fmuladd.f64(double %252, double %253, double %259)
%261 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 1
%262 = load double, double* %261
%263 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 1
%264 = load double, double* %263
%265 = fneg double %262
%266 = call double @llvm.fmuladd.f64(double %265, double %264, double %260)
%267 = call double @sqrt(double %266)
%268 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 2
store double %267, double* %268
ret i32 0
}

LLVM IR (-O2):
@M_PI = global double 0.000000e+00
define i32 @matMakeO2T(double %0, double %1, double %2, double %3, double %4, double %5, double* nocapture writeonly %6) {
%8 = load double, double* @M_PI
%9 = fdiv double %8, 1.800000e+02
%10 = fmul double %9, %3
%11 = fmul double %9, %4
%12 = fmul double %9, %5
%13 = fcmp oeq double %10, 0xFFF0000000000000
%14 = fcmp oeq double %10, 0x7FF0000000000000
%15 = or i1 %14, %13
br i1 %15, label %16, label %18
16: ; preds = %7
%17 = tail call double @sin(double %10)
br label %18
18: ; preds = %7, %16
%19 = tail call double @sin(double %11)
%20 = tail call double @sin(double %12)
%21 = tail call double @cos(double %10)
%22 = tail call double @cos(double %11)
%23 = tail call double @cos(double %12)
%24 = fmul double %0, %1
%25 = fmul double %24, %2
%26 = fneg double %21
%27 = tail call double @llvm.fmuladd.f64(double %26, double %21, double 1.000000e+00)
%28 = fneg double %22
%29 = tail call double @llvm.fmuladd.f64(double %28, double %22, double %27)
%30 = fneg double %23
%31 = tail call double @llvm.fmuladd.f64(double %30, double %23, double %29)
%32 = fmul double %21, 2.000000e+00
%33 = fmul double %32, %22
%34 = tail call double @llvm.fmuladd.f64(double %33, double %23, double %31)
%35 = tail call double @sqrt(double %34)
%36 = fmul double %25, %35
%37 = fdiv double 1.000000e+00, %36
%38 = fmul double %1, %2
%39 = fmul double %24, %20
%40 = fmul double %39, %37
%41 = tail call double @llvm.fmuladd.f64(double %22, double %23, double %26)
%42 = fmul double %19, %20
%43 = fdiv double %41, %42
store double %0, double* %6
%44 = getelementptr inbounds double, double* %6, i64 1
%45 = fmul double %23, %1
%46 = getelementptr inbounds double, double* %6, i64 4
%47 = bitcast double* %44 to i8*
store double %45, double* %46
%48 = fmul double %20, %1
%49 = getelementptr inbounds double, double* %6, i64 5
store double %48, double* %49
%50 = getelementptr inbounds double, double* %6, i64 6
%51 = getelementptr inbounds double, double* %6, i64 8
%52 = bitcast double* %50 to i8*
%53 = fneg double %2
%54 = fmul double %19, %53
%55 = insertelement <2 x double> poison, double %22, i64 0
%56 = insertelement <2 x double> %55, double %54, i64 1
%57 = insertelement <2 x double> poison, double %2, i64 0
%58 = insertelement <2 x double> %57, double %43, i64 1
%59 = fmul <2 x double> %56, %58
%60 = bitcast double* %51 to <2 x double>*
store <2 x double> %59, <2 x double>* %60
%61 = fdiv double 1.000000e+00, %40
%62 = getelementptr inbounds double, double* %6, i64 10
store double %61, double* %62
%63 = getelementptr inbounds double, double* %6, i64 11
%64 = getelementptr inbounds double, double* %6, i64 15
%65 = bitcast double* %63 to i8*
store double 1.000000e+00, double* %64
%66 = fneg double %38
%67 = fmul double %23, %66
%68 = fmul double %22, %67
%69 = tail call double @llvm.fmuladd.f64(double %38, double %21, double %68)
%70 = fdiv double %69, %48
%71 = extractelement <2 x double> %59, i64 0
%72 = fneg double %71
%73 = fmul double %71, %72
%74 = tail call double @llvm.fmuladd.f64(double %2, double %2, double %73)
%75 = fneg double %70
%76 = tail call double @llvm.fmuladd.f64(double %75, double %70, double %74)
%77 = fcmp olt double %76, 0.000000e+00
br i1 %77, label %78, label %80
78: ; preds = %18
%79 = tail call double @sqrt(double %76)
br label %80
80: ; preds = %18, %78
ret i32 0
}