#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 FABS(double);
double exp(double);
void coreGTRGAMMAINVAR(double propInvar, double *frequencies,
                       double *gammaRates, double *EIGN, double *sumtable,
                       volatile double *ext_dlnLdlz,
                       volatile double *ext_d2lnLdlz2, int *iptr, int *wrptr,
                       int upper, double lz) {
  double *sum, diagptable[64];
  int i, j, k;
  double dlnLdlz = 0;
  double d2lnLdlz2 = 0;
  double ki, kisqr;
  double freqs[4];
  double scaler = 0.25 * (1.0 - propInvar);
  double tmp_1;
  double inv_Li, dlnLidlz, d2lnLidlz2;
  freqs[0] = frequencies[0] * propInvar;
  freqs[1] = frequencies[1] * propInvar;
  freqs[2] = frequencies[2] * propInvar;
  freqs[3] = frequencies[3] * propInvar;
  for (i = 0; i < 4; i++) {
    ki = gammaRates[i];
    kisqr = ki * ki;
    diagptable[i * 16] = exp(EIGN[0] * ki * lz);
    diagptable[i * 16 + 1] = exp(EIGN[1] * ki * lz);
    diagptable[i * 16 + 2] = exp(EIGN[2] * ki * lz);
    diagptable[i * 16 + 3] = EIGN[0] * ki;
    diagptable[i * 16 + 4] = EIGN[0] * EIGN[0] * kisqr;
    diagptable[i * 16 + 5] = EIGN[1] * ki;
    diagptable[i * 16 + 6] = EIGN[1] * EIGN[1] * kisqr;
    diagptable[i * 16 + 7] = EIGN[2] * ki;
    diagptable[i * 16 + 8] = EIGN[2] * EIGN[2] * kisqr;
  }
  for (i = 0; i < upper; i++) {
    sum = &(sumtable[i * 16]);
    inv_Li = 0.0;
    dlnLidlz = 0.0;
    d2lnLidlz2 = 0.0;
    for (j = 0; j < 4; j++) {
      inv_Li += sum[4 * j];
      for (k = 0; k < 3; k++) {
        tmp_1 = diagptable[16 * j + k] * sum[4 * j + k + 1];
        inv_Li += tmp_1;
        dlnLidlz += tmp_1 * diagptable[16 * j + k * 2 + 3];
        d2lnLidlz2 += tmp_1 * diagptable[16 * j + k * 2 + 4];
      }
    }
    inv_Li = FABS(inv_Li) * scaler;
    if (iptr[i] < 4)
      inv_Li += freqs[iptr[i]];
    inv_Li = 1.0 / inv_Li;
    dlnLidlz *= inv_Li;
    d2lnLidlz2 *= inv_Li;
    dlnLidlz *= scaler;
    d2lnLidlz2 *= scaler;
    dlnLdlz += wrptr[i] * dlnLidlz;
    d2lnLdlz2 += wrptr[i] * (d2lnLidlz2 - dlnLidlz * dlnLidlz);
  }
  *ext_dlnLdlz = dlnLdlz;
  *ext_d2lnLdlz2 = d2lnLdlz2;
}


x86 assembly (clang, -O0):
<coreGTRGAMMAINVAR>:
push   %rbp
mov    %rsp,%rbp
sub    $0x2d0,%rsp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
movsd  %xmm1,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x268(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x270(%rbp)
ff
movsd  0xe92(%rip),%xmm1
subsd  -0x8(%rbp),%xmm1
movsd  0xe7d(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x2a8(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x2a0(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x298(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x290(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x288(%rbp)
ff
movl   $0x0,-0x254(%rbp)
cmpl   $0x4,-0x254(%rbp)
jge    13db <coreGTRGAMMAINVAR+0x2bb>
mov    -0x18(%rbp),%rax
movslq -0x254(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x278(%rbp)
ff
movsd  -0x278(%rbp),%xmm0
ff
mulsd  -0x278(%rbp),%xmm0
ff
movsd  %xmm0,-0x280(%rbp)
ff
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  -0x278(%rbp),%xmm0
ff
mulsd  -0x40(%rbp),%xmm0
call   1040 <exp@plt>
mov    -0x254(%rbp),%eax
shl    $0x4,%eax
cltq
movsd  %xmm0,-0x250(%rbp,%rax,8)
ff ff
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
mulsd  -0x278(%rbp),%xmm0
ff
mulsd  -0x40(%rbp),%xmm0
call   1040 <exp@plt>
mov    -0x254(%rbp),%eax
shl    $0x4,%eax
add    $0x1,%eax
cltq
movsd  %xmm0,-0x250(%rbp,%rax,8)
ff ff
mov    -0x20(%rbp),%rax
movsd  0x10(%rax),%xmm0
mulsd  -0x278(%rbp),%xmm0
ff
mulsd  -0x40(%rbp),%xmm0
call   1040 <exp@plt>
mov    -0x254(%rbp),%eax
shl    $0x4,%eax
add    $0x2,%eax
cltq
movsd  %xmm0,-0x250(%rbp,%rax,8)
ff ff
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  -0x278(%rbp),%xmm0
ff
mov    -0x254(%rbp),%eax
shl    $0x4,%eax
add    $0x3,%eax
cltq
movsd  %xmm0,-0x250(%rbp,%rax,8)
ff ff
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
mulsd  (%rax),%xmm0
mulsd  -0x280(%rbp),%xmm0
ff
mov    -0x254(%rbp),%eax
shl    $0x4,%eax
add    $0x4,%eax
cltq
movsd  %xmm0,-0x250(%rbp,%rax,8)
ff ff
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
mulsd  -0x278(%rbp),%xmm0
ff
mov    -0x254(%rbp),%eax
shl    $0x4,%eax
add    $0x5,%eax
cltq
movsd  %xmm0,-0x250(%rbp,%rax,8)
ff ff
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x20(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mulsd  -0x280(%rbp),%xmm0
ff
mov    -0x254(%rbp),%eax
shl    $0x4,%eax
add    $0x6,%eax
cltq
movsd  %xmm0,-0x250(%rbp,%rax,8)
ff ff
mov    -0x20(%rbp),%rax
movsd  0x10(%rax),%xmm0
mulsd  -0x278(%rbp),%xmm0
ff
mov    -0x254(%rbp),%eax
shl    $0x4,%eax
add    $0x7,%eax
cltq
movsd  %xmm0,-0x250(%rbp,%rax,8)
ff ff
mov    -0x20(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x20(%rbp),%rax
mulsd  0x10(%rax),%xmm0
mulsd  -0x280(%rbp),%xmm0
ff
mov    -0x254(%rbp),%eax
shl    $0x4,%eax
add    $0x8,%eax
cltq
movsd  %xmm0,-0x250(%rbp,%rax,8)
ff ff
mov    -0x254(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x254(%rbp)
jmp    11f0 <coreGTRGAMMAINVAR+0xd0>
movl   $0x0,-0x254(%rbp)
mov    -0x254(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    1701 <coreGTRGAMMAINVAR+0x5e1>
mov    -0x28(%rbp),%rax
mov    -0x254(%rbp),%ecx
shl    $0x4,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2b8(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2c0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2c8(%rbp)
ff
movl   $0x0,-0x258(%rbp)
cmpl   $0x4,-0x258(%rbp)
jge    158c <coreGTRGAMMAINVAR+0x46c>
mov    -0x48(%rbp),%rax
mov    -0x258(%rbp),%ecx
shl    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x2b8(%rbp),%xmm0
ff
movsd  %xmm0,-0x2b8(%rbp)
ff
movl   $0x0,-0x25c(%rbp)
cmpl   $0x3,-0x25c(%rbp)
jge    1573 <coreGTRGAMMAINVAR+0x453>
mov    -0x258(%rbp),%eax
shl    $0x4,%eax
add    -0x25c(%rbp),%eax
cltq
movsd  -0x250(%rbp,%rax,8),%xmm0
ff ff
mov    -0x48(%rbp),%rax
mov    -0x258(%rbp),%ecx
shl    $0x2,%ecx
add    -0x25c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x2b0(%rbp)
ff
movsd  -0x2b0(%rbp),%xmm0
ff
addsd  -0x2b8(%rbp),%xmm0
ff
movsd  %xmm0,-0x2b8(%rbp)
ff
movsd  -0x2b0(%rbp),%xmm0
ff
mov    -0x258(%rbp),%eax
shl    $0x4,%eax
mov    -0x25c(%rbp),%ecx
shl    $0x1,%ecx
add    %ecx,%eax
add    $0x3,%eax
cltq
movsd  -0x250(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x2c0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x2c0(%rbp)
ff
movsd  -0x2b0(%rbp),%xmm0
ff
mov    -0x258(%rbp),%eax
shl    $0x4,%eax
mov    -0x25c(%rbp),%ecx
shl    $0x1,%ecx
add    %ecx,%eax
add    $0x4,%eax
cltq
movsd  -0x250(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x2c8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x2c8(%rbp)
ff
mov    -0x25c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x25c(%rbp)
jmp    1476 <coreGTRGAMMAINVAR+0x356>
jmp    1578 <coreGTRGAMMAINVAR+0x458>
mov    -0x258(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x258(%rbp)
jmp    143a <coreGTRGAMMAINVAR+0x31a>
movsd  -0x2b8(%rbp),%xmm0
ff
call   1030 <FABS@plt>
mulsd  -0x2a8(%rbp),%xmm0
ff
movsd  %xmm0,-0x2b8(%rbp)
ff
mov    0x10(%rbp),%rax
movslq -0x254(%rbp),%rcx
cmpl   $0x4,(%rax,%rcx,4)
jge    15e6 <coreGTRGAMMAINVAR+0x4c6>
mov    0x10(%rbp),%rax
movslq -0x254(%rbp),%rcx
movslq (%rax,%rcx,4),%rax
movsd  -0x2a0(%rbp,%rax,8),%xmm0
ff ff
addsd  -0x2b8(%rbp),%xmm0
ff
movsd  %xmm0,-0x2b8(%rbp)
ff
movsd  0xa1a(%rip),%xmm0
divsd  -0x2b8(%rbp),%xmm0
ff
movsd  %xmm0,-0x2b8(%rbp)
ff
movsd  -0x2b8(%rbp),%xmm0
ff
mulsd  -0x2c0(%rbp),%xmm0
ff
movsd  %xmm0,-0x2c0(%rbp)
ff
movsd  -0x2b8(%rbp),%xmm0
ff
mulsd  -0x2c8(%rbp),%xmm0
ff
movsd  %xmm0,-0x2c8(%rbp)
ff
movsd  -0x2a8(%rbp),%xmm0
ff
mulsd  -0x2c0(%rbp),%xmm0
ff
movsd  %xmm0,-0x2c0(%rbp)
ff
movsd  -0x2a8(%rbp),%xmm0
ff
mulsd  -0x2c8(%rbp),%xmm0
ff
movsd  %xmm0,-0x2c8(%rbp)
ff
mov    0x18(%rbp),%rax
movslq -0x254(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x2c0(%rbp),%xmm2
ff
movsd  -0x268(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x268(%rbp)
ff
mov    0x18(%rbp),%rax
movslq -0x254(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x2c8(%rbp),%xmm1
ff
movsd  -0x2c0(%rbp),%xmm2
ff
movsd  -0x2c0(%rbp),%xmm3
ff
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  -0x270(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x270(%rbp)
ff
mov    -0x254(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x254(%rbp)
jmp    13e5 <coreGTRGAMMAINVAR+0x2c5>
movsd  -0x268(%rbp),%xmm0
ff
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x270(%rbp),%xmm0
ff
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x2d0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<coreGTRGAMMAINVAR>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x408,%rsp
movsd  %xmm1,0x10(%rsp)
mov    %r9,0x48(%rsp)
mov    %r8,0x40(%rsp)
mov    %rcx,%r12
mov    %rdx,%rbx
mov    %rsi,%r13
mov    0x440(%rsp),%r14
movsd  0xea6(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
movupd (%rdi),%xmm1
unpcklpd %xmm0,%xmm0
mulpd  %xmm0,%xmm1
movapd %xmm1,0x1e0(%rsp)
movupd 0x10(%rdi),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0x1f0(%rsp)
movsd  (%rdx),%xmm0
lea    0x200(%rsp),%rbp
xor    %r15d,%r15d
nopw   0x0(%rax,%rax,1)
movsd  0x0(%r13,%r15,8),%xmm1
movapd %xmm1,0x20(%rsp)
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
movapd %xmm2,0x30(%rsp)
mulsd  %xmm1,%xmm0
mulsd  0x10(%rsp),%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0x0(%rbp)
movsd  0x8(%rbx),%xmm0
mulsd  0x20(%rsp),%xmm0
mulsd  0x10(%rsp),%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0x8(%rbp)
movsd  0x10(%rbx),%xmm0
mulsd  0x20(%rsp),%xmm0
mulsd  0x10(%rsp),%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0x10(%rbp)
movsd  (%rbx),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
movapd 0x20(%rsp),%xmm3
unpcklpd 0x30(%rsp),%xmm3
movapd %xmm0,%xmm2
unpcklpd %xmm1,%xmm2
mulpd  %xmm3,%xmm2
movupd %xmm2,0x18(%rbp)
movsd  0x8(%rbx),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
unpcklpd %xmm2,%xmm1
mulpd  %xmm3,%xmm1
movupd %xmm1,0x28(%rbp)
movsd  0x10(%rbx),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
unpcklpd %xmm2,%xmm1
mulpd  %xmm3,%xmm1
movupd %xmm1,0x38(%rbp)
add    $0x1,%r15
sub    $0xffffffffffffff80,%rbp
cmp    $0x4,%r15
jne    11a0 <coreGTRGAMMAINVAR+0x80>
mov    0x450(%rsp),%eax
test   %eax,%eax
jle    1894 <coreGTRGAMMAINVAR+0x774>
movsd  0x8(%rsp),%xmm3
mulsd  0xd77(%rip),%xmm3
mov    %eax,%ebx
movsd  0x2b0(%rsp),%xmm0
movsd  %xmm0,0x138(%rsp)
movsd  0x2b8(%rsp),%xmm0
movsd  %xmm0,0x130(%rsp)
movsd  0x2c0(%rsp),%xmm0
movsd  %xmm0,0x128(%rsp)
movsd  0x300(%rsp),%xmm0
movsd  %xmm0,0x120(%rsp)
movsd  0x318(%rsp),%xmm0
movsd  %xmm0,0x118(%rsp)
movsd  0x320(%rsp),%xmm0
movsd  %xmm0,0x110(%rsp)
movsd  0x328(%rsp),%xmm0
movsd  %xmm0,0x108(%rsp)
movsd  0x330(%rsp),%xmm0
movsd  %xmm0,0x100(%rsp)
movsd  0x200(%rsp),%xmm0
movsd  %xmm0,0xf8(%rsp)
movsd  0x208(%rsp),%xmm0
movsd  %xmm0,0xf0(%rsp)
movsd  0x218(%rsp),%xmm0
movsd  %xmm0,0xe8(%rsp)
movsd  0x220(%rsp),%xmm0
movsd  %xmm0,0xe0(%rsp)
movsd  0x228(%rsp),%xmm0
movsd  %xmm0,0xd8(%rsp)
movsd  0x230(%rsp),%xmm0
movsd  %xmm0,0xd0(%rsp)
movsd  0x210(%rsp),%xmm0
movsd  %xmm0,0xc8(%rsp)
movsd  0x238(%rsp),%xmm0
movsd  %xmm0,0xc0(%rsp)
add    $0x70,%r12
xor    %ebp,%ebp
movsd  0x240(%rsp),%xmm0
movsd  %xmm0,0xb8(%rsp)
movsd  0x280(%rsp),%xmm0
movsd  %xmm0,0xb0(%rsp)
movsd  0x298(%rsp),%xmm0
movsd  %xmm0,0xa8(%rsp)
movsd  0x2a0(%rsp),%xmm0
movsd  %xmm0,0xa0(%rsp)
movups 0x288(%rsp),%xmm0
movaps %xmm0,0x190(%rsp)
movsd  0x2a8(%rsp),%xmm0
movsd  %xmm0,0x98(%rsp)
movups 0x308(%rsp),%xmm0
movaps %xmm0,0x180(%rsp)
movsd  0x338(%rsp),%xmm0
movsd  %xmm0,0x90(%rsp)
movsd  0x340(%rsp),%xmm0
movsd  %xmm0,0x88(%rsp)
movsd  0x380(%rsp),%xmm0
movsd  %xmm0,0x80(%rsp)
movsd  0x398(%rsp),%xmm0
movsd  %xmm0,0x78(%rsp)
movsd  0x3a0(%rsp),%xmm0
movsd  %xmm0,0x70(%rsp)
movsd  0x3a8(%rsp),%xmm0
movsd  %xmm0,0x68(%rsp)
movsd  0x3b0(%rsp),%xmm0
movsd  %xmm0,0x60(%rsp)
movups 0x388(%rsp),%xmm0
movaps %xmm0,0x170(%rsp)
movsd  0x3b8(%rsp),%xmm0
movsd  %xmm0,0x58(%rsp)
movsd  0x3c0(%rsp),%xmm0
movsd  %xmm0,0x50(%rsp)
xorpd  %xmm0,%xmm0
xorpd  %xmm4,%xmm4
mov    0x448(%rsp),%r15
movsd  %xmm3,0x8(%rsp)
jmp    170c <coreGTRGAMMAINVAR+0x5ec>
xchg   %ax,%ax
movsd  0x140(%rsp),%xmm5
movapd %xmm5,%xmm1
mulsd  0xe8(%rsp),%xmm1
xorpd  %xmm4,%xmm4
addsd  %xmm4,%xmm1
movsd  0x148(%rsp),%xmm6
movapd %xmm6,%xmm2
mulsd  0xd8(%rsp),%xmm2
addsd  %xmm1,%xmm2
mulsd  0xe0(%rsp),%xmm5
addsd  %xmm4,%xmm5
mulsd  0xd0(%rsp),%xmm6
addsd  %xmm5,%xmm6
movsd  0x150(%rsp),%xmm5
movapd %xmm5,%xmm1
mulsd  0xc0(%rsp),%xmm1
addsd  %xmm2,%xmm1
mulsd  0xb8(%rsp),%xmm5
addsd  %xmm6,%xmm5
movsd  0x158(%rsp),%xmm6
movapd %xmm6,%xmm2
mulsd  0xa8(%rsp),%xmm2
addsd  %xmm1,%xmm2
mulsd  0xa0(%rsp),%xmm6
addsd  %xmm5,%xmm6
movapd 0x1b0(%rsp),%xmm5
movapd %xmm5,%xmm1
mulsd  0x98(%rsp),%xmm1
addsd  %xmm2,%xmm1
mulsd  0x138(%rsp),%xmm5
addsd  %xmm6,%xmm5
movapd 0x1a0(%rsp),%xmm7
movapd %xmm7,%xmm2
mulsd  0x130(%rsp),%xmm2
addsd  %xmm1,%xmm2
mulsd  0x128(%rsp),%xmm7
addsd  %xmm5,%xmm7
movsd  0x160(%rsp),%xmm6
movapd %xmm6,%xmm1
mulsd  0x118(%rsp),%xmm1
addsd  %xmm2,%xmm1
mulsd  0x110(%rsp),%xmm6
addsd  %xmm7,%xmm6
movapd 0x30(%rsp),%xmm5
movapd %xmm5,%xmm2
mulsd  0x108(%rsp),%xmm2
addsd  %xmm1,%xmm2
mulsd  0x100(%rsp),%xmm5
addsd  %xmm6,%xmm5
movapd 0x1c0(%rsp),%xmm7
movapd %xmm7,%xmm1
mulsd  0x90(%rsp),%xmm1
addsd  %xmm2,%xmm1
mulsd  0x88(%rsp),%xmm7
addsd  %xmm5,%xmm7
movsd  0x168(%rsp),%xmm6
movapd %xmm6,%xmm2
mulsd  0x78(%rsp),%xmm2
addsd  %xmm1,%xmm2
mulsd  0x70(%rsp),%xmm6
addsd  %xmm7,%xmm6
movapd 0x10(%rsp),%xmm5
movapd %xmm5,%xmm1
mulsd  0x68(%rsp),%xmm1
addsd  %xmm2,%xmm1
mulsd  0x60(%rsp),%xmm5
addsd  %xmm6,%xmm5
movapd 0x1d0(%rsp),%xmm6
movapd %xmm6,%xmm2
mulsd  0x58(%rsp),%xmm2
addsd  %xmm1,%xmm2
mulsd  0x50(%rsp),%xmm6
addsd  %xmm5,%xmm6
movsd  0x948(%rip),%xmm1
divsd  %xmm0,%xmm1
mulsd  %xmm1,%xmm2
mulsd  %xmm6,%xmm1
mulsd  %xmm3,%xmm2
mulsd  %xmm3,%xmm1
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
subsd  %xmm0,%xmm1
unpcklpd %xmm2,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sdl (%r15,%rbp,4),%xmm0
unpcklpd %xmm0,%xmm0
mulpd  %xmm0,%xmm1
movapd %xmm1,%xmm0
movapd 0x20(%rsp),%xmm2
addpd  %xmm2,%xmm0
add    $0x1,%rbp
sub    $0xffffffffffffff80,%r12
cmp    %rbp,%rbx
je     188e <coreGTRGAMMAINVAR+0x76e>
movapd %xmm0,0x20(%rsp)
movsd  -0x70(%r12),%xmm0
addsd  %xmm4,%xmm0
movsd  -0x68(%r12),%xmm1
mulsd  0xf8(%rsp),%xmm1
movsd  %xmm1,0x140(%rsp)
addsd  %xmm1,%xmm0
movsd  -0x60(%r12),%xmm1
mulsd  0xf0(%rsp),%xmm1
movsd  %xmm1,0x148(%rsp)
addsd  %xmm1,%xmm0
movsd  -0x58(%r12),%xmm1
mulsd  0xc8(%rsp),%xmm1
movsd  %xmm1,0x150(%rsp)
addsd  %xmm1,%xmm0
addsd  -0x50(%r12),%xmm0
movsd  -0x48(%r12),%xmm1
mulsd  0xb0(%rsp),%xmm1
movsd  %xmm1,0x158(%rsp)
addsd  %xmm1,%xmm0
movupd -0x40(%r12),%xmm3
mulpd  0x190(%rsp),%xmm3
movupd -0x20(%r12),%xmm2
mulpd  0x180(%rsp),%xmm2
movupd (%r12),%xmm1
mulpd  0x170(%rsp),%xmm1
addsd  %xmm3,%xmm0
movapd %xmm3,0x1b0(%rsp)
unpckhpd %xmm3,%xmm3
movapd %xmm3,0x1a0(%rsp)
addsd  %xmm3,%xmm0
addsd  -0x30(%r12),%xmm0
movsd  -0x28(%r12),%xmm3
mulsd  0x120(%rsp),%xmm3
movsd  %xmm3,0x160(%rsp)
addsd  %xmm3,%xmm0
addsd  %xmm2,%xmm0
movapd %xmm2,0x30(%rsp)
unpckhpd %xmm2,%xmm2
movapd %xmm2,0x1c0(%rsp)
addsd  %xmm2,%xmm0
addsd  -0x10(%r12),%xmm0
movsd  -0x8(%r12),%xmm2
mulsd  0x80(%rsp),%xmm2
movsd  %xmm2,0x168(%rsp)
addsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movapd %xmm1,0x10(%rsp)
unpckhpd %xmm1,%xmm1
movapd %xmm1,0x1d0(%rsp)
addsd  %xmm1,%xmm0
call   1030 <FABS@plt>
movsd  0x8(%rsp),%xmm3
mulsd  %xmm3,%xmm0
movslq (%r14,%rbp,4),%rax
cmp    $0x3,%rax
jg     14f0 <coreGTRGAMMAINVAR+0x3d0>
addsd  0x1e0(%rsp,%rax,8),%xmm0
jmp    14f0 <coreGTRGAMMAINVAR+0x3d0>
addpd  %xmm2,%xmm1
jmp    1898 <coreGTRGAMMAINVAR+0x778>
xorpd  %xmm1,%xmm1
mov    0x40(%rsp),%rax
movhpd %xmm1,(%rax)
mov    0x48(%rsp),%rax
movlpd %xmm1,(%rax)
add    $0x408,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @coreGTRGAMMAINVAR(double %0, double* %1, double* %2, double* %3, double* %4, double* %5, double* %6, i32* %7, i32* %8, i32 %9, double %10) {
%12 = alloca double
%13 = alloca double*
%14 = alloca double*
%15 = alloca double*
%16 = alloca double*
%17 = alloca double*
%18 = alloca double*
%19 = alloca i32*
%20 = alloca i32*
%21 = alloca i32
%22 = alloca double
%23 = alloca double*
%24 = alloca [64 x double]
%25 = alloca i32
%26 = alloca i32
%27 = alloca i32
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca [4 x double]
%33 = alloca double
%34 = alloca double
%35 = alloca double
%36 = alloca double
%37 = alloca double
store double %0, double* %12
store double* %1, double** %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 i32* %7, i32** %19
store i32* %8, i32** %20
store i32 %9, i32* %21
store double %10, double* %22
store double 0.000000e+00, double* %28
store double 0.000000e+00, double* %29
%38 = load double, double* %12
%39 = fsub double 1.000000e+00, %38
%40 = fmul double 2.500000e-01, %39
store double %40, double* %33
%41 = load double*, double** %13
%42 = getelementptr inbounds double, double* %41, i64 0
%43 = load double, double* %42
%44 = load double, double* %12
%45 = fmul double %43, %44
%46 = getelementptr inbounds [4 x double], [4 x double]* %32, i64 0, i64 0
store double %45, double* %46
%47 = load double*, double** %13
%48 = getelementptr inbounds double, double* %47, i64 1
%49 = load double, double* %48
%50 = load double, double* %12
%51 = fmul double %49, %50
%52 = getelementptr inbounds [4 x double], [4 x double]* %32, i64 0, i64 1
store double %51, double* %52
%53 = load double*, double** %13
%54 = getelementptr inbounds double, double* %53, i64 2
%55 = load double, double* %54
%56 = load double, double* %12
%57 = fmul double %55, %56
%58 = getelementptr inbounds [4 x double], [4 x double]* %32, i64 0, i64 2
store double %57, double* %58
%59 = load double*, double** %13
%60 = getelementptr inbounds double, double* %59, i64 3
%61 = load double, double* %60
%62 = load double, double* %12
%63 = fmul double %61, %62
%64 = getelementptr inbounds [4 x double], [4 x double]* %32, i64 0, i64 3
store double %63, double* %64
store i32 0, i32* %25
br label %65
65: ; preds = %187, %11
%66 = load i32, i32* %25
%67 = icmp slt i32 %66, 4
br i1 %67, label %68, label %190
68: ; preds = %65
%69 = load double*, double** %14
%70 = load i32, i32* %25
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds double, double* %69, i64 %71
%73 = load double, double* %72
store double %73, double* %30
%74 = load double, double* %30
%75 = load double, double* %30
%76 = fmul double %74, %75
store double %76, double* %31
%77 = load double*, double** %15
%78 = getelementptr inbounds double, double* %77, i64 0
%79 = load double, double* %78
%80 = load double, double* %30
%81 = fmul double %79, %80
%82 = load double, double* %22
%83 = fmul double %81, %82
%84 = call double @exp(double %83)
%85 = load i32, i32* %25
%86 = mul nsw i32 %85, 16
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %87
store double %84, double* %88
%89 = load double*, double** %15
%90 = getelementptr inbounds double, double* %89, i64 1
%91 = load double, double* %90
%92 = load double, double* %30
%93 = fmul double %91, %92
%94 = load double, double* %22
%95 = fmul double %93, %94
%96 = call double @exp(double %95)
%97 = load i32, i32* %25
%98 = mul nsw i32 %97, 16
%99 = add nsw i32 %98, 1
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %100
store double %96, double* %101
%102 = load double*, double** %15
%103 = getelementptr inbounds double, double* %102, i64 2
%104 = load double, double* %103
%105 = load double, double* %30
%106 = fmul double %104, %105
%107 = load double, double* %22
%108 = fmul double %106, %107
%109 = call double @exp(double %108)
%110 = load i32, i32* %25
%111 = mul nsw i32 %110, 16
%112 = add nsw i32 %111, 2
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %113
store double %109, double* %114
%115 = load double*, double** %15
%116 = getelementptr inbounds double, double* %115, i64 0
%117 = load double, double* %116
%118 = load double, double* %30
%119 = fmul double %117, %118
%120 = load i32, i32* %25
%121 = mul nsw i32 %120, 16
%122 = add nsw i32 %121, 3
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %123
store double %119, double* %124
%125 = load double*, double** %15
%126 = getelementptr inbounds double, double* %125, i64 0
%127 = load double, double* %126
%128 = load double*, double** %15
%129 = getelementptr inbounds double, double* %128, i64 0
%130 = load double, double* %129
%131 = fmul double %127, %130
%132 = load double, double* %31
%133 = fmul double %131, %132
%134 = load i32, i32* %25
%135 = mul nsw i32 %134, 16
%136 = add nsw i32 %135, 4
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %137
store double %133, double* %138
%139 = load double*, double** %15
%140 = getelementptr inbounds double, double* %139, i64 1
%141 = load double, double* %140
%142 = load double, double* %30
%143 = fmul double %141, %142
%144 = load i32, i32* %25
%145 = mul nsw i32 %144, 16
%146 = add nsw i32 %145, 5
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %147
store double %143, double* %148
%149 = load double*, double** %15
%150 = getelementptr inbounds double, double* %149, i64 1
%151 = load double, double* %150
%152 = load double*, double** %15
%153 = getelementptr inbounds double, double* %152, i64 1
%154 = load double, double* %153
%155 = fmul double %151, %154
%156 = load double, double* %31
%157 = fmul double %155, %156
%158 = load i32, i32* %25
%159 = mul nsw i32 %158, 16
%160 = add nsw i32 %159, 6
%161 = sext i32 %160 to i64
%162 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %161
store double %157, double* %162
%163 = load double*, double** %15
%164 = getelementptr inbounds double, double* %163, i64 2
%165 = load double, double* %164
%166 = load double, double* %30
%167 = fmul double %165, %166
%168 = load i32, i32* %25
%169 = mul nsw i32 %168, 16
%170 = add nsw i32 %169, 7
%171 = sext i32 %170 to i64
%172 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %171
store double %167, double* %172
%173 = load double*, double** %15
%174 = getelementptr inbounds double, double* %173, i64 2
%175 = load double, double* %174
%176 = load double*, double** %15
%177 = getelementptr inbounds double, double* %176, i64 2
%178 = load double, double* %177
%179 = fmul double %175, %178
%180 = load double, double* %31
%181 = fmul double %179, %180
%182 = load i32, i32* %25
%183 = mul nsw i32 %182, 16
%184 = add nsw i32 %183, 8
%185 = sext i32 %184 to i64
%186 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %185
store double %181, double* %186
br label %187
187: ; preds = %68
%188 = load i32, i32* %25
%189 = add nsw i32 %188, 1
store i32 %189, i32* %25
br label %65
190: ; preds = %65
store i32 0, i32* %25
br label %191
191: ; preds = %327, %190
%192 = load i32, i32* %25
%193 = load i32, i32* %21
%194 = icmp slt i32 %192, %193
br i1 %194, label %195, label %330
195: ; preds = %191
%196 = load double*, double** %16
%197 = load i32, i32* %25
%198 = mul nsw i32 %197, 16
%199 = sext i32 %198 to i64
%200 = getelementptr inbounds double, double* %196, i64 %199
store double* %200, double** %23
store double 0.000000e+00, double* %35
store double 0.000000e+00, double* %36
store double 0.000000e+00, double* %37
store i32 0, i32* %26
br label %201
201: ; preds = %265, %195
%202 = load i32, i32* %26
%203 = icmp slt i32 %202, 4
br i1 %203, label %204, label %268
204: ; preds = %201
%205 = load double*, double** %23
%206 = load i32, i32* %26
%207 = mul nsw i32 4, %206
%208 = sext i32 %207 to i64
%209 = getelementptr inbounds double, double* %205, i64 %208
%210 = load double, double* %209
%211 = load double, double* %35
%212 = fadd double %211, %210
store double %212, double* %35
store i32 0, i32* %27
br label %213
213: ; preds = %261, %204
%214 = load i32, i32* %27
%215 = icmp slt i32 %214, 3
br i1 %215, label %216, label %264
216: ; preds = %213
%217 = load i32, i32* %26
%218 = mul nsw i32 16, %217
%219 = load i32, i32* %27
%220 = add nsw i32 %218, %219
%221 = sext i32 %220 to i64
%222 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %221
%223 = load double, double* %222
%224 = load double*, double** %23
%225 = load i32, i32* %26
%226 = mul nsw i32 4, %225
%227 = load i32, i32* %27
%228 = add nsw i32 %226, %227
%229 = add nsw i32 %228, 1
%230 = sext i32 %229 to i64
%231 = getelementptr inbounds double, double* %224, i64 %230
%232 = load double, double* %231
%233 = fmul double %223, %232
store double %233, double* %34
%234 = load double, double* %34
%235 = load double, double* %35
%236 = fadd double %235, %234
store double %236, double* %35
%237 = load double, double* %34
%238 = load i32, i32* %26
%239 = mul nsw i32 16, %238
%240 = load i32, i32* %27
%241 = mul nsw i32 %240, 2
%242 = add nsw i32 %239, %241
%243 = add nsw i32 %242, 3
%244 = sext i32 %243 to i64
%245 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %244
%246 = load double, double* %245
%247 = load double, double* %36
%248 = call double @llvm.fmuladd.f64(double %237, double %246, double %247)
store double %248, double* %36
%249 = load double, double* %34
%250 = load i32, i32* %26
%251 = mul nsw i32 16, %250
%252 = load i32, i32* %27
%253 = mul nsw i32 %252, 2
%254 = add nsw i32 %251, %253
%255 = add nsw i32 %254, 4
%256 = sext i32 %255 to i64
%257 = getelementptr inbounds [64 x double], [64 x double]* %24, i64 0, i64 %256
%258 = load double, double* %257
%259 = load double, double* %37
%260 = call double @llvm.fmuladd.f64(double %249, double %258, double %259)
store double %260, double* %37
br label %261
261: ; preds = %216
%262 = load i32, i32* %27
%263 = add nsw i32 %262, 1
store i32 %263, i32* %27
br label %213
264: ; preds = %213
br label %265
265: ; preds = %264
%266 = load i32, i32* %26
%267 = add nsw i32 %266, 1
store i32 %267, i32* %26
br label %201
268: ; preds = %201
%269 = load double, double* %35
%270 = call double @FABS(double %269)
%271 = load double, double* %33
%272 = fmul double %270, %271
store double %272, double* %35
%273 = load i32*, i32** %19
%274 = load i32, i32* %25
%275 = sext i32 %274 to i64
%276 = getelementptr inbounds i32, i32* %273, i64 %275
%277 = load i32, i32* %276
%278 = icmp slt i32 %277, 4
br i1 %278, label %279, label %290
279: ; preds = %268
%280 = load i32*, i32** %19
%281 = load i32, i32* %25
%282 = sext i32 %281 to i64
%283 = getelementptr inbounds i32, i32* %280, i64 %282
%284 = load i32, i32* %283
%285 = sext i32 %284 to i64
%286 = getelementptr inbounds [4 x double], [4 x double]* %32, i64 0, i64 %285
%287 = load double, double* %286
%288 = load double, double* %35
%289 = fadd double %288, %287
store double %289, double* %35
br label %290
290: ; preds = %279, %268
%291 = load double, double* %35
%292 = fdiv double 1.000000e+00, %291
store double %292, double* %35
%293 = load double, double* %35
%294 = load double, double* %36
%295 = fmul double %294, %293
store double %295, double* %36
%296 = load double, double* %35
%297 = load double, double* %37
%298 = fmul double %297, %296
store double %298, double* %37
%299 = load double, double* %33
%300 = load double, double* %36
%301 = fmul double %300, %299
store double %301, double* %36
%302 = load double, double* %33
%303 = load double, double* %37
%304 = fmul double %303, %302
store double %304, double* %37
%305 = load i32*, i32** %20
%306 = load i32, i32* %25
%307 = sext i32 %306 to i64
%308 = getelementptr inbounds i32, i32* %305, i64 %307
%309 = load i32, i32* %308
%310 = sitofp i32 %309 to double
%311 = load double, double* %36
%312 = load double, double* %28
%313 = call double @llvm.fmuladd.f64(double %310, double %311, double %312)
store double %313, double* %28
%314 = load i32*, i32** %20
%315 = load i32, i32* %25
%316 = sext i32 %315 to i64
%317 = getelementptr inbounds i32, i32* %314, i64 %316
%318 = load i32, i32* %317
%319 = sitofp i32 %318 to double
%320 = load double, double* %37
%321 = load double, double* %36
%322 = load double, double* %36
%323 = fneg double %321
%324 = call double @llvm.fmuladd.f64(double %323, double %322, double %320)
%325 = load double, double* %29
%326 = call double @llvm.fmuladd.f64(double %319, double %324, double %325)
store double %326, double* %29
br label %327
327: ; preds = %290
%328 = load i32, i32* %25
%329 = add nsw i32 %328, 1
store i32 %329, i32* %25
br label %191
330: ; preds = %191
%331 = load double, double* %28
%332 = load double*, double** %17
store volatile double %331, double* %332
%333 = load double, double* %29
%334 = load double*, double** %18
store volatile double %333, double* %334
ret void
}

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