#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ exp(double);
int /*<<< orphan*/ free(double *);
int log(double);
scalar_t__ malloc(int);
void evaluateGTRGAMMAPROTMULT_VECTOR(int *ex2, int *modelptr, double *x2,
                                     double *v, double *EIGN,
                                     double *gammaRates, double *tipVector,
                                     double *pz, char *tipX1, int lower, int n,
                                     int numberOfModels, int multiBranch) {
  double z, lz = 0.0, term, ki;
  int i, j;
  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 < n; i++) {
    model = modelptr[i];
    left = &(tipVector[460 * model + 20 * tipX1[i]]);
    diagptable = &diagptable_start[76 * model];
    term = 0;
    for (j = 0; j < 4; j++) {
      right = &(x2[80 * i + 20 * j]);
      term += left[0] * right[0];
      term += left[1] * right[1] * *diagptable++;
      term += left[2] * right[2] * *diagptable++;
      term += left[3] * right[3] * *diagptable++;
      term += left[4] * right[4] * *diagptable++;
      term += left[5] * right[5] * *diagptable++;
      term += left[6] * right[6] * *diagptable++;
      term += left[7] * right[7] * *diagptable++;
      term += left[8] * right[8] * *diagptable++;
      term += left[9] * right[9] * *diagptable++;
      term += left[10] * right[10] * *diagptable++;
      term += left[11] * right[11] * *diagptable++;
      term += left[12] * right[12] * *diagptable++;
      term += left[13] * right[13] * *diagptable++;
      term += left[14] * right[14] * *diagptable++;
      term += left[15] * right[15] * *diagptable++;
      term += left[16] * right[16] * *diagptable++;
      term += left[17] * right[17] * *diagptable++;
      term += left[18] * right[18] * *diagptable++;
      term += left[19] * right[19] * *diagptable++;
    }
    term =
        log(0.25 * term) +
        ex2[i] *
            log((
                1.0 /
                115792089237316195423570985008687907853269984665640564039457584007913129639936.0));
    v[i] = term;
  }
  free(diagptable_start);
}


x86 assembly (clang, -O0):
<evaluateGTRGAMMAPROTMULT_VECTOR>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    0x40(%rbp),%eax
mov    0x38(%rbp),%eax
mov    0x30(%rbp),%eax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
cmpl   $0x0,0x40(%rbp)
jne    11cd <evaluateGTRGAMMAPROTMULT_VECTOR+0x8d>
mov    0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe5e(%rip),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    11ba <evaluateGTRGAMMAPROTMULT_VECTOR+0x7a>
movsd  0xe4b(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1060 <log@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x40(%rbp)
movslq 0x38(%rbp),%rax
imul   $0x260,%rax,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x68(%rbp)
mov    %rax,-0x60(%rbp)
movl   $0x0,-0x7c(%rbp)
mov    -0x7c(%rbp),%eax
cmp    0x38(%rbp),%eax
jge    1308 <evaluateGTRGAMMAPROTMULT_VECTOR+0x1c8>
cmpl   $0x0,0x40(%rbp)
je     1249 <evaluateGTRGAMMAPROTMULT_VECTOR+0x109>
mov    0x18(%rbp),%rax
movslq -0x7c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xde2(%rip),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    1236 <evaluateGTRGAMMAPROTMULT_VECTOR+0xf6>
movsd  0xdcf(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1060 <log@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x68(%rbp),%rax
imul   $0x4c,-0x7c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
movl   $0x0,-0x54(%rbp)
cmpl   $0x4,-0x54(%rbp)
jge    12f5 <evaluateGTRGAMMAPROTMULT_VECTOR+0x1b5>
mov    -0x30(%rbp),%rax
mov    -0x7c(%rbp),%ecx
shl    $0x2,%ecx
add    -0x54(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x58(%rbp)
cmpl   $0x13,-0x58(%rbp)
jge    12e2 <evaluateGTRGAMMAPROTMULT_VECTOR+0x1a2>
mov    -0x28(%rbp),%rax
imul   $0x13,-0x7c(%rbp),%ecx
add    -0x58(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x50(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
call   1050 <exp@plt>
cvtsi2sd %eax,%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  %xmm0,(%rax)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1291 <evaluateGTRGAMMAPROTMULT_VECTOR+0x151>
jmp    12e7 <evaluateGTRGAMMAPROTMULT_VECTOR+0x1a7>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1266 <evaluateGTRGAMMAPROTMULT_VECTOR+0x126>
jmp    12fa <evaluateGTRGAMMAPROTMULT_VECTOR+0x1ba>
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    11ee <evaluateGTRGAMMAPROTMULT_VECTOR+0xae>
mov    0x28(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%eax
cmp    0x30(%rbp),%eax
jge    1866 <evaluateGTRGAMMAPROTMULT_VECTOR+0x726>
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x7c(%rbp)
mov    0x10(%rbp),%rax
imul   $0x1cc,-0x7c(%rbp),%ecx
mov    0x20(%rbp),%rdx
movslq -0x54(%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,-0x70(%rbp)
mov    -0x68(%rbp),%rax
imul   $0x4c,-0x7c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x58(%rbp)
cmpl   $0x4,-0x58(%rbp)
jge    17fa <evaluateGTRGAMMAPROTMULT_VECTOR+0x6ba>
mov    -0x18(%rbp),%rax
imul   $0x50,-0x54(%rbp),%ecx
imul   $0x14,-0x58(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x78(%rbp)
mov    -0x70(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x78(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x10(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x18(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x20(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x20(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x28(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x30(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x30(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x38(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x38(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x40(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x40(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x48(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x48(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x50(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x50(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x58(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x58(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x60(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x60(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x68(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x68(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x70(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x70(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x78(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x78(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x80(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x80(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x88(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x88(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x90(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x90(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x70(%rbp),%rax
movsd  0x98(%rax),%xmm0
mov    -0x78(%rbp),%rax
mulsd  0x98(%rax),%xmm0
mov    -0x60(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x60(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1377 <evaluateGTRGAMMAPROTMULT_VECTOR+0x237>
movsd  0x80e(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
call   1060 <log@plt>
mov    %eax,-0x80(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x84(%rbp)
movsd  0x7e0(%rip),%xmm0
call   1060 <log@plt>
mov    -0x84(%rbp),%ecx
mov    %eax,%edx
mov    -0x80(%rbp),%eax
imul   %edx,%ecx
add    %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    130e <evaluateGTRGAMMAPROTMULT_VECTOR+0x1ce>
mov    -0x68(%rbp),%rdi
call   1030 <free@plt>
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<evaluateGTRGAMMAPROTMULT_VECTOR>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %r9,%r13
mov    %r8,0x38(%rsp)
mov    %rcx,0x80(%rsp)
mov    %rdx,%rbx
mov    %rsi,0x78(%rsp)
mov    %rdi,0x70(%rsp)
mov    0xf0(%rsp),%ecx
mov    0xe8(%rsp),%ebp
mov    0xe0(%rsp),%eax
mov    %eax,0x54(%rsp)
mov    0xd8(%rsp),%eax
mov    %eax,0x50(%rsp)
mov    0xd0(%rsp),%rax
mov    %rax,0x68(%rsp)
mov    0xc8(%rsp),%rax
mov    %rax,0x58(%rsp)
mov    0xc0(%rsp),%rax
mov    %rax,0x60(%rsp)
xorpd  %xmm0,%xmm0
mov    %ecx,0x30(%rsp)
test   %ecx,%ecx
jne    11e2 <evaluateGTRGAMMAPROTMULT_VECTOR+0xa2>
movsd  0xe33(%rip),%xmm0
mov    0x58(%rsp),%rax
maxsd  (%rax),%xmm0
call   1060 <log@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x10(%rsp)
imul   $0x260,%ebp,%edi
call   1040 <malloc@plt>
movsd  0x10(%rsp),%xmm1
mov    %rax,0x18(%rsp)
test   %ebp,%ebp
jle    13e1 <evaluateGTRGAMMAPROTMULT_VECTOR+0x2a1>
mov    %ebp,%eax
mov    %rax,0x20(%rsp)
xor    %ebp,%ebp
xor    %r12d,%r12d
mov    %r13,0x28(%rsp)
nopw   0x0(%rax,%rax,1)
mov    %ebp,%r14d
and    $0x3fffffff,%r14d
shl    $0x5,%r14
add    0x18(%rsp),%r14
cmpl   $0x0,0x30(%rsp)
je     1259 <evaluateGTRGAMMAPROTMULT_VECTOR+0x119>
movsd  0xdbe(%rip),%xmm0
mov    0x58(%rsp),%rax
maxsd  (%rax,%r12,8),%xmm0
call   1060 <log@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
mov    %rbp,0x48(%rsp)
mov    %ebp,%eax
mov    0x38(%rsp),%rcx
lea    (%rcx,%rax,8),%r15
mov    %r12,0x40(%rsp)
lea    0x0(,%r12,4),%r12d
movsd  0x0(%r13,%r12,8),%xmm0
movsd  %xmm0,0x8(%rsp)
xor    %ebp,%ebp
movsd  %xmm1,0x10(%rsp)
nopl   0x0(%rax,%rax,1)
movsd  (%r15,%rbp,1),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x10(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%r14,%rbp,1)
add    $0x8,%rbp
cmp    $0x98,%rbp
jne    1290 <evaluateGTRGAMMAPROTMULT_VECTOR+0x150>
shl    $0x3,%r12
mov    %r12,%rax
or     $0x8,%rax
movsd  0x0(%r13,%rax,1),%xmm0
movsd  %xmm0,0x8(%rsp)
add    %rbp,%r14
mov    $0x13,%r13d
mov    %r15,%rbp
nopl   0x0(%rax)
movsd  0x0(%rbp),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x10(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%r14)
add    $0x8,%r14
add    $0x8,%rbp
add    $0xffffffffffffffff,%r13
jne    12f0 <evaluateGTRGAMMAPROTMULT_VECTOR+0x1b0>
mov    %r12,%rax
or     $0x10,%rax
mov    0x28(%rsp),%r13
movsd  0x0(%r13,%rax,1),%xmm0
movsd  %xmm0,0x8(%rsp)
xor    %ebp,%ebp
nop
movsd  (%r15,%rbp,8),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x10(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%r14)
add    $0x8,%r14
add    $0x1,%rbp
cmp    $0x13,%rbp
jne    1340 <evaluateGTRGAMMAPROTMULT_VECTOR+0x200>
or     $0x18,%r12
movsd  0x0(%r13,%r12,1),%xmm0
movsd  %xmm0,0x8(%rsp)
xor    %ebp,%ebp
nopl   0x0(%rax,%rax,1)
movsd  (%r15,%rbp,8),%xmm0
mulsd  0x8(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  0x10(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%r14)
add    $0x8,%r14
add    $0x1,%rbp
cmp    $0x13,%rbp
jne    1390 <evaluateGTRGAMMAPROTMULT_VECTOR+0x250>
mov    0x40(%rsp),%r12
add    $0x1,%r12
mov    0x48(%rsp),%rbp
add    $0x13,%ebp
cmp    0x20(%rsp),%r12
jne    1220 <evaluateGTRGAMMAPROTMULT_VECTOR+0xe0>
mov    0x54(%rsp),%eax
mov    0x50(%rsp),%ecx
cmp    %eax,%ecx
jge    1721 <evaluateGTRGAMMAPROTMULT_VECTOR+0x5e1>
movslq %ecx,%r14
movslq %eax,%r15
lea    (%r14,%r14,4),%rax
shl    $0x7,%rax
add    %rax,%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x78(%rsp),%rax
movslq (%rax,%r14,4),%rax
imul   $0x1cc,%rax,%rcx
mov    0x68(%rsp),%rdx
movsbq (%rdx,%r14,1),%rdx
lea    (%rdx,%rdx,4),%rdx
lea    (%rcx,%rdx,4),%rdx
imul   $0x260,%rax,%rax
add    0x18(%rsp),%rax
mov    0x60(%rsp),%rsi
movsd  (%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  0x8(%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x10(%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x48(%rsp)
movsd  0x18(%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x40(%rsp)
movsd  0x20(%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x38(%rsp)
movsd  0x28(%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x30(%rsp)
movsd  0x30(%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x28(%rsp)
movsd  0x38(%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x20(%rsp)
movsd  0x40(%rsi,%rdx,8),%xmm14
movsd  0x48(%rsi,%rdx,8),%xmm15
movsd  0x50(%rsi,%rdx,8),%xmm4
movsd  0x58(%rsi,%rdx,8),%xmm5
movsd  0x60(%rsi,%rdx,8),%xmm6
movsd  0x68(%rsi,%rdx,8),%xmm7
movsd  0x70(%rsi,%rdx,8),%xmm1
movsd  0x78(%rsi,%rdx,8),%xmm8
mov    $0x98,%ecx
movsd  0x80(%rsi,%rdx,8),%xmm9
movsd  0x88(%rsi,%rdx,8),%xmm10
movsd  0x90(%rsi,%rdx,8),%xmm11
movsd  0x98(%rsi,%rdx,8),%xmm12
xorpd  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
movsd  -0x98(%rbx,%rcx,1),%xmm13
ff ff ff
mulsd  0x10(%rsp),%xmm13
addsd  %xmm0,%xmm13
movsd  -0x90(%rbx,%rcx,1),%xmm0
ff ff
mulsd  0x8(%rsp),%xmm0
mulsd  (%rax),%xmm0
addsd  %xmm13,%xmm0
movsd  -0x88(%rbx,%rcx,1),%xmm2
ff ff
mulsd  0x48(%rsp),%xmm2
mulsd  0x8(%rax),%xmm2
addsd  %xmm0,%xmm2
movsd  -0x80(%rbx,%rcx,1),%xmm0
mulsd  0x40(%rsp),%xmm0
mulsd  0x10(%rax),%xmm0
movsd  -0x78(%rbx,%rcx,1),%xmm3
mulsd  0x38(%rsp),%xmm3
mulsd  0x18(%rax),%xmm3
addsd  %xmm2,%xmm0
addsd  %xmm0,%xmm3
movsd  -0x70(%rbx,%rcx,1),%xmm0
mulsd  0x30(%rsp),%xmm0
mulsd  0x20(%rax),%xmm0
addsd  %xmm3,%xmm0
movsd  -0x68(%rbx,%rcx,1),%xmm2
mulsd  0x28(%rsp),%xmm2
mulsd  0x28(%rax),%xmm2
addsd  %xmm0,%xmm2
movsd  -0x60(%rbx,%rcx,1),%xmm0
mulsd  0x20(%rsp),%xmm0
mulsd  0x30(%rax),%xmm0
movsd  -0x58(%rbx,%rcx,1),%xmm3
mulsd  %xmm14,%xmm3
mulsd  0x38(%rax),%xmm3
addsd  %xmm2,%xmm0
addsd  %xmm0,%xmm3
movsd  -0x50(%rbx,%rcx,1),%xmm0
mulsd  %xmm15,%xmm0
mulsd  0x40(%rax),%xmm0
addsd  %xmm3,%xmm0
movsd  -0x48(%rbx,%rcx,1),%xmm2
mulsd  %xmm4,%xmm2
mulsd  0x48(%rax),%xmm2
addsd  %xmm0,%xmm2
movsd  -0x40(%rbx,%rcx,1),%xmm0
mulsd  %xmm5,%xmm0
mulsd  0x50(%rax),%xmm0
movsd  -0x38(%rbx,%rcx,1),%xmm3
mulsd  %xmm6,%xmm3
mulsd  0x58(%rax),%xmm3
addsd  %xmm2,%xmm0
addsd  %xmm0,%xmm3
movsd  -0x30(%rbx,%rcx,1),%xmm0
mulsd  %xmm7,%xmm0
mulsd  0x60(%rax),%xmm0
addsd  %xmm3,%xmm0
movsd  -0x28(%rbx,%rcx,1),%xmm2
mulsd  %xmm1,%xmm2
mulsd  0x68(%rax),%xmm2
addsd  %xmm0,%xmm2
movsd  -0x20(%rbx,%rcx,1),%xmm0
mulsd  %xmm8,%xmm0
mulsd  0x70(%rax),%xmm0
movsd  -0x18(%rbx,%rcx,1),%xmm3
mulsd  %xmm9,%xmm3
mulsd  0x78(%rax),%xmm3
addsd  %xmm2,%xmm0
addsd  %xmm0,%xmm3
movsd  -0x10(%rbx,%rcx,1),%xmm0
mulsd  %xmm10,%xmm0
mulsd  0x80(%rax),%xmm0
addsd  %xmm3,%xmm0
movsd  -0x8(%rbx,%rcx,1),%xmm2
mulsd  %xmm11,%xmm2
mulsd  0x88(%rax),%xmm2
addsd  %xmm0,%xmm2
movsd  (%rbx,%rcx,1),%xmm0
mulsd  %xmm12,%xmm0
mulsd  0x90(%rax),%xmm0
addsd  %xmm2,%xmm0
add    $0xa0,%rcx
add    $0x98,%rax
cmp    $0x318,%rcx
jne    1510 <evaluateGTRGAMMAPROTMULT_VECTOR+0x3d0>
mulsd  0x934(%rip),%xmm0
call   1060 <log@plt>
mov    %eax,%r12d
mov    0x70(%rsp),%rax
mov    (%rax,%r14,4),%ebp
movsd  0x923(%rip),%xmm0
call   1060 <log@plt>
imul   %ebp,%eax
add    %r12d,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    0x80(%rsp),%rax
movsd  %xmm0,(%rax,%r14,8)
add    $0x1,%r14
add    $0x280,%rbx
cmp    %r15,%r14
jne    1410 <evaluateGTRGAMMAPROTMULT_VECTOR+0x2d0>
mov    0x18(%rsp),%rdi
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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