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


x86 assembly (clang, -O0):
<coreGTRGAMMAPROTMULTINVAR>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
mov    0x40(%rbp),%rax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%eax
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)
movslq 0x20(%rbp),%rax
imul   $0x720,%rax,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x48(%rbp)
mov    %rax,-0x40(%rbp)
imul   $0x14,0x20(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x88(%rbp)
movslq 0x20(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x80(%rbp)
movl   $0x0,-0x74(%rbp)
mov    -0x74(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    1270 <coreGTRGAMMAPROTMULTINVAR+0x140>
mov    0x30(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  0xe16(%rip),%xmm1
subsd  (%rax,%rcx,8),%xmm1
movsd  0xe11(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x80(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x4c(%rbp)
cmpl   $0x14,-0x4c(%rbp)
jge    125d <coreGTRGAMMAPROTMULTINVAR+0x12d>
mov    0x38(%rbp),%rax
imul   $0x14,-0x74(%rbp),%ecx
add    -0x4c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x30(%rbp),%rax
movslq -0x74(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x88(%rbp),%rax
imul   $0x14,-0x74(%rbp),%ecx
add    -0x4c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    120f <coreGTRGAMMAPROTMULTINVAR+0xdf>
jmp    1262 <coreGTRGAMMAPROTMULTINVAR+0x132>
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    11ce <coreGTRGAMMAPROTMULTINVAR+0x9e>
movl   $0x0,-0x74(%rbp)
mov    -0x74(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    13bd <coreGTRGAMMAPROTMULTINVAR+0x28d>
mov    -0x48(%rbp),%rax
imul   $0xe4,-0x74(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x4c(%rbp)
cmpl   $0x4,-0x4c(%rbp)
jge    13aa <coreGTRGAMMAPROTMULTINVAR+0x27a>
mov    -0x8(%rbp),%rax
mov    -0x74(%rbp),%ecx
shl    $0x2,%ecx
add    -0x4c(%rbp),%ecx
movslq %ecx,%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    1397 <coreGTRGAMMAPROTMULTINVAR+0x267>
mov    -0x10(%rbp),%rax
imul   $0x13,-0x74(%rbp),%ecx
add    -0x50(%rbp),%ecx
movslq %ecx,%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
imul   $0x13,-0x74(%rbp),%ecx
add    -0x50(%rbp),%ecx
movslq %ecx,%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
imul   $0x13,-0x74(%rbp),%ecx
add    -0x50(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
imul   $0x13,-0x74(%rbp),%ecx
add    -0x50(%rbp),%ecx
movslq %ecx,%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    12dd <coreGTRGAMMAPROTMULTINVAR+0x1ad>
jmp    139c <coreGTRGAMMAPROTMULTINVAR+0x26c>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    12a3 <coreGTRGAMMAPROTMULTINVAR+0x173>
jmp    13af <coreGTRGAMMAPROTMULTINVAR+0x27f>
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    1277 <coreGTRGAMMAPROTMULTINVAR+0x147>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    192f <coreGTRGAMMAPROTMULTINVAR+0x7ff>
mov    0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x74(%rbp)
mov    -0x48(%rbp),%rax
imul   $0xe4,-0x74(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%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,-0x98(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movl   $0x1,-0x50(%rbp)
cmpl   $0x14,-0x50(%rbp)
jge    1513 <coreGTRGAMMAPROTMULTINVAR+0x3e3>
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,-0x90(%rbp)
ff
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0xa0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0xa8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1457 <coreGTRGAMMAPROTMULTINVAR+0x327>
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movl   $0x1,-0x50(%rbp)
cmpl   $0x14,-0x50(%rbp)
jge    15f9 <coreGTRGAMMAPROTMULTINVAR+0x4c9>
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,-0x90(%rbp)
ff
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0xa0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0xa8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    153d <coreGTRGAMMAPROTMULTINVAR+0x40d>
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movl   $0x1,-0x50(%rbp)
cmpl   $0x14,-0x50(%rbp)
jge    16df <coreGTRGAMMAPROTMULTINVAR+0x5af>
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,-0x90(%rbp)
ff
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0xa0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0xa8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1623 <coreGTRGAMMAPROTMULTINVAR+0x4f3>
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movl   $0x1,-0x50(%rbp)
cmpl   $0x14,-0x50(%rbp)
jge    17c5 <coreGTRGAMMAPROTMULTINVAR+0x695>
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,-0x90(%rbp)
ff
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0xa0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0xa8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1709 <coreGTRGAMMAPROTMULTINVAR+0x5d9>
mov    -0x80(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
mov    0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
cmpl   $0x14,(%rax,%rcx,4)
jge    1822 <coreGTRGAMMAPROTMULTINVAR+0x6f2>
mov    -0x88(%rbp),%rax
imul   $0x14,-0x74(%rbp),%ecx
mov    0x40(%rbp),%rdx
movslq -0x4c(%rbp),%rsi
add    (%rdx,%rsi,4),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0x7d6(%rip),%xmm0
divsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
mulsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x80(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x80(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0xa0(%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  -0xa8(%rbp),%xmm1
ff
movsd  -0xa0(%rbp),%xmm2
ff
movsd  -0xa0(%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    13c3 <coreGTRGAMMAPROTMULTINVAR+0x293>
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>
mov    -0x88(%rbp),%rdi
call   1030 <free@plt>
mov    -0x80(%rbp),%rdi
call   1030 <free@plt>
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<coreGTRGAMMAPROTMULTINVAR>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
movsd  %xmm0,0x18(%rsp)
mov    %r9,0x40(%rsp)
mov    %r8d,0x34(%rsp)
mov    %ecx,0x30(%rsp)
mov    %rdx,0x38(%rsp)
mov    %rsi,%rbp
mov    %rdi,%r14
mov    0xe0(%rsp),%rax
mov    %rax,0x68(%rsp)
mov    0xd8(%rsp),%r12
mov    0xd0(%rsp),%rbx
mov    0xc8(%rsp),%rax
mov    %rax,0x60(%rsp)
mov    0xc0(%rsp),%r13d
mov    0xb8(%rsp),%rax
mov    %rax,0x50(%rsp)
mov    0xb0(%rsp),%rax
mov    %rax,0x48(%rsp)
imul   $0x720,%r13d,%edi
call   1040 <malloc@plt>
mov    %rax,0x20(%rsp)
mov    %r13d,%eax
shl    $0x5,%eax
lea    (%rax,%rax,4),%edi
call   1040 <malloc@plt>
mov    %rax,%r15
lea    0x0(,%r13,8),%edi
call   1040 <malloc@plt>
mov    %rax,0x28(%rsp)
test   %r13d,%r13d
jle    162c <coreGTRGAMMAPROTMULTINVAR+0x4fc>
mov    %r13d,%esi
mov    $0x13,%eax
xor    %ecx,%ecx
movsd  0xe08(%rip),%xmm0
movsd  0xe08(%rip),%xmm1
mov    0x28(%rsp),%rdx
cs nopw 0x0(%rax,%rax,1)
nop
movapd %xmm0,%xmm2
subsd  (%rbx,%rcx,8),%xmm2
mulsd  %xmm1,%xmm2
movsd  %xmm2,(%rdx,%rcx,8)
movsd  -0x98(%r12,%rax,8),%xmm2
ff ff ff
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x98(%r15,%rax,8)
ff ff ff
movsd  -0x90(%r12,%rax,8),%xmm2
ff ff ff
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x90(%r15,%rax,8)
ff ff ff
movsd  -0x88(%r12,%rax,8),%xmm2
ff ff ff
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x88(%r15,%rax,8)
ff ff ff
movsd  -0x80(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x80(%r15,%rax,8)
movsd  -0x78(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x78(%r15,%rax,8)
movsd  -0x70(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x70(%r15,%rax,8)
movsd  -0x68(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x68(%r15,%rax,8)
movsd  -0x60(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x60(%r15,%rax,8)
movsd  -0x58(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x58(%r15,%rax,8)
movsd  -0x50(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x50(%r15,%rax,8)
movsd  -0x48(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x48(%r15,%rax,8)
movsd  -0x40(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x40(%r15,%rax,8)
movsd  -0x38(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x38(%r15,%rax,8)
movsd  -0x30(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x30(%r15,%rax,8)
movsd  -0x28(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x28(%r15,%rax,8)
movsd  -0x20(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x20(%r15,%rax,8)
movsd  -0x18(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x18(%r15,%rax,8)
movsd  -0x10(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x10(%r15,%rax,8)
movsd  -0x8(%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,-0x8(%r15,%rax,8)
movsd  (%r12,%rax,8),%xmm2
mulsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,(%r15,%rax,8)
add    $0x1,%rcx
add    $0x14,%rax
cmp    %rcx,%rsi
jne    1210 <coreGTRGAMMAPROTMULTINVAR+0xe0>
test   %r13d,%r13d
jle    162c <coreGTRGAMMAPROTMULTINVAR+0x4fc>
xor    %edx,%edx
mov    %rsi,0x58(%rsp)
nop
imul   $0xe4,%edx,%eax
mov    0x20(%rsp),%rcx
lea    (%rcx,%rax,8),%r13
lea    0x0(,%rdx,4),%r12
mov    %rdx,0x70(%rsp)
mov    %rdx,%rax
shl    $0x5,%rax
mov    %r14,%rbx
movsd  (%r14,%rax,1),%xmm0
movsd  %xmm0,0x8(%rsp)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
xor    %r14d,%r14d
movsd  0x18(%rsp),%xmm1
cs nopw 0x0(%rax,%rax,1)
nop
movsd  0x0(%rbp,%r14,8),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x18(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x0(%r13)
movsd  0x0(%rbp,%r14,8),%xmm0
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%r13)
movsd  0x0(%rbp,%r14,8),%xmm0
mulsd  %xmm0,%xmm0
mulsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%r13)
add    $0x18,%r13
add    $0x1,%r14
cmp    $0x13,%r14
jne    1420 <coreGTRGAMMAPROTMULTINVAR+0x2f0>
shl    $0x3,%r12
mov    %r12,%rax
or     $0x8,%rax
mov    %rbx,%r14
movsd  (%rbx,%rax,1),%xmm0
movsd  %xmm0,0x8(%rsp)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rbx,8),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x18(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x0(%r13)
movsd  0x0(%rbp,%rbx,8),%xmm0
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%r13)
movsd  0x0(%rbp,%rbx,8),%xmm0
mulsd  %xmm0,%xmm0
mulsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%r13)
add    $0x18,%r13
add    $0x1,%rbx
cmp    $0x13,%rbx
jne    14b0 <coreGTRGAMMAPROTMULTINVAR+0x380>
mov    %r12,%rax
or     $0x10,%rax
movsd  (%r14,%rax,1),%xmm0
movsd  %xmm0,0x8(%rsp)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
xor    %ebx,%ebx
nopl   (%rax)
movsd  0x0(%rbp,%rbx,8),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x18(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x0(%r13)
movsd  0x0(%rbp,%rbx,8),%xmm0
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%r13)
movsd  0x0(%rbp,%rbx,8),%xmm0
mulsd  %xmm0,%xmm0
mulsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%r13)
add    $0x18,%r13
add    $0x1,%rbx
cmp    $0x13,%rbx
jne    1530 <coreGTRGAMMAPROTMULTINVAR+0x400>
or     $0x18,%r12
movsd  (%r14,%r12,1),%xmm0
movsd  %xmm0,0x8(%rsp)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rbx,8),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x18(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x0(%r13)
movsd  0x0(%rbp,%rbx,8),%xmm0
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%r13)
movsd  0x0(%rbp,%rbx,8),%xmm0
mulsd  %xmm0,%xmm0
mulsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%r13)
add    $0x1,%rbx
add    $0x18,%r13
cmp    $0x13,%rbx
jne    15b0 <coreGTRGAMMAPROTMULTINVAR+0x480>
mov    0x70(%rsp),%rdx
add    $0x1,%rdx
add    $0x98,%rbp
mov    0x58(%rsp),%rsi
cmp    %rsi,%rdx
jne    13d0 <coreGTRGAMMAPROTMULTINVAR+0x2a0>
xorpd  %xmm3,%xmm3
mov    0x34(%rsp),%ecx
mov    0x30(%rsp),%edx
cmp    %ecx,%edx
jge    1914 <coreGTRGAMMAPROTMULTINVAR+0x7e4>
mov    0x20(%rsp),%rax
lea    0x1c8(%rax),%r8
lea    0x390(%rax),%r9
lea    0x558(%rax),%r10
movslq %edx,%r12
movslq %ecx,%r11
lea    (%r12,%r12,4),%rbp
shl    $0x7,%rbp
add    0x38(%rsp),%rbp
xorpd  %xmm4,%xmm4
xorpd  %xmm0,%xmm0
movsd  0x983(%rip),%xmm1
mov    0x40(%rsp),%r13
jmp    16ee <coreGTRGAMMAPROTMULTINVAR+0x5be>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movapd %xmm1,%xmm6
divsd  %xmm3,%xmm6
movapd %xmm5,%xmm3
mulsd  %xmm6,%xmm3
unpckhpd %xmm5,%xmm5
mulsd  %xmm6,%xmm5
mulsd  %xmm4,%xmm3
mulsd  %xmm4,%xmm5
xorps  %xmm4,%xmm4
cvtsi2sdl 0x0(%r13,%r12,4),%xmm4
movapd %xmm3,%xmm6
mulsd  %xmm3,%xmm6
subsd  %xmm6,%xmm5
unpcklpd %xmm4,%xmm4
unpcklpd %xmm5,%xmm3
mulpd  %xmm4,%xmm3
movapd %xmm3,%xmm4
addpd  %xmm2,%xmm4
add    $0x1,%r12
add    $0x280,%rbp
cmp    %r11,%r12
je     1910 <coreGTRGAMMAPROTMULTINVAR+0x7e0>
movapd %xmm4,%xmm2
mov    0x60(%rsp),%rax
movslq (%rax,%r12,4),%r14
imul   $0xe4,%r14,%rbx
mov    0x20(%rsp),%rax
lea    (%rax,%rbx,8),%rdx
lea    (%r12,%r12,4),%rax
shl    $0x7,%rax
mov    0x38(%rsp),%rcx
movsd  (%rcx,%rax,1),%xmm3
addsd  %xmm0,%xmm3
xor    %ecx,%ecx
xorpd  %xmm5,%xmm5
nopw   0x0(%rax,%rax,1)
movsd  (%rdx),%xmm4
mulsd  0x8(%rbp,%rcx,1),%xmm4
addsd  %xmm4,%xmm3
movupd 0x8(%rdx),%xmm6
unpcklpd %xmm4,%xmm4
mulpd  %xmm6,%xmm4
addpd  %xmm5,%xmm4
cmp    $0x90,%ecx
je     1790 <coreGTRGAMMAPROTMULTINVAR+0x660>
movsd  0x18(%rdx),%xmm5
mulsd  0x10(%rbp,%rcx,1),%xmm5
addsd  %xmm5,%xmm3
movupd 0x20(%rdx),%xmm6
unpcklpd %xmm5,%xmm5
mulpd  %xmm6,%xmm5
addpd  %xmm4,%xmm5
add    $0x10,%rcx
add    $0x30,%rdx
jmp    1730 <coreGTRGAMMAPROTMULTINVAR+0x600>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
addsd  0x10(%rbp,%rcx,1),%xmm3
lea    (%r8,%rbx,8),%rdi
lea    (%rcx,%rbp,1),%rax
add    $0x18,%rax
lea    (%r9,%rbx,8),%rdx
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rsi,8),%xmm5
mulsd  (%rax),%xmm5
addsd  %xmm5,%xmm3
movupd 0x8(%rdi,%rsi,8),%xmm6
unpcklpd %xmm5,%xmm5
mulpd  %xmm6,%xmm5
cmp    $0x36,%esi
je     1800 <coreGTRGAMMAPROTMULTINVAR+0x6d0>
addpd  %xmm4,%xmm5
movsd  0x18(%rdi,%rsi,8),%xmm4
mulsd  0x8(%rax),%xmm4
addsd  %xmm4,%xmm3
movupd 0x20(%rdi,%rsi,8),%xmm6
unpcklpd %xmm4,%xmm4
mulpd  %xmm6,%xmm4
addpd  %xmm5,%xmm4
add    $0x6,%rsi
add    $0x10,%rax
jmp    17b0 <coreGTRGAMMAPROTMULTINVAR+0x680>
nop
addsd  0xb0(%rbp,%rcx,1),%xmm3
addpd  %xmm4,%xmm5
lea    (%rcx,%rbp,1),%rax
add    $0xb8,%rax
lea    (%r10,%rbx,8),%rbx
xor    %esi,%esi
nopl   (%rax)
movsd  (%rdx,%rsi,8),%xmm4
mulsd  (%rax),%xmm4
addsd  %xmm4,%xmm3
movupd 0x8(%rdx,%rsi,8),%xmm6
unpcklpd %xmm4,%xmm4
mulpd  %xmm6,%xmm4
addpd  %xmm5,%xmm4
cmp    $0x36,%esi
je     1870 <coreGTRGAMMAPROTMULTINVAR+0x740>
movsd  0x18(%rdx,%rsi,8),%xmm5
mulsd  0x8(%rax),%xmm5
addsd  %xmm5,%xmm3
movupd 0x20(%rdx,%rsi,8),%xmm6
unpcklpd %xmm5,%xmm5
mulpd  %xmm6,%xmm5
addpd  %xmm4,%xmm5
add    $0x6,%rsi
add    $0x10,%rax
jmp    1820 <coreGTRGAMMAPROTMULTINVAR+0x6f0>
nop
addsd  0x150(%rbp,%rcx,1),%xmm3
lea    (%rcx,%rbp,1),%rax
add    $0x158,%rax
xor    %ecx,%ecx
mov    0x28(%rsp),%rdx
nopw   0x0(%rax,%rax,1)
movsd  (%rbx,%rcx,8),%xmm5
mulsd  (%rax),%xmm5
addsd  %xmm5,%xmm3
movupd 0x8(%rbx,%rcx,8),%xmm6
unpcklpd %xmm5,%xmm5
mulpd  %xmm6,%xmm5
addpd  %xmm4,%xmm5
cmp    $0x36,%ecx
je     18e0 <coreGTRGAMMAPROTMULTINVAR+0x7b0>
movsd  0x18(%rbx,%rcx,8),%xmm4
mulsd  0x8(%rax),%xmm4
addsd  %xmm4,%xmm3
movupd 0x20(%rbx,%rcx,8),%xmm6
unpcklpd %xmm4,%xmm4
mulpd  %xmm6,%xmm4
addpd  %xmm5,%xmm4
add    $0x6,%rcx
add    $0x10,%rax
jmp    1890 <coreGTRGAMMAPROTMULTINVAR+0x760>
nop
movsd  (%rdx,%r14,8),%xmm4
mulsd  %xmm4,%xmm3
mov    0x68(%rsp),%rax
mov    (%rax,%r12,4),%eax
cmp    $0x13,%eax
jg     1690 <coreGTRGAMMAPROTMULTINVAR+0x560>
lea    (%r14,%r14,4),%ecx
lea    (%rax,%rcx,4),%eax
cltq
addsd  (%r15,%rax,8),%xmm3
jmp    1690 <coreGTRGAMMAPROTMULTINVAR+0x560>
addpd  %xmm2,%xmm3
mov    0x48(%rsp),%rax
movlpd %xmm3,(%rax)
mov    0x50(%rsp),%rax
movhpd %xmm3,(%rax)
mov    0x20(%rsp),%rdi
call   1030 <free@plt>
mov    %r15,%rdi
call   1030 <free@plt>
mov    0x28(%rsp),%rdi
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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