#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 */
typedef double d2lnLdlz2;

/* Variables and functions */
double FABS(double);
double exp(double);
double coreGammaAsc(double *gammaRates, double *EIGN, double *sumtable,
                    int upper, volatile double *ext_dlnLdlz,
                    volatile double *ext_d2lnLdlz2, double lz,
                    const int numStates) {
  double diagptable[1024], lh = 0.0, dlnLdlz = 0.0, d2lnLdlz2 = 0.0, ki, kisqr;
  int i, j, l;
  const int gammaStates = 4 * numStates;
  for (i = 0; i < 4; i++) {
    ki = gammaRates[i];
    kisqr = ki * ki;
    for (l = 1; l < numStates; l++) {
      diagptable[i * gammaStates + l * 4] = exp(EIGN[l - 1] * ki * lz);
      diagptable[i * gammaStates + l * 4 + 1] = EIGN[l - 1] * ki;
      diagptable[i * gammaStates + l * 4 + 2] =
          EIGN[l - 1] * EIGN[l - 1] * kisqr;
    }
  }
  for (i = 0; i < upper; i++) {
    double *sum = &sumtable[i * gammaStates], tmp, inv_Li = 0.0, dlnLidlz = 0.0,
           d2lnLidlz2 = 0.0;
    for (j = 0; j < 4; j++) {
      inv_Li += sum[j * numStates];
      for (l = 1; l < numStates; l++) {
        inv_Li += (tmp = diagptable[j * gammaStates + l * 4] *
                         sum[j * numStates + l]);
        dlnLidlz += tmp * diagptable[j * gammaStates + l * 4 + 1];
        d2lnLidlz2 += tmp * diagptable[j * gammaStates + l * 4 + 2];
      }
    }
    inv_Li = 0.25 * FABS(inv_Li);
    dlnLidlz *= 0.25;
    d2lnLidlz2 *= 0.25;
    lh += inv_Li;
    dlnLdlz += dlnLidlz;
    d2lnLdlz2 += d2lnLidlz2;
  }
  *ext_dlnLdlz = (dlnLdlz / (lh - 1.0));
  *ext_d2lnLdlz2 = (((lh - 1.0) * (d2lnLdlz2) - (dlnLdlz * dlnLdlz)) /
                    ((lh - 1.0) * (lh - 1.0)));
  return lh;
}


x86 assembly (clang, -O0):
<coreGammaAsc>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20a0,%rsp
mov    0x10(%rbp),%eax
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,-0x2048(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2050(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2058(%rbp)
ff
mov    0x10(%rbp),%eax
shl    $0x2,%eax
mov    %eax,-0x2078(%rbp)
movl   $0x0,-0x206c(%rbp)
cmpl   $0x4,-0x206c(%rbp)
jge    12e9 <coreGammaAsc+0x1c9>
mov    -0x8(%rbp),%rax
movslq -0x206c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x2060(%rbp)
ff
movsd  -0x2060(%rbp),%xmm0
ff
mulsd  -0x2060(%rbp),%xmm0
ff
movsd  %xmm0,-0x2068(%rbp)
ff
movl   $0x1,-0x2074(%rbp)
mov    -0x2074(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    12d0 <coreGammaAsc+0x1b0>
mov    -0x10(%rbp),%rax
mov    -0x2074(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x2060(%rbp),%xmm0
ff
mulsd  -0x38(%rbp),%xmm0
call   1040 <exp@plt>
mov    -0x206c(%rbp),%eax
imul   -0x2078(%rbp),%eax
mov    -0x2074(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
cltq
movsd  %xmm0,-0x2040(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%rax
mov    -0x2074(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x2060(%rbp),%xmm0
ff
mov    -0x206c(%rbp),%eax
imul   -0x2078(%rbp),%eax
mov    -0x2074(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x1,%eax
cltq
movsd  %xmm0,-0x2040(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%rax
mov    -0x2074(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2074(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x2068(%rbp),%xmm0
ff
mov    -0x206c(%rbp),%eax
imul   -0x2078(%rbp),%eax
mov    -0x2074(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x2,%eax
cltq
movsd  %xmm0,-0x2040(%rbp,%rax,8)
ff ff
mov    -0x2074(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2074(%rbp)
jmp    11c8 <coreGammaAsc+0xa8>
jmp    12d5 <coreGammaAsc+0x1b5>
mov    -0x206c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x206c(%rbp)
jmp    1181 <coreGammaAsc+0x61>
movl   $0x0,-0x206c(%rbp)
mov    -0x206c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1561 <coreGammaAsc+0x441>
mov    -0x18(%rbp),%rax
mov    -0x206c(%rbp),%ecx
imul   -0x2078(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x2080(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2090(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2098(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20a0(%rbp)
ff
movl   $0x0,-0x2070(%rbp)
cmpl   $0x4,-0x2070(%rbp)
jge    14b1 <coreGammaAsc+0x391>
mov    -0x2080(%rbp),%rax
mov    -0x2070(%rbp),%ecx
imul   0x10(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x2090(%rbp),%xmm0
ff
movsd  %xmm0,-0x2090(%rbp)
ff
movl   $0x1,-0x2074(%rbp)
mov    -0x2074(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    1498 <coreGammaAsc+0x378>
mov    -0x2070(%rbp),%eax
imul   -0x2078(%rbp),%eax
mov    -0x2074(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
cltq
movsd  -0x2040(%rbp,%rax,8),%xmm0
ff ff
mov    -0x2080(%rbp),%rax
mov    -0x2070(%rbp),%ecx
imul   0x10(%rbp),%ecx
add    -0x2074(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x2088(%rbp)
ff
addsd  -0x2090(%rbp),%xmm0
ff
movsd  %xmm0,-0x2090(%rbp)
ff
movsd  -0x2088(%rbp),%xmm0
ff
mov    -0x2070(%rbp),%eax
imul   -0x2078(%rbp),%eax
mov    -0x2074(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x1,%eax
cltq
movsd  -0x2040(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x2098(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x2098(%rbp)
ff
movsd  -0x2088(%rbp),%xmm0
ff
mov    -0x2070(%rbp),%eax
imul   -0x2078(%rbp),%eax
mov    -0x2074(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
add    $0x2,%eax
cltq
movsd  -0x2040(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x20a0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20a0(%rbp)
ff
mov    -0x2074(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2074(%rbp)
jmp    138f <coreGammaAsc+0x26f>
jmp    149d <coreGammaAsc+0x37d>
mov    -0x2070(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2070(%rbp)
jmp    134f <coreGammaAsc+0x22f>
movsd  -0x2090(%rbp),%xmm0
ff
call   1030 <FABS@plt>
movaps %xmm0,%xmm1
movsd  0xb3f(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x2090(%rbp)
ff
movsd  0xb2b(%rip),%xmm0
mulsd  -0x2098(%rbp),%xmm0
ff
movsd  %xmm0,-0x2098(%rbp)
ff
movsd  0xb13(%rip),%xmm0
mulsd  -0x20a0(%rbp),%xmm0
ff
movsd  %xmm0,-0x20a0(%rbp)
ff
movsd  -0x2090(%rbp),%xmm0
ff
addsd  -0x2048(%rbp),%xmm0
ff
movsd  %xmm0,-0x2048(%rbp)
ff
movsd  -0x2098(%rbp),%xmm0
ff
addsd  -0x2050(%rbp),%xmm0
ff
movsd  %xmm0,-0x2050(%rbp)
ff
movsd  -0x20a0(%rbp),%xmm0
ff
addsd  -0x2058(%rbp),%xmm0
ff
movsd  %xmm0,-0x2058(%rbp)
ff
mov    -0x206c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x206c(%rbp)
jmp    12f3 <coreGammaAsc+0x1d3>
movsd  -0x2050(%rbp),%xmm0
ff
movsd  -0x2048(%rbp),%xmm1
ff
movsd  0xa87(%rip),%xmm2
subsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x2048(%rbp),%xmm0
ff
movsd  0xa67(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  -0x2058(%rbp),%xmm2
ff
movsd  -0x2050(%rbp),%xmm1
ff
mulsd  -0x2050(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x2048(%rbp),%xmm1
ff
movsd  0xa1c(%rip),%xmm2
subsd  %xmm2,%xmm1
movsd  -0x2048(%rbp),%xmm2
ff
movsd  0xa08(%rip),%xmm3
subsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x2048(%rbp),%xmm0
ff
add    $0x20a0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<coreGammaAsc>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x20c8,%rsp
mov    %r9,0x98(%rsp)
mov    %r8,0x90(%rsp)
mov    %ecx,%r12d
mov    %rdx,0x80(%rsp)
mov    0x2100(%rsp),%ecx
lea    0x0(,%rcx,4),%eax
cltq
mov    %rax,0x38(%rsp)
mov    %ecx,%eax
mov    %rax,0x8(%rsp)
cmp    $0x1,%ecx
jle    1433 <coreGammaAsc+0x313>
movapd %xmm0,%xmm1
mov    %rsi,%r13
mov    %rdi,%r14
movsd  (%rdi),%xmm2
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movapd %xmm2,0x10(%rsp)
unpcklpd %xmm0,%xmm2
movapd %xmm2,0x20(%rsp)
mov    0x8(%rsp),%rax
lea    -0x4(,%rax,4),%rbx
ff
xor    %ebp,%ebp
movsd  %xmm1,0x50(%rsp)
nopl   0x0(%rax)
movsd  0x0(%r13,%rbp,2),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,0xe0(%rsp,%rbp,8)
movsd  0x0(%r13,%rbp,2),%xmm0
lea    0x5(%rbp),%eax
cltq
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
unpcklpd %xmm1,%xmm0
movsd  0x50(%rsp),%xmm1
mulpd  0x20(%rsp),%xmm0
movupd %xmm0,0xc0(%rsp,%rax,8)
add    $0x4,%rbp
cmp    %rbp,%rbx
jne    11b0 <coreGammaAsc+0x90>
mov    0x2100(%rsp),%eax
cmp    $0x2,%eax
jl     1433 <coreGammaAsc+0x313>
mov    %r14,0x40(%rsp)
movsd  0x8(%r14),%xmm0
movsd  %xmm0,0x10(%rsp)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x20(%rsp)
mov    0x8(%rsp),%rax
lea    -0x8(,%rax,8),%r14
ff
mov    0x38(%rsp),%rax
lea    (%rsp,%rax,8),%r15
add    $0xe0,%r15
lea    0x5(%rax),%ebx
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0x0(%r13,%rbp,1),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,(%r15,%rbp,4)
movsd  0x0(%r13,%rbp,1),%xmm0
movsd  0x10(%rsp),%xmm1
mulsd  %xmm0,%xmm1
movslq %ebx,%rbx
movsd  %xmm1,0xc0(%rsp,%rbx,8)
movsd  0x50(%rsp),%xmm1
mulsd  %xmm0,%xmm0
mulsd  0x20(%rsp),%xmm0
lea    0x1(%rbx),%eax
cltq
movsd  %xmm0,0xc0(%rsp,%rax,8)
add    $0x8,%rbp
add    $0x4,%ebx
cmp    %rbp,%r14
jne    1260 <coreGammaAsc+0x140>
mov    0x2100(%rsp),%eax
cmp    $0x2,%eax
mov    0x40(%rsp),%rax
jl     1433 <coreGammaAsc+0x313>
movsd  0x10(%rax),%xmm0
movsd  %xmm0,0x10(%rsp)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x20(%rsp)
mov    0x38(%rsp),%rax
lea    (%rax,%rax,1),%rbx
shl    $0x4,%rax
lea    0xc0(%rsp),%rcx
lea    (%rax,%rcx,1),%r15
add    $0x20,%r15
add    $0x6,%ebx
xor    %ebp,%ebp
nop
movsd  0x0(%r13,%rbp,1),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,(%r15,%rbp,4)
movsd  0x0(%r13,%rbp,1),%xmm0
movsd  0x10(%rsp),%xmm1
mulsd  %xmm0,%xmm1
lea    -0x1(%rbx),%eax
cltq
movsd  %xmm1,0xc0(%rsp,%rax,8)
movsd  0x50(%rsp),%xmm1
mulsd  %xmm0,%xmm0
mulsd  0x20(%rsp),%xmm0
movslq %ebx,%rbx
movsd  %xmm0,0xc0(%rsp,%rbx,8)
add    $0x8,%rbp
add    $0x4,%ebx
cmp    %rbp,%r14
jne    1310 <coreGammaAsc+0x1f0>
mov    0x2100(%rsp),%eax
cmp    $0x2,%eax
mov    0x40(%rsp),%rax
jl     1433 <coreGammaAsc+0x313>
movsd  0x18(%rax),%xmm0
movsd  %xmm0,0x10(%rsp)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x20(%rsp)
mov    0x38(%rsp),%rax
lea    (%rax,%rax,2),%rax
mov    0x8(%rsp),%rcx
lea    -0x8(,%rcx,8),%r14
ff
lea    (%rsp,%rax,8),%r15
add    $0xe0,%r15
lea    0x6(%rax),%ebx
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  0x0(%r13,%rbp,1),%xmm0
mulsd  0x10(%rsp),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,(%r15,%rbp,4)
movsd  0x0(%r13,%rbp,1),%xmm0
movsd  0x10(%rsp),%xmm1
mulsd  %xmm0,%xmm1
lea    -0x1(%rbx),%eax
cltq
movsd  %xmm1,0xc0(%rsp,%rax,8)
movsd  0x50(%rsp),%xmm1
mulsd  %xmm0,%xmm0
mulsd  0x20(%rsp),%xmm0
movslq %ebx,%rbx
movsd  %xmm0,0xc0(%rsp,%rbx,8)
add    $0x8,%rbp
add    $0x4,%ebx
cmp    %rbp,%r14
jne    13d0 <coreGammaAsc+0x2b0>
test   %r12d,%r12d
jle    17c8 <coreGammaAsc+0x6a8>
mov    0x2100(%rsp),%eax
movslq %eax,%rbx
mov    %r12d,%eax
mov    %rax,0x40(%rsp)
lea    (%rbx,%rbx,1),%rax
mov    %rax,0x70(%rsp)
lea    (%rbx,%rbx,2),%r8
mov    0x38(%rsp),%rcx
lea    (%rcx,%rcx,1),%rsi
lea    (%rcx,%rcx,2),%rbp
lea    0x0(,%rcx,8),%rdi
mov    %rdi,0xb8(%rsp)
mov    0x8(%rsp),%rdi
lea    0x0(,%rdi,4),%r14
shl    $0x3,%rdi
mov    0x80(%rsp),%r13
lea    0x0(,%rbx,8),%r15
add    %r13,%r15
lea    0xc0(%rsp,%rcx,8),%rax
lea    0x5(%rcx),%edx
mov    %edx,0x4c(%rsp)
mov    %rbx,0x78(%rsp)
shl    $0x4,%rbx
add    %r13,%rbx
shl    $0x4,%rcx
add    %rsp,%rcx
add    $0xc0,%rcx
mov    %rcx,0xa0(%rsp)
add    $0x6,%esi
mov    %rsi,0xa8(%rsp)
mov    %r8,0x68(%rsp)
lea    0x0(,%r8,8),%r12
add    %r13,%r12
lea    0xc0(%rsp,%rbp,8),%rcx
mov    %rcx,0x88(%rsp)
lea    0x6(%rbp),%ecx
mov    %ecx,0x48(%rsp)
xorpd  %xmm4,%xmm4
xorpd  %xmm5,%xmm5
xor    %ebp,%ebp
mov    %rdi,0x8(%rsp)
mov    %rax,0xb0(%rsp)
cs nopw 0x0(%rax,%rax,1)
mov    %rbp,%rcx
imul   0x38(%rsp),%rcx
mov    0x80(%rsp),%rdx
lea    (%rdx,%rcx,8),%r8
movsd  (%rdx,%rcx,8),%xmm0
addsd  0xace(%rip),%xmm0
mov    0x2100(%rsp),%ecx
cmp    $0x2,%ecx
movapd %xmm4,0x50(%rsp)
movsd  %xmm5,0x10(%rsp)
jl     1730 <coreGammaAsc+0x610>
xorpd  %xmm3,%xmm3
mov    $0x4,%ecx
nopw   0x0(%rax,%rax,1)
movapd %xmm3,%xmm2
movsd  0xc0(%rsp,%rcx,8),%xmm1
mulsd  0x0(%r13,%rcx,2),%xmm1
addsd  %xmm1,%xmm0
lea    0x1(%rcx),%edx
movslq %edx,%rdx
movupd 0xc0(%rsp,%rdx,8),%xmm3
unpcklpd %xmm1,%xmm1
mulpd  %xmm3,%xmm1
movapd %xmm1,%xmm3
addpd  %xmm2,%xmm3
add    $0x4,%rcx
cmp    %rcx,%r14
jne    1570 <coreGammaAsc+0x450>
addpd  %xmm2,%xmm1
mov    0x78(%rsp),%rcx
addsd  (%r8,%rcx,8),%xmm0
mov    0x2100(%rsp),%ecx
cmp    $0x2,%ecx
jl     173f <coreGammaAsc+0x61f>
mov    $0x8,%ecx
mov    0x4c(%rsp),%edx
nopl   0x0(%rax,%rax,1)
movapd %xmm1,%xmm3
movsd  (%rax,%rcx,4),%xmm1
mulsd  (%r15,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
movslq %edx,%rdx
lea    0x1(%rdx),%esi
movslq %esi,%rsi
unpcklpd %xmm1,%xmm1
movsd  0xc0(%rsp,%rdx,8),%xmm2
movhpd 0xc0(%rsp,%rsi,8),%xmm2
mulpd  %xmm1,%xmm2
movapd %xmm2,%xmm1
addpd  %xmm3,%xmm1
add    $0x8,%rcx
add    $0x4,%edx
cmp    %rcx,%rdi
jne    15e0 <coreGammaAsc+0x4c0>
addpd  %xmm3,%xmm2
mov    0x70(%rsp),%rcx
addsd  (%r8,%rcx,8),%xmm0
mov    0x2100(%rsp),%ecx
cmp    $0x2,%ecx
jl     174e <coreGammaAsc+0x62e>
mov    $0x8,%ecx
mov    0xa8(%rsp),%rdx
mov    0xa0(%rsp),%rax
xchg   %ax,%ax
movapd %xmm2,%xmm1
movsd  (%rax,%rcx,4),%xmm2
mulsd  (%rbx,%rcx,1),%xmm2
addsd  %xmm2,%xmm0
lea    -0x1(%rdx),%esi
movslq %esi,%rsi
movslq %edx,%rdx
unpcklpd %xmm2,%xmm2
movsd  0xc0(%rsp,%rsi,8),%xmm3
movhpd 0xc0(%rsp,%rdx,8),%xmm3
mulpd  %xmm2,%xmm3
movapd %xmm3,%xmm2
addpd  %xmm1,%xmm2
add    $0x8,%rcx
add    $0x4,%edx
cmp    %rcx,%rdi
jne    1660 <coreGammaAsc+0x540>
addpd  %xmm1,%xmm3
mov    0x68(%rsp),%rcx
addsd  (%r8,%rcx,8),%xmm0
mov    0x2100(%rsp),%eax
cmp    $0x2,%eax
jl     17c0 <coreGammaAsc+0x6a0>
mov    $0x8,%eax
mov    0x48(%rsp),%ecx
mov    0x88(%rsp),%rsi
nopl   0x0(%rax)
movapd %xmm3,%xmm1
movsd  (%rsi,%rax,4),%xmm3
mulsd  (%r12,%rax,1),%xmm3
addsd  %xmm3,%xmm0
lea    -0x1(%rcx),%edx
movslq %edx,%rdx
movslq %ecx,%rcx
unpcklpd %xmm3,%xmm3
movsd  0xc0(%rsp,%rdx,8),%xmm2
movhpd 0xc0(%rsp,%rcx,8),%xmm2
mulpd  %xmm3,%xmm2
movapd %xmm2,%xmm3
addpd  %xmm1,%xmm3
add    $0x8,%rax
add    $0x4,%ecx
cmp    %rax,%rdi
jne    16e0 <coreGammaAsc+0x5c0>
addpd  %xmm1,%xmm2
jmp    1759 <coreGammaAsc+0x639>
mov    0x78(%rsp),%rcx
addsd  (%r8,%rcx,8),%xmm0
xorpd  %xmm1,%xmm1
mov    0x70(%rsp),%rcx
addsd  (%r8,%rcx,8),%xmm0
movapd %xmm1,%xmm2
mov    0x68(%rsp),%rcx
addsd  (%r8,%rcx,8),%xmm0
movapd %xmm2,0x20(%rsp)
call   1030 <FABS@plt>
mulsd  0x894(%rip),%xmm0
movsd  0x10(%rsp),%xmm5
addsd  %xmm0,%xmm5
movapd 0x20(%rsp),%xmm0
mulpd  0x89c(%rip),%xmm0
movapd 0x50(%rsp),%xmm4
addpd  %xmm0,%xmm4
add    $0x1,%rbp
mov    0xb8(%rsp),%rax
add    %rax,%r13
add    %rax,%r15
add    %rax,%rbx
add    %rax,%r12
cmp    0x40(%rsp),%rbp
mov    0x8(%rsp),%rdi
mov    0xb0(%rsp),%rax
jne    1520 <coreGammaAsc+0x400>
jmp    17d0 <coreGammaAsc+0x6b0>
movapd %xmm3,0x20(%rsp)
jmp    175f <coreGammaAsc+0x63f>
xorpd  %xmm4,%xmm4
xorpd  %xmm5,%xmm5
movsd  0x830(%rip),%xmm0
addsd  %xmm5,%xmm0
movapd %xmm4,%xmm1
divsd  %xmm0,%xmm1
mov    0x90(%rsp),%rax
movsd  %xmm1,(%rax)
movapd %xmm4,%xmm1
mulsd  %xmm4,%xmm1
unpckhpd %xmm4,%xmm4
mulsd  %xmm0,%xmm4
subsd  %xmm1,%xmm4
mulsd  %xmm0,%xmm0
divsd  %xmm0,%xmm4
mov    0x98(%rsp),%rax
movsd  %xmm4,(%rax)
movapd %xmm5,%xmm0
add    $0x20c8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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