void coreGTRGAMMAMULT(int lower, int upper, double *sumtable,
                      double *ext_dlnLdlz, double *ext_d2lnLdlz2, double *EIGN,
                      double *gammaRates, double lz, int *wrptr,
                      int numberOfModels, int *modelptr) {
  double *sum, *diagptable, *diagptable_start;
  int i;
  double dlnLdlz = 0;
  double d2lnLdlz2 = 0;
  double tmp_1, tmp_2, tmp_3;
  double inv_Li, dlnLidlz, d2lnLidlz2;
  double ki, kisqr;
  int model;
  diagptable = diagptable_start =
      (double *)malloc(sizeof(double) * 36 * numberOfModels);
  for (model = 0; model < numberOfModels; model++) {
    diagptable = &diagptable_start[model * 36];
    for (i = 0; i < 4; i++) {
      ki = gammaRates[model * 4 + i];
      kisqr = ki * ki;
      *diagptable++ = exp(EIGN[model * 3] * ki * lz);
      *diagptable++ = exp(EIGN[model * 3 + 1] * ki * lz);
      *diagptable++ = exp(EIGN[model * 3 + 2] * ki * lz);
      *diagptable++ = EIGN[model * 3] * ki;
      *diagptable++ = EIGN[model * 3] * EIGN[model * 3] * kisqr;
      *diagptable++ = EIGN[model * 3 + 1] * ki;
      *diagptable++ = EIGN[model * 3 + 1] * EIGN[model * 3 + 1] * kisqr;
      *diagptable++ = EIGN[model * 3 + 2] * ki;
      *diagptable++ = EIGN[model * 3 + 2] * EIGN[model * 3 + 2] * kisqr;
    }
  }
  for (i = lower; i < upper; i++) {
    diagptable = &diagptable_start[36 * modelptr[i]];
    sum = &(sumtable[i * 16]);
    inv_Li = *sum++;
    inv_Li += (tmp_1 = *diagptable++ * *sum++);
    inv_Li += (tmp_2 = *diagptable++ * *sum++);
    inv_Li += (tmp_3 = *diagptable++ * *sum++);
    dlnLidlz = tmp_1 * *diagptable++;
    d2lnLidlz2 = tmp_1 * *diagptable++;
    dlnLidlz += tmp_2 * *diagptable++;
    d2lnLidlz2 += tmp_2 * *diagptable++;
    dlnLidlz += tmp_3 * *diagptable++;
    d2lnLidlz2 += tmp_3 * *diagptable++;
    inv_Li += *sum++;
    inv_Li += (tmp_1 = *sum++ * *diagptable++);
    inv_Li += (tmp_2 = *sum++ * *diagptable++);
    inv_Li += (tmp_3 = *sum++ * *diagptable++);
    dlnLidlz += tmp_1 * *diagptable++;
    d2lnLidlz2 += tmp_1 * *diagptable++;
    dlnLidlz += tmp_2 * *diagptable++;
    d2lnLidlz2 += tmp_2 * *diagptable++;
    dlnLidlz += tmp_3 * *diagptable++;
    d2lnLidlz2 += tmp_3 * *diagptable++;
    inv_Li += *sum++;
    inv_Li += (tmp_1 = *sum++ * *diagptable++);
    inv_Li += (tmp_2 = *sum++ * *diagptable++);
    inv_Li += (tmp_3 = *sum++ * *diagptable++);
    dlnLidlz += tmp_1 * *diagptable++;
    d2lnLidlz2 += tmp_1 * *diagptable++;
    dlnLidlz += tmp_2 * *diagptable++;
    d2lnLidlz2 += tmp_2 * *diagptable++;
    dlnLidlz += tmp_3 * *diagptable++;
    d2lnLidlz2 += tmp_3 * *diagptable++;
    inv_Li += *sum++;
    inv_Li += (tmp_1 = *sum++ * *diagptable++);
    inv_Li += (tmp_2 = *sum++ * *diagptable++);
    inv_Li += (tmp_3 = *sum++ * *diagptable++);
    dlnLidlz += tmp_1 * *diagptable++;
    d2lnLidlz2 += tmp_1 * *diagptable++;
    dlnLidlz += tmp_2 * *diagptable++;
    d2lnLidlz2 += tmp_2 * *diagptable++;
    dlnLidlz += tmp_3 * *diagptable++;
    d2lnLidlz2 += tmp_3 * *diagptable++;
    inv_Li = 1.0 / inv_Li;
    dlnLidlz *= inv_Li;
    d2lnLidlz2 *= inv_Li;
    dlnLdlz += wrptr[i] * dlnLidlz;
    d2lnLdlz2 += wrptr[i] * (d2lnLidlz2 - dlnLidlz * dlnLidlz);
  }
  *ext_dlnLdlz = dlnLdlz;
  *ext_d2lnLdlz2 = d2lnLdlz2;
  free(diagptable_start);
}


x86 assembly (clang, -O0):
<coreGTRGAMMAMULT>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
movslq 0x20(%rbp),%rax
imul   $0x120,%rax,%rdi
call   1040 <malloc@plt>
mov    %rax,-0x48(%rbp)
mov    %rax,-0x40(%rbp)
movl   $0x0,-0xa4(%rbp)
mov    -0xa4(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    143c <coreGTRGAMMAMULT+0x30c>
mov    -0x48(%rbp),%rax
imul   $0x24,-0xa4(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x4c(%rbp)
cmpl   $0x4,-0x4c(%rbp)
jge    1423 <coreGTRGAMMAMULT+0x2f3>
mov    0x10(%rbp),%rax
mov    -0xa4(%rbp),%ecx
shl    $0x2,%ecx
add    -0x4c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
mulsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
mulsd  -0x30(%rbp),%xmm0
call   1050 <exp@plt>
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
mulsd  -0x30(%rbp),%xmm0
call   1050 <exp@plt>
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
mulsd  -0x30(%rbp),%xmm0
call   1050 <exp@plt>
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
imul   $0x3,-0xa4(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11c6 <coreGTRGAMMAMULT+0x96>
jmp    1428 <coreGTRGAMMAMULT+0x2f8>
mov    -0xa4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa4(%rbp)
jmp    1197 <coreGTRGAMMAMULT+0x67>
mov    -0x4(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1c96 <coreGTRGAMMAMULT+0xb66>
mov    -0x48(%rbp),%rax
mov    0x28(%rbp),%rcx
movslq -0x4c(%rbp),%rdx
imul   $0x24,(%rcx,%rdx,4),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
shl    $0x4,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x68(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x68(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x68(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x68(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  0x429(%rip),%xmm0
divsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
mov    0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x90(%rbp),%xmm1
ff
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x88(%rbp),%xmm3
ff
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1442 <coreGTRGAMMAMULT+0x312>
movsd  -0x58(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x60(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<coreGTRGAMMAMULT>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
movsd  %xmm0,(%rsp)
mov    %r9,%r14
mov    %r8,0x48(%rsp)
mov    %rcx,0x40(%rsp)
mov    %rdx,0x30(%rsp)
mov    %esi,0xc(%rsp)
mov    %rdi,0x38(%rsp)
movslq 0xc0(%rsp),%rbp
mov    0xc8(%rsp),%rax
mov    %rax,0x28(%rsp)
mov    0xb8(%rsp),%rax
mov    %rax,0x20(%rsp)
mov    0xb0(%rsp),%r15
mov    %rbp,%rax
shl    $0x5,%rax
lea    (%rax,%rax,8),%rdi
call   1040 <malloc@plt>
mov    %rax,%rbx
test   %rbp,%rbp
jle    12cd <coreGTRGAMMAMULT+0x19d>
mov    %ebp,%eax
mov    %rax,0x50(%rsp)
xor    %ecx,%ecx
nopl   (%rax)
lea    0x0(,%rcx,4),%eax
lea    (%rax,%rax,8),%eax
lea    (%rbx,%rax,8),%r12
mov    %rcx,0x58(%rsp)
lea    (%rcx,%rcx,2),%r13
movsd  (%r14,%r13,8),%xmm0
xor    %ebp,%ebp
nop
movsd  (%r15,%rbp,8),%xmm1
movapd %xmm1,0x10(%rsp)
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
movapd %xmm2,0x60(%rsp)
mulsd  %xmm1,%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,(%r12)
movsd  0x8(%r14,%r13,8),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x8(%r12)
movsd  0x10(%r14,%r13,8),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x10(%r12)
movsd  (%r14,%r13,8),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
movapd 0x10(%rsp),%xmm3
unpcklpd 0x60(%rsp),%xmm3
movapd %xmm0,%xmm2
unpcklpd %xmm1,%xmm2
mulpd  %xmm3,%xmm2
movupd %xmm2,0x18(%r12)
movsd  0x8(%r14,%r13,8),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
unpcklpd %xmm2,%xmm1
mulpd  %xmm3,%xmm1
movupd %xmm1,0x28(%r12)
movsd  0x10(%r14,%r13,8),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
unpcklpd %xmm2,%xmm1
mulpd  %xmm3,%xmm1
movupd %xmm1,0x38(%r12)
add    $0x1,%rbp
add    $0x48,%r12
cmp    $0x4,%rbp
jne    11d0 <coreGTRGAMMAMULT+0xa0>
mov    0x58(%rsp),%rcx
add    $0x1,%rcx
add    $0x20,%r15
cmp    0x50(%rsp),%rcx
jne    11b0 <coreGTRGAMMAMULT+0x80>
xorpd  %xmm4,%xmm4
mov    0x38(%rsp),%rdi
mov    0xc(%rsp),%ecx
cmp    %ecx,%edi
jge    15ba <coreGTRGAMMAMULT+0x48a>
movslq %edi,%rax
movslq %ecx,%r8
shl    $0x4,%edi
xorpd  %xmm9,%xmm9
movsd  0xd07(%rip),%xmm8
mov    0x30(%rsp),%rsi
mov    0x28(%rsp),%rbp
mov    0x20(%rsp),%rcx
nopl   0x0(%rax,%rax,1)
movslq 0x0(%rbp,%rax,4),%rdx
lea    (%rdx,%rdx,8),%rdx
shl    $0x5,%rdx
movslq %edi,%rdi
movsd  (%rbx,%rdx,1),%xmm7
mulsd  0x8(%rsi,%rdi,8),%xmm7
movsd  0x8(%rbx,%rdx,1),%xmm6
mulsd  0x10(%rsi,%rdi,8),%xmm6
movsd  (%rsi,%rdi,8),%xmm2
movsd  0x10(%rbx,%rdx,1),%xmm3
mulsd  0x18(%rsi,%rdi,8),%xmm3
addsd  %xmm7,%xmm2
movsd  0x18(%rbx,%rdx,1),%xmm0
mulsd  %xmm7,%xmm0
mulsd  0x20(%rbx,%rdx,1),%xmm7
addsd  %xmm6,%xmm2
movsd  0x28(%rbx,%rdx,1),%xmm1
mulsd  %xmm6,%xmm1
addsd  %xmm0,%xmm1
mulsd  0x30(%rbx,%rdx,1),%xmm6
movsd  0x38(%rbx,%rdx,1),%xmm4
addsd  %xmm3,%xmm2
mulsd  %xmm3,%xmm4
movsd  0x28(%rsi,%rdi,8),%xmm0
mulsd  0x40(%rbx,%rdx,1),%xmm3
mulsd  0x48(%rbx,%rdx,1),%xmm0
addsd  0x20(%rsi,%rdi,8),%xmm2
movsd  0x30(%rsi,%rdi,8),%xmm5
mulsd  0x50(%rbx,%rdx,1),%xmm5
addsd  %xmm1,%xmm4
movsd  0x38(%rsi,%rdi,8),%xmm10
mulsd  0x58(%rbx,%rdx,1),%xmm10
addsd  %xmm7,%xmm6
addsd  %xmm0,%xmm2
movsd  0x60(%rbx,%rdx,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm6,%xmm3
addsd  %xmm4,%xmm1
mulsd  0x68(%rbx,%rdx,1),%xmm0
addsd  %xmm5,%xmm2
movsd  0x70(%rbx,%rdx,1),%xmm4
mulsd  %xmm5,%xmm4
mulsd  0x78(%rbx,%rdx,1),%xmm5
addsd  %xmm1,%xmm4
movsd  0x80(%rbx,%rdx,1),%xmm1
mulsd  %xmm10,%xmm1
addsd  %xmm10,%xmm2
addsd  %xmm4,%xmm1
mulsd  0x88(%rbx,%rdx,1),%xmm10
addsd  %xmm3,%xmm0
addsd  0x40(%rsi,%rdi,8),%xmm2
movsd  0x48(%rsi,%rdi,8),%xmm3
addsd  %xmm0,%xmm5
mulsd  0x90(%rbx,%rdx,1),%xmm3
movupd 0x50(%rsi,%rdi,8),%xmm0
addsd  %xmm3,%xmm2
movupd 0x98(%rbx,%rdx,1),%xmm4
mulpd  %xmm0,%xmm4
movsd  0xb8(%rbx,%rdx,1),%xmm0
addsd  %xmm5,%xmm10
mulsd  %xmm4,%xmm0
movsd  0xc0(%rbx,%rdx,1),%xmm11
mulsd  %xmm4,%xmm11
addsd  %xmm4,%xmm2
movsd  0xa8(%rbx,%rdx,1),%xmm5
mulsd  %xmm3,%xmm5
addsd  %xmm1,%xmm5
unpckhpd %xmm4,%xmm4
movsd  0xc8(%rbx,%rdx,1),%xmm12
mulsd  %xmm4,%xmm12
movsd  0xd0(%rbx,%rdx,1),%xmm13
mulsd  %xmm4,%xmm13
movapd %xmm4,%xmm1
mulsd  0xb0(%rbx,%rdx,1),%xmm3
addsd  %xmm5,%xmm0
addsd  %xmm2,%xmm1
addsd  0x60(%rsi,%rdi,8),%xmm1
movsd  0x68(%rsi,%rdi,8),%xmm2
addsd  %xmm0,%xmm12
mulsd  0xd8(%rbx,%rdx,1),%xmm2
movupd 0x70(%rsi,%rdi,8),%xmm0
addsd  %xmm10,%xmm3
movupd 0xe0(%rbx,%rdx,1),%xmm4
mulpd  %xmm0,%xmm4
movsd  0x100(%rbx,%rdx,1),%xmm0
addsd  %xmm2,%xmm1
mulsd  %xmm4,%xmm0
movsd  0x108(%rbx,%rdx,1),%xmm5
mulsd  %xmm4,%xmm5
addsd  %xmm4,%xmm1
movapd %xmm4,%xmm6
unpckhpd %xmm4,%xmm6
movsd  0x110(%rbx,%rdx,1),%xmm7
movsd  0x118(%rbx,%rdx,1),%xmm4
mulsd  %xmm6,%xmm7
mulsd  %xmm6,%xmm4
addsd  %xmm1,%xmm6
movapd %xmm8,%xmm1
divsd  %xmm6,%xmm1
addsd  %xmm3,%xmm11
movsd  0xf0(%rbx,%rdx,1),%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm12,%xmm3
addsd  %xmm11,%xmm13
addsd  %xmm3,%xmm0
mulsd  0xf8(%rbx,%rdx,1),%xmm2
addsd  %xmm13,%xmm2
addsd  %xmm2,%xmm5
addsd  %xmm0,%xmm7
mulsd  %xmm1,%xmm7
addsd  %xmm5,%xmm4
mulsd  %xmm1,%xmm4
xorps  %xmm0,%xmm0
cvtsi2sdl (%rcx,%rax,4),%xmm0
movapd %xmm7,%xmm1
mulsd  %xmm7,%xmm1
subsd  %xmm1,%xmm4
unpcklpd %xmm0,%xmm0
unpcklpd %xmm7,%xmm4
mulpd  %xmm0,%xmm4
movapd %xmm9,%xmm0
movapd %xmm4,%xmm9
addpd  %xmm0,%xmm9
add    $0x1,%rax
add    $0x10,%edi
cmp    %rax,%r8
jne    1310 <coreGTRGAMMAMULT+0x1e0>
addpd  %xmm0,%xmm4
mov    0x40(%rsp),%rax
movhpd %xmm4,(%rax)
mov    0x48(%rsp),%rax
movlpd %xmm4,(%rax)
mov    %rbx,%rdi
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @coreGTRGAMMAMULT(i32 %0, i32 %1, double* %2, double* %3, double* %4, double* %5, double* %6, double %7, i32* %8, i32 %9, i32* %10) {
%12 = alloca i32
%13 = alloca i32
%14 = alloca double*
%15 = alloca double*
%16 = alloca double*
%17 = alloca double*
%18 = alloca double*
%19 = alloca double
%20 = alloca i32*
%21 = alloca i32
%22 = alloca i32*
%23 = alloca double*
%24 = alloca double*
%25 = alloca double*
%26 = alloca i32
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca double
%33 = alloca double
%34 = alloca double
%35 = alloca double
%36 = alloca double
%37 = alloca i32
store i32 %0, i32* %12
store i32 %1, i32* %13
store double* %2, double** %14
store double* %3, double** %15
store double* %4, double** %16
store double* %5, double** %17
store double* %6, double** %18
store double %7, double* %19
store i32* %8, i32** %20
store i32 %9, i32* %21
store i32* %10, i32** %22
store double 0.000000e+00, double* %27
store double 0.000000e+00, double* %28
%38 = load i32, i32* %21
%39 = sext i32 %38 to i64
%40 = mul i64 288, %39
%41 = call i8* @malloc(i64 %40)
%42 = bitcast i8* %41 to double*
store double* %42, double** %25
store double* %42, double** %24
store i32 0, i32* %37
br label %43
43: ; preds = %200, %11
%44 = load i32, i32* %37
%45 = load i32, i32* %21
%46 = icmp slt i32 %44, %45
br i1 %46, label %47, label %203
47: ; preds = %43
%48 = load double*, double** %25
%49 = load i32, i32* %37
%50 = mul nsw i32 %49, 36
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %48, i64 %51
store double* %52, double** %24
store i32 0, i32* %26
br label %53
53: ; preds = %196, %47
%54 = load i32, i32* %26
%55 = icmp slt i32 %54, 4
br i1 %55, label %56, label %199
56: ; preds = %53
%57 = load double*, double** %18
%58 = load i32, i32* %37
%59 = mul nsw i32 %58, 4
%60 = load i32, i32* %26
%61 = add nsw i32 %59, %60
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %57, i64 %62
%64 = load double, double* %63
store double %64, double* %35
%65 = load double, double* %35
%66 = load double, double* %35
%67 = fmul double %65, %66
store double %67, double* %36
%68 = load double*, double** %17
%69 = load i32, i32* %37
%70 = mul nsw i32 %69, 3
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds double, double* %68, i64 %71
%73 = load double, double* %72
%74 = load double, double* %35
%75 = fmul double %73, %74
%76 = load double, double* %19
%77 = fmul double %75, %76
%78 = call double @exp(double %77)
%79 = load double*, double** %24
%80 = getelementptr inbounds double, double* %79, i32 1
store double* %80, double** %24
store double %78, double* %79
%81 = load double*, double** %17
%82 = load i32, i32* %37
%83 = mul nsw i32 %82, 3
%84 = add nsw i32 %83, 1
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %81, i64 %85
%87 = load double, double* %86
%88 = load double, double* %35
%89 = fmul double %87, %88
%90 = load double, double* %19
%91 = fmul double %89, %90
%92 = call double @exp(double %91)
%93 = load double*, double** %24
%94 = getelementptr inbounds double, double* %93, i32 1
store double* %94, double** %24
store double %92, double* %93
%95 = load double*, double** %17
%96 = load i32, i32* %37
%97 = mul nsw i32 %96, 3
%98 = add nsw i32 %97, 2
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double, double* %95, i64 %99
%101 = load double, double* %100
%102 = load double, double* %35
%103 = fmul double %101, %102
%104 = load double, double* %19
%105 = fmul double %103, %104
%106 = call double @exp(double %105)
%107 = load double*, double** %24
%108 = getelementptr inbounds double, double* %107, i32 1
store double* %108, double** %24
store double %106, double* %107
%109 = load double*, double** %17
%110 = load i32, i32* %37
%111 = mul nsw i32 %110, 3
%112 = sext i32 %111 to i64
%113 = getelementptr inbounds double, double* %109, i64 %112
%114 = load double, double* %113
%115 = load double, double* %35
%116 = fmul double %114, %115
%117 = load double*, double** %24
%118 = getelementptr inbounds double, double* %117, i32 1
store double* %118, double** %24
store double %116, double* %117
%119 = load double*, double** %17
%120 = load i32, i32* %37
%121 = mul nsw i32 %120, 3
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds double, double* %119, i64 %122
%124 = load double, double* %123
%125 = load double*, double** %17
%126 = load i32, i32* %37
%127 = mul nsw i32 %126, 3
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds double, double* %125, i64 %128
%130 = load double, double* %129
%131 = fmul double %124, %130
%132 = load double, double* %36
%133 = fmul double %131, %132
%134 = load double*, double** %24
%135 = getelementptr inbounds double, double* %134, i32 1
store double* %135, double** %24
store double %133, double* %134
%136 = load double*, double** %17
%137 = load i32, i32* %37
%138 = mul nsw i32 %137, 3
%139 = add nsw i32 %138, 1
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds double, double* %136, i64 %140
%142 = load double, double* %141
%143 = load double, double* %35
%144 = fmul double %142, %143
%145 = load double*, double** %24
%146 = getelementptr inbounds double, double* %145, i32 1
store double* %146, double** %24
store double %144, double* %145
%147 = load double*, double** %17
%148 = load i32, i32* %37
%149 = mul nsw i32 %148, 3
%150 = add nsw i32 %149, 1
%151 = sext i32 %150 to i64
%152 = getelementptr inbounds double, double* %147, i64 %151
%153 = load double, double* %152
%154 = load double*, double** %17
%155 = load i32, i32* %37
%156 = mul nsw i32 %155, 3
%157 = add nsw i32 %156, 1
%158 = sext i32 %157 to i64
%159 = getelementptr inbounds double, double* %154, i64 %158
%160 = load double, double* %159
%161 = fmul double %153, %160
%162 = load double, double* %36
%163 = fmul double %161, %162
%164 = load double*, double** %24
%165 = getelementptr inbounds double, double* %164, i32 1
store double* %165, double** %24
store double %163, double* %164
%166 = load double*, double** %17
%167 = load i32, i32* %37
%168 = mul nsw i32 %167, 3
%169 = add nsw i32 %168, 2
%170 = sext i32 %169 to i64
%171 = getelementptr inbounds double, double* %166, i64 %170
%172 = load double, double* %171
%173 = load double, double* %35
%174 = fmul double %172, %173
%175 = load double*, double** %24
%176 = getelementptr inbounds double, double* %175, i32 1
store double* %176, double** %24
store double %174, double* %175
%177 = load double*, double** %17
%178 = load i32, i32* %37
%179 = mul nsw i32 %178, 3
%180 = add nsw i32 %179, 2
%181 = sext i32 %180 to i64
%182 = getelementptr inbounds double, double* %177, i64 %181
%183 = load double, double* %182
%184 = load double*, double** %17
%185 = load i32, i32* %37
%186 = mul nsw i32 %185, 3
%187 = add nsw i32 %186, 2
%188 = sext i32 %187 to i64
%189 = getelementptr inbounds double, double* %184, i64 %188
%190 = load double, double* %189
%191 = fmul double %183, %190
%192 = load double, double* %36
%193 = fmul double %191, %192
%194 = load double*, double** %24
%195 = getelementptr inbounds double, double* %194, i32 1
store double* %195, double** %24
store double %193, double* %194
br label %196
196: ; preds = %56
%197 = load i32, i32* %26
%198 = add nsw i32 %197, 1
store i32 %198, i32* %26
br label %53
199: ; preds = %53
br label %200
200: ; preds = %199
%201 = load i32, i32* %37
%202 = add nsw i32 %201, 1
store i32 %202, i32* %37
br label %43
203: ; preds = %43
%204 = load i32, i32* %12
store i32 %204, i32* %26
br label %205
205: ; preds = %522, %203
%206 = load i32, i32* %26
%207 = load i32, i32* %13
%208 = icmp slt i32 %206, %207
br i1 %208, label %209, label %525
209: ; preds = %205
%210 = load double*, double** %25
%211 = load i32*, i32** %22
%212 = load i32, i32* %26
%213 = sext i32 %212 to i64
%214 = getelementptr inbounds i32, i32* %211, i64 %213
%215 = load i32, i32* %214
%216 = mul nsw i32 36, %215
%217 = sext i32 %216 to i64
%218 = getelementptr inbounds double, double* %210, i64 %217
store double* %218, double** %24
%219 = load double*, double** %14
%220 = load i32, i32* %26
%221 = mul nsw i32 %220, 16
%222 = sext i32 %221 to i64
%223 = getelementptr inbounds double, double* %219, i64 %222
store double* %223, double** %23
%224 = load double*, double** %23
%225 = getelementptr inbounds double, double* %224, i32 1
store double* %225, double** %23
%226 = load double, double* %224
store double %226, double* %32
%227 = load double*, double** %24
%228 = getelementptr inbounds double, double* %227, i32 1
store double* %228, double** %24
%229 = load double, double* %227
%230 = load double*, double** %23
%231 = getelementptr inbounds double, double* %230, i32 1
store double* %231, double** %23
%232 = load double, double* %230
%233 = fmul double %229, %232
store double %233, double* %29
%234 = load double, double* %32
%235 = fadd double %234, %233
store double %235, double* %32
%236 = load double*, double** %24
%237 = getelementptr inbounds double, double* %236, i32 1
store double* %237, double** %24
%238 = load double, double* %236
%239 = load double*, double** %23
%240 = getelementptr inbounds double, double* %239, i32 1
store double* %240, double** %23
%241 = load double, double* %239
%242 = fmul double %238, %241
store double %242, double* %30
%243 = load double, double* %32
%244 = fadd double %243, %242
store double %244, double* %32
%245 = load double*, double** %24
%246 = getelementptr inbounds double, double* %245, i32 1
store double* %246, double** %24
%247 = load double, double* %245
%248 = load double*, double** %23
%249 = getelementptr inbounds double, double* %248, i32 1
store double* %249, double** %23
%250 = load double, double* %248
%251 = fmul double %247, %250
store double %251, double* %31
%252 = load double, double* %32
%253 = fadd double %252, %251
store double %253, double* %32
%254 = load double, double* %29
%255 = load double*, double** %24
%256 = getelementptr inbounds double, double* %255, i32 1
store double* %256, double** %24
%257 = load double, double* %255
%258 = fmul double %254, %257
store double %258, double* %33
%259 = load double, double* %29
%260 = load double*, double** %24
%261 = getelementptr inbounds double, double* %260, i32 1
store double* %261, double** %24
%262 = load double, double* %260
%263 = fmul double %259, %262
store double %263, double* %34
%264 = load double, double* %30
%265 = load double*, double** %24
%266 = getelementptr inbounds double, double* %265, i32 1
store double* %266, double** %24
%267 = load double, double* %265
%268 = load double, double* %33
%269 = call double @llvm.fmuladd.f64(double %264, double %267, double %268)
store double %269, double* %33
%270 = load double, double* %30
%271 = load double*, double** %24
%272 = getelementptr inbounds double, double* %271, i32 1
store double* %272, double** %24
%273 = load double, double* %271
%274 = load double, double* %34
%275 = call double @llvm.fmuladd.f64(double %270, double %273, double %274)
store double %275, double* %34
%276 = load double, double* %31
%277 = load double*, double** %24
%278 = getelementptr inbounds double, double* %277, i32 1
store double* %278, double** %24
%279 = load double, double* %277
%280 = load double, double* %33
%281 = call double @llvm.fmuladd.f64(double %276, double %279, double %280)
store double %281, double* %33
%282 = load double, double* %31
%283 = load double*, double** %24
%284 = getelementptr inbounds double, double* %283, i32 1
store double* %284, double** %24
%285 = load double, double* %283
%286 = load double, double* %34
%287 = call double @llvm.fmuladd.f64(double %282, double %285, double %286)
store double %287, double* %34
%288 = load double*, double** %23
%289 = getelementptr inbounds double, double* %288, i32 1
store double* %289, double** %23
%290 = load double, double* %288
%291 = load double, double* %32
%292 = fadd double %291, %290
store double %292, double* %32
%293 = load double*, double** %23
%294 = getelementptr inbounds double, double* %293, i32 1
store double* %294, double** %23
%295 = load double, double* %293
%296 = load double*, double** %24
%297 = getelementptr inbounds double, double* %296, i32 1
store double* %297, double** %24
%298 = load double, double* %296
%299 = fmul double %295, %298
store double %299, double* %29
%300 = load double, double* %32
%301 = fadd double %300, %299
store double %301, double* %32
%302 = load double*, double** %23
%303 = getelementptr inbounds double, double* %302, i32 1
store double* %303, double** %23
%304 = load double, double* %302
%305 = load double*, double** %24
%306 = getelementptr inbounds double, double* %305, i32 1
store double* %306, double** %24
%307 = load double, double* %305
%308 = fmul double %304, %307
store double %308, double* %30
%309 = load double, double* %32
%310 = fadd double %309, %308
store double %310, double* %32
%311 = load double*, double** %23
%312 = getelementptr inbounds double, double* %311, i32 1
store double* %312, double** %23
%313 = load double, double* %311
%314 = load double*, double** %24
%315 = getelementptr inbounds double, double* %314, i32 1
store double* %315, double** %24
%316 = load double, double* %314
%317 = fmul double %313, %316
store double %317, double* %31
%318 = load double, double* %32
%319 = fadd double %318, %317
store double %319, double* %32
%320 = load double, double* %29
%321 = load double*, double** %24
%322 = getelementptr inbounds double, double* %321, i32 1
store double* %322, double** %24
%323 = load double, double* %321
%324 = load double, double* %33
%325 = call double @llvm.fmuladd.f64(double %320, double %323, double %324)
store double %325, double* %33
%326 = load double, double* %29
%327 = load double*, double** %24
%328 = getelementptr inbounds double, double* %327, i32 1
store double* %328, double** %24
%329 = load double, double* %327
%330 = load double, double* %34
%331 = call double @llvm.fmuladd.f64(double %326, double %329, double %330)
store double %331, double* %34
%332 = load double, double* %30
%333 = load double*, double** %24
%334 = getelementptr inbounds double, double* %333, i32 1
store double* %334, double** %24
%335 = load double, double* %333
%336 = load double, double* %33
%337 = call double @llvm.fmuladd.f64(double %332, double %335, double %336)
store double %337, double* %33
%338 = load double, double* %30
%339 = load double*, double** %24
%340 = getelementptr inbounds double, double* %339, i32 1
store double* %340, double** %24
%341 = load double, double* %339
%342 = load double, double* %34
%343 = call double @llvm.fmuladd.f64(double %338, double %341, double %342)
store double %343, double* %34
%344 = load double, double* %31
%345 = load double*, double** %24
%346 = getelementptr inbounds double, double* %345, i32 1
store double* %346, double** %24
%347 = load double, double* %345
%348 = load double, double* %33
%349 = call double @llvm.fmuladd.f64(double %344, double %347, double %348)
store double %349, double* %33
%350 = load double, double* %31
%351 = load double*, double** %24
%352 = getelementptr inbounds double, double* %351, i32 1
store double* %352, double** %24
%353 = load double, double* %351
%354 = load double, double* %34
%355 = call double @llvm.fmuladd.f64(double %350, double %353, double %354)
store double %355, double* %34
%356 = load double*, double** %23
%357 = getelementptr inbounds double, double* %356, i32 1
store double* %357, double** %23
%358 = load double, double* %356
%359 = load double, double* %32
%360 = fadd double %359, %358
store double %360, double* %32
%361 = load double*, double** %23
%362 = getelementptr inbounds double, double* %361, i32 1
store double* %362, double** %23
%363 = load double, double* %361
%364 = load double*, double** %24
%365 = getelementptr inbounds double, double* %364, i32 1
store double* %365, double** %24
%366 = load double, double* %364
%367 = fmul double %363, %366
store double %367, double* %29
%368 = load double, double* %32
%369 = fadd double %368, %367
store double %369, double* %32
%370 = load double*, double** %23
%371 = getelementptr inbounds double, double* %370, i32 1
store double* %371, double** %23
%372 = load double, double* %370
%373 = load double*, double** %24
%374 = getelementptr inbounds double, double* %373, i32 1
store double* %374, double** %24
%375 = load double, double* %373
%376 = fmul double %372, %375
store double %376, double* %30
%377 = load double, double* %32
%378 = fadd double %377, %376
store double %378, double* %32
%379 = load double*, double** %23
%380 = getelementptr inbounds double, double* %379, i32 1
store double* %380, double** %23
%381 = load double, double* %379
%382 = load double*, double** %24
%383 = getelementptr inbounds double, double* %382, i32 1
store double* %383, double** %24
%384 = load double, double* %382
%385 = fmul double %381, %384
store double %385, double* %31
%386 = load double, double* %32
%387 = fadd double %386, %385
store double %387, double* %32
%388 = load double, double* %29
%389 = load double*, double** %24
%390 = getelementptr inbounds double, double* %389, i32 1
store double* %390, double** %24
%391 = load double, double* %389
%392 = load double, double* %33
%393 = call double @llvm.fmuladd.f64(double %388, double %391, double %392)
store double %393, double* %33
%394 = load double, double* %29
%395 = load double*, double** %24
%396 = getelementptr inbounds double, double* %395, i32 1
store double* %396, double** %24
%397 = load double, double* %395
%398 = load double, double* %34
%399 = call double @llvm.fmuladd.f64(double %394, double %397, double %398)
store double %399, double* %34
%400 = load double, double* %30
%401 = load double*, double** %24
%402 = getelementptr inbounds double, double* %401, i32 1
store double* %402, double** %24
%403 = load double, double* %401
%404 = load double, double* %33
%405 = call double @llvm.fmuladd.f64(double %400, double %403, double %404)
store double %405, double* %33
%406 = load double, double* %30
%407 = load double*, double** %24
%408 = getelementptr inbounds double, double* %407, i32 1
store double* %408, double** %24
%409 = load double, double* %407
%410 = load double, double* %34
%411 = call double @llvm.fmuladd.f64(double %406, double %409, double %410)
store double %411, double* %34
%412 = load double, double* %31
%413 = load double*, double** %24
%414 = getelementptr inbounds double, double* %413, i32 1
store double* %414, double** %24
%415 = load double, double* %413
%416 = load double, double* %33
%417 = call double @llvm.fmuladd.f64(double %412, double %415, double %416)
store double %417, double* %33
%418 = load double, double* %31
%419 = load double*, double** %24
%420 = getelementptr inbounds double, double* %419, i32 1
store double* %420, double** %24
%421 = load double, double* %419
%422 = load double, double* %34
%423 = call double @llvm.fmuladd.f64(double %418, double %421, double %422)
store double %423, double* %34
%424 = load double*, double** %23
%425 = getelementptr inbounds double, double* %424, i32 1
store double* %425, double** %23
%426 = load double, double* %424
%427 = load double, double* %32
%428 = fadd double %427, %426
store double %428, double* %32
%429 = load double*, double** %23
%430 = getelementptr inbounds double, double* %429, i32 1
store double* %430, double** %23
%431 = load double, double* %429
%432 = load double*, double** %24
%433 = getelementptr inbounds double, double* %432, i32 1
store double* %433, double** %24
%434 = load double, double* %432
%435 = fmul double %431, %434
store double %435, double* %29
%436 = load double, double* %32
%437 = fadd double %436, %435
store double %437, double* %32
%438 = load double*, double** %23
%439 = getelementptr inbounds double, double* %438, i32 1
store double* %439, double** %23
%440 = load double, double* %438
%441 = load double*, double** %24
%442 = getelementptr inbounds double, double* %441, i32 1
store double* %442, double** %24
%443 = load double, double* %441
%444 = fmul double %440, %443
store double %444, double* %30
%445 = load double, double* %32
%446 = fadd double %445, %444
store double %446, double* %32
%447 = load double*, double** %23
%448 = getelementptr inbounds double, double* %447, i32 1
store double* %448, double** %23
%449 = load double, double* %447
%450 = load double*, double** %24
%451 = getelementptr inbounds double, double* %450, i32 1
store double* %451, double** %24
%452 = load double, double* %450
%453 = fmul double %449, %452
store double %453, double* %31
%454 = load double, double* %32
%455 = fadd double %454, %453
store double %455, double* %32
%456 = load double, double* %29
%457 = load double*, double** %24
%458 = getelementptr inbounds double, double* %457, i32 1
store double* %458, double** %24
%459 = load double, double* %457
%460 = load double, double* %33
%461 = call double @llvm.fmuladd.f64(double %456, double %459, double %460)
store double %461, double* %33
%462 = load double, double* %29
%463 = load double*, double** %24
%464 = getelementptr inbounds double, double* %463, i32 1
store double* %464, double** %24
%465 = load double, double* %463
%466 = load double, double* %34
%467 = call double @llvm.fmuladd.f64(double %462, double %465, double %466)
store double %467, double* %34
%468 = load double, double* %30
%469 = load double*, double** %24
%470 = getelementptr inbounds double, double* %469, i32 1
store double* %470, double** %24
%471 = load double, double* %469
%472 = load double, double* %33
%473 = call double @llvm.fmuladd.f64(double %468, double %471, double %472)
store double %473, double* %33
%474 = load double, double* %30
%475 = load double*, double** %24
%476 = getelementptr inbounds double, double* %475, i32 1
store double* %476, double** %24
%477 = load double, double* %475
%478 = load double, double* %34
%479 = call double @llvm.fmuladd.f64(double %474, double %477, double %478)
store double %479, double* %34
%480 = load double, double* %31
%481 = load double*, double** %24
%482 = getelementptr inbounds double, double* %481, i32 1
store double* %482, double** %24
%483 = load double, double* %481
%484 = load double, double* %33
%485 = call double @llvm.fmuladd.f64(double %480, double %483, double %484)
store double %485, double* %33
%486 = load double, double* %31
%487 = load double*, double** %24
%488 = getelementptr inbounds double, double* %487, i32 1
store double* %488, double** %24
%489 = load double, double* %487
%490 = load double, double* %34
%491 = call double @llvm.fmuladd.f64(double %486, double %489, double %490)
store double %491, double* %34
%492 = load double, double* %32
%493 = fdiv double 1.000000e+00, %492
store double %493, double* %32
%494 = load double, double* %32
%495 = load double, double* %33
%496 = fmul double %495, %494
store double %496, double* %33
%497 = load double, double* %32
%498 = load double, double* %34
%499 = fmul double %498, %497
store double %499, double* %34
%500 = load i32*, i32** %20
%501 = load i32, i32* %26
%502 = sext i32 %501 to i64
%503 = getelementptr inbounds i32, i32* %500, i64 %502
%504 = load i32, i32* %503
%505 = sitofp i32 %504 to double
%506 = load double, double* %33
%507 = load double, double* %27
%508 = call double @llvm.fmuladd.f64(double %505, double %506, double %507)
store double %508, double* %27
%509 = load i32*, i32** %20
%510 = load i32, i32* %26
%511 = sext i32 %510 to i64
%512 = getelementptr inbounds i32, i32* %509, i64 %511
%513 = load i32, i32* %512
%514 = sitofp i32 %513 to double
%515 = load double, double* %34
%516 = load double, double* %33
%517 = load double, double* %33
%518 = fneg double %516
%519 = call double @llvm.fmuladd.f64(double %518, double %517, double %515)
%520 = load double, double* %28
%521 = call double @llvm.fmuladd.f64(double %514, double %519, double %520)
store double %521, double* %28
br label %522
522: ; preds = %209
%523 = load i32, i32* %26
%524 = add nsw i32 %523, 1
store i32 %524, i32* %26
br label %205
525: ; preds = %205
%526 = load double, double* %27
%527 = load double*, double** %15
store double %526, double* %527
%528 = load double, double* %28
%529 = load double*, double** %16
store double %528, double* %529
%530 = load double*, double** %25
%531 = bitcast double* %530 to i8*
ret void
}

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