#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 */
int /*<<< orphan*/ exp(double);
int /*<<< orphan*/ free(double *);
scalar_t__ malloc(int);
void coreGTRGAMMAPROTINVAR(double *gammaRates, double *EIGN, double *sumtable,
                           int lower, int upper, int *wrptr,
                           double *ext_dlnLdlz, double *ext_d2lnLdlz2,
                           double lz, double *frequencies, double *invariants,
                           int *iptr) {
  double *sum, *diagptable, *diagptable_start;
  int i, l;
  double dlnLdlz = 0;
  double d2lnLdlz2 = 0;
  double ki, kisqr;
  double freqs[20];
  double scaler = 0.25 * (1.0 - invariants[0]);
  double tmp;
  double inv_Li, dlnLidlz, d2lnLidlz2;
  for (i = 0; i < 20; i++)
    freqs[i] = frequencies[i] * invariants[0];
  diagptable = diagptable_start = (double *)malloc(sizeof(double) * 228);
  for (i = 0; i < 4; i++) {
    ki = gammaRates[i];
    kisqr = ki * ki;
    for (l = 0; l < 19; l++) {
      *diagptable++ = exp(EIGN[l] * ki * lz);
      *diagptable++ = EIGN[l] * ki;
      *diagptable++ = EIGN[l] * EIGN[l] * kisqr;
    }
  }
  for (i = lower; i < upper; i++) {
    diagptable = diagptable_start;
    sum = &sumtable[i * 80];
    inv_Li = 0.0;
    dlnLidlz = 0.0;
    d2lnLidlz2 = 0.0;
    inv_Li += *sum++;
    for (l = 1; l < 20; l++) {
      inv_Li += (tmp = *diagptable++ * *sum++);
      dlnLidlz += tmp * *diagptable++;
      d2lnLidlz2 += tmp * *diagptable++;
    }
    inv_Li += *sum++;
    for (l = 1; l < 20; l++) {
      inv_Li += (tmp = *diagptable++ * *sum++);
      dlnLidlz += tmp * *diagptable++;
      d2lnLidlz2 += tmp * *diagptable++;
    }
    inv_Li += *sum++;
    for (l = 1; l < 20; l++) {
      inv_Li += (tmp = *diagptable++ * *sum++);
      dlnLidlz += tmp * *diagptable++;
      d2lnLidlz2 += tmp * *diagptable++;
    }
    inv_Li += *sum++;
    for (l = 1; l < 20; l++) {
      inv_Li += (tmp = *diagptable++ * *sum++);
      dlnLidlz += tmp * *diagptable++;
      d2lnLidlz2 += tmp * *diagptable++;
    }
    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;
  free(diagptable_start);
}


x86 assembly (clang, -O0):
<coreGTRGAMMAPROTINVAR>:
push   %rbp
mov    %rsp,%rbp
sub    $0x140,%rsp
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
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    %r8d,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    0x28(%rbp),%rax
movsd  0xe81(%rip),%xmm1
subsd  (%rax),%xmm1
movsd  0xe6d(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
movl   $0x0,-0x4c(%rbp)
cmpl   $0x14,-0x4c(%rbp)
jge    11e0 <coreGTRGAMMAPROTINVAR+0xb0>
mov    0x20(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x28(%rbp),%rax
mulsd  (%rax),%xmm0
movslq -0x4c(%rbp),%rax
movsd  %xmm0,-0x110(%rbp,%rax,8)
ff ff
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11a6 <coreGTRGAMMAPROTINVAR+0x76>
mov    $0x720,%edi
call   1040 <malloc@plt>
mov    %rax,-0x48(%rbp)
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x4c(%rbp)
cmpl   $0x4,-0x4c(%rbp)
jge    12e0 <coreGTRGAMMAPROTINVAR+0x1b0>
mov    -0x8(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
mulsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x50(%rbp)
cmpl   $0x13,-0x50(%rbp)
jge    12cd <coreGTRGAMMAPROTINVAR+0x19d>
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x68(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
call   1050 <exp@plt>
cvtsi2sd %eax,%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x68(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x70(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    122b <coreGTRGAMMAPROTINVAR+0xfb>
jmp    12d2 <coreGTRGAMMAPROTINVAR+0x1a2>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11f9 <coreGTRGAMMAPROTINVAR+0xc9>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    181b <coreGTRGAMMAPROTINVAR+0x6eb>
mov    -0x48(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x18(%rbp),%rax
imul   $0x50,-0x4c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x128(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x138(%rbp)
ff
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
addsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x128(%rbp)
ff
movl   $0x1,-0x50(%rbp)
cmpl   $0x14,-0x50(%rbp)
jge    1417 <coreGTRGAMMAPROTINVAR+0x2e7>
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
addsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x120(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x130(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0x120(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x138(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x138(%rbp)
ff
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    135b <coreGTRGAMMAPROTINVAR+0x22b>
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
addsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x128(%rbp)
ff
movl   $0x1,-0x50(%rbp)
cmpl   $0x14,-0x50(%rbp)
jge    14fd <coreGTRGAMMAPROTINVAR+0x3cd>
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
addsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x120(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x130(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0x120(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x138(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x138(%rbp)
ff
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1441 <coreGTRGAMMAPROTINVAR+0x311>
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
addsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x128(%rbp)
ff
movl   $0x1,-0x50(%rbp)
cmpl   $0x14,-0x50(%rbp)
jge    15e3 <coreGTRGAMMAPROTINVAR+0x4b3>
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
addsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x120(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x130(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0x120(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x138(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x138(%rbp)
ff
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1527 <coreGTRGAMMAPROTINVAR+0x3f7>
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
addsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x128(%rbp)
ff
movl   $0x1,-0x50(%rbp)
cmpl   $0x14,-0x50(%rbp)
jge    16c9 <coreGTRGAMMAPROTINVAR+0x599>
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
addsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x120(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x130(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0x120(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x138(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x138(%rbp)
ff
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    160d <coreGTRGAMMAPROTINVAR+0x4dd>
movsd  -0x118(%rbp),%xmm0
ff
mulsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x128(%rbp)
ff
mov    0x30(%rbp),%rax
movslq -0x4c(%rbp),%rcx
cmpl   $0x14,(%rax,%rcx,4)
jge    1718 <coreGTRGAMMAPROTINVAR+0x5e8>
mov    0x30(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movslq (%rax,%rcx,4),%rax
movsd  -0x110(%rbp,%rax,8),%xmm0
ff ff
addsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x128(%rbp)
ff
movsd  0x8e8(%rip),%xmm0
divsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x128(%rbp),%xmm0
ff
mulsd  -0x130(%rbp),%xmm0
ff
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0x128(%rbp),%xmm0
ff
mulsd  -0x138(%rbp),%xmm0
ff
movsd  %xmm0,-0x138(%rbp)
ff
movsd  -0x118(%rbp),%xmm0
ff
mulsd  -0x130(%rbp),%xmm0
ff
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0x118(%rbp),%xmm0
ff
mulsd  -0x138(%rbp),%xmm0
ff
movsd  %xmm0,-0x138(%rbp)
ff
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x130(%rbp),%xmm2
ff
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x138(%rbp),%xmm1
ff
movsd  -0x130(%rbp),%xmm2
ff
movsd  -0x130(%rbp),%xmm3
ff
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    12e6 <coreGTRGAMMAPROTINVAR+0x1b6>
movsd  -0x58(%rbp),%xmm0
mov    0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x60(%rbp),%xmm0
mov    0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
add    $0x140,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<coreGTRGAMMAPROTINVAR>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xf8,%rsp
movsd  %xmm0,0x10(%rsp)
mov    %r9,0x30(%rsp)
mov    %r8d,0x1c(%rsp)
mov    %ecx,0x18(%rsp)
mov    %rdx,%r15
mov    %rsi,%r13
mov    %rdi,%rbp
mov    0x150(%rsp),%rax
mov    %rax,0x28(%rsp)
mov    0x138(%rsp),%rax
mov    %rax,0x48(%rsp)
mov    0x130(%rsp),%rax
mov    %rax,0x40(%rsp)
mov    0x140(%rsp),%rax
mov    0x148(%rsp),%rcx
movsd  (%rcx),%xmm0
movsd  0xe5f(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x20(%rsp)
movupd (%rax),%xmm1
unpcklpd %xmm0,%xmm0
mulpd  %xmm0,%xmm1
movapd %xmm1,0x50(%rsp)
movupd 0x10(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0x60(%rsp)
movupd 0x20(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0x70(%rsp)
movupd 0x30(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0x80(%rsp)
movupd 0x40(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0x90(%rsp)
movupd 0x50(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0xa0(%rsp)
movupd 0x60(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0xb0(%rsp)
movupd 0x70(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0xc0(%rsp)
movupd 0x80(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0xd0(%rsp)
movupd 0x90(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0xe0(%rsp)
mov    $0x720,%edi
call   1040 <malloc@plt>
mov    %rax,%rbx
mov    %rbp,0x38(%rsp)
movsd  0x0(%rbp),%xmm0
movsd  %xmm0,(%rsp)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x8(%rsp)
xor    %r14d,%r14d
mov    %rax,%r12
nopl   0x0(%rax,%rax,1)
movsd  0x0(%r13,%r14,8),%xmm0
mulsd  (%rsp),%xmm0
mulsd  0x10(%rsp),%xmm0
call   1050 <exp@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%r12)
movsd  0x0(%r13,%r14,8),%xmm0
mulsd  (%rsp),%xmm0
movsd  %xmm0,0x8(%r12)
movsd  0x0(%r13,%r14,8),%xmm0
mulsd  %xmm0,%xmm0
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x10(%r12)
add    $0x18,%r12
add    $0x1,%r14
cmp    $0x13,%r14
jne    1290 <coreGTRGAMMAPROTINVAR+0x160>
mov    0x38(%rsp),%r14
movsd  0x8(%r14),%xmm0
movsd  %xmm0,(%rsp)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x8(%rsp)
xor    %ebp,%ebp
nopl   0x0(%rax)
movsd  0x0(%r13,%rbp,8),%xmm0
mulsd  (%rsp),%xmm0
mulsd  0x10(%rsp),%xmm0
call   1050 <exp@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%r12)
movsd  0x0(%r13,%rbp,8),%xmm0
mulsd  (%rsp),%xmm0
movsd  %xmm0,0x8(%r12)
movsd  0x0(%r13,%rbp,8),%xmm0
mulsd  %xmm0,%xmm0
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x10(%r12)
add    $0x18,%r12
add    $0x1,%rbp
cmp    $0x13,%rbp
jne    1310 <coreGTRGAMMAPROTINVAR+0x1e0>
movsd  0x10(%r14),%xmm0
movsd  %xmm0,(%rsp)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x8(%rsp)
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  0x0(%r13,%rbp,8),%xmm0
mulsd  (%rsp),%xmm0
mulsd  0x10(%rsp),%xmm0
call   1050 <exp@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%r12)
movsd  0x0(%r13,%rbp,8),%xmm0
mulsd  (%rsp),%xmm0
movsd  %xmm0,0x8(%r12)
movsd  0x0(%r13,%rbp,8),%xmm0
mulsd  %xmm0,%xmm0
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x10(%r12)
add    $0x18,%r12
add    $0x1,%rbp
cmp    $0x13,%rbp
jne    1390 <coreGTRGAMMAPROTINVAR+0x260>
movsd  0x18(%r14),%xmm0
movsd  %xmm0,(%rsp)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x8(%rsp)
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  0x0(%r13,%rbp,8),%xmm0
mulsd  (%rsp),%xmm0
mulsd  0x10(%rsp),%xmm0
call   1050 <exp@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%r12)
movsd  0x0(%r13,%rbp,8),%xmm0
mulsd  (%rsp),%xmm0
movsd  %xmm0,0x8(%r12)
movsd  0x0(%r13,%rbp,8),%xmm0
mulsd  %xmm0,%xmm0
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x10(%r12)
add    $0x1,%rbp
add    $0x18,%r12
cmp    $0x13,%rbp
jne    1410 <coreGTRGAMMAPROTINVAR+0x2e0>
xorpd  %xmm2,%xmm2
mov    0x1c(%rsp),%ecx
mov    0x18(%rsp),%eax
cmp    %ecx,%eax
jge    1731 <coreGTRGAMMAPROTINVAR+0x601>
movsd  0x20(%rsp),%xmm7
mulsd  0xb79(%rip),%xmm7
cltq
movslq %ecx,%r8
lea    (%rax,%rax,4),%rdx
shl    $0x7,%rdx
add    %r15,%rdx
xorpd  %xmm3,%xmm3
xorpd  %xmm0,%xmm0
mov    0x30(%rsp),%r9
mov    0x28(%rsp),%rcx
movsd  0xb47(%rip),%xmm6
jmp    151d <coreGTRGAMMAPROTINVAR+0x3ed>
nopl   0x0(%rax,%rax,1)
movapd %xmm6,%xmm3
divsd  %xmm2,%xmm3
movapd %xmm4,%xmm2
mulsd  %xmm3,%xmm2
unpckhpd %xmm4,%xmm4
mulsd  %xmm3,%xmm4
mulsd  %xmm7,%xmm2
mulsd  %xmm7,%xmm4
xorps  %xmm3,%xmm3
cvtsi2sdl (%r9,%rax,4),%xmm3
movapd %xmm2,%xmm5
mulsd  %xmm2,%xmm5
subsd  %xmm5,%xmm4
unpcklpd %xmm3,%xmm3
unpcklpd %xmm4,%xmm2
mulpd  %xmm3,%xmm2
movapd %xmm2,%xmm3
addpd  %xmm1,%xmm3
add    $0x1,%rax
add    $0x280,%rdx
cmp    %r8,%rax
je     172d <coreGTRGAMMAPROTINVAR+0x5fd>
movapd %xmm3,%xmm1
lea    (%rax,%rax,4),%rsi
shl    $0x7,%rsi
movsd  (%r15,%rsi,1),%xmm2
addsd  %xmm0,%xmm2
xor    %esi,%esi
mov    %rbx,%rdi
xorpd  %xmm4,%xmm4
nopl   0x0(%rax)
movsd  (%rdi),%xmm3
mulsd  0x8(%rdx,%rsi,1),%xmm3
addsd  %xmm3,%xmm2
movupd 0x8(%rdi),%xmm5
unpcklpd %xmm3,%xmm3
mulpd  %xmm5,%xmm3
addpd  %xmm4,%xmm3
cmp    $0x90,%esi
je     15a0 <coreGTRGAMMAPROTINVAR+0x470>
movsd  0x18(%rdi),%xmm4
mulsd  0x10(%rdx,%rsi,1),%xmm4
addsd  %xmm4,%xmm2
movupd 0x20(%rdi),%xmm5
unpcklpd %xmm4,%xmm4
mulpd  %xmm5,%xmm4
addpd  %xmm3,%xmm4
add    $0x10,%rsi
add    $0x30,%rdi
jmp    1540 <coreGTRGAMMAPROTINVAR+0x410>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
addsd  0x10(%rdx,%rsi,1),%xmm2
lea    (%rdx,%rsi,1),%rdi
add    $0x18,%rdi
xor    %ebp,%ebp
movsd  0x1c8(%rbx,%rbp,8),%xmm4
mulsd  (%rdi),%xmm4
addsd  %xmm4,%xmm2
movupd 0x1d0(%rbx,%rbp,8),%xmm5
unpcklpd %xmm4,%xmm4
mulpd  %xmm5,%xmm4
addpd  %xmm3,%xmm4
cmp    $0x36,%ebp
je     1610 <coreGTRGAMMAPROTINVAR+0x4e0>
movsd  0x1e0(%rbx,%rbp,8),%xmm3
mulsd  0x8(%rdi),%xmm3
addsd  %xmm3,%xmm2
movupd 0x1e8(%rbx,%rbp,8),%xmm5
unpcklpd %xmm3,%xmm3
mulpd  %xmm5,%xmm3
addpd  %xmm4,%xmm3
add    $0x6,%rbp
add    $0x10,%rdi
jmp    15b0 <coreGTRGAMMAPROTINVAR+0x480>
nopl   0x0(%rax)
addsd  0xb0(%rdx,%rsi,1),%xmm2
lea    (%rdx,%rsi,1),%rdi
add    $0xb8,%rdi
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
movsd  0x390(%rbx,%rbp,8),%xmm3
mulsd  (%rdi),%xmm3
addsd  %xmm3,%xmm2
movupd 0x398(%rbx,%rbp,8),%xmm5
unpcklpd %xmm3,%xmm3
mulpd  %xmm5,%xmm3
addpd  %xmm4,%xmm3
cmp    $0x36,%ebp
je     1690 <coreGTRGAMMAPROTINVAR+0x560>
movsd  0x3a8(%rbx,%rbp,8),%xmm4
mulsd  0x8(%rdi),%xmm4
addsd  %xmm4,%xmm2
movupd 0x3b0(%rbx,%rbp,8),%xmm5
unpcklpd %xmm4,%xmm4
mulpd  %xmm5,%xmm4
addpd  %xmm3,%xmm4
add    $0x6,%rbp
add    $0x10,%rdi
jmp    1630 <coreGTRGAMMAPROTINVAR+0x500>
nopl   0x0(%rax)
addsd  0x150(%rdx,%rsi,1),%xmm2
add    %rdx,%rsi
add    $0x158,%rsi
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nop
movsd  0x558(%rbx,%rdi,8),%xmm4
mulsd  (%rsi),%xmm4
addsd  %xmm4,%xmm2
movupd 0x560(%rbx,%rdi,8),%xmm5
unpcklpd %xmm4,%xmm4
mulpd  %xmm5,%xmm4
addpd  %xmm3,%xmm4
cmp    $0x36,%edi
je     1710 <coreGTRGAMMAPROTINVAR+0x5e0>
movsd  0x570(%rbx,%rdi,8),%xmm3
mulsd  0x8(%rsi),%xmm3
addsd  %xmm3,%xmm2
movupd 0x578(%rbx,%rdi,8),%xmm5
unpcklpd %xmm3,%xmm3
mulpd  %xmm5,%xmm3
addpd  %xmm4,%xmm3
add    $0x6,%rdi
add    $0x10,%rsi
jmp    16b0 <coreGTRGAMMAPROTINVAR+0x580>
nopl   0x0(%rax)
mulsd  %xmm7,%xmm2
movslq (%rcx,%rax,4),%rsi
cmp    $0x13,%rsi
jg     14c0 <coreGTRGAMMAPROTINVAR+0x390>
addsd  0x50(%rsp,%rsi,8),%xmm2
jmp    14c0 <coreGTRGAMMAPROTINVAR+0x390>
addpd  %xmm1,%xmm2
mov    0x40(%rsp),%rax
movlpd %xmm2,(%rax)
mov    0x48(%rsp),%rax
movhpd %xmm2,(%rax)
mov    %rbx,%rdi
add    $0xf8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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