#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 coreGTRGAMMAPROTINVAR(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[512];
  int i, l, j;
  double dlnLdlz = 0;
  double d2lnLdlz2 = 0;
  double ki, kisqr;
  double freqs[20];
  double scaler = 0.25 * (1.0 - propInvar);
  double tmp;
  double inv_Li, dlnLidlz, d2lnLidlz2;
  for (i = 0; i < 20; i++)
    freqs[i] = frequencies[i] * propInvar;
  for (i = 0; i < 4; i++) {
    ki = gammaRates[i];
    kisqr = ki * ki;
    for (l = 1; l < 20; l++) {
      diagptable[i * 128 + l * 4] = exp(EIGN[l - 1] * ki * lz);
      diagptable[i * 128 + l * 4 + 1] = EIGN[l - 1] * ki;
      diagptable[i * 128 + l * 4 + 2] = EIGN[l - 1] * EIGN[l - 1] * kisqr;
    }
  }
  for (i = 0; i < upper; i++) {
    sum = &sumtable[i * 80];
    inv_Li = 0.0;
    dlnLidlz = 0.0;
    d2lnLidlz2 = 0.0;
    for (j = 0; j < 4; j++) {
      inv_Li += sum[j * 20];
      for (l = 1; l < 20; l++) {
        inv_Li += (tmp = diagptable[j * 128 + l * 4] * sum[j * 20 + l]);
        dlnLidlz += tmp * diagptable[j * 128 + l * 4 + 1];
        d2lnLidlz2 += tmp * diagptable[j * 128 + l * 4 + 2];
      }
    }
    inv_Li = FABS(inv_Li) * scaler;
    if (iptr[i] < 20)
      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):
<coreGTRGAMMAPROTINVAR>:
push   %rbp
mov    %rsp,%rbp
sub    $0x1150,%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,-0x1068(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1070(%rbp)
ff
movsd  0xe92(%rip),%xmm1
subsd  -0x40(%rbp),%xmm1
movsd  0xe7d(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x1128(%rbp)
ff
movl   $0x0,-0x1054(%rbp)
cmpl   $0x14,-0x1054(%rbp)
jge    11df <coreGTRGAMMAPROTINVAR+0xbf>
mov    0x18(%rbp),%rax
movslq -0x1054(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x40(%rbp),%xmm0
movslq -0x1054(%rbp),%rax
movsd  %xmm0,-0x1120(%rbp,%rax,8)
ff ff
mov    -0x1054(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1054(%rbp)
jmp    1199 <coreGTRGAMMAPROTINVAR+0x79>
movl   $0x0,-0x1054(%rbp)
cmpl   $0x4,-0x1054(%rbp)
jge    1343 <coreGTRGAMMAPROTINVAR+0x223>
mov    -0x8(%rbp),%rax
movslq -0x1054(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x1078(%rbp)
ff
movsd  -0x1078(%rbp),%xmm0
ff
mulsd  -0x1078(%rbp),%xmm0
ff
movsd  %xmm0,-0x1080(%rbp)
ff
movl   $0x1,-0x1058(%rbp)
cmpl   $0x14,-0x1058(%rbp)
jge    132a <coreGTRGAMMAPROTINVAR+0x20a>
mov    -0x10(%rbp),%rax
mov    -0x1058(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x1078(%rbp),%xmm0
ff
mulsd  -0x38(%rbp),%xmm0
call   1040 <exp@plt>
mov    -0x1054(%rbp),%eax
shl    $0x7,%eax
mov    -0x1058(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
cltq
movsd  %xmm0,-0x1050(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%rax
mov    -0x1058(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x1078(%rbp),%xmm0
ff
mov    -0x1054(%rbp),%eax
shl    $0x7,%eax
mov    -0x1058(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x1,%eax
cltq
movsd  %xmm0,-0x1050(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%rax
mov    -0x1058(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1058(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x1080(%rbp),%xmm0
ff
mov    -0x1054(%rbp),%eax
shl    $0x7,%eax
mov    -0x1058(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x2,%eax
cltq
movsd  %xmm0,-0x1050(%rbp,%rax,8)
ff ff
mov    -0x1058(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1058(%rbp)
jmp    1230 <coreGTRGAMMAPROTINVAR+0x110>
jmp    132f <coreGTRGAMMAPROTINVAR+0x20f>
mov    -0x1054(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1054(%rbp)
jmp    11e9 <coreGTRGAMMAPROTINVAR+0xc9>
movl   $0x0,-0x1054(%rbp)
mov    -0x1054(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    165d <coreGTRGAMMAPROTINVAR+0x53d>
mov    -0x18(%rbp),%rax
imul   $0x50,-0x1054(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1138(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1140(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1148(%rbp)
ff
movl   $0x0,-0x105c(%rbp)
cmpl   $0x4,-0x105c(%rbp)
jge    14e8 <coreGTRGAMMAPROTINVAR+0x3c8>
mov    -0x48(%rbp),%rax
imul   $0x14,-0x105c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x1138(%rbp),%xmm0
ff
movsd  %xmm0,-0x1138(%rbp)
ff
movl   $0x1,-0x1058(%rbp)
cmpl   $0x14,-0x1058(%rbp)
jge    14cf <coreGTRGAMMAPROTINVAR+0x3af>
mov    -0x105c(%rbp),%eax
shl    $0x7,%eax
mov    -0x1058(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
cltq
movsd  -0x1050(%rbp,%rax,8),%xmm0
ff ff
mov    -0x48(%rbp),%rax
imul   $0x14,-0x105c(%rbp),%ecx
add    -0x1058(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x1130(%rbp)
ff
addsd  -0x1138(%rbp),%xmm0
ff
movsd  %xmm0,-0x1138(%rbp)
ff
movsd  -0x1130(%rbp),%xmm0
ff
mov    -0x105c(%rbp),%eax
shl    $0x7,%eax
mov    -0x1058(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x1,%eax
cltq
movsd  -0x1050(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x1140(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1140(%rbp)
ff
movsd  -0x1130(%rbp),%xmm0
ff
mov    -0x105c(%rbp),%eax
shl    $0x7,%eax
mov    -0x1058(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x2,%eax
cltq
movsd  -0x1050(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x1148(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1148(%rbp)
ff
mov    -0x1058(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1058(%rbp)
jmp    13da <coreGTRGAMMAPROTINVAR+0x2ba>
jmp    14d4 <coreGTRGAMMAPROTINVAR+0x3b4>
mov    -0x105c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x105c(%rbp)
jmp    13a0 <coreGTRGAMMAPROTINVAR+0x280>
movsd  -0x1138(%rbp),%xmm0
ff
call   1030 <FABS@plt>
mulsd  -0x1128(%rbp),%xmm0
ff
movsd  %xmm0,-0x1138(%rbp)
ff
mov    0x20(%rbp),%rax
movslq -0x1054(%rbp),%rcx
cmpl   $0x14,(%rax,%rcx,4)
jge    1542 <coreGTRGAMMAPROTINVAR+0x422>
mov    0x20(%rbp),%rax
movslq -0x1054(%rbp),%rcx
movslq (%rax,%rcx,4),%rax
movsd  -0x1120(%rbp,%rax,8),%xmm0
ff ff
addsd  -0x1138(%rbp),%xmm0
ff
movsd  %xmm0,-0x1138(%rbp)
ff
movsd  0xabe(%rip),%xmm0
divsd  -0x1138(%rbp),%xmm0
ff
movsd  %xmm0,-0x1138(%rbp)
ff
movsd  -0x1138(%rbp),%xmm0
ff
mulsd  -0x1140(%rbp),%xmm0
ff
movsd  %xmm0,-0x1140(%rbp)
ff
movsd  -0x1138(%rbp),%xmm0
ff
mulsd  -0x1148(%rbp),%xmm0
ff
movsd  %xmm0,-0x1148(%rbp)
ff
movsd  -0x1128(%rbp),%xmm0
ff
mulsd  -0x1140(%rbp),%xmm0
ff
movsd  %xmm0,-0x1140(%rbp)
ff
movsd  -0x1128(%rbp),%xmm0
ff
mulsd  -0x1148(%rbp),%xmm0
ff
movsd  %xmm0,-0x1148(%rbp)
ff
mov    -0x28(%rbp),%rax
movslq -0x1054(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x1140(%rbp),%xmm2
ff
movsd  -0x1068(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1068(%rbp)
ff
mov    -0x28(%rbp),%rax
movslq -0x1054(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x1148(%rbp),%xmm1
ff
movsd  -0x1140(%rbp),%xmm2
ff
movsd  -0x1140(%rbp),%xmm3
ff
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  -0x1070(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1070(%rbp)
ff
mov    -0x1054(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1054(%rbp)
jmp    134d <coreGTRGAMMAPROTINVAR+0x22d>
movsd  -0x1068(%rbp),%xmm0
ff
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x1070(%rbp),%xmm0
ff
mov    0x10(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x1150,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<coreGTRGAMMAPROTINVAR>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x10f8,%rsp
movapd %xmm0,%xmm2
mov    %r9,%r15
mov    %r8,0x48(%rsp)
mov    %ecx,%r12d
mov    %rdx,%r14
mov    %rsi,%rbp
mov    %rdi,%r13
mov    0x1138(%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)
movupd 0x80(%rax),%xmm0
mulpd  %xmm1,%xmm0
movapd %xmm0,0xd0(%rsp)
movupd 0x90(%rax),%xmm0
mulpd  %xmm1,%xmm0
movapd %xmm0,0xe0(%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)
nopw   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm2,%xmm0
call   1040 <exp@plt>
movsd  0x30(%rsp),%xmm2
movsd  %xmm0,0x110(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
mulpd  0x10(%rsp),%xmm0
movupd %xmm0,0x118(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    1240 <coreGTRGAMMAPROTINVAR+0x120>
movsd  0x8(%r13),%xmm1
movapd %xmm1,%xmm0
mulsd  %xmm1,%xmm0
movapd %xmm1,(%rsp)
unpcklpd %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
xor    %ebx,%ebx
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm2,%xmm0
call   1040 <exp@plt>
movsd  0x30(%rsp),%xmm2
movsd  %xmm0,0x510(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
mulpd  0x10(%rsp),%xmm0
movupd %xmm0,0x518(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    12b0 <coreGTRGAMMAPROTINVAR+0x190>
movsd  0x10(%r13),%xmm1
movapd %xmm1,%xmm0
mulsd  %xmm1,%xmm0
movapd %xmm1,(%rsp)
unpcklpd %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
xor    %ebx,%ebx
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm2,%xmm0
call   1040 <exp@plt>
movsd  0x30(%rsp),%xmm2
movsd  %xmm0,0x910(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
mulpd  0x10(%rsp),%xmm0
movupd %xmm0,0x918(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    1320 <coreGTRGAMMAPROTINVAR+0x200>
movsd  0x18(%r13),%xmm1
movapd %xmm1,%xmm0
mulsd  %xmm1,%xmm0
movapd %xmm1,(%rsp)
unpcklpd %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
xor    %ebx,%ebx
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm2,%xmm0
call   1040 <exp@plt>
movsd  0x30(%rsp),%xmm2
movsd  %xmm0,0xd10(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
mulpd  0x10(%rsp),%xmm0
movupd %xmm0,0xd18(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    1390 <coreGTRGAMMAPROTINVAR+0x270>
test   %r12d,%r12d
jle    16e0 <coreGTRGAMMAPROTINVAR+0x5c0>
mov    %r15,0x40(%rsp)
movsd  0x28(%rsp),%xmm3
mulsd  0xc0b(%rip),%xmm3
mov    %r12d,%eax
mov    %rax,0x10(%rsp)
lea    0x10(%r14),%r13
lea    0xb0(%r14),%rbx
lea    0x150(%r14),%rbp
lea    0x1f0(%r14),%r15
xorpd  %xmm0,%xmm0
xor    %r12d,%r12d
movsd  %xmm3,0x28(%rsp)
jmp    14b8 <coreGTRGAMMAPROTINVAR+0x398>
movapd (%rsp),%xmm4
movsd  0xbc3(%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    $0x280,%r13
add    $0x280,%rbx
add    $0x280,%rbp
add    $0x280,%r15
cmp    0x10(%rsp),%r12
je     16d5 <coreGTRGAMMAPROTINVAR+0x5b5>
movapd %xmm0,0x30(%rsp)
mov    %r12,%rax
shl    $0x4,%rax
lea    (%rax,%rax,4),%rax
movsd  (%r14,%rax,8),%xmm0
addsd  0xb39(%rip),%xmm0
xor    %ecx,%ecx
xorpd  %xmm2,%xmm2
nopl   (%rax)
movsd  0x110(%rsp,%rcx,4),%xmm1
mulsd  -0x8(%r13,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0x118(%rsp,%rcx,4),%xmm3
unpcklpd %xmm1,%xmm1
mulpd  %xmm3,%xmm1
addpd  %xmm2,%xmm1
cmp    $0x90,%rcx
je     1550 <coreGTRGAMMAPROTINVAR+0x430>
movsd  0x130(%rsp,%rcx,4),%xmm2
mulsd  0x0(%r13,%rcx,1),%xmm2
addsd  %xmm2,%xmm0
movupd 0x138(%rsp,%rcx,4),%xmm3
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
add    $0x10,%rcx
jmp    14e0 <coreGTRGAMMAPROTINVAR+0x3c0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
addsd  0xa0(%r14,%rax,8),%xmm0
xor    %ecx,%ecx
nopl   0x0(%rax)
movsd  0x510(%rsp,%rcx,4),%xmm2
mulsd  -0x8(%rbx,%rcx,1),%xmm2
addsd  %xmm2,%xmm0
movupd 0x518(%rsp,%rcx,4),%xmm3
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
cmp    $0x90,%rcx
je     15c0 <coreGTRGAMMAPROTINVAR+0x4a0>
movsd  0x530(%rsp,%rcx,4),%xmm1
mulsd  (%rbx,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0x538(%rsp,%rcx,4),%xmm3
unpcklpd %xmm1,%xmm1
mulpd  %xmm3,%xmm1
addpd  %xmm2,%xmm1
add    $0x10,%rcx
jmp    1560 <coreGTRGAMMAPROTINVAR+0x440>
xchg   %ax,%ax
addsd  0x140(%r14,%rax,8),%xmm0
xor    %ecx,%ecx
nopl   0x0(%rax)
movsd  0x910(%rsp,%rcx,4),%xmm1
mulsd  -0x8(%rbp,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0x918(%rsp,%rcx,4),%xmm3
unpcklpd %xmm1,%xmm1
mulpd  %xmm3,%xmm1
addpd  %xmm2,%xmm1
cmp    $0x90,%rcx
je     1630 <coreGTRGAMMAPROTINVAR+0x510>
movsd  0x930(%rsp,%rcx,4),%xmm2
mulsd  0x0(%rbp,%rcx,1),%xmm2
addsd  %xmm2,%xmm0
movupd 0x938(%rsp,%rcx,4),%xmm3
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
add    $0x10,%rcx
jmp    15d0 <coreGTRGAMMAPROTINVAR+0x4b0>
nop
addsd  0x1e0(%r14,%rax,8),%xmm0
xor    %eax,%eax
nopl   0x0(%rax)
movsd  0xd10(%rsp,%rax,4),%xmm3
mulsd  -0x8(%r15,%rax,1),%xmm3
addsd  %xmm3,%xmm0
movupd 0xd18(%rsp,%rax,4),%xmm2
unpcklpd %xmm3,%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm1,%xmm3
cmp    $0x90,%rax
je     16a0 <coreGTRGAMMAPROTINVAR+0x580>
movsd  0xd30(%rsp,%rax,4),%xmm1
mulsd  (%r15,%rax,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0xd38(%rsp,%rax,4),%xmm2
unpcklpd %xmm1,%xmm1
mulpd  %xmm2,%xmm1
addpd  %xmm3,%xmm1
add    $0x10,%rax
jmp    1640 <coreGTRGAMMAPROTINVAR+0x520>
nop
movapd %xmm3,(%rsp)
call   1030 <FABS@plt>
movsd  0x28(%rsp),%xmm3
mulsd  %xmm3,%xmm0
mov    0x1140(%rsp),%rax
movslq (%rax,%r12,4),%rax
cmp    $0x13,%rax
jg     1430 <coreGTRGAMMAPROTINVAR+0x310>
addsd  0x50(%rsp,%rax,8),%xmm0
jmp    1430 <coreGTRGAMMAPROTINVAR+0x310>
addpd  %xmm1,%xmm4
mov    0x40(%rsp),%r15
jmp    16e4 <coreGTRGAMMAPROTINVAR+0x5c4>
xorpd  %xmm4,%xmm4
movhpd %xmm4,(%r15)
mov    0x1130(%rsp),%rax
movlpd %xmm4,(%rax)
add    $0x10f8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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