#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 */
double FABS(double);
double exp(double);
void coreGTRGAMMAPROT(double *gammaRates, double *EIGN, double *sumtable,
                      int upper, int *wrptr, volatile double *ext_dlnLdlz,
                      volatile double *ext_d2lnLdlz2, double lz) {
  double *sum, diagptable[512];
  int i, j, l;
  double dlnLdlz = 0;
  double d2lnLdlz2 = 0;
  double ki, kisqr;
  double tmp;
  double inv_Li, dlnLidlz, d2lnLidlz2;
  for (i = 0; i < 4; i++) {
    ki = gammaRates[i];
    kisqr = ki * ki;
    for (l = 1; l < 20; l++) {
      diagptable[i * 128 + l * 4] = exp(EIGN[l - 1] * ki * lz);
      diagptable[i * 128 + l * 4 + 1] = EIGN[l - 1] * ki;
      diagptable[i * 128 + l * 4 + 2] = EIGN[l - 1] * EIGN[l - 1] * kisqr;
    }
  }
  for (i = 0; i < upper; i++) {
    sum = &sumtable[i * 80];
    inv_Li = 0.0;
    dlnLidlz = 0.0;
    d2lnLidlz2 = 0.0;
    for (j = 0; j < 4; j++) {
      inv_Li += sum[j * 20];
      for (l = 1; l < 20; l++) {
        inv_Li += (tmp = diagptable[j * 128 + l * 4] * sum[j * 20 + l]);
        dlnLidlz += tmp * diagptable[j * 128 + l * 4 + 1];
        d2lnLidlz2 += tmp * diagptable[j * 128 + l * 4 + 2];
      }
    }
    inv_Li = 1.0 / FABS(inv_Li);
    dlnLidlz *= inv_Li;
    d2lnLidlz2 *= inv_Li;
    dlnLdlz += wrptr[i] * dlnLidlz;
    d2lnLdlz2 += wrptr[i] * (d2lnLidlz2 - dlnLidlz * dlnLidlz);
  }
  *ext_dlnLdlz = dlnLdlz;
  *ext_d2lnLdlz2 = d2lnLdlz2;
}


x86 assembly (clang, -O0):
<coreGTRGAMMAPROT>:
push   %rbp
mov    %rsp,%rbp
sub    $0x1090,%rsp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1058(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1060(%rbp)
ff
movl   $0x0,-0x1044(%rbp)
cmpl   $0x4,-0x1044(%rbp)
jge    12c5 <coreGTRGAMMAPROT+0x1a5>
mov    -0x8(%rbp),%rax
movslq -0x1044(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x1068(%rbp)
ff
movsd  -0x1068(%rbp),%xmm0
ff
mulsd  -0x1068(%rbp),%xmm0
ff
movsd  %xmm0,-0x1070(%rbp)
ff
movl   $0x1,-0x104c(%rbp)
cmpl   $0x14,-0x104c(%rbp)
jge    12ac <coreGTRGAMMAPROT+0x18c>
mov    -0x10(%rbp),%rax
mov    -0x104c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x1068(%rbp),%xmm0
ff
mulsd  -0x38(%rbp),%xmm0
call   1040 <exp@plt>
mov    -0x1044(%rbp),%eax
shl    $0x7,%eax
mov    -0x104c(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
cltq
movsd  %xmm0,-0x1040(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%rax
mov    -0x104c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x1068(%rbp),%xmm0
ff
mov    -0x1044(%rbp),%eax
shl    $0x7,%eax
mov    -0x104c(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x1,%eax
cltq
movsd  %xmm0,-0x1040(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%rax
mov    -0x104c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x104c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x1070(%rbp),%xmm0
ff
mov    -0x1044(%rbp),%eax
shl    $0x7,%eax
mov    -0x104c(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x2,%eax
cltq
movsd  %xmm0,-0x1040(%rbp,%rax,8)
ff ff
mov    -0x104c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x104c(%rbp)
jmp    11b2 <coreGTRGAMMAPROT+0x92>
jmp    12b1 <coreGTRGAMMAPROT+0x191>
mov    -0x1044(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1044(%rbp)
jmp    116b <coreGTRGAMMAPROT+0x4b>
movl   $0x0,-0x1044(%rbp)
mov    -0x1044(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1561 <coreGTRGAMMAPROT+0x441>
mov    -0x18(%rbp),%rax
imul   $0x50,-0x1044(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1080(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1088(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1090(%rbp)
ff
movl   $0x0,-0x1048(%rbp)
cmpl   $0x4,-0x1048(%rbp)
jge    146a <coreGTRGAMMAPROT+0x34a>
mov    -0x40(%rbp),%rax
imul   $0x14,-0x1048(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x1080(%rbp),%xmm0
ff
movsd  %xmm0,-0x1080(%rbp)
ff
movl   $0x1,-0x104c(%rbp)
cmpl   $0x14,-0x104c(%rbp)
jge    1451 <coreGTRGAMMAPROT+0x331>
mov    -0x1048(%rbp),%eax
shl    $0x7,%eax
mov    -0x104c(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
cltq
movsd  -0x1040(%rbp,%rax,8),%xmm0
ff ff
mov    -0x40(%rbp),%rax
imul   $0x14,-0x1048(%rbp),%ecx
add    -0x104c(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x1078(%rbp)
ff
addsd  -0x1080(%rbp),%xmm0
ff
movsd  %xmm0,-0x1080(%rbp)
ff
movsd  -0x1078(%rbp),%xmm0
ff
mov    -0x1048(%rbp),%eax
shl    $0x7,%eax
mov    -0x104c(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x1,%eax
cltq
movsd  -0x1040(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x1088(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1088(%rbp)
ff
movsd  -0x1078(%rbp),%xmm0
ff
mov    -0x1048(%rbp),%eax
shl    $0x7,%eax
mov    -0x104c(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x2,%eax
cltq
movsd  -0x1040(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x1090(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1090(%rbp)
ff
mov    -0x104c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x104c(%rbp)
jmp    135c <coreGTRGAMMAPROT+0x23c>
jmp    1456 <coreGTRGAMMAPROT+0x336>
mov    -0x1048(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1048(%rbp)
jmp    1322 <coreGTRGAMMAPROT+0x202>
movsd  -0x1080(%rbp),%xmm0
ff
call   1030 <FABS@plt>
movaps %xmm0,%xmm1
movsd  0xb7e(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x1080(%rbp)
ff
movsd  -0x1080(%rbp),%xmm0
ff
mulsd  -0x1088(%rbp),%xmm0
ff
movsd  %xmm0,-0x1088(%rbp)
ff
movsd  -0x1080(%rbp),%xmm0
ff
mulsd  -0x1090(%rbp),%xmm0
ff
movsd  %xmm0,-0x1090(%rbp)
ff
mov    -0x28(%rbp),%rax
movslq -0x1044(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x1088(%rbp),%xmm2
ff
movsd  -0x1058(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1058(%rbp)
ff
mov    -0x28(%rbp),%rax
movslq -0x1044(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x1090(%rbp),%xmm1
ff
movsd  -0x1088(%rbp),%xmm2
ff
movsd  -0x1088(%rbp),%xmm3
ff
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  -0x1060(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1060(%rbp)
ff
mov    -0x1044(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1044(%rbp)
jmp    12cf <coreGTRGAMMAPROT+0x1af>
movsd  -0x1058(%rbp),%xmm0
ff
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x1060(%rbp),%xmm0
ff
mov    0x10(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x1090,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<coreGTRGAMMAPROT>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x1048,%rsp
movapd %xmm0,%xmm1
mov    %r9,%r15
mov    %r8,0x38(%rsp)
mov    %ecx,%r12d
mov    %rdx,%r14
mov    %rsi,%rbp
mov    %rdi,%r13
movsd  (%rdi),%xmm2
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movapd %xmm2,0x10(%rsp)
unpcklpd %xmm0,%xmm2
movapd %xmm2,(%rsp)
xor    %ebx,%ebx
movsd  %xmm1,0x20(%rsp)
nopl   0x0(%rax)
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0x60(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
movsd  0x20(%rsp),%xmm1
mulpd  (%rsp),%xmm0
movupd %xmm0,0x68(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    1170 <coreGTRGAMMAPROT+0x50>
movsd  0x8(%r13),%xmm2
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movapd %xmm2,0x10(%rsp)
unpcklpd %xmm0,%xmm2
movapd %xmm2,(%rsp)
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0x460(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
movsd  0x20(%rsp),%xmm1
mulpd  (%rsp),%xmm0
movupd %xmm0,0x468(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    11e0 <coreGTRGAMMAPROT+0xc0>
movsd  0x10(%r13),%xmm2
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movapd %xmm2,0x10(%rsp)
unpcklpd %xmm0,%xmm2
movapd %xmm2,(%rsp)
xor    %ebx,%ebx
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0x860(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
movsd  0x20(%rsp),%xmm1
mulpd  (%rsp),%xmm0
movupd %xmm0,0x868(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    1250 <coreGTRGAMMAPROT+0x130>
movsd  0x18(%r13),%xmm2
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movapd %xmm2,0x10(%rsp)
unpcklpd %xmm0,%xmm2
movapd %xmm2,(%rsp)
xor    %ebx,%ebx
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0xc60(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
movsd  0x20(%rsp),%xmm1
mulpd  (%rsp),%xmm0
movupd %xmm0,0xc68(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    12c0 <coreGTRGAMMAPROT+0x1a0>
test   %r12d,%r12d
jle    15be <coreGTRGAMMAPROT+0x49e>
mov    %r15,0x30(%rsp)
mov    %r12d,%eax
mov    %rax,(%rsp)
lea    0x10(%r14),%r13
lea    0xb0(%r14),%rbx
lea    0x150(%r14),%rbp
lea    0x1f0(%r14),%r15
xorpd  %xmm1,%xmm1
xor    %r12d,%r12d
jmp    13db <coreGTRGAMMAPROT+0x2bb>
nopl   0x0(%rax,%rax,1)
movapd %xmm3,0x10(%rsp)
call   1030 <FABS@plt>
movsd  0xc9d(%rip),%xmm1
divsd  %xmm0,%xmm1
movapd 0x10(%rsp),%xmm3
movapd %xmm3,%xmm0
mulsd  %xmm1,%xmm0
unpckhpd %xmm3,%xmm3
mulsd  %xmm1,%xmm3
mov    0x38(%rsp),%rax
xorps  %xmm1,%xmm1
cvtsi2sdl (%rax,%r12,4),%xmm1
movapd %xmm0,%xmm2
mulsd  %xmm0,%xmm2
subsd  %xmm2,%xmm3
unpcklpd %xmm1,%xmm1
unpcklpd %xmm3,%xmm0
mulpd  %xmm1,%xmm0
movapd %xmm0,%xmm1
movapd 0x20(%rsp),%xmm2
addpd  %xmm2,%xmm1
add    $0x1,%r12
add    $0x280,%r13
add    $0x280,%rbx
add    $0x280,%rbp
add    $0x280,%r15
cmp    (%rsp),%r12
je     15b3 <coreGTRGAMMAPROT+0x493>
movapd %xmm1,0x20(%rsp)
mov    %r12,%rax
shl    $0x4,%rax
lea    (%rax,%rax,4),%rax
movsd  (%r14,%rax,8),%xmm0
addsd  0xc0e(%rip),%xmm0
xor    %ecx,%ecx
xorpd  %xmm2,%xmm2
movsd  0x60(%rsp,%rcx,4),%xmm1
mulsd  -0x8(%r13,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0x68(%rsp,%rcx,4),%xmm3
unpcklpd %xmm1,%xmm1
mulpd  %xmm3,%xmm1
addpd  %xmm2,%xmm1
cmp    $0x90,%rcx
je     1460 <coreGTRGAMMAPROT+0x340>
movsd  0x80(%rsp,%rcx,4),%xmm2
mulsd  0x0(%r13,%rcx,1),%xmm2
addsd  %xmm2,%xmm0
movupd 0x88(%rsp,%rcx,4),%xmm3
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
add    $0x10,%rcx
jmp    1400 <coreGTRGAMMAPROT+0x2e0>
nopl   0x0(%rax,%rax,1)
addsd  0xa0(%r14,%rax,8),%xmm0
xor    %ecx,%ecx
nopl   0x0(%rax)
movsd  0x460(%rsp,%rcx,4),%xmm2
mulsd  -0x8(%rbx,%rcx,1),%xmm2
addsd  %xmm2,%xmm0
movupd 0x468(%rsp,%rcx,4),%xmm3
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
cmp    $0x90,%rcx
je     14d0 <coreGTRGAMMAPROT+0x3b0>
movsd  0x480(%rsp,%rcx,4),%xmm1
mulsd  (%rbx,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0x488(%rsp,%rcx,4),%xmm3
unpcklpd %xmm1,%xmm1
mulpd  %xmm3,%xmm1
addpd  %xmm2,%xmm1
add    $0x10,%rcx
jmp    1470 <coreGTRGAMMAPROT+0x350>
xchg   %ax,%ax
addsd  0x140(%r14,%rax,8),%xmm0
xor    %ecx,%ecx
nopl   0x0(%rax)
movsd  0x860(%rsp,%rcx,4),%xmm1
mulsd  -0x8(%rbp,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0x868(%rsp,%rcx,4),%xmm3
unpcklpd %xmm1,%xmm1
mulpd  %xmm3,%xmm1
addpd  %xmm2,%xmm1
cmp    $0x90,%rcx
je     1540 <coreGTRGAMMAPROT+0x420>
movsd  0x880(%rsp,%rcx,4),%xmm2
mulsd  0x0(%rbp,%rcx,1),%xmm2
addsd  %xmm2,%xmm0
movupd 0x888(%rsp,%rcx,4),%xmm3
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
add    $0x10,%rcx
jmp    14e0 <coreGTRGAMMAPROT+0x3c0>
nop
addsd  0x1e0(%r14,%rax,8),%xmm0
xor    %eax,%eax
nopl   0x0(%rax)
movsd  0xc60(%rsp,%rax,4),%xmm3
mulsd  -0x8(%r15,%rax,1),%xmm3
addsd  %xmm3,%xmm0
movupd 0xc68(%rsp,%rax,4),%xmm2
unpcklpd %xmm3,%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm1,%xmm3
cmp    $0x90,%rax
je     1350 <coreGTRGAMMAPROT+0x230>
movsd  0xc80(%rsp,%rax,4),%xmm1
mulsd  (%r15,%rax,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0xc88(%rsp,%rax,4),%xmm2
unpcklpd %xmm1,%xmm1
mulpd  %xmm2,%xmm1
addpd  %xmm3,%xmm1
add    $0x10,%rax
jmp    1550 <coreGTRGAMMAPROT+0x430>
addpd  %xmm2,%xmm0
mov    0x30(%rsp),%r15
jmp    15c2 <coreGTRGAMMAPROT+0x4a2>
xorpd  %xmm0,%xmm0
movlpd %xmm0,(%r15)
mov    0x1080(%rsp),%rax
movhpd %xmm0,(%rax)
add    $0x1048,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @coreGTRGAMMAPROT(double* %0, double* %1, double* %2, i32 %3, i32* %4, double* %5, double* %6, double %7) {
%9 = alloca double*
%10 = alloca double*
%11 = alloca double*
%12 = alloca i32
%13 = alloca i32*
%14 = alloca double*
%15 = alloca double*
%16 = alloca double
%17 = alloca double*
%18 = alloca [512 x double]
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca double
store double* %0, double** %9
store double* %1, double** %10
store double* %2, double** %11
store i32 %3, i32* %12
store i32* %4, i32** %13
store double* %5, double** %14
store double* %6, double** %15
store double %7, double* %16
store double 0.000000e+00, double* %22
store double 0.000000e+00, double* %23
store i32 0, i32* %19
br label %30
30: ; preds = %107, %8
%31 = load i32, i32* %19
%32 = icmp slt i32 %31, 4
br i1 %32, label %33, label %110
33: ; preds = %30
%34 = load double*, double** %9
%35 = load i32, i32* %19
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
%38 = load double, double* %37
store double %38, double* %24
%39 = load double, double* %24
%40 = load double, double* %24
%41 = fmul double %39, %40
store double %41, double* %25
store i32 1, i32* %21
br label %42
42: ; preds = %103, %33
%43 = load i32, i32* %21
%44 = icmp slt i32 %43, 20
br i1 %44, label %45, label %106
45: ; preds = %42
%46 = load double*, double** %10
%47 = load i32, i32* %21
%48 = sub nsw i32 %47, 1
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %46, i64 %49
%51 = load double, double* %50
%52 = load double, double* %24
%53 = fmul double %51, %52
%54 = load double, double* %16
%55 = fmul double %53, %54
%56 = call double @exp(double %55)
%57 = load i32, i32* %19
%58 = mul nsw i32 %57, 128
%59 = load i32, i32* %21
%60 = mul nsw i32 %59, 4
%61 = add nsw i32 %58, %60
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds [512 x double], [512 x double]* %18, i64 0, i64 %62
store double %56, double* %63
%64 = load double*, double** %10
%65 = load i32, i32* %21
%66 = sub nsw i32 %65, 1
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %64, i64 %67
%69 = load double, double* %68
%70 = load double, double* %24
%71 = fmul double %69, %70
%72 = load i32, i32* %19
%73 = mul nsw i32 %72, 128
%74 = load i32, i32* %21
%75 = mul nsw i32 %74, 4
%76 = add nsw i32 %73, %75
%77 = add nsw i32 %76, 1
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds [512 x double], [512 x double]* %18, i64 0, i64 %78
store double %71, double* %79
%80 = load double*, double** %10
%81 = load i32, i32* %21
%82 = sub nsw i32 %81, 1
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds double, double* %80, i64 %83
%85 = load double, double* %84
%86 = load double*, double** %10
%87 = load i32, i32* %21
%88 = sub nsw i32 %87, 1
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %86, i64 %89
%91 = load double, double* %90
%92 = fmul double %85, %91
%93 = load double, double* %25
%94 = fmul double %92, %93
%95 = load i32, i32* %19
%96 = mul nsw i32 %95, 128
%97 = load i32, i32* %21
%98 = mul nsw i32 %97, 4
%99 = add nsw i32 %96, %98
%100 = add nsw i32 %99, 2
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds [512 x double], [512 x double]* %18, i64 0, i64 %101
store double %94, double* %102
br label %103
103: ; preds = %45
%104 = load i32, i32* %21
%105 = add nsw i32 %104, 1
store i32 %105, i32* %21
br label %42
106: ; preds = %42
br label %107
107: ; preds = %106
%108 = load i32, i32* %19
%109 = add nsw i32 %108, 1
store i32 %109, i32* %19
br label %30
110: ; preds = %30
store i32 0, i32* %19
br label %111
111: ; preds = %219, %110
%112 = load i32, i32* %19
%113 = load i32, i32* %12
%114 = icmp slt i32 %112, %113
br i1 %114, label %115, label %222
115: ; preds = %111
%116 = load double*, double** %11
%117 = load i32, i32* %19
%118 = mul nsw i32 %117, 80
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds double, double* %116, i64 %119
store double* %120, double** %17
store double 0.000000e+00, double* %27
store double 0.000000e+00, double* %28
store double 0.000000e+00, double* %29
store i32 0, i32* %20
br label %121
121: ; preds = %184, %115
%122 = load i32, i32* %20
%123 = icmp slt i32 %122, 4
br i1 %123, label %124, label %187
124: ; preds = %121
%125 = load double*, double** %17
%126 = load i32, i32* %20
%127 = mul nsw i32 %126, 20
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds double, double* %125, i64 %128
%130 = load double, double* %129
%131 = load double, double* %27
%132 = fadd double %131, %130
store double %132, double* %27
store i32 1, i32* %21
br label %133
133: ; preds = %180, %124
%134 = load i32, i32* %21
%135 = icmp slt i32 %134, 20
br i1 %135, label %136, label %183
136: ; preds = %133
%137 = load i32, i32* %20
%138 = mul nsw i32 %137, 128
%139 = load i32, i32* %21
%140 = mul nsw i32 %139, 4
%141 = add nsw i32 %138, %140
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds [512 x double], [512 x double]* %18, i64 0, i64 %142
%144 = load double, double* %143
%145 = load double*, double** %17
%146 = load i32, i32* %20
%147 = mul nsw i32 %146, 20
%148 = load i32, i32* %21
%149 = add nsw i32 %147, %148
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds double, double* %145, i64 %150
%152 = load double, double* %151
%153 = fmul double %144, %152
store double %153, double* %26
%154 = load double, double* %27
%155 = fadd double %154, %153
store double %155, double* %27
%156 = load double, double* %26
%157 = load i32, i32* %20
%158 = mul nsw i32 %157, 128
%159 = load i32, i32* %21
%160 = mul nsw i32 %159, 4
%161 = add nsw i32 %158, %160
%162 = add nsw i32 %161, 1
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds [512 x double], [512 x double]* %18, i64 0, i64 %163
%165 = load double, double* %164
%166 = load double, double* %28
%167 = call double @llvm.fmuladd.f64(double %156, double %165, double %166)
store double %167, double* %28
%168 = load double, double* %26
%169 = load i32, i32* %20
%170 = mul nsw i32 %169, 128
%171 = load i32, i32* %21
%172 = mul nsw i32 %171, 4
%173 = add nsw i32 %170, %172
%174 = add nsw i32 %173, 2
%175 = sext i32 %174 to i64
%176 = getelementptr inbounds [512 x double], [512 x double]* %18, i64 0, i64 %175
%177 = load double, double* %176
%178 = load double, double* %29
%179 = call double @llvm.fmuladd.f64(double %168, double %177, double %178)
store double %179, double* %29
br label %180
180: ; preds = %136
%181 = load i32, i32* %21
%182 = add nsw i32 %181, 1
store i32 %182, i32* %21
br label %133
183: ; preds = %133
br label %184
184: ; preds = %183
%185 = load i32, i32* %20
%186 = add nsw i32 %185, 1
store i32 %186, i32* %20
br label %121
187: ; preds = %121
%188 = load double, double* %27
%189 = call double @FABS(double %188)
%190 = fdiv double 1.000000e+00, %189
store double %190, double* %27
%191 = load double, double* %27
%192 = load double, double* %28
%193 = fmul double %192, %191
store double %193, double* %28
%194 = load double, double* %27
%195 = load double, double* %29
%196 = fmul double %195, %194
store double %196, double* %29
%197 = load i32*, i32** %13
%198 = load i32, i32* %19
%199 = sext i32 %198 to i64
%200 = getelementptr inbounds i32, i32* %197, i64 %199
%201 = load i32, i32* %200
%202 = sitofp i32 %201 to double
%203 = load double, double* %28
%204 = load double, double* %22
%205 = call double @llvm.fmuladd.f64(double %202, double %203, double %204)
store double %205, double* %22
%206 = load i32*, i32** %13
%207 = load i32, i32* %19
%208 = sext i32 %207 to i64
%209 = getelementptr inbounds i32, i32* %206, i64 %208
%210 = load i32, i32* %209
%211 = sitofp i32 %210 to double
%212 = load double, double* %29
%213 = load double, double* %28
%214 = load double, double* %28
%215 = fneg double %213
%216 = call double @llvm.fmuladd.f64(double %215, double %214, double %212)
%217 = load double, double* %23
%218 = call double @llvm.fmuladd.f64(double %211, double %216, double %217)
store double %218, double* %23
br label %219
219: ; preds = %187
%220 = load i32, i32* %19
%221 = add nsw i32 %220, 1
store i32 %221, i32* %19
br label %111
222: ; preds = %111
%223 = load double, double* %22
%224 = load double*, double** %14
store volatile double %223, double* %224
%225 = load double, double* %23
%226 = load double*, double** %15
store volatile double %225, double* %226
ret void
}

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