void categorizeGTRGAMMAPROTMULT(char *tipX1, double *x2, int *ex2, double *pz,
                                double *gammaRates, double *EIGN,
                                double *tipVector, int *rateCategory, int lower,
                                int upper, int *modelptr, int NumberOfModels,
                                int multiBranch) {
  double z, lz = 0.0, term[4], ki, min, allSum;
  int i, j, index = 0;
  double *diagptable, *diagptable_start;
  double *left, *right;
  int model;
  if (!multiBranch) {
    z = pz[0];
    if (z < 1.0E-15)
      z = 1.0E-15;
    lz = log(z);
  }
  diagptable = diagptable_start =
      (double *)malloc(sizeof(double) * 4 * 19 * NumberOfModels);
  for (model = 0; model < NumberOfModels; model++) {
    if (multiBranch) {
      z = pz[model];
      if (z < 1.0E-15)
        z = 1.0E-15;
      lz = log(z);
    }
    diagptable = &diagptable_start[4 * 19 * model];
    for (i = 0; i < 4; i++) {
      ki = gammaRates[model * 4 + i];
      for (j = 0; j < 19; j++)
        *diagptable++ = exp(EIGN[model * 19 + j] * ki * lz);
    }
  }
  for (i = lower; i < upper; i++) {
    model = modelptr[i];
    left = &(tipVector[460 * model + 20 * tipX1[i]]);
    diagptable = &diagptable_start[76 * model];
    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++;
      term[j] =
          log(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):
<categorizeGTRGAMMAPROTMULT>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    0x40(%rbp),%eax
mov    0x38(%rbp),%eax
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)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x84(%rbp)
cmpl   $0x0,0x40(%rbp)
jne    11d3 <categorizeGTRGAMMAPROTMULT+0x93>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe54(%rip),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    11c4 <categorizeGTRGAMMAPROTMULT+0x84>
movsd  0xe41(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1060 <log@plt>
movsd  %xmm0,-0x40(%rbp)
movslq 0x38(%rbp),%rax
imul   $0x260,%rax,%rdi
call   1040 <malloc@plt>
mov    %rax,-0x98(%rbp)
mov    %rax,-0x90(%rbp)
movl   $0x0,-0xac(%rbp)
mov    -0xac(%rbp),%eax
cmp    0x38(%rbp),%eax
jge    132e <categorizeGTRGAMMAPROTMULT+0x1ee>
cmpl   $0x0,0x40(%rbp)
je     1258 <categorizeGTRGAMMAPROTMULT+0x118>
mov    -0x20(%rbp),%rax
movslq -0xac(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xdcf(%rip),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    1249 <categorizeGTRGAMMAPROTMULT+0x109>
movsd  0xdbc(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1060 <log@plt>
movsd  %xmm0,-0x40(%rbp)
mov    -0x98(%rbp),%rax
imul   $0x4c,-0xac(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x90(%rbp)
movl   $0x0,-0x7c(%rbp)
cmpl   $0x4,-0x7c(%rbp)
jge    1315 <categorizeGTRGAMMAPROTMULT+0x1d5>
mov    -0x28(%rbp),%rax
mov    -0xac(%rbp),%ecx
shl    $0x2,%ecx
add    -0x7c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x80(%rbp)
cmpl   $0x13,-0x80(%rbp)
jge    1302 <categorizeGTRGAMMAPROTMULT+0x1c2>
mov    -0x30(%rbp),%rax
imul   $0x13,-0xac(%rbp),%ecx
add    -0x80(%rbp),%ecx
movslq %ecx,%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    12ac <categorizeGTRGAMMAPROTMULT+0x16c>
jmp    1307 <categorizeGTRGAMMAPROTMULT+0x1c7>
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    127e <categorizeGTRGAMMAPROTMULT+0x13e>
jmp    131a <categorizeGTRGAMMAPROTMULT+0x1da>
mov    -0xac(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xac(%rbp)
jmp    11fb <categorizeGTRGAMMAPROTMULT+0xbb>
mov    0x20(%rbp),%eax
mov    %eax,-0x7c(%rbp)
mov    -0x7c(%rbp),%eax
cmp    0x28(%rbp),%eax
jge    1b59 <categorizeGTRGAMMAPROTMULT+0xa19>
mov    0x30(%rbp),%rax
movslq -0x7c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0xac(%rbp)
mov    0x10(%rbp),%rax
imul   $0x1cc,-0xac(%rbp),%ecx
mov    -0x8(%rbp),%rdx
movslq -0x7c(%rbp),%rsi
movsbl (%rdx,%rsi,1),%edx
imul   $0x14,%edx,%edx
add    %edx,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0xa0(%rbp)
mov    -0x98(%rbp),%rax
imul   $0x4c,-0xac(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x90(%rbp)
movl   $0x0,-0x80(%rbp)
cmpl   $0x4,-0x80(%rbp)
jge    19fa <categorizeGTRGAMMAPROTMULT+0x8ba>
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)
movslq -0x80(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm0
call   1060 <log@plt>
movsd  %xmm0,-0xc0(%rbp)
ff
mov    -0x18(%rbp),%rax
movslq -0x7c(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  0x65e(%rip),%xmm0
call   1060 <log@plt>
movsd  -0xc0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xb8(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movslq -0x80(%rbp),%rax
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0x80(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x80(%rbp)
jmp    13a7 <categorizeGTRGAMMAPROTMULT+0x267>
movsd  -0x60(%rbp),%xmm1
addsd  -0x58(%rbp),%xmm1
addsd  -0x50(%rbp),%xmm1
addsd  -0x48(%rbp),%xmm1
movsd  0x5fa(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0x5f1(%rip),%xmm0
divsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0x5d7(%rip),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  0x5c0(%rip),%xmm0
mulsd  -0x58(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0x5a9(%rip),%xmm0
mulsd  -0x50(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0x592(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0x573(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x60(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1ac2 <categorizeGTRGAMMAPROTMULT+0x982>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x84(%rbp)
movsd  -0x58(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1aea <categorizeGTRGAMMAPROTMULT+0x9aa>
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x1,-0x84(%rbp)
movsd  -0x50(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1b12 <categorizeGTRGAMMAPROTMULT+0x9d2>
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x2,-0x84(%rbp)
movsd  -0x48(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1b3a <categorizeGTRGAMMAPROTMULT+0x9fa>
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    1334 <categorizeGTRGAMMAPROTMULT+0x1f4>
mov    -0x98(%rbp),%rdi
call   1030 <free@plt>
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<categorizeGTRGAMMAPROTMULT>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
mov    %r9,0x78(%rsp)
mov    %r8,%r14
mov    %rcx,0x30(%rsp)
mov    %rdx,%rbp
mov    %rsi,0x38(%rsp)
mov    %rdi,0x58(%rsp)
mov    0x110(%rsp),%ecx
mov    0x108(%rsp),%ebx
mov    0x100(%rsp),%rax
mov    %rax,0x50(%rsp)
mov    0xf8(%rsp),%eax
mov    %eax,0x28(%rsp)
mov    0xf0(%rsp),%eax
mov    %eax,0x24(%rsp)
mov    0xe8(%rsp),%rax
mov    %rax,0x48(%rsp)
mov    0xe0(%rsp),%rax
mov    %rax,0x40(%rsp)
xorpd  %xmm0,%xmm0
mov    %ecx,0x2c(%rsp)
test   %ecx,%ecx
jne    11d8 <categorizeGTRGAMMAPROTMULT+0x98>
movsd  0xe36(%rip),%xmm0
mov    0x30(%rsp),%rax
maxsd  (%rax),%xmm0
call   1060 <log@plt>
movsd  %xmm0,0x8(%rsp)
mov    %rbp,0x60(%rsp)
movslq %ebx,%rbp
imul   $0x260,%rbp,%rdi
call   1040 <malloc@plt>
movsd  0x8(%rsp),%xmm1
mov    %rax,0x10(%rsp)
test   %ebp,%ebp
jle    13a3 <categorizeGTRGAMMAPROTMULT+0x263>
mov    %ebx,%eax
mov    %rax,0x68(%rsp)
xor    %ebp,%ebp
xor    %r13d,%r13d
mov    %r14,0x70(%rsp)
cs nopw 0x0(%rax,%rax,1)
mov    %ebp,%ebx
and    $0x3fffffff,%ebx
shl    $0x5,%rbx
add    0x10(%rsp),%rbx
cmpl   $0x0,0x2c(%rsp)
je     1254 <categorizeGTRGAMMAPROTMULT+0x114>
movsd  0xdc0(%rip),%xmm0
mov    0x30(%rsp),%rax
maxsd  (%rax,%r13,8),%xmm0
call   1060 <log@plt>
movapd %xmm0,%xmm1
mov    %rbp,0x18(%rsp)
mov    %ebp,%eax
mov    0x78(%rsp),%rcx
lea    (%rcx,%rax,8),%r15
lea    0x0(,%r13,4),%r12d
movsd  (%r14,%r12,8),%xmm0
movsd  %xmm0,(%rsp)
xor    %ebp,%ebp
movsd  %xmm1,0x8(%rsp)
nop
movsd  (%r15,%rbp,1),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x8(%rsp),%xmm1
movsd  %xmm0,(%rbx,%rbp,1)
add    $0x8,%rbp
cmp    $0x98,%rbp
jne    1280 <categorizeGTRGAMMAPROTMULT+0x140>
shl    $0x3,%r12
mov    %r12,%rax
or     $0x8,%rax
movsd  (%r14,%rax,1),%xmm0
movsd  %xmm0,(%rsp)
add    %rbp,%rbx
mov    $0x13,%ebp
mov    %r15,%r14
nopl   (%rax)
movsd  (%r14),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x8(%rsp),%xmm1
movsd  %xmm0,(%rbx)
add    $0x8,%rbx
add    $0x8,%r14
add    $0xffffffffffffffff,%rbp
jne    12d0 <categorizeGTRGAMMAPROTMULT+0x190>
mov    %r12,%rax
or     $0x10,%rax
mov    0x70(%rsp),%r14
movsd  (%r14,%rax,1),%xmm0
movsd  %xmm0,(%rsp)
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%r15,%rbp,8),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x8(%rsp),%xmm1
movsd  %xmm0,(%rbx)
add    $0x8,%rbx
add    $0x1,%rbp
cmp    $0x13,%rbp
jne    1320 <categorizeGTRGAMMAPROTMULT+0x1e0>
or     $0x18,%r12
movsd  (%r14,%r12,1),%xmm0
movsd  %xmm0,(%rsp)
xor    %ebp,%ebp
nopl   (%rax)
movsd  (%r15,%rbp,8),%xmm0
mulsd  (%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x8(%rsp),%xmm1
movsd  %xmm0,(%rbx)
add    $0x8,%rbx
add    $0x1,%rbp
cmp    $0x13,%rbp
jne    1360 <categorizeGTRGAMMAPROTMULT+0x220>
add    $0x1,%r13
mov    0x18(%rsp),%rbp
add    $0x13,%ebp
cmp    0x68(%rsp),%r13
jne    1220 <categorizeGTRGAMMAPROTMULT+0xe0>
mov    0x28(%rsp),%eax
mov    0x24(%rsp),%ecx
cmp    %eax,%ecx
jge    172c <categorizeGTRGAMMAPROTMULT+0x5ec>
movslq %ecx,%rbp
cltq
mov    %rax,0x18(%rsp)
lea    0x0(,%rbp,4),%rax
add    %rbp,%rax
shl    $0x7,%rax
mov    0x38(%rsp),%rcx
lea    (%rax,%rcx,1),%r13
add    $0x98,%r13
xorpd  %xmm4,%xmm4
xor    %ecx,%ecx
mov    0x60(%rsp),%r14
nopw   0x0(%rax,%rax,1)
mov    %ecx,(%rsp)
mov    0x50(%rsp),%rax
movslq (%rax,%rbp,4),%rax
imul   $0x1cc,%rax,%rcx
mov    0x58(%rsp),%rdx
movsbq (%rdx,%rbp,1),%rdx
lea    (%rdx,%rdx,4),%rdx
lea    (%rcx,%rdx,4),%rcx
mov    0x40(%rsp),%rdx
lea    (%rdx,%rcx,8),%rbx
imul   $0x260,%rax,%r15
add    0x10(%rsp),%r15
mov    %r13,0x8(%rsp)
xor    %r12d,%r12d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rbx),%xmm0
movsd  0x8(%rbx),%xmm1
mulsd  -0x98(%r13),%xmm0
ff ff
mulsd  -0x90(%r13),%xmm1
ff ff
addsd  %xmm4,%xmm0
mulsd  (%r15),%xmm1
movsd  0x10(%rbx),%xmm2
mulsd  -0x88(%r13),%xmm2
ff ff
mulsd  0x8(%r15),%xmm2
addsd  %xmm0,%xmm1
movsd  0x18(%rbx),%xmm0
mulsd  -0x80(%r13),%xmm0
addsd  %xmm1,%xmm2
mulsd  0x10(%r15),%xmm0
movsd  0x20(%rbx),%xmm1
mulsd  -0x78(%r13),%xmm1
mulsd  0x18(%r15),%xmm1
addsd  %xmm2,%xmm0
movsd  0x28(%rbx),%xmm2
mulsd  -0x70(%r13),%xmm2
addsd  %xmm0,%xmm1
mulsd  0x20(%r15),%xmm2
movsd  0x30(%rbx),%xmm0
mulsd  -0x68(%r13),%xmm0
mulsd  0x28(%r15),%xmm0
addsd  %xmm1,%xmm2
movsd  0x38(%rbx),%xmm1
mulsd  -0x60(%r13),%xmm1
addsd  %xmm2,%xmm0
mulsd  0x30(%r15),%xmm1
movsd  0x40(%rbx),%xmm2
mulsd  -0x58(%r13),%xmm2
mulsd  0x38(%r15),%xmm2
addsd  %xmm0,%xmm1
movsd  0x48(%rbx),%xmm0
mulsd  -0x50(%r13),%xmm0
addsd  %xmm1,%xmm2
mulsd  0x40(%r15),%xmm0
movsd  0x50(%rbx),%xmm1
mulsd  -0x48(%r13),%xmm1
mulsd  0x48(%r15),%xmm1
addsd  %xmm2,%xmm0
movsd  0x58(%rbx),%xmm2
mulsd  -0x40(%r13),%xmm2
addsd  %xmm0,%xmm1
mulsd  0x50(%r15),%xmm2
movsd  0x60(%rbx),%xmm0
mulsd  -0x38(%r13),%xmm0
mulsd  0x58(%r15),%xmm0
addsd  %xmm1,%xmm2
movsd  0x68(%rbx),%xmm1
mulsd  -0x30(%r13),%xmm1
addsd  %xmm2,%xmm0
mulsd  0x60(%r15),%xmm1
movsd  0x70(%rbx),%xmm2
mulsd  -0x28(%r13),%xmm2
mulsd  0x68(%r15),%xmm2
addsd  %xmm0,%xmm1
movsd  0x78(%rbx),%xmm0
mulsd  -0x20(%r13),%xmm0
addsd  %xmm1,%xmm2
mulsd  0x70(%r15),%xmm0
movsd  0x80(%rbx),%xmm1
mulsd  -0x18(%r13),%xmm1
mulsd  0x78(%r15),%xmm1
addsd  %xmm2,%xmm0
movsd  0x88(%rbx),%xmm2
mulsd  -0x10(%r13),%xmm2
addsd  %xmm0,%xmm1
mulsd  0x80(%r15),%xmm2
movsd  0x90(%rbx),%xmm3
mulsd  -0x8(%r13),%xmm3
mulsd  0x88(%r15),%xmm3
addsd  %xmm1,%xmm2
movsd  0x98(%rbx),%xmm0
mulsd  0x0(%r13),%xmm0
addsd  %xmm2,%xmm3
mulsd  0x90(%r15),%xmm0
addsd  %xmm3,%xmm0
call   1060 <log@plt>
movsd  0xa01(%rip),%xmm2
xorpd  %xmm4,%xmm4
xorps  %xmm1,%xmm1
cvtsi2sdl (%r14,%rbp,4),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x80(%rsp,%r12,8)
add    $0x1,%r12
add    $0xa0,%r13
add    $0x98,%r15
cmp    $0x4,%r12
jne    1440 <categorizeGTRGAMMAPROTMULT+0x300>
movapd 0x80(%rsp),%xmm1
movapd 0x90(%rsp),%xmm0
movapd %xmm1,%xmm2
unpckhpd %xmm1,%xmm2
addsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
movapd %xmm0,%xmm3
unpckhpd %xmm0,%xmm3
addsd  %xmm2,%xmm3
mulsd  0x998(%rip),%xmm3
movsd  0x998(%rip),%xmm2
divsd  %xmm3,%xmm2
movapd 0x9a4(%rip),%xmm5
mulpd  %xmm5,%xmm1
unpcklpd %xmm2,%xmm2
mulpd  %xmm2,%xmm1
movsd  0x980(%rip),%xmm6
ucomisd %xmm1,%xmm6
movapd %xmm1,0x80(%rsp)
mov    $0x0,%eax
mov    (%rsp),%ecx
cmova  %eax,%ecx
movapd %xmm1,%xmm3
minsd  %xmm6,%xmm3
unpckhpd %xmm1,%xmm1
ucomisd %xmm1,%xmm3
mulpd  %xmm5,%xmm0
mov    $0x1,%eax
cmova  %eax,%ecx
mulpd  %xmm2,%xmm0
minsd  %xmm3,%xmm1
ucomisd %xmm0,%xmm1
movapd %xmm0,0x90(%rsp)
mov    $0x2,%eax
cmova  %eax,%ecx
movapd %xmm0,%xmm2
minsd  %xmm1,%xmm2
unpckhpd %xmm0,%xmm0
ucomisd %xmm0,%xmm2
mov    $0x3,%eax
cmova  %eax,%ecx
mov    0x48(%rsp),%rax
mov    %ecx,(%rax,%rbp,4)
add    $0x1,%rbp
mov    0x8(%rsp),%r13
add    $0x280,%r13
cmp    0x18(%rsp),%rbp
jne    13f0 <categorizeGTRGAMMAPROTMULT+0x2b0>
mov    0x10(%rsp),%rdi
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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