#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 coreGTRGAMMASECONDARYINVAR(double *gammaRates, double *EIGN,
                                double *sumtable, int upper, int *wrptr,
                                volatile double *ext_dlnLdlz,
                                volatile double *ext_d2lnLdlz2, double lz,
                                double *frequencies, double propInvar,
                                int *iptr) {
  double *sum, diagptable[256];
  int i, l, j;
  double dlnLdlz = 0;
  double d2lnLdlz2 = 0;
  double ki, kisqr;
  double freqs[16];
  double scaler = 0.25 * (1.0 - propInvar);
  double tmp;
  double inv_Li, dlnLidlz, d2lnLidlz2;
  for (i = 0; i < 16; i++)
    freqs[i] = frequencies[i] * propInvar;
  for (i = 0; i < 4; i++) {
    ki = gammaRates[i];
    kisqr = ki * ki;
    for (l = 1; l < 16; l++) {
      diagptable[i * 64 + l * 4] = exp(EIGN[l - 1] * ki * lz);
      diagptable[i * 64 + l * 4 + 1] = EIGN[l - 1] * ki;
      diagptable[i * 64 + l * 4 + 2] = EIGN[l - 1] * EIGN[l - 1] * kisqr;
    }
  }
  for (i = 0; i < upper; i++) {
    sum = &sumtable[i * 64];
    inv_Li = 0.0;
    dlnLidlz = 0.0;
    d2lnLidlz2 = 0.0;
    for (j = 0; j < 4; j++) {
      inv_Li += sum[j * 16];
      for (l = 1; l < 16; l++) {
        inv_Li += (tmp = diagptable[j * 64 + l * 4] * sum[j * 16 + l]);
        dlnLidlz += tmp * diagptable[j * 64 + l * 4 + 1];
        d2lnLidlz2 += tmp * diagptable[j * 64 + l * 4 + 2];
      }
    }
    inv_Li = FABS(inv_Li) * scaler;
    if (iptr[i] < 16)
      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):
<coreGTRGAMMASECONDARYINVAR>:
push   %rbp
mov    %rsp,%rbp
sub    $0x930,%rsp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movsd  %xmm0,-0x38(%rbp)
movsd  %xmm1,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x868(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x870(%rbp)
ff
movsd  0xe92(%rip),%xmm1
subsd  -0x40(%rbp),%xmm1
movsd  0xe7d(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x908(%rbp)
ff
movl   $0x0,-0x854(%rbp)
cmpl   $0x10,-0x854(%rbp)
jge    11df <coreGTRGAMMASECONDARYINVAR+0xbf>
mov    0x18(%rbp),%rax
movslq -0x854(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x40(%rbp),%xmm0
movslq -0x854(%rbp),%rax
movsd  %xmm0,-0x900(%rbp,%rax,8)
ff ff
mov    -0x854(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x854(%rbp)
jmp    1199 <coreGTRGAMMASECONDARYINVAR+0x79>
movl   $0x0,-0x854(%rbp)
cmpl   $0x4,-0x854(%rbp)
jge    1343 <coreGTRGAMMASECONDARYINVAR+0x223>
mov    -0x8(%rbp),%rax
movslq -0x854(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x878(%rbp)
ff
movsd  -0x878(%rbp),%xmm0
ff
mulsd  -0x878(%rbp),%xmm0
ff
movsd  %xmm0,-0x880(%rbp)
ff
movl   $0x1,-0x858(%rbp)
cmpl   $0x10,-0x858(%rbp)
jge    132a <coreGTRGAMMASECONDARYINVAR+0x20a>
mov    -0x10(%rbp),%rax
mov    -0x858(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x878(%rbp),%xmm0
ff
mulsd  -0x38(%rbp),%xmm0
call   1040 <exp@plt>
mov    -0x854(%rbp),%eax
shl    $0x6,%eax
mov    -0x858(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
cltq
movsd  %xmm0,-0x850(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%rax
mov    -0x858(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x878(%rbp),%xmm0
ff
mov    -0x854(%rbp),%eax
shl    $0x6,%eax
mov    -0x858(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x1,%eax
cltq
movsd  %xmm0,-0x850(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%rax
mov    -0x858(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x858(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x880(%rbp),%xmm0
ff
mov    -0x854(%rbp),%eax
shl    $0x6,%eax
mov    -0x858(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x2,%eax
cltq
movsd  %xmm0,-0x850(%rbp,%rax,8)
ff ff
mov    -0x858(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x858(%rbp)
jmp    1230 <coreGTRGAMMASECONDARYINVAR+0x110>
jmp    132f <coreGTRGAMMASECONDARYINVAR+0x20f>
mov    -0x854(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x854(%rbp)
jmp    11e9 <coreGTRGAMMASECONDARYINVAR+0xc9>
movl   $0x0,-0x854(%rbp)
mov    -0x854(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1663 <coreGTRGAMMASECONDARYINVAR+0x543>
mov    -0x18(%rbp),%rax
mov    -0x854(%rbp),%ecx
shl    $0x6,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x918(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x920(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x928(%rbp)
ff
movl   $0x0,-0x85c(%rbp)
cmpl   $0x4,-0x85c(%rbp)
jge    14ee <coreGTRGAMMASECONDARYINVAR+0x3ce>
mov    -0x48(%rbp),%rax
mov    -0x85c(%rbp),%ecx
shl    $0x4,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x918(%rbp),%xmm0
ff
movsd  %xmm0,-0x918(%rbp)
ff
movl   $0x1,-0x858(%rbp)
cmpl   $0x10,-0x858(%rbp)
jge    14d5 <coreGTRGAMMASECONDARYINVAR+0x3b5>
mov    -0x85c(%rbp),%eax
shl    $0x6,%eax
mov    -0x858(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
cltq
movsd  -0x850(%rbp,%rax,8),%xmm0
ff ff
mov    -0x48(%rbp),%rax
mov    -0x85c(%rbp),%ecx
shl    $0x4,%ecx
add    -0x858(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x910(%rbp)
ff
addsd  -0x918(%rbp),%xmm0
ff
movsd  %xmm0,-0x918(%rbp)
ff
movsd  -0x910(%rbp),%xmm0
ff
mov    -0x85c(%rbp),%eax
shl    $0x6,%eax
mov    -0x858(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x1,%eax
cltq
movsd  -0x850(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x920(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x920(%rbp)
ff
movsd  -0x910(%rbp),%xmm0
ff
mov    -0x85c(%rbp),%eax
shl    $0x6,%eax
mov    -0x858(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x2,%eax
cltq
movsd  -0x850(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x928(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x928(%rbp)
ff
mov    -0x858(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x858(%rbp)
jmp    13de <coreGTRGAMMASECONDARYINVAR+0x2be>
jmp    14da <coreGTRGAMMASECONDARYINVAR+0x3ba>
mov    -0x85c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x85c(%rbp)
jmp    13a2 <coreGTRGAMMASECONDARYINVAR+0x282>
movsd  -0x918(%rbp),%xmm0
ff
call   1030 <FABS@plt>
mulsd  -0x908(%rbp),%xmm0
ff
movsd  %xmm0,-0x918(%rbp)
ff
mov    0x20(%rbp),%rax
movslq -0x854(%rbp),%rcx
cmpl   $0x10,(%rax,%rcx,4)
jge    1548 <coreGTRGAMMASECONDARYINVAR+0x428>
mov    0x20(%rbp),%rax
movslq -0x854(%rbp),%rcx
movslq (%rax,%rcx,4),%rax
movsd  -0x900(%rbp,%rax,8),%xmm0
ff ff
addsd  -0x918(%rbp),%xmm0
ff
movsd  %xmm0,-0x918(%rbp)
ff
movsd  0xab8(%rip),%xmm0
divsd  -0x918(%rbp),%xmm0
ff
movsd  %xmm0,-0x918(%rbp)
ff
movsd  -0x918(%rbp),%xmm0
ff
mulsd  -0x920(%rbp),%xmm0
ff
movsd  %xmm0,-0x920(%rbp)
ff
movsd  -0x918(%rbp),%xmm0
ff
mulsd  -0x928(%rbp),%xmm0
ff
movsd  %xmm0,-0x928(%rbp)
ff
movsd  -0x908(%rbp),%xmm0
ff
mulsd  -0x920(%rbp),%xmm0
ff
movsd  %xmm0,-0x920(%rbp)
ff
movsd  -0x908(%rbp),%xmm0
ff
mulsd  -0x928(%rbp),%xmm0
ff
movsd  %xmm0,-0x928(%rbp)
ff
mov    -0x28(%rbp),%rax
movslq -0x854(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x920(%rbp),%xmm2
ff
movsd  -0x868(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x868(%rbp)
ff
mov    -0x28(%rbp),%rax
movslq -0x854(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x928(%rbp),%xmm1
ff
movsd  -0x920(%rbp),%xmm2
ff
movsd  -0x920(%rbp),%xmm3
ff
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  -0x870(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x870(%rbp)
ff
mov    -0x854(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x854(%rbp)
jmp    134d <coreGTRGAMMASECONDARYINVAR+0x22d>
movsd  -0x868(%rbp),%xmm0
ff
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x870(%rbp),%xmm0
ff
mov    0x10(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x930,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<coreGTRGAMMASECONDARYINVAR>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x8d8,%rsp
movapd %xmm0,%xmm2
mov    %r9,%r15
mov    %r8,0x48(%rsp)
mov    %ecx,%r12d
mov    %rdx,%r14
mov    %rsi,%rbp
mov    %rdi,%r13
mov    0x918(%rsp),%rax
movsd  0xea7(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x28(%rsp)
movupd (%rax),%xmm0
unpcklpd %xmm1,%xmm1
mulpd  %xmm1,%xmm0
movapd %xmm0,0x50(%rsp)
movupd 0x10(%rax),%xmm0
mulpd  %xmm1,%xmm0
movapd %xmm0,0x60(%rsp)
movupd 0x20(%rax),%xmm0
mulpd  %xmm1,%xmm0
movapd %xmm0,0x70(%rsp)
movupd 0x30(%rax),%xmm0
mulpd  %xmm1,%xmm0
movapd %xmm0,0x80(%rsp)
movupd 0x40(%rax),%xmm0
mulpd  %xmm1,%xmm0
movapd %xmm0,0x90(%rsp)
movupd 0x50(%rax),%xmm0
mulpd  %xmm1,%xmm0
movapd %xmm0,0xa0(%rsp)
movupd 0x60(%rax),%xmm0
mulpd  %xmm1,%xmm0
movapd %xmm0,0xb0(%rsp)
movupd 0x70(%rax),%xmm0
mulpd  %xmm1,%xmm0
movapd %xmm0,0xc0(%rsp)
movsd  (%rdi),%xmm1
movapd %xmm1,%xmm0
mulsd  %xmm1,%xmm0
movapd %xmm1,(%rsp)
unpcklpd %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
xor    %ebx,%ebx
movsd  %xmm2,0x30(%rsp)
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm2,%xmm0
call   1040 <exp@plt>
movsd  0x30(%rsp),%xmm2
movsd  %xmm0,0xf0(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
mulpd  0x10(%rsp),%xmm0
movupd %xmm0,0xf8(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x78,%rbx
jne    1210 <coreGTRGAMMASECONDARYINVAR+0xf0>
movsd  0x8(%r13),%xmm1
movapd %xmm1,%xmm0
mulsd  %xmm1,%xmm0
movapd %xmm1,(%rsp)
unpcklpd %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
xor    %ebx,%ebx
nopl   (%rax)
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm2,%xmm0
call   1040 <exp@plt>
movsd  0x30(%rsp),%xmm2
movsd  %xmm0,0x2f0(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
mulpd  0x10(%rsp),%xmm0
movupd %xmm0,0x2f8(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x78,%rbx
jne    1280 <coreGTRGAMMASECONDARYINVAR+0x160>
movsd  0x10(%r13),%xmm1
movapd %xmm1,%xmm0
mulsd  %xmm1,%xmm0
movapd %xmm1,(%rsp)
unpcklpd %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
xor    %ebx,%ebx
nopl   (%rax)
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm2,%xmm0
call   1040 <exp@plt>
movsd  0x30(%rsp),%xmm2
movsd  %xmm0,0x4f0(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
mulpd  0x10(%rsp),%xmm0
movupd %xmm0,0x4f8(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x78,%rbx
jne    12f0 <coreGTRGAMMASECONDARYINVAR+0x1d0>
movsd  0x18(%r13),%xmm1
movapd %xmm1,%xmm0
mulsd  %xmm1,%xmm0
movapd %xmm1,(%rsp)
unpcklpd %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
xor    %ebx,%ebx
nopl   (%rax)
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm2,%xmm0
call   1040 <exp@plt>
movsd  0x30(%rsp),%xmm2
movsd  %xmm0,0x6f0(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
mulpd  0x10(%rsp),%xmm0
movupd %xmm0,0x6f8(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x78,%rbx
jne    1360 <coreGTRGAMMASECONDARYINVAR+0x240>
test   %r12d,%r12d
jle    16a0 <coreGTRGAMMASECONDARYINVAR+0x580>
mov    %r15,0x40(%rsp)
movsd  0x28(%rsp),%xmm3
mulsd  0xc3e(%rip),%xmm3
mov    %r12d,%eax
mov    %rax,0x10(%rsp)
lea    0x10(%r14),%r13
lea    0x90(%r14),%rbx
lea    0x110(%r14),%rbp
lea    0x190(%r14),%r15
xorpd  %xmm0,%xmm0
xor    %r12d,%r12d
movsd  %xmm3,0x28(%rsp)
jmp    1488 <coreGTRGAMMASECONDARYINVAR+0x368>
nopl   (%rax)
movapd (%rsp),%xmm4
movsd  0xbf3(%rip),%xmm1
divsd  %xmm0,%xmm1
movapd %xmm4,%xmm0
mulsd  %xmm1,%xmm0
unpckhpd %xmm4,%xmm4
mulsd  %xmm1,%xmm4
mulsd  %xmm3,%xmm0
mulsd  %xmm3,%xmm4
mov    0x48(%rsp),%rax
xorps  %xmm1,%xmm1
cvtsi2sdl (%rax,%r12,4),%xmm1
movapd %xmm0,%xmm2
mulsd  %xmm0,%xmm2
subsd  %xmm2,%xmm4
unpcklpd %xmm1,%xmm1
unpcklpd %xmm0,%xmm4
mulpd  %xmm1,%xmm4
movapd %xmm4,%xmm0
movapd 0x30(%rsp),%xmm1
addpd  %xmm1,%xmm0
add    $0x1,%r12
add    $0x200,%r13
add    $0x200,%rbx
add    $0x200,%rbp
add    $0x200,%r15
cmp    0x10(%rsp),%r12
je     1695 <coreGTRGAMMASECONDARYINVAR+0x575>
movapd %xmm0,0x30(%rsp)
mov    %r12,%rax
shl    $0x6,%rax
mov    %r12,%rcx
shl    $0x9,%rcx
movsd  (%r14,%rcx,1),%xmm0
addsd  0xb66(%rip),%xmm0
xor    %ecx,%ecx
xorpd  %xmm2,%xmm2
movsd  0xf0(%rsp,%rcx,4),%xmm1
mulsd  -0x8(%r13,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0xf8(%rsp,%rcx,4),%xmm3
unpcklpd %xmm1,%xmm1
mulpd  %xmm3,%xmm1
addpd  %xmm2,%xmm1
cmp    $0x70,%rcx
je     1510 <coreGTRGAMMASECONDARYINVAR+0x3f0>
movsd  0x110(%rsp,%rcx,4),%xmm2
mulsd  0x0(%r13,%rcx,1),%xmm2
addsd  %xmm2,%xmm0
movupd 0x118(%rsp,%rcx,4),%xmm3
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
add    $0x10,%rcx
jmp    14b0 <coreGTRGAMMASECONDARYINVAR+0x390>
xchg   %ax,%ax
addsd  0x80(%r14,%rax,8),%xmm0
xor    %ecx,%ecx
nopl   0x0(%rax)
movsd  0x2f0(%rsp,%rcx,4),%xmm2
mulsd  -0x8(%rbx,%rcx,1),%xmm2
addsd  %xmm2,%xmm0
movupd 0x2f8(%rsp,%rcx,4),%xmm3
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
cmp    $0x70,%rcx
je     1580 <coreGTRGAMMASECONDARYINVAR+0x460>
movsd  0x310(%rsp,%rcx,4),%xmm1
mulsd  (%rbx,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0x318(%rsp,%rcx,4),%xmm3
unpcklpd %xmm1,%xmm1
mulpd  %xmm3,%xmm1
addpd  %xmm2,%xmm1
add    $0x10,%rcx
jmp    1520 <coreGTRGAMMASECONDARYINVAR+0x400>
nopl   0x0(%rax,%rax,1)
addsd  0x100(%r14,%rax,8),%xmm0
xor    %ecx,%ecx
nopl   0x0(%rax)
movsd  0x4f0(%rsp,%rcx,4),%xmm1
mulsd  -0x8(%rbp,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0x4f8(%rsp,%rcx,4),%xmm3
unpcklpd %xmm1,%xmm1
mulpd  %xmm3,%xmm1
addpd  %xmm2,%xmm1
cmp    $0x70,%rcx
je     15f0 <coreGTRGAMMASECONDARYINVAR+0x4d0>
movsd  0x510(%rsp,%rcx,4),%xmm2
mulsd  0x0(%rbp,%rcx,1),%xmm2
addsd  %xmm2,%xmm0
movupd 0x518(%rsp,%rcx,4),%xmm3
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
add    $0x10,%rcx
jmp    1590 <coreGTRGAMMASECONDARYINVAR+0x470>
nopl   0x0(%rax)
addsd  0x180(%r14,%rax,8),%xmm0
xor    %eax,%eax
nopl   0x0(%rax)
movsd  0x6f0(%rsp,%rax,4),%xmm3
mulsd  -0x8(%r15,%rax,1),%xmm3
addsd  %xmm3,%xmm0
movupd 0x6f8(%rsp,%rax,4),%xmm2
unpcklpd %xmm3,%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm1,%xmm3
cmp    $0x70,%rax
je     1660 <coreGTRGAMMASECONDARYINVAR+0x540>
movsd  0x710(%rsp,%rax,4),%xmm1
mulsd  (%r15,%rax,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0x718(%rsp,%rax,4),%xmm2
unpcklpd %xmm1,%xmm1
mulpd  %xmm2,%xmm1
addpd  %xmm3,%xmm1
add    $0x10,%rax
jmp    1600 <coreGTRGAMMASECONDARYINVAR+0x4e0>
nopl   (%rax)
movapd %xmm3,(%rsp)
call   1030 <FABS@plt>
movsd  0x28(%rsp),%xmm3
mulsd  %xmm3,%xmm0
mov    0x920(%rsp),%rax
movslq (%rax,%r12,4),%rax
cmp    $0xf,%rax
jg     1400 <coreGTRGAMMASECONDARYINVAR+0x2e0>
addsd  0x50(%rsp,%rax,8),%xmm0
jmp    1400 <coreGTRGAMMASECONDARYINVAR+0x2e0>
addpd  %xmm1,%xmm4
mov    0x40(%rsp),%r15
jmp    16a4 <coreGTRGAMMASECONDARYINVAR+0x584>
xorpd  %xmm4,%xmm4
movhpd %xmm4,(%r15)
mov    0x910(%rsp),%rax
movlpd %xmm4,(%rax)
add    $0x8d8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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