#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(int);
int matMakeT2O(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;
  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] = 1.0 / a;
  rm[1] = 0.0;
  rm[2] = 0.0;
  rm[3] = 0.0;
  rm[4] = -cc / (a * sc);
  rm[5] = 1.0 / (b * sc);
  rm[6] = 0.0;
  rm[7] = 0.0;
  rm[8] = ax * cbx;
  rm[9] = bx * cax;
  rm[10] = cx;
  rm[11] = 0.0;
  rm[12] = 0.0;
  rm[13] = 0.0;
  rm[14] = 0.0;
  rm[15] = 1.0;
  return 0;
}


x86 assembly (clang, -O0):
<matMakeT2O>:
push   %rbp
mov    %rsp,%rbp
sub    $0xd0,%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,-0xc8(%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
cvttsd2si %xmm0,%edi
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0xc8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  0xd05(%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  0xb46(%rip),%xmm0
divsd  -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  -0x80(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x8(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,0x20(%rax)
movsd  -0x10(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
movsd  0xad0(%rip),%xmm0
divsd  %xmm1,%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  -0x40(%rbp),%xmm0
mulsd  -0xa8(%rbp),%xmm0
ff
mov    -0x38(%rbp),%rax
movsd  %xmm0,0x40(%rax)
movsd  -0x48(%rbp),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
mov    -0x38(%rbp),%rax
movsd  %xmm0,0x48(%rax)
movsd  -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  0xa35(%rip),%xmm0
movsd  %xmm0,0x78(%rax)
xor    %eax,%eax
add    $0xd0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<matMakeT2O>:
push   %rbx
sub    $0xd0,%rsp
mov    0x2ea9(%rip),%rax
movsd  (%rax),%xmm6
divsd  0xeb5(%rip),%xmm6
mov    %rdi,%rbx
movaps %xmm2,0x50(%rsp)
movaps %xmm1,0x30(%rsp)
movaps %xmm0,0x40(%rsp)
movapd %xmm3,%xmm0
mulsd  %xmm6,%xmm0
movsd  %xmm0,0x20(%rsp)
mulsd  %xmm6,%xmm4
movsd  %xmm4,(%rsp)
mulsd  %xmm5,%xmm6
movsd  %xmm6,0x10(%rsp)
call   1030 <sin@plt>
movapd %xmm0,0xc0(%rsp)
movsd  (%rsp),%xmm0
call   1030 <sin@plt>
movaps %xmm0,0xa0(%rsp)
movsd  0x10(%rsp),%xmm0
call   1030 <sin@plt>
movaps %xmm0,0xb0(%rsp)
movsd  0x20(%rsp),%xmm0
call   1040 <cos@plt>
movaps %xmm0,0x20(%rsp)
movsd  (%rsp),%xmm0
call   1040 <cos@plt>
movaps %xmm0,(%rsp)
movq   0x10(%rsp),%xmm0
call   1040 <cos@plt>
movapd 0x40(%rsp),%xmm1
mulsd  0x30(%rsp),%xmm1
movapd %xmm1,0x80(%rsp)
mulsd  0x50(%rsp),%xmm1
movapd %xmm1,0x70(%rsp)
movdqa (%rsp),%xmm4
movdqa %xmm4,%xmm5
movdqa 0x20(%rsp),%xmm3
punpcklqdq %xmm3,%xmm5
movdqa %xmm5,0x90(%rsp)
movdqa 0xdf1(%rip),%xmm1
pxor   %xmm1,%xmm5
movdqa %xmm1,%xmm2
pshufd $0xee,%xmm5,%xmm1
mulsd  %xmm3,%xmm1
addsd  0xdd0(%rip),%xmm1
mulsd  %xmm4,%xmm5
addsd  %xmm1,%xmm5
pxor   %xmm0,%xmm2
movdqa %xmm2,0x10(%rsp)
movdqa %xmm0,%xmm1
movdqa %xmm0,0x60(%rsp)
mulsd  %xmm0,%xmm1
subsd  %xmm1,%xmm5
movapd %xmm3,%xmm1
addsd  %xmm3,%xmm1
mulsd  %xmm4,%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm5,%xmm1
cvttsd2si %xmm1,%edi
call   1050 <sqrt@plt>
mulsd  0x70(%rsp),%xmm0
movsd  0xd80(%rip),%xmm8
movapd %xmm8,%xmm1
divsd  %xmm0,%xmm1
movapd 0x30(%rsp),%xmm2
movapd %xmm2,%xmm0
movapd 0x40(%rsp),%xmm3
unpcklpd %xmm3,%xmm0
movapd 0x50(%rsp),%xmm4
unpcklpd %xmm4,%xmm4
mulpd  %xmm0,%xmm4
movapd 0xc0(%rsp),%xmm0
unpcklpd 0xa0(%rsp),%xmm0
mulpd  %xmm0,%xmm4
movapd %xmm0,%xmm5
movapd 0xb0(%rsp),%xmm0
movapd 0x80(%rsp),%xmm7
mulsd  %xmm0,%xmm7
mulsd  %xmm1,%xmm7
unpcklpd %xmm1,%xmm1
mulpd  %xmm4,%xmm1
movapd 0x20(%rsp),%xmm4
unpcklpd (%rsp),%xmm4
movapd 0x60(%rsp),%xmm6
unpcklpd %xmm6,%xmm6
mulpd  %xmm4,%xmm6
subpd  0x90(%rsp),%xmm6
unpcklpd %xmm0,%xmm0
mulpd  %xmm0,%xmm5
movapd %xmm0,%xmm4
divpd  %xmm5,%xmm6
mulpd  %xmm1,%xmm6
movapd %xmm3,%xmm1
movapd %xmm8,%xmm0
divsd  %xmm3,%xmm0
movsd  %xmm0,(%rbx)
xorpd  %xmm0,%xmm0
movupd %xmm0,0x8(%rbx)
movq   $0x0,0x18(%rbx)
unpcklpd %xmm2,%xmm1
mulpd  %xmm4,%xmm1
movapd 0x10(%rsp),%xmm2
shufpd $0x2,0xcc2(%rip),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,0x20(%rbx)
movupd %xmm0,0x30(%rbx)
movupd %xmm6,0x40(%rbx)
movsd  %xmm7,0x50(%rbx)
movupd %xmm0,0x58(%rbx)
movupd %xmm0,0x68(%rbx)
movabs $0x3ff0000000000000,%rax
mov    %rax,0x78(%rbx)
xor    %eax,%eax
add    $0xd0,%rsp
pop    %rbx
ret

LLVM IR (-O0):
@M_PI = global double 0.000000e+00
define i32 @matMakeT2O(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
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
%32 = load double, double* @M_PI
%33 = fdiv double %32, 1.800000e+02
%34 = load double, double* %11
%35 = fmul double %34, %33
store double %35, double* %11
%36 = load double, double* @M_PI
%37 = fdiv double %36, 1.800000e+02
%38 = load double, double* %12
%39 = fmul double %38, %37
store double %39, double* %12
%40 = load double, double* @M_PI
%41 = fdiv double %40, 1.800000e+02
%42 = load double, double* %13
%43 = fmul double %42, %41
store double %43, double* %13
%44 = load double, double* %11
%45 = call double @sin(double %44)
store double %45, double* %18
%46 = load double, double* %12
%47 = call double @sin(double %46)
store double %47, double* %19
%48 = load double, double* %13
%49 = call double @sin(double %48)
store double %49, double* %20
%50 = load double, double* %11
%51 = call double @cos(double %50)
store double %51, double* %21
%52 = load double, double* %12
%53 = call double @cos(double %52)
store double %53, double* %22
%54 = load double, double* %13
%55 = call double @cos(double %54)
store double %55, double* %23
%56 = load double, double* %8
%57 = load double, double* %9
%58 = fmul double %56, %57
%59 = load double, double* %10
%60 = fmul double %58, %59
%61 = load double, double* %21
%62 = load double, double* %21
%63 = fneg double %61
%64 = call double @llvm.fmuladd.f64(double %63, double %62, double 1.000000e+00)
%65 = load double, double* %22
%66 = load double, double* %22
%67 = fneg double %65
%68 = call double @llvm.fmuladd.f64(double %67, double %66, double %64)
%69 = load double, double* %23
%70 = load double, double* %23
%71 = fneg double %69
%72 = call double @llvm.fmuladd.f64(double %71, double %70, double %68)
%73 = load double, double* %21
%74 = fmul double 2.000000e+00, %73
%75 = load double, double* %22
%76 = fmul double %74, %75
%77 = load double, double* %23
%78 = call double @llvm.fmuladd.f64(double %76, double %77, double %72)
%79 = fptosi double %78 to i32
%80 = call double @sqrt(i32 %79)
%81 = fmul double %60, %80
store double %81, double* %30
%82 = load double, double* %30
%83 = fdiv double 1.000000e+00, %82
store double %83, double* %31
%84 = load double, double* %9
%85 = load double, double* %10
%86 = fmul double %84, %85
%87 = load double, double* %18
%88 = fmul double %86, %87
%89 = load double, double* %31
%90 = fmul double %88, %89
store double %90, double* %15
%91 = load double, double* %8
%92 = load double, double* %10
%93 = fmul double %91, %92
%94 = load double, double* %19
%95 = fmul double %93, %94
%96 = load double, double* %31
%97 = fmul double %95, %96
store double %97, double* %16
%98 = load double, double* %8
%99 = load double, double* %9
%100 = fmul double %98, %99
%101 = load double, double* %20
%102 = fmul double %100, %101
%103 = load double, double* %31
%104 = fmul double %102, %103
store double %104, double* %17
%105 = load double, double* %30
%106 = load double, double* %8
%107 = load double, double* %9
%108 = fmul double %106, %107
%109 = load double, double* %10
%110 = fmul double %108, %109
%111 = load double, double* %19
%112 = fmul double %110, %111
%113 = load double, double* %20
%114 = fmul double %112, %113
%115 = fdiv double %105, %114
store double %115, double* %24
%116 = load double, double* %30
%117 = load double, double* %8
%118 = load double, double* %9
%119 = fmul double %117, %118
%120 = load double, double* %10
%121 = fmul double %119, %120
%122 = load double, double* %18
%123 = fmul double %121, %122
%124 = load double, double* %20
%125 = fmul double %123, %124
%126 = fdiv double %116, %125
store double %126, double* %25
%127 = load double, double* %30
%128 = load double, double* %8
%129 = load double, double* %9
%130 = fmul double %128, %129
%131 = load double, double* %10
%132 = fmul double %130, %131
%133 = load double, double* %18
%134 = fmul double %132, %133
%135 = load double, double* %19
%136 = fmul double %134, %135
%137 = fdiv double %127, %136
store double %137, double* %26
%138 = load double, double* %22
%139 = load double, double* %23
%140 = load double, double* %21
%141 = fneg double %140
%142 = call double @llvm.fmuladd.f64(double %138, double %139, double %141)
%143 = load double, double* %19
%144 = load double, double* %20
%145 = fmul double %143, %144
%146 = fdiv double %142, %145
store double %146, double* %27
%147 = load double, double* %21
%148 = load double, double* %23
%149 = load double, double* %22
%150 = fneg double %149
%151 = call double @llvm.fmuladd.f64(double %147, double %148, double %150)
%152 = load double, double* %18
%153 = load double, double* %20
%154 = fmul double %152, %153
%155 = fdiv double %151, %154
store double %155, double* %28
%156 = load double, double* %21
%157 = load double, double* %22
%158 = load double, double* %23
%159 = fneg double %158
%160 = call double @llvm.fmuladd.f64(double %156, double %157, double %159)
%161 = load double, double* %18
%162 = load double, double* %19
%163 = fmul double %161, %162
%164 = fdiv double %160, %163
store double %164, double* %29
%165 = load double, double* %8
%166 = fdiv double 1.000000e+00, %165
%167 = load double*, double** %14
%168 = getelementptr inbounds double, double* %167, i64 0
store double %166, double* %168
%169 = load double*, double** %14
%170 = getelementptr inbounds double, double* %169, i64 1
store double 0.000000e+00, double* %170
%171 = load double*, double** %14
%172 = getelementptr inbounds double, double* %171, i64 2
store double 0.000000e+00, double* %172
%173 = load double*, double** %14
%174 = getelementptr inbounds double, double* %173, i64 3
store double 0.000000e+00, double* %174
%175 = load double, double* %23
%176 = fneg double %175
%177 = load double, double* %8
%178 = load double, double* %20
%179 = fmul double %177, %178
%180 = fdiv double %176, %179
%181 = load double*, double** %14
%182 = getelementptr inbounds double, double* %181, i64 4
store double %180, double* %182
%183 = load double, double* %9
%184 = load double, double* %20
%185 = fmul double %183, %184
%186 = fdiv double 1.000000e+00, %185
%187 = load double*, double** %14
%188 = getelementptr inbounds double, double* %187, i64 5
store double %186, double* %188
%189 = load double*, double** %14
%190 = getelementptr inbounds double, double* %189, i64 6
store double 0.000000e+00, double* %190
%191 = load double*, double** %14
%192 = getelementptr inbounds double, double* %191, i64 7
store double 0.000000e+00, double* %192
%193 = load double, double* %15
%194 = load double, double* %28
%195 = fmul double %193, %194
%196 = load double*, double** %14
%197 = getelementptr inbounds double, double* %196, i64 8
store double %195, double* %197
%198 = load double, double* %16
%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 = load double*, double** %14
%205 = getelementptr inbounds double, double* %204, i64 10
store double %203, double* %205
%206 = load double*, double** %14
%207 = getelementptr inbounds double, double* %206, i64 11
store double 0.000000e+00, double* %207
%208 = load double*, double** %14
%209 = getelementptr inbounds double, double* %208, i64 12
store double 0.000000e+00, double* %209
%210 = load double*, double** %14
%211 = getelementptr inbounds double, double* %210, i64 13
store double 0.000000e+00, double* %211
%212 = load double*, double** %14
%213 = getelementptr inbounds double, double* %212, i64 14
store double 0.000000e+00, double* %213
%214 = load double*, double** %14
%215 = getelementptr inbounds double, double* %214, i64 15
store double 1.000000e+00, double* %215
ret i32 0
}

LLVM IR (-O2):
@M_PI = global double 0.000000e+00
define i32 @matMakeT2O(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 = tail call double @sin(double %10)
%14 = tail call double @sin(double %11)
%15 = tail call double @sin(double %12)
%16 = tail call double @cos(double %10)
%17 = tail call double @cos(double %11)
%18 = tail call double @cos(double %12)
%19 = fmul double %0, %1
%20 = fmul double %19, %2
%21 = insertelement <2 x double> poison, double %17, i64 0
%22 = insertelement <2 x double> %21, double %16, i64 1
%23 = fneg <2 x double> %22
%24 = extractelement <2 x double> %23, i64 1
%25 = tail call double @llvm.fmuladd.f64(double %24, double %16, double 1.000000e+00)
%26 = extractelement <2 x double> %23, i64 0
%27 = tail call double @llvm.fmuladd.f64(double %26, double %17, double %25)
%28 = fneg double %18
%29 = tail call double @llvm.fmuladd.f64(double %28, double %18, double %27)
%30 = fmul double %16, 2.000000e+00
%31 = fmul double %30, %17
%32 = tail call double @llvm.fmuladd.f64(double %31, double %18, double %29)
%33 = fptosi double %32 to i32
%34 = tail call double @sqrt(i32 %33)
%35 = fmul double %20, %34
%36 = fdiv double 1.000000e+00, %35
%37 = insertelement <2 x double> poison, double %1, i64 0
%38 = insertelement <2 x double> %37, double %0, i64 1
%39 = insertelement <2 x double> poison, double %2, i64 0
%40 = shufflevector <2 x double> %39, <2 x double> poison, <2 x i32> zeroinitializer
%41 = fmul <2 x double> %38, %40
%42 = insertelement <2 x double> poison, double %13, i64 0
%43 = insertelement <2 x double> %42, double %14, i64 1
%44 = fmul <2 x double> %41, %43
%45 = insertelement <2 x double> poison, double %36, i64 0
%46 = shufflevector <2 x double> %45, <2 x double> poison, <2 x i32> zeroinitializer
%47 = fmul <2 x double> %44, %46
%48 = fmul double %19, %15
%49 = fmul double %48, %36
%50 = insertelement <2 x double> poison, double %16, i64 0
%51 = insertelement <2 x double> %50, double %17, i64 1
%52 = insertelement <2 x double> poison, double %18, i64 0
%53 = shufflevector <2 x double> %52, <2 x double> poison, <2 x i32> zeroinitializer
%54 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %51, <2 x double> %53, <2 x double> %23)
%55 = insertelement <2 x double> poison, double %15, i64 0
%56 = shufflevector <2 x double> %55, <2 x double> poison, <2 x i32> zeroinitializer
%57 = fmul <2 x double> %43, %56
%58 = fdiv <2 x double> %54, %57
%59 = fdiv double 1.000000e+00, %0
store double %59, double* %6
%60 = getelementptr inbounds double, double* %6, i64 1
%61 = getelementptr inbounds double, double* %6, i64 4
%62 = bitcast double* %60 to i8*
%63 = insertelement <2 x double> poison, double %0, i64 0
%64 = insertelement <2 x double> %63, double %1, i64 1
%65 = fmul <2 x double> %56, %64
%66 = insertelement <2 x double> <double poison, double 1.000000e+00>, double %28, i64 0
%67 = fdiv <2 x double> %66, %65
%68 = bitcast double* %61 to <2 x double>*
store <2 x double> %67, <2 x double>* %68
%69 = getelementptr inbounds double, double* %6, i64 6
%70 = getelementptr inbounds double, double* %6, i64 8
%71 = bitcast double* %69 to i8*
%72 = fmul <2 x double> %58, %47
%73 = bitcast double* %70 to <2 x double>*
store <2 x double> %72, <2 x double>* %73
%74 = getelementptr inbounds double, double* %6, i64 10
store double %49, double* %74
%75 = getelementptr inbounds double, double* %6, i64 11
%76 = getelementptr inbounds double, double* %6, i64 15
%77 = bitcast double* %75 to i8*
store double 1.000000e+00, double* %76
ret i32 0
}