#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_LG4(double *gammaRates, double *EIGN[4], double *sumtable,
                          int upper, int *wrptr, volatile double *ext_dlnLdlz,
                          volatile double *ext_d2lnLdlz2, double lz,
                          double *weights) {
  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[i][l - 1] * ki * lz);
      diagptable[i * 128 + l * 4 + 1] = EIGN[i][l - 1] * ki;
      diagptable[i * 128 + l * 4 + 2] = EIGN[i][l - 1] * EIGN[i][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++) {
      double a1 = 0.0, a2 = 0.0, a3 = 0.0;
      a1 += sum[j * 20];
      for (l = 1; l < 20; l++) {
        a1 += (tmp = diagptable[j * 128 + l * 4] * sum[j * 20 + l]);
        a2 += tmp * diagptable[j * 128 + l * 4 + 1];
        a3 += tmp * diagptable[j * 128 + l * 4 + 2];
      }
      inv_Li += weights[j] * a1;
      dlnLidlz += weights[j] * a2;
      d2lnLidlz2 += weights[j] * a3;
    }
    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_LG4>:
push   %rbp
mov    %rsp,%rbp
sub    $0x10b0,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %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    12f5 <coreGTRGAMMAPROT_LG4+0x1d5>
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    12dc <coreGTRGAMMAPROT_LG4+0x1bc>
mov    -0x10(%rbp),%rax
movslq -0x1044(%rbp),%rcx
mov    (%rax,%rcx,8),%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
movslq -0x1044(%rbp),%rcx
mov    (%rax,%rcx,8),%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
movslq -0x1044(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x104c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1044(%rbp),%rcx
mov    (%rax,%rcx,8),%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    11b6 <coreGTRGAMMAPROT_LG4+0x96>
jmp    12e1 <coreGTRGAMMAPROT_LG4+0x1c1>
mov    -0x1044(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1044(%rbp)
jmp    116f <coreGTRGAMMAPROT_LG4+0x4f>
movl   $0x0,-0x1044(%rbp)
mov    -0x1044(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    163d <coreGTRGAMMAPROT_LG4+0x51d>
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    1546 <coreGTRGAMMAPROT_LG4+0x426>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1098(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x10a0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x10a8(%rbp)
ff
mov    -0x40(%rbp),%rax
imul   $0x14,-0x1048(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x1098(%rbp),%xmm0
ff
movsd  %xmm0,-0x1098(%rbp)
ff
movl   $0x1,-0x104c(%rbp)
cmpl   $0x14,-0x104c(%rbp)
jge    14a2 <coreGTRGAMMAPROT_LG4+0x382>
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  -0x1098(%rbp),%xmm0
ff
movsd  %xmm0,-0x1098(%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  -0x10a0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x10a0(%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  -0x10a8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x10a8(%rbp)
ff
mov    -0x104c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x104c(%rbp)
jmp    13ad <coreGTRGAMMAPROT_LG4+0x28d>
mov    0x18(%rbp),%rax
movslq -0x1048(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x1098(%rbp),%xmm2
ff
movsd  -0x1080(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1080(%rbp)
ff
mov    0x18(%rbp),%rax
movslq -0x1048(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x10a0(%rbp),%xmm2
ff
movsd  -0x1088(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1088(%rbp)
ff
mov    0x18(%rbp),%rax
movslq -0x1048(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x10a8(%rbp),%xmm2
ff
movsd  -0x1090(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1090(%rbp)
ff
mov    -0x1048(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1048(%rbp)
jmp    1352 <coreGTRGAMMAPROT_LG4+0x232>
movsd  -0x1080(%rbp),%xmm0
ff
call   1030 <FABS@plt>
movaps %xmm0,%xmm1
movsd  0xaa2(%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    12ff <coreGTRGAMMAPROT_LG4+0x1df>
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    $0x10b0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<coreGTRGAMMAPROT_LG4>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x1058,%rsp
movapd %xmm0,%xmm1
mov    %r9,%r15
mov    %r8,0x48(%rsp)
mov    %ecx,%r12d
mov    %rdx,%r13
mov    %rsi,%rbp
mov    %rdi,%r14
movsd  (%rdi),%xmm2
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mov    (%rsi),%rax
movapd %xmm2,0x10(%rsp)
unpcklpd %xmm0,%xmm2
movapd %xmm2,(%rsp)
xor    %ebx,%ebx
movsd  %xmm1,0x20(%rsp)
nop
movsd  (%rax,%rbx,1),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0x70(%rsp,%rbx,4)
mov    0x0(%rbp),%rax
movsd  (%rax,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
movsd  0x20(%rsp),%xmm1
mulpd  (%rsp),%xmm0
movupd %xmm0,0x78(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    1170 <coreGTRGAMMAPROT_LG4+0x50>
movsd  0x8(%r14),%xmm2
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mov    0x8(%rbp),%rax
movapd %xmm2,0x10(%rsp)
unpcklpd %xmm0,%xmm2
movapd %xmm2,(%rsp)
xor    %ebx,%ebx
movsd  (%rax,%rbx,1),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0x470(%rsp,%rbx,4)
mov    0x8(%rbp),%rax
movsd  (%rax,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
movsd  0x20(%rsp),%xmm1
mulpd  (%rsp),%xmm0
movupd %xmm0,0x478(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    11e0 <coreGTRGAMMAPROT_LG4+0xc0>
movsd  0x10(%r14),%xmm2
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mov    0x10(%rbp),%rax
movapd %xmm2,0x10(%rsp)
unpcklpd %xmm0,%xmm2
movapd %xmm2,(%rsp)
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rax,%rbx,1),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0x870(%rsp,%rbx,4)
mov    0x10(%rbp),%rax
movsd  (%rax,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
movsd  0x20(%rsp),%xmm1
mulpd  (%rsp),%xmm0
movupd %xmm0,0x878(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    1260 <coreGTRGAMMAPROT_LG4+0x140>
movsd  0x18(%r14),%xmm2
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mov    0x18(%rbp),%rax
movapd %xmm2,0x10(%rsp)
unpcklpd %xmm0,%xmm2
movapd %xmm2,(%rsp)
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rax,%rbx,1),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0xc70(%rsp,%rbx,4)
mov    0x18(%rbp),%rax
movsd  (%rax,%rbx,1),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
movsd  0x20(%rsp),%xmm1
mulpd  (%rsp),%xmm0
movupd %xmm0,0xc78(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    $0x98,%rbx
jne    12e0 <coreGTRGAMMAPROT_LG4+0x1c0>
test   %r12d,%r12d
jle    16be <coreGTRGAMMAPROT_LG4+0x59e>
mov    %r15,0x38(%rsp)
mov    %r12d,%eax
mov    %rax,0x40(%rsp)
lea    0x10(%r13),%rbp
lea    0xb0(%r13),%rbx
lea    0x150(%r13),%r15
lea    0x1f0(%r13),%r12
xorpd  %xmm0,%xmm0
xor    %r14d,%r14d
xorpd  %xmm10,%xmm10
jmp    144f <coreGTRGAMMAPROT_LG4+0x32f>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mulsd  %xmm8,%xmm4
addsd  %xmm1,%xmm4
unpckhpd %xmm5,%xmm5
mulsd  %xmm8,%xmm5
addsd  %xmm2,%xmm5
movsd  0x18(%rdx),%xmm1
mulsd  %xmm1,%xmm0
addsd  %xmm4,%xmm0
movapd %xmm1,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm11,%xmm2
movapd %xmm2,0x10(%rsp)
unpckhpd %xmm3,%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm5,%xmm3
movapd %xmm3,(%rsp)
call   1030 <FABS@plt>
xorpd  %xmm10,%xmm10
movsd  0xc27(%rip),%xmm1
divsd  %xmm0,%xmm1
movapd 0x10(%rsp),%xmm3
mulsd  %xmm1,%xmm3
mulsd  (%rsp),%xmm1
mov    0x48(%rsp),%rax
xorps  %xmm0,%xmm0
cvtsi2sdl (%rax,%r14,4),%xmm0
movapd %xmm3,%xmm2
mulsd  %xmm3,%xmm2
subsd  %xmm2,%xmm1
movapd %xmm3,%xmm2
unpcklpd %xmm0,%xmm0
unpcklpd %xmm1,%xmm2
mulpd  %xmm0,%xmm2
movapd %xmm2,%xmm0
movapd 0x20(%rsp),%xmm1
addpd  %xmm1,%xmm0
add    $0x1,%r14
add    $0x280,%rbp
add    $0x280,%rbx
add    $0x280,%r15
add    $0x280,%r12
cmp    0x40(%rsp),%r14
je     16b3 <coreGTRGAMMAPROT_LG4+0x593>
movapd %xmm0,0x20(%rsp)
mov    %r14,%rax
shl    $0x4,%rax
lea    (%rax,%rax,4),%rax
movsd  0x0(%r13,%rax,8),%xmm0
addsd  %xmm10,%xmm0
xor    %ecx,%ecx
xorpd  %xmm1,%xmm1
mov    0x1098(%rsp),%rdx
nopw   0x0(%rax,%rax,1)
movsd  0x70(%rsp,%rcx,4),%xmm3
mulsd  -0x8(%rbp,%rcx,1),%xmm3
addsd  %xmm3,%xmm0
movupd 0x78(%rsp,%rcx,4),%xmm2
unpcklpd %xmm3,%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm1,%xmm3
cmp    $0x90,%rcx
je     14e0 <coreGTRGAMMAPROT_LG4+0x3c0>
movsd  0x90(%rsp,%rcx,4),%xmm1
mulsd  0x0(%rbp,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
movupd 0x98(%rsp,%rcx,4),%xmm2
unpcklpd %xmm1,%xmm1
mulpd  %xmm2,%xmm1
addpd  %xmm3,%xmm1
add    $0x10,%rcx
jmp    1480 <coreGTRGAMMAPROT_LG4+0x360>
nopl   0x0(%rax)
movsd  (%rdx),%xmm5
movapd %xmm5,%xmm4
mulsd  %xmm3,%xmm4
addsd  %xmm10,%xmm4
movsd  0xa0(%r13,%rax,8),%xmm1
addsd  %xmm10,%xmm1
xorpd  %xmm6,%xmm6
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movsd  0x470(%rsp,%rcx,4),%xmm2
mulsd  -0x8(%rbx,%rcx,1),%xmm2
addsd  %xmm2,%xmm1
movupd 0x478(%rsp,%rcx,4),%xmm7
unpcklpd %xmm2,%xmm2
mulpd  %xmm7,%xmm2
addpd  %xmm6,%xmm2
cmp    $0x90,%rcx
je     1570 <coreGTRGAMMAPROT_LG4+0x450>
movsd  0x490(%rsp,%rcx,4),%xmm6
mulsd  (%rbx,%rcx,1),%xmm6
addsd  %xmm6,%xmm1
movupd 0x498(%rsp,%rcx,4),%xmm7
unpcklpd %xmm6,%xmm6
mulpd  %xmm7,%xmm6
addpd  %xmm2,%xmm6
add    $0x10,%rcx
jmp    1510 <coreGTRGAMMAPROT_LG4+0x3f0>
xchg   %ax,%ax
mulsd  %xmm5,%xmm0
addsd  %xmm10,%xmm0
unpckhpd %xmm3,%xmm3
mulsd  %xmm5,%xmm3
addsd  %xmm10,%xmm3
movsd  0x8(%rdx),%xmm8
movapd %xmm8,%xmm9
mulsd  %xmm2,%xmm9
addsd  %xmm4,%xmm9
movsd  0x140(%r13,%rax,8),%xmm4
addsd  %xmm10,%xmm4
xorpd  %xmm7,%xmm7
xor    %ecx,%ecx
movsd  0x870(%rsp,%rcx,4),%xmm5
mulsd  -0x8(%r15,%rcx,1),%xmm5
addsd  %xmm5,%xmm4
movupd 0x878(%rsp,%rcx,4),%xmm6
unpcklpd %xmm5,%xmm5
mulpd  %xmm6,%xmm5
addpd  %xmm7,%xmm5
cmp    $0x90,%rcx
je     1610 <coreGTRGAMMAPROT_LG4+0x4f0>
movsd  0x890(%rsp,%rcx,4),%xmm7
mulsd  (%r15,%rcx,1),%xmm7
addsd  %xmm7,%xmm4
movupd 0x898(%rsp,%rcx,4),%xmm6
unpcklpd %xmm7,%xmm7
mulpd  %xmm6,%xmm7
addpd  %xmm5,%xmm7
add    $0x10,%rcx
jmp    15b0 <coreGTRGAMMAPROT_LG4+0x490>
mulsd  %xmm8,%xmm1
addsd  %xmm0,%xmm1
unpckhpd %xmm2,%xmm2
mulsd  %xmm8,%xmm2
addsd  %xmm3,%xmm2
movsd  0x10(%rdx),%xmm8
movapd %xmm8,%xmm11
mulsd  %xmm5,%xmm11
addsd  %xmm9,%xmm11
movsd  0x1e0(%r13,%rax,8),%xmm0
addsd  %xmm10,%xmm0
xorpd  %xmm6,%xmm6
xor    %eax,%eax
movsd  0xc70(%rsp,%rax,4),%xmm3
mulsd  -0x8(%r12,%rax,1),%xmm3
addsd  %xmm3,%xmm0
movupd 0xc78(%rsp,%rax,4),%xmm7
unpcklpd %xmm3,%xmm3
mulpd  %xmm7,%xmm3
addpd  %xmm6,%xmm3
cmp    $0x90,%rax
je     1380 <coreGTRGAMMAPROT_LG4+0x260>
movsd  0xc90(%rsp,%rax,4),%xmm6
mulsd  (%r12,%rax,1),%xmm6
addsd  %xmm6,%xmm0
movupd 0xc98(%rsp,%rax,4),%xmm7
unpcklpd %xmm6,%xmm6
mulpd  %xmm7,%xmm6
addpd  %xmm3,%xmm6
add    $0x10,%rax
jmp    1650 <coreGTRGAMMAPROT_LG4+0x530>
addpd  %xmm1,%xmm2
mov    0x38(%rsp),%r15
jmp    16c2 <coreGTRGAMMAPROT_LG4+0x5a2>
xorpd  %xmm2,%xmm2
movlpd %xmm2,(%r15)
mov    0x1090(%rsp),%rax
movhpd %xmm2,(%rax)
add    $0x1058,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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