void categorizeGTRGAMMAPROTINVAR(char *tipX1, double *x2, int *ex2, double *pz,
                                 double *gammaRates, double *EIGN,
                                 double *tipVector, int *rateCategory,
                                 int lower, int upper, int *iptr,
                                 double *invariants, double *frequencies) {
  double z, lz, term[4], ki, min, allSum;
  int i, j, index = 0;
  double *diagptable, *diagptable_start;
  double *left, *right;
  double scaler = (1.0 - invariants[0]);
  double freqs[20];
  for (i = 0; i < 20; i++)
    freqs[i] = frequencies[i] * invariants[0];
  z = pz[0];
  if (z < 1.0E-15)
    z = 1.0E-15;
  lz = log(z);
  diagptable = diagptable_start = (double *)malloc(sizeof(double) * 4 * 19);
  for (i = 0; i < 4; i++) {
    ki = gammaRates[i];
    for (j = 0; j < 19; j++)
      *diagptable++ = exp(EIGN[j] * ki * lz);
  }
  for (i = lower; i < upper; i++) {
    left = &(tipVector[20 * tipX1[i]]);
    diagptable = diagptable_start;
    for (j = 0; j < 4; j++) {
      right = &(x2[80 * i + 20 * j]);
      term[j] = 0;
      term[j] += left[0] * right[0];
      term[j] += left[1] * right[1] * *diagptable++;
      term[j] += left[2] * right[2] * *diagptable++;
      term[j] += left[3] * right[3] * *diagptable++;
      term[j] += left[4] * right[4] * *diagptable++;
      term[j] += left[5] * right[5] * *diagptable++;
      term[j] += left[6] * right[6] * *diagptable++;
      term[j] += left[7] * right[7] * *diagptable++;
      term[j] += left[8] * right[8] * *diagptable++;
      term[j] += left[9] * right[9] * *diagptable++;
      term[j] += left[10] * right[10] * *diagptable++;
      term[j] += left[11] * right[11] * *diagptable++;
      term[j] += left[12] * right[12] * *diagptable++;
      term[j] += left[13] * right[13] * *diagptable++;
      term[j] += left[14] * right[14] * *diagptable++;
      term[j] += left[15] * right[15] * *diagptable++;
      term[j] += left[16] * right[16] * *diagptable++;
      term[j] += left[17] * right[17] * *diagptable++;
      term[j] += left[18] * right[18] * *diagptable++;
      term[j] += left[19] * right[19] * *diagptable++;
      if (iptr[i] < 20)
        term[j] = log(
            (scaler * term[j] + freqs[iptr[i]]) *
            pow((1.0 /
                 115792089237316195423570985008687907853269984665640564039457584007913129639936.0),
                ex2[i]));
      else
        term[j] =
            log(scaler * term[j]) +
            ex2[i] *
                log((
                    1.0 /
                    115792089237316195423570985008687907853269984665640564039457584007913129639936.0));
    }
    allSum = 0.25 * (term[0] + term[1] + term[2] + term[3]);
    allSum = 1.0 / allSum;
    term[0] = 0.25 * term[0] * allSum;
    term[1] = 0.25 * term[1] * allSum;
    term[2] = 0.25 * term[2] * allSum;
    term[3] = 0.25 * term[3] * allSum;
    min = 1.0E300;
    if (term[0] < min) {
      min = term[0];
      index = 0;
    }
    if (term[1] < min) {
      min = term[1];
      index = 1;
    }
    if (term[2] < min) {
      min = term[2];
      index = 2;
    }
    if (term[3] < min) {
      min = term[3];
      index = 3;
    }
    rateCategory[i] = index;
  }
  free(diagptable_start);
  return;
}


x86 assembly (clang, -O0):
<categorizeGTRGAMMAPROTINVAR>:
push   %rbp
mov    %rsp,%rbp
sub    $0x170,%rsp
mov    0x40(%rbp),%rax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%eax
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    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x84(%rbp)
mov    0x38(%rbp),%rax
movsd  0xe5d(%rip),%xmm0
subsd  (%rax),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movl   $0x0,-0x7c(%rbp)
cmpl   $0x14,-0x7c(%rbp)
jge    11f0 <categorizeGTRGAMMAPROTINVAR+0xa0>
mov    0x40(%rbp),%rax
movslq -0x7c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x38(%rbp),%rax
mulsd  (%rax),%xmm0
movslq -0x7c(%rbp),%rax
movsd  %xmm0,-0x150(%rbp,%rax,8)
ff ff
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    11b6 <categorizeGTRGAMMAPROTINVAR+0x66>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe03(%rip),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    121d <categorizeGTRGAMMAPROTINVAR+0xcd>
movsd  0xdf0(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1060 <log@plt>
movsd  %xmm0,-0x40(%rbp)
mov    $0x260,%edi
call   1040 <malloc@plt>
mov    %rax,-0x98(%rbp)
mov    %rax,-0x90(%rbp)
movl   $0x0,-0x7c(%rbp)
cmpl   $0x4,-0x7c(%rbp)
jge    12ce <categorizeGTRGAMMAPROTINVAR+0x17e>
mov    -0x28(%rbp),%rax
movslq -0x7c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x80(%rbp)
cmpl   $0x13,-0x80(%rbp)
jge    12bb <categorizeGTRGAMMAPROTINVAR+0x16b>
mov    -0x30(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x68(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
call   1050 <exp@plt>
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  %xmm0,(%rax)
mov    -0x80(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x80(%rbp)
jmp    126e <categorizeGTRGAMMAPROTINVAR+0x11e>
jmp    12c0 <categorizeGTRGAMMAPROTINVAR+0x170>
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    124b <categorizeGTRGAMMAPROTINVAR+0xfb>
mov    0x20(%rbp),%eax
mov    %eax,-0x7c(%rbp)
mov    -0x7c(%rbp),%eax
cmp    0x28(%rbp),%eax
jge    1b5e <categorizeGTRGAMMAPROTINVAR+0xa0e>
mov    0x10(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x7c(%rbp),%rdx
movsbl (%rcx,%rdx,1),%ecx
imul   $0x14,%ecx,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0xa0(%rbp)
mov    -0x98(%rbp),%rax
mov    %rax,-0x90(%rbp)
movl   $0x0,-0x80(%rbp)
cmpl   $0x4,-0x80(%rbp)
jge    19ff <categorizeGTRGAMMAPROTINVAR+0x8af>
mov    -0x10(%rbp),%rax
imul   $0x50,-0x7c(%rbp),%ecx
imul   $0x14,-0x80(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0xa8(%rbp)
movslq -0x80(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0xa8(%rbp),%rax
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x10(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x18(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x20(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x20(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x28(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x30(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x30(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x38(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x38(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x40(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x40(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x48(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x48(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x50(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x50(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x58(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x58(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x60(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x60(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x68(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x68(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x70(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x70(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x78(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x78(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x80(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x80(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x88(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x88(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x90(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x90(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xa0(%rbp),%rax
movsd  0x98(%rax),%xmm0
mov    -0xa8(%rbp),%rax
mulsd  0x98(%rax),%xmm0
mov    -0x90(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x90(%rbp)
movsd  (%rax),%xmm2
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    0x30(%rbp),%rax
movslq -0x7c(%rbp),%rcx
cmpl   $0x14,(%rax,%rcx,4)
jge    1986 <categorizeGTRGAMMAPROTINVAR+0x836>
movsd  -0xb0(%rbp),%xmm0
ff
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm2
mov    0x30(%rbp),%rax
movslq -0x7c(%rbp),%rcx
movslq (%rax,%rcx,4),%rax
movsd  -0x150(%rbp,%rax,8),%xmm1
ff ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x158(%rbp)
ff
mov    -0x18(%rbp),%rax
movslq -0x7c(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm1
movsd  0x6c2(%rip),%xmm0
call   1070 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0x158(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
call   1060 <log@plt>
movslq -0x80(%rbp),%rax
movsd  %xmm0,-0x60(%rbp,%rax,8)
jmp    19ec <categorizeGTRGAMMAPROTINVAR+0x89c>
movsd  -0xb0(%rbp),%xmm0
ff
movslq -0x80(%rbp),%rax
mulsd  -0x60(%rbp,%rax,8),%xmm0
call   1060 <log@plt>
movsd  %xmm0,-0x168(%rbp)
ff
mov    -0x18(%rbp),%rax
movslq -0x7c(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  %xmm0,-0x160(%rbp)
ff
movsd  0x65e(%rip),%xmm0
call   1060 <log@plt>
movsd  -0x168(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x160(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movslq -0x80(%rbp),%rax
movsd  %xmm0,-0x60(%rbp,%rax,8)
jmp    19f1 <categorizeGTRGAMMAPROTINVAR+0x8a1>
mov    -0x80(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x80(%rbp)
jmp    1319 <categorizeGTRGAMMAPROTINVAR+0x1c9>
movsd  -0x60(%rbp),%xmm1
addsd  -0x58(%rbp),%xmm1
addsd  -0x50(%rbp),%xmm1
addsd  -0x48(%rbp),%xmm1
movsd  0x5fd(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0x5d4(%rip),%xmm0
divsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0x5da(%rip),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  0x5c3(%rip),%xmm0
mulsd  -0x58(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0x5ac(%rip),%xmm0
mulsd  -0x50(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0x595(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0x576(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x60(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1ac7 <categorizeGTRGAMMAPROTINVAR+0x977>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x84(%rbp)
movsd  -0x58(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1aef <categorizeGTRGAMMAPROTINVAR+0x99f>
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x1,-0x84(%rbp)
movsd  -0x50(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1b17 <categorizeGTRGAMMAPROTINVAR+0x9c7>
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x2,-0x84(%rbp)
movsd  -0x48(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1b3f <categorizeGTRGAMMAPROTINVAR+0x9ef>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x3,-0x84(%rbp)
mov    -0x84(%rbp),%edx
mov    0x18(%rbp),%rax
movslq -0x7c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    12d4 <categorizeGTRGAMMAPROTINVAR+0x184>
mov    -0x98(%rbp),%rdi
call   1030 <free@plt>
add    $0x170,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<categorizeGTRGAMMAPROTINVAR>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x128,%rsp
mov    %r9,%rbx
mov    %r8,%r14
mov    %rdx,0x18(%rsp)
mov    %rsi,0x10(%rsp)
mov    %rdi,0x50(%rsp)
mov    0x180(%rsp),%rax
mov    %rax,0x58(%rsp)
mov    0x178(%rsp),%r12d
mov    0x170(%rsp),%r15d
mov    0x168(%rsp),%rax
mov    %rax,0x48(%rsp)
mov    0x160(%rsp),%rax
mov    %rax,0x40(%rsp)
mov    0x190(%rsp),%rax
mov    0x188(%rsp),%rdx
movsd  (%rdx),%xmm0
movupd (%rax),%xmm1
movapd %xmm0,0x20(%rsp)
unpcklpd %xmm0,%xmm0
mulpd  %xmm0,%xmm1
movapd %xmm1,0x80(%rsp)
movupd 0x10(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0x90(%rsp)
movupd 0x20(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0xa0(%rsp)
movupd 0x30(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0xb0(%rsp)
movupd 0x40(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0xc0(%rsp)
movupd 0x50(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0xd0(%rsp)
movupd 0x60(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0xe0(%rsp)
movupd 0x70(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0xf0(%rsp)
movupd 0x80(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0x100(%rsp)
movupd 0x90(%rax),%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,0x110(%rsp)
movsd  0xd74(%rip),%xmm0
maxsd  (%rcx),%xmm0
call   1060 <log@plt>
movsd  %xmm0,(%rsp)
mov    $0x260,%edi
call   1040 <malloc@plt>
mov    %rax,%r13
xor    %ebp,%ebp
nopl   0x0(%rax)
movsd  (%r14),%xmm1
movsd  %xmm1,0x8(%rsp)
movsd  (%rbx),%xmm0
mulsd  %xmm1,%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x0(%r13,%rbp,8)
movsd  0x8(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x8(%r13,%rbp,8)
movsd  0x10(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x10(%r13,%rbp,8)
movsd  0x18(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x18(%r13,%rbp,8)
movsd  0x20(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x20(%r13,%rbp,8)
movsd  0x28(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x28(%r13,%rbp,8)
movsd  0x30(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x30(%r13,%rbp,8)
movsd  0x38(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x38(%r13,%rbp,8)
movsd  0x40(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x40(%r13,%rbp,8)
movsd  0x48(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x48(%r13,%rbp,8)
movsd  0x50(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x50(%r13,%rbp,8)
movsd  0x58(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x58(%r13,%rbp,8)
movsd  0x60(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x60(%r13,%rbp,8)
movsd  0x68(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x68(%r13,%rbp,8)
movsd  0x70(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x70(%r13,%rbp,8)
movsd  0x78(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x78(%r13,%rbp,8)
movsd  0x80(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x80(%r13,%rbp,8)
movsd  0x88(%rbx),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x88(%r13,%rbp,8)
movsd  0x8(%rsp),%xmm0
mulsd  0x90(%rbx),%xmm0
mulsd  (%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x90(%r13,%rbp,8)
add    $0x8,%r14
add    $0x13,%rbp
cmp    $0x4c,%rbp
jne    12b0 <categorizeGTRGAMMAPROTINVAR+0x160>
cmp    %r12d,%r15d
jge    18d8 <categorizeGTRGAMMAPROTINVAR+0x788>
movsd  0xb07(%rip),%xmm0
subsd  0x20(%rsp),%xmm0
movsd  %xmm0,(%rsp)
movslq %r15d,%rbx
movslq %r12d,%rax
mov    %rax,0x20(%rsp)
lea    (%rbx,%rbx,4),%rax
shl    $0x7,%rax
mov    0x10(%rsp),%rcx
lea    (%rax,%rcx,1),%r12
add    $0x98,%r12
xor    %ebp,%ebp
mov    %r13,0x38(%rsp)
jmp    1620 <categorizeGTRGAMMAPROTINVAR+0x4d0>
nopl   0x0(%rax,%rax,1)
movapd 0x60(%rsp),%xmm1
movapd 0x70(%rsp),%xmm0
movapd %xmm1,%xmm2
unpckhpd %xmm1,%xmm2
addsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
movapd %xmm0,%xmm3
unpckhpd %xmm0,%xmm3
addsd  %xmm2,%xmm3
mulsd  0xab0(%rip),%xmm3
movsd  0xa90(%rip),%xmm2
divsd  %xmm3,%xmm2
movapd 0xabc(%rip),%xmm4
mulpd  %xmm4,%xmm1
unpcklpd %xmm2,%xmm2
mulpd  %xmm2,%xmm1
movsd  0xa90(%rip),%xmm5
ucomisd %xmm1,%xmm5
movapd %xmm1,0x60(%rsp)
mov    $0x0,%eax
cmova  %eax,%ebp
movapd %xmm1,%xmm3
minsd  %xmm5,%xmm3
unpckhpd %xmm1,%xmm1
ucomisd %xmm1,%xmm3
mulpd  %xmm4,%xmm0
mov    $0x1,%eax
cmova  %eax,%ebp
mulpd  %xmm2,%xmm0
minsd  %xmm3,%xmm1
ucomisd %xmm0,%xmm1
movapd %xmm0,0x70(%rsp)
mov    $0x2,%eax
cmova  %eax,%ebp
movapd %xmm0,%xmm2
minsd  %xmm1,%xmm2
unpckhpd %xmm0,%xmm0
ucomisd %xmm0,%xmm2
mov    $0x3,%eax
cmova  %eax,%ebp
mov    0x48(%rsp),%rax
mov    %ebp,(%rax,%rbx,4)
add    $0x1,%rbx
mov    0x10(%rsp),%r12
add    $0x280,%r12
cmp    0x20(%rsp),%rbx
mov    0x38(%rsp),%r13
je     18d8 <categorizeGTRGAMMAPROTINVAR+0x788>
mov    0x50(%rsp),%rax
movsbq (%rax,%rbx,1),%rax
lea    (%rax,%rax,4),%r14
shl    $0x5,%r14
add    0x40(%rsp),%r14
mov    %r12,0x10(%rsp)
xor    %r15d,%r15d
jmp    169a <categorizeGTRGAMMAPROTINVAR+0x54a>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mulsd  (%rsp),%xmm0
call   1060 <log@plt>
movapd %xmm0,%xmm1
mov    0x18(%rsp),%rax
xorps  %xmm0,%xmm0
cvtsi2sdl (%rax,%rbx,4),%xmm0
mulsd  0x99d(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x60(%rsp,%r15,8)
add    $0x1,%r15
add    $0xa0,%r12
add    $0x98,%r13
cmp    $0x4,%r15
je     1540 <categorizeGTRGAMMAPROTINVAR+0x3f0>
movsd  (%r14),%xmm0
movsd  0x8(%r14),%xmm1
mulsd  -0x98(%r12),%xmm0
ff ff ff
addsd  0x979(%rip),%xmm0
mulsd  -0x90(%r12),%xmm1
ff ff ff
mulsd  0x0(%r13),%xmm1
addsd  %xmm0,%xmm1
movsd  0x10(%r14),%xmm0
mulsd  -0x88(%r12),%xmm0
ff ff ff
mulsd  0x8(%r13),%xmm0
addsd  %xmm1,%xmm0
movsd  0x18(%r14),%xmm1
mulsd  -0x80(%r12),%xmm1
mulsd  0x10(%r13),%xmm1
addsd  %xmm0,%xmm1
movsd  0x20(%r14),%xmm0
mulsd  -0x78(%r12),%xmm0
mulsd  0x18(%r13),%xmm0
addsd  %xmm1,%xmm0
movsd  0x28(%r14),%xmm1
mulsd  -0x70(%r12),%xmm1
mulsd  0x20(%r13),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x60(%rsp,%r15,8)
movsd  0x30(%r14),%xmm0
movsd  0x38(%r14),%xmm2
mulsd  -0x68(%r12),%xmm0
mulsd  0x28(%r13),%xmm0
mulsd  -0x60(%r12),%xmm2
mulsd  0x30(%r13),%xmm2
addsd  %xmm1,%xmm0
movsd  0x40(%r14),%xmm1
mulsd  -0x58(%r12),%xmm1
addsd  %xmm0,%xmm2
mulsd  0x38(%r13),%xmm1
movsd  0x48(%r14),%xmm0
mulsd  -0x50(%r12),%xmm0
mulsd  0x40(%r13),%xmm0
addsd  %xmm2,%xmm1
movsd  0x50(%r14),%xmm2
mulsd  -0x48(%r12),%xmm2
addsd  %xmm1,%xmm0
mulsd  0x48(%r13),%xmm2
movsd  0x58(%r14),%xmm1
mulsd  -0x40(%r12),%xmm1
mulsd  0x50(%r13),%xmm1
addsd  %xmm0,%xmm2
movsd  0x60(%r14),%xmm0
mulsd  -0x38(%r12),%xmm0
addsd  %xmm2,%xmm1
mulsd  0x58(%r13),%xmm0
movsd  0x68(%r14),%xmm2
mulsd  -0x30(%r12),%xmm2
mulsd  0x60(%r13),%xmm2
addsd  %xmm1,%xmm0
movsd  0x70(%r14),%xmm1
mulsd  -0x28(%r12),%xmm1
addsd  %xmm0,%xmm2
mulsd  0x68(%r13),%xmm1
movsd  0x78(%r14),%xmm0
mulsd  -0x20(%r12),%xmm0
mulsd  0x70(%r13),%xmm0
addsd  %xmm2,%xmm1
movsd  0x80(%r14),%xmm2
mulsd  -0x18(%r12),%xmm2
addsd  %xmm1,%xmm0
mulsd  0x78(%r13),%xmm2
movsd  0x88(%r14),%xmm1
mulsd  -0x10(%r12),%xmm1
mulsd  0x80(%r13),%xmm1
addsd  %xmm0,%xmm2
movsd  0x90(%r14),%xmm3
mulsd  -0x8(%r12),%xmm3
addsd  %xmm2,%xmm1
mulsd  0x88(%r13),%xmm3
movsd  0x98(%r14),%xmm0
mulsd  (%r12),%xmm0
mulsd  0x90(%r13),%xmm0
addsd  %xmm1,%xmm3
addsd  %xmm3,%xmm0
mov    0x58(%rsp),%rax
movslq (%rax,%rbx,4),%rax
cmp    $0x13,%rax
jg     1650 <categorizeGTRGAMMAPROTINVAR+0x500>
mulsd  (%rsp),%xmm0
addsd  0x80(%rsp,%rax,8),%xmm0
movsd  %xmm0,0x8(%rsp)
mov    0x18(%rsp),%rax
xorps  %xmm1,%xmm1
cvtsi2sdl (%rax,%rbx,4),%xmm1
movsd  0x755(%rip),%xmm0
call   1070 <pow@plt>
mulsd  0x8(%rsp),%xmm0
call   1060 <log@plt>
jmp    1677 <categorizeGTRGAMMAPROTINVAR+0x527>
mov    %r13,%rdi
add    $0x128,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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