#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 coreCatAsc(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, l;
  ki = 1.0;
  kisqr = 1.0;
  for (l = 1; l < numStates; l++) {
    diagptable[l * 4] = exp(EIGN[l - 1] * ki * lz);
    diagptable[l * 4 + 1] = EIGN[l - 1] * ki;
    diagptable[l * 4 + 2] = EIGN[l - 1] * EIGN[l - 1] * kisqr;
  }
  for (i = 0; i < upper; i++) {
    double *sum = &sumtable[i * numStates], tmp, inv_Li = 0.0, dlnLidlz = 0.0,
           d2lnLidlz2 = 0.0;
    inv_Li += sum[0];
    for (l = 1; l < numStates; l++) {
      inv_Li += (tmp = diagptable[l * 4] * sum[l]);
      dlnLidlz += tmp * diagptable[l * 4 + 1];
      d2lnLidlz2 += tmp * diagptable[l * 4 + 2];
    }
    inv_Li = FABS(inv_Li);
    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):
<coreCatAsc>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20a0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
mov    %r9d,-0x34(%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
movsd  0xe90(%rip),%xmm0
movsd  %xmm0,-0x2060(%rbp)
ff
movsd  0xe80(%rip),%xmm0
movsd  %xmm0,-0x2068(%rbp)
ff
movl   $0x1,-0x2070(%rbp)
mov    -0x2070(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    126d <coreCatAsc+0x14d>
mov    -0x8(%rbp),%rax
mov    -0x2070(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x2060(%rbp),%xmm0
ff
mulsd  -0x30(%rbp),%xmm0
call   1040 <exp@plt>
mov    -0x2070(%rbp),%eax
shl    $0x2,%eax
cltq
movsd  %xmm0,-0x2040(%rbp,%rax,8)
ff ff
mov    -0x8(%rbp),%rax
mov    -0x2070(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x2060(%rbp),%xmm0
ff
mov    -0x2070(%rbp),%eax
shl    $0x2,%eax
add    $0x1,%eax
cltq
movsd  %xmm0,-0x2040(%rbp,%rax,8)
ff ff
mov    -0x8(%rbp),%rax
mov    -0x2070(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2070(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x2068(%rbp),%xmm0
ff
mov    -0x2070(%rbp),%eax
shl    $0x2,%eax
add    $0x2,%eax
cltq
movsd  %xmm0,-0x2040(%rbp,%rax,8)
ff ff
mov    -0x2070(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2070(%rbp)
jmp    1192 <coreCatAsc+0x72>
movl   $0x0,-0x206c(%rbp)
mov    -0x206c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    142c <coreCatAsc+0x30c>
mov    -0x10(%rbp),%rax
mov    -0x206c(%rbp),%ecx
imul   -0x34(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x2078(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2088(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2090(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2098(%rbp)
ff
mov    -0x2078(%rbp),%rax
movsd  (%rax),%xmm0
addsd  -0x2088(%rbp),%xmm0
ff
movsd  %xmm0,-0x2088(%rbp)
ff
movl   $0x1,-0x2070(%rbp)
mov    -0x2070(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    13bb <coreCatAsc+0x29b>
mov    -0x2070(%rbp),%eax
shl    $0x2,%eax
cltq
movsd  -0x2040(%rbp,%rax,8),%xmm0
ff ff
mov    -0x2078(%rbp),%rax
movslq -0x2070(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x2080(%rbp)
ff
addsd  -0x2088(%rbp),%xmm0
ff
movsd  %xmm0,-0x2088(%rbp)
ff
movsd  -0x2080(%rbp),%xmm0
ff
mov    -0x2070(%rbp),%eax
shl    $0x2,%eax
add    $0x1,%eax
cltq
movsd  -0x2040(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x2090(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x2090(%rbp)
ff
movsd  -0x2080(%rbp),%xmm0
ff
mov    -0x2070(%rbp),%eax
shl    $0x2,%eax
add    $0x2,%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
mov    -0x2070(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2070(%rbp)
jmp    12eb <coreCatAsc+0x1cb>
movsd  -0x2088(%rbp),%xmm0
ff
call   1030 <FABS@plt>
movsd  %xmm0,-0x2088(%rbp)
ff
movsd  -0x2088(%rbp),%xmm0
ff
addsd  -0x2048(%rbp),%xmm0
ff
movsd  %xmm0,-0x2048(%rbp)
ff
movsd  -0x2090(%rbp),%xmm0
ff
addsd  -0x2050(%rbp),%xmm0
ff
movsd  %xmm0,-0x2050(%rbp)
ff
movsd  -0x2098(%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    1277 <coreCatAsc+0x157>
movsd  -0x2050(%rbp),%xmm0
ff
movsd  -0x2048(%rbp),%xmm1
ff
movsd  0xbbc(%rip),%xmm2
subsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x2048(%rbp),%xmm0
ff
movsd  0xb9c(%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  0xb51(%rip),%xmm2
subsd  %xmm2,%xmm1
movsd  -0x2048(%rbp),%xmm2
ff
movsd  0xb3d(%rip),%xmm3
subsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x2048(%rbp),%xmm0
ff
add    $0x20a0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<coreCatAsc>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x2058,%rsp
mov    %r8,0x30(%rsp)
mov    %rcx,0x28(%rsp)
mov    %edx,%r12d
mov    %rsi,%r14
mov    %r9d,%r13d
mov    %r9d,0xc(%rsp)
cmp    $0x2,%r9d
jl     11ad <coreCatAsc+0x8d>
movapd %xmm0,%xmm1
mov    %rdi,%rbp
lea    -0x8(,%r13,8),%r15
ff
xor    %ebx,%ebx
movsd  %xmm0,0x10(%rsp)
cs nopw 0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rbx,1),%xmm0
mulsd  %xmm1,%xmm0
call   1040 <exp@plt>
movsd  0x10(%rsp),%xmm1
movsd  %xmm0,0x70(%rsp,%rbx,4)
movsd  0x0(%rbp,%rbx,1),%xmm0
movsd  %xmm0,0x78(%rsp,%rbx,4)
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x80(%rsp,%rbx,4)
add    $0x8,%rbx
cmp    %rbx,%r15
jne    1170 <coreCatAsc+0x50>
test   %r12d,%r12d
jle    131a <coreCatAsc+0x1fa>
movslq 0xc(%rsp),%rbp
mov    %r12d,%eax
mov    %rax,0x40(%rsp)
add    $0xffffffffffffffff,%r13
lea    0x8(%r14),%r15
lea    0x0(,%rbp,8),%rbx
mov    %r13,%rax
mov    %r13,0x38(%rsp)
and    $0xfffffffffffffffe,%r13
neg    %r13
xorpd  %xmm0,%xmm0
movapd %xmm0,0x10(%rsp)
xor    %r12d,%r12d
xorpd  %xmm3,%xmm3
jmp    1235 <coreCatAsc+0x115>
cs nopw 0x0(%rax,%rax,1)
nop
xorpd  %xmm2,%xmm2
movapd 0x10(%rsp),%xmm1
addpd  %xmm2,%xmm1
movapd %xmm1,0x10(%rsp)
call   1030 <FABS@plt>
movsd  0x48(%rsp),%xmm3
addsd  %xmm0,%xmm3
add    $0x1,%r12
add    %rbx,%r15
cmp    0x40(%rsp),%r12
je     1328 <coreCatAsc+0x208>
mov    %r12,%rax
imul   %rbp,%rax
movsd  (%r14,%rax,8),%xmm0
addsd  0xdbe(%rip),%xmm0
cmpl   $0x2,0xc(%rsp)
movsd  %xmm3,0x48(%rsp)
jl     1200 <coreCatAsc+0xe0>
xorpd  %xmm1,%xmm1
mov    $0x1,%ecx
jne    1270 <coreCatAsc+0x150>
testb  $0x1,0x38(%rsp)
je     1204 <coreCatAsc+0xe4>
jmp    12e9 <coreCatAsc+0x1c9>
nopl   0x0(%rax,%rax,1)
lea    0x98(%rsp),%rdx
nopl   0x0(%rax,%rax,1)
movsd  -0x28(%rdx),%xmm3
mulsd  -0x8(%r15,%rcx,8),%xmm3
addsd  %xmm3,%xmm0
movupd -0x20(%rdx),%xmm2
unpcklpd %xmm3,%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm1,%xmm3
movsd  -0x8(%rdx),%xmm2
mulsd  (%r15,%rcx,8),%xmm2
addsd  %xmm2,%xmm0
movupd (%rdx),%xmm1
unpcklpd %xmm2,%xmm2
mulpd  %xmm1,%xmm2
movapd %xmm2,%xmm1
addpd  %xmm3,%xmm1
add    $0x40,%rdx
lea    (%rcx,%r13,1),%rsi
add    $0x2,%rsi
add    $0x2,%rcx
cmp    $0x1,%rsi
jne    1280 <coreCatAsc+0x160>
addpd  %xmm3,%xmm2
testb  $0x1,0x38(%rsp)
je     1204 <coreCatAsc+0xe4>
lea    (%r14,%rax,8),%rax
mov    %rcx,%rdx
shl    $0x5,%rdx
movsd  0x50(%rsp,%rdx,1),%xmm2
mulsd  (%rax,%rcx,8),%xmm2
addsd  %xmm2,%xmm0
movupd 0x58(%rsp,%rdx,1),%xmm3
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
jmp    1204 <coreCatAsc+0xe4>
xorpd  %xmm0,%xmm0
movapd %xmm0,0x10(%rsp)
xorpd  %xmm3,%xmm3
movsd  0xcd0(%rip),%xmm0
addsd  %xmm3,%xmm0
movapd 0x10(%rsp),%xmm2
movapd %xmm2,%xmm1
divsd  %xmm0,%xmm1
mov    0x28(%rsp),%rax
movsd  %xmm1,(%rax)
movapd %xmm2,%xmm1
mulsd  %xmm2,%xmm1
unpckhpd %xmm2,%xmm2
mulsd  %xmm0,%xmm2
subsd  %xmm1,%xmm2
mulsd  %xmm0,%xmm0
divsd  %xmm0,%xmm2
mov    0x30(%rsp),%rax
movsd  %xmm2,(%rax)
movapd %xmm3,%xmm0
add    $0x2058,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @coreCatAsc(double* nocapture %0, double* nocapture %1, i32 %2, double* %3, double* %4, double %5, i32 %6) {
%8 = alloca [1024 x double]
%9 = bitcast [1024 x double]* %8 to i8*
%10 = icmp sgt i32 %6, 1
br i1 %10, label %11, label %13
11: ; preds = %7
%12 = zext i32 %6 to i64
br label %24
13: ; preds = %24, %7
%14 = icmp sgt i32 %2, 0
br i1 %14, label %15, label %118
15: ; preds = %13
%16 = sext i32 %6 to i64
%17 = zext i32 %2 to i64
%18 = zext i32 %6 to i64
%19 = add nsw i64 %18, -1
%20 = and i64 %19, 1
%21 = icmp eq i32 %6, 2
%22 = and i64 %19, -2
%23 = icmp eq i64 %20, 0
br label %41
24: ; preds = %11, %24
%25 = phi i64 [ 1, %11 ], [ %39, %24 ]
%26 = add nsw i64 %25, -1
%27 = getelementptr inbounds double, double* %0, i64 %26
%28 = load double, double* %27
%29 = fmul double %28, %5
%30 = tail call double @exp(double %29)
%31 = shl nsw i64 %25, 2
%32 = getelementptr inbounds [1024 x double], [1024 x double]* %8, i64 0, i64 %31
store double %30, double* %32
%33 = load double, double* %27
%34 = or i64 %31, 1
%35 = getelementptr inbounds [1024 x double], [1024 x double]* %8, i64 0, i64 %34
store double %33, double* %35
%36 = fmul double %33, %33
%37 = or i64 %31, 2
%38 = getelementptr inbounds [1024 x double], [1024 x double]* %8, i64 0, i64 %37
store double %36, double* %38
%39 = add nuw nsw i64 %25, 1
%40 = icmp eq i64 %39, %12
br i1 %40, label %13, label %24
41: ; preds = %15, %110
%42 = phi i64 [ 0, %15 ], [ %116, %110 ]
%43 = phi double [ 0.000000e+00, %15 ], [ %114, %110 ]
%44 = phi <2 x double> [ zeroinitializer, %15 ], [ %115, %110 ]
%45 = mul nsw i64 %42, %16
%46 = getelementptr inbounds double, double* %1, i64 %45
%47 = load double, double* %46
%48 = fadd double %47, 0.000000e+00
br i1 %10, label %49, label %110
49: ; preds = %41
br i1 %21, label %89, label %50
50: ; preds = %49, %50
%51 = phi i64 [ %84, %50 ], [ 1, %49 ]
%52 = phi double [ %76, %50 ], [ %48, %49 ]
%53 = phi <2 x double> [ %83, %50 ], [ zeroinitializer, %49 ]
%54 = phi i64 [ %85, %50 ], [ 0, %49 ]
%55 = shl nsw i64 %51, 2
%56 = getelementptr inbounds [1024 x double], [1024 x double]* %8, i64 0, i64 %55
%57 = load double, double* %56
%58 = getelementptr inbounds double, double* %46, i64 %51
%59 = load double, double* %58
%60 = fmul double %57, %59
%61 = fadd double %52, %60
%62 = or i64 %55, 1
%63 = getelementptr inbounds [1024 x double], [1024 x double]* %8, i64 0, i64 %62
%64 = bitcast double* %63 to <2 x double>*
%65 = load <2 x double>, <2 x double>* %64
%66 = insertelement <2 x double> poison, double %60, i64 0
%67 = shufflevector <2 x double> %66, <2 x double> poison, <2 x i32> zeroinitializer
%68 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %67, <2 x double> %65, <2 x double> %53)
%69 = add nuw nsw i64 %51, 1
%70 = shl nsw i64 %69, 2
%71 = getelementptr inbounds [1024 x double], [1024 x double]* %8, i64 0, i64 %70
%72 = load double, double* %71
%73 = getelementptr inbounds double, double* %46, i64 %69
%74 = load double, double* %73
%75 = fmul double %72, %74
%76 = fadd double %61, %75
%77 = or i64 %70, 1
%78 = getelementptr inbounds [1024 x double], [1024 x double]* %8, i64 0, i64 %77
%79 = bitcast double* %78 to <2 x double>*
%80 = load <2 x double>, <2 x double>* %79
%81 = insertelement <2 x double> poison, double %75, i64 0
%82 = shufflevector <2 x double> %81, <2 x double> poison, <2 x i32> zeroinitializer
%83 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %82, <2 x double> %80, <2 x double> %68)
%84 = add nuw nsw i64 %51, 2
%85 = add i64 %54, 2
%86 = icmp eq i64 %85, %22
br i1 %86, label %87, label %50
87: ; preds = %50
%88 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %82, <2 x double> %80, <2 x double> %68)
br label %89
89: ; preds = %87, %49
%90 = phi double [ undef, %49 ], [ %76, %87 ]
%91 = phi <2 x double> [ undef, %49 ], [ %88, %87 ]
%92 = phi i64 [ 1, %49 ], [ %84, %87 ]
%93 = phi double [ %48, %49 ], [ %76, %87 ]
%94 = phi <2 x double> [ zeroinitializer, %49 ], [ %83, %87 ]
br i1 %23, label %110, label %95
95: ; preds = %89
%96 = shl nsw i64 %92, 2
%97 = getelementptr inbounds [1024 x double], [1024 x double]* %8, i64 0, i64 %96
%98 = load double, double* %97
%99 = getelementptr inbounds double, double* %46, i64 %92
%100 = load double, double* %99
%101 = fmul double %98, %100
%102 = fadd double %93, %101
%103 = or i64 %96, 1
%104 = getelementptr inbounds [1024 x double], [1024 x double]* %8, i64 0, i64 %103
%105 = bitcast double* %104 to <2 x double>*
%106 = load <2 x double>, <2 x double>* %105
%107 = insertelement <2 x double> poison, double %101, i64 0
%108 = shufflevector <2 x double> %107, <2 x double> poison, <2 x i32> zeroinitializer
%109 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %108, <2 x double> %106, <2 x double> %94)
br label %110
110: ; preds = %95, %89, %41
%111 = phi double [ %48, %41 ], [ %90, %89 ], [ %102, %95 ]
%112 = phi <2 x double> [ zeroinitializer, %41 ], [ %91, %89 ], [ %109, %95 ]
%113 = tail call double @FABS(double %111)
%114 = fadd double %43, %113
%115 = fadd <2 x double> %44, %112
%116 = add nuw nsw i64 %42, 1
%117 = icmp eq i64 %116, %17
br i1 %117, label %118, label %41
118: ; preds = %110, %13
%119 = phi double [ 0.000000e+00, %13 ], [ %114, %110 ]
%120 = phi <2 x double> [ zeroinitializer, %13 ], [ %115, %110 ]
%121 = fadd double %119, -1.000000e+00
%122 = extractelement <2 x double> %120, i64 0
%123 = fdiv double %122, %121
store volatile double %123, double* %3
%124 = fneg double %122
%125 = fmul double %122, %124
%126 = extractelement <2 x double> %120, i64 1
%127 = tail call double @llvm.fmuladd.f64(double %121, double %126, double %125)
%128 = fmul double %121, %121
%129 = fdiv double %127, %128
store volatile double %129, double* %4
ret double %119
}