double normpart(double *n, double *T, double *fracshift, int g) {
  double tmp1[g];
  double tmp2[g];
  int i, j;
  for (i = 0; i < g; i++)
    tmp1[i] = n[i] + fracshift[i];
  double sum;
  for (i = 0; i < g; i++) {
    sum = 0;
    for (j = 0; j < g; j++)
      sum += T[i * g + j] * tmp1[j];
    tmp2[i] = sum;
  }
  double norm = 0;
  for (i = 0; i < g; i++)
    norm += tmp2[i] * tmp2[i];
  return -3.14159265358979323846 * norm;
}


x86 assembly (clang, -O0):
<normpart>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    %rsp,%rcx
mov    %rcx,-0x28(%rbp)
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x60(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x30(%rbp)
mov    -0x1c(%rbp),%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x58(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11a6 <normpart+0xa6>
mov    -0x60(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x3c(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
mov    -0x18(%rbp),%rcx
movslq -0x3c(%rbp),%rdx
addsd  (%rcx,%rdx,8),%xmm0
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1165 <normpart+0x65>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1237 <normpart+0x137>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1217 <normpart+0x117>
mov    -0x60(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    -0x3c(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    -0x40(%rbp),%edx
movslq %edx,%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11c8 <normpart+0xc8>
mov    -0x58(%rbp),%rax
movsd  -0x48(%rbp),%xmm0
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11ad <normpart+0xad>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1288 <normpart+0x188>
mov    -0x58(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1246 <normpart+0x146>
movsd  0xd70(%rip),%xmm0
mulsd  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<normpart>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %ecx,%r9d
mov    %rsp,%r13
lea    0xf(,%r9,8),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%r13
mov    %r13,%rsp
mov    %rsp,%r8
sub    %rax,%r8
mov    %r8,%rsp
movsd  0xec9(%rip),%xmm0
test   %ecx,%ecx
jle    13de <normpart+0x2de>
cmp    $0x4,%ecx
jae    114c <normpart+0x4c>
xor    %r10d,%r10d
jmp    1220 <normpart+0x120>
mov    %r9d,%r10d
and    $0xfffffffc,%r10d
lea    -0x4(%r10),%rax
mov    %rax,%r11
shr    $0x2,%r11
add    $0x1,%r11
test   %rax,%rax
je     1348 <normpart+0x248>
mov    %r11,%r14
and    $0xfffffffffffffffe,%r14
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd (%rdi,%rbx,8),%xmm1
movupd 0x10(%rdi,%rbx,8),%xmm2
movupd (%rdx,%rbx,8),%xmm3
addpd  %xmm1,%xmm3
movupd 0x10(%rdx,%rbx,8),%xmm1
addpd  %xmm2,%xmm1
movapd %xmm3,0x0(%r13,%rbx,8)
movapd %xmm1,0x10(%r13,%rbx,8)
movupd 0x20(%rdi,%rbx,8),%xmm1
movupd 0x30(%rdi,%rbx,8),%xmm2
movupd 0x20(%rdx,%rbx,8),%xmm3
addpd  %xmm1,%xmm3
movupd 0x30(%rdx,%rbx,8),%xmm1
addpd  %xmm2,%xmm1
movapd %xmm3,0x20(%r13,%rbx,8)
movapd %xmm1,0x30(%r13,%rbx,8)
add    $0x8,%rbx
add    $0xfffffffffffffffe,%r14
jne    1180 <normpart+0x80>
test   $0x1,%r11b
je     1216 <normpart+0x116>
movupd (%rdi,%rbx,8),%xmm1
movupd 0x10(%rdi,%rbx,8),%xmm2
movupd (%rdx,%rbx,8),%xmm3
addpd  %xmm1,%xmm3
movupd 0x10(%rdx,%rbx,8),%xmm1
addpd  %xmm2,%xmm1
movapd %xmm3,0x0(%r13,%rbx,8)
movapd %xmm1,0x10(%r13,%rbx,8)
cmp    %r9,%r10
je     123c <normpart+0x13c>
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%r10,8),%xmm1
addsd  (%rdx,%r10,8),%xmm1
movsd  %xmm1,0x0(%r13,%r10,8)
add    $0x1,%r10
cmp    %r10,%r9
jne    1220 <normpart+0x120>
test   %ecx,%ecx
jle    13de <normpart+0x2de>
lea    -0x1(%r9),%r10
mov    %r9d,%r12d
and    $0x3,%r12d
mov    %r9d,%r11d
and    $0xfffffffc,%r11d
lea    0x18(%rsi),%rdx
lea    0x0(,%r9,8),%r14
xor    %r15d,%r15d
jmp    1289 <normpart+0x189>
nopw   0x0(%rax,%rax,1)
movsd  %xmm1,(%r8,%r15,8)
add    $0x1,%r15
add    %r14,%rdx
add    %r14,%rsi
cmp    %r9,%r15
je     132c <normpart+0x22c>
xorpd  %xmm1,%xmm1
xor    %ebx,%ebx
cmp    $0x3,%r10
jb     12ec <normpart+0x1ec>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  -0x18(%rdx,%rbx,8),%xmm2
movsd  -0x10(%rdx,%rbx,8),%xmm3
mulsd  0x0(%r13,%rbx,8),%xmm2
mulsd  0x8(%r13,%rbx,8),%xmm3
addsd  %xmm1,%xmm2
movsd  -0x8(%rdx,%rbx,8),%xmm4
mulsd  0x10(%r13,%rbx,8),%xmm4
addsd  %xmm2,%xmm3
movsd  (%rdx,%rbx,8),%xmm1
mulsd  0x18(%r13,%rbx,8),%xmm1
addsd  %xmm3,%xmm4
addsd  %xmm4,%xmm1
add    $0x4,%rbx
cmp    %rbx,%r11
jne    12a0 <normpart+0x1a0>
test   %r12,%r12
je     1270 <normpart+0x170>
lea    0x0(,%rbx,8),%rdi
add    %r13,%rdi
lea    (%rsi,%rbx,8),%rbx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
movsd  (%rbx,%rax,8),%xmm2
mulsd  (%rdi,%rax,8),%xmm2
addsd  %xmm2,%xmm1
add    $0x1,%rax
cmp    %rax,%r12
jne    1310 <normpart+0x210>
jmp    1270 <normpart+0x170>
test   %ecx,%ecx
jle    13de <normpart+0x2de>
mov    %r9d,%eax
and    $0x3,%eax
cmp    $0x3,%r10
jae    1359 <normpart+0x259>
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    13b4 <normpart+0x2b4>
xor    %ebx,%ebx
test   $0x1,%r11b
jne    11ea <normpart+0xea>
jmp    1216 <normpart+0x116>
and    $0xfffffffc,%r9d
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%r8,%rcx,8),%xmm1
movsd  0x8(%r8,%rcx,8),%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
movsd  0x10(%r8,%rcx,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
movsd  0x18(%r8,%rcx,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
add    $0x4,%rcx
cmp    %rcx,%r9
jne    1370 <normpart+0x270>
test   %rax,%rax
je     13d6 <normpart+0x2d6>
lea    (%r8,%rcx,8),%rcx
xor    %edx,%edx
nop
movsd  (%rcx,%rdx,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rdx
cmp    %rdx,%rax
jne    13c0 <normpart+0x2c0>
mulsd  0xc2a(%rip),%xmm0
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @normpart(double* %0, double* %1, double* %2, i32 %3) {
%5 = alloca double*
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca i8*
%10 = alloca i64
%11 = alloca i64
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
%15 = alloca double
store double* %0, double** %5
store double* %1, double** %6
store double* %2, double** %7
store i32 %3, i32* %8
%16 = load i32, i32* %8
%17 = zext i32 %16 to i64
%18 = call i8* @llvm.stacksave()
store i8* %18, i8** %9
%19 = alloca double, i64 %17
store i64 %17, i64* %10
%20 = load i32, i32* %8
%21 = zext i32 %20 to i64
%22 = alloca double, i64 %21
store i64 %21, i64* %11
store i32 0, i32* %12
br label %23
23: ; preds = %42, %4
%24 = load i32, i32* %12
%25 = load i32, i32* %8
%26 = icmp slt i32 %24, %25
br i1 %26, label %27, label %45
27: ; preds = %23
%28 = load double*, double** %5
%29 = load i32, i32* %12
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %28, i64 %30
%32 = load double, double* %31
%33 = load double*, double** %7
%34 = load i32, i32* %12
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %33, i64 %35
%37 = load double, double* %36
%38 = fadd double %32, %37
%39 = load i32, i32* %12
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %19, i64 %40
store double %38, double* %41
br label %42
42: ; preds = %27
%43 = load i32, i32* %12
%44 = add nsw i32 %43, 1
store i32 %44, i32* %12
br label %23
45: ; preds = %23
store i32 0, i32* %12
br label %46
46: ; preds = %79, %45
%47 = load i32, i32* %12
%48 = load i32, i32* %8
%49 = icmp slt i32 %47, %48
br i1 %49, label %50, label %82
50: ; preds = %46
store double 0.000000e+00, double* %14
store i32 0, i32* %13
br label %51
51: ; preds = %71, %50
%52 = load i32, i32* %13
%53 = load i32, i32* %8
%54 = icmp slt i32 %52, %53
br i1 %54, label %55, label %74
55: ; preds = %51
%56 = load double*, double** %6
%57 = load i32, i32* %12
%58 = load i32, i32* %8
%59 = mul nsw i32 %57, %58
%60 = load i32, i32* %13
%61 = add nsw i32 %59, %60
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %56, i64 %62
%64 = load double, double* %63
%65 = load i32, i32* %13
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double, double* %19, i64 %66
%68 = load double, double* %67
%69 = load double, double* %14
%70 = call double @llvm.fmuladd.f64(double %64, double %68, double %69)
store double %70, double* %14
br label %71
71: ; preds = %55
%72 = load i32, i32* %13
%73 = add nsw i32 %72, 1
store i32 %73, i32* %13
br label %51
74: ; preds = %51
%75 = load double, double* %14
%76 = load i32, i32* %12
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %22, i64 %77
store double %75, double* %78
br label %79
79: ; preds = %74
%80 = load i32, i32* %12
%81 = add nsw i32 %80, 1
store i32 %81, i32* %12
br label %46
82: ; preds = %46
store double 0.000000e+00, double* %15
store i32 0, i32* %12
br label %83
83: ; preds = %98, %82
%84 = load i32, i32* %12
%85 = load i32, i32* %8
%86 = icmp slt i32 %84, %85
br i1 %86, label %87, label %101
87: ; preds = %83
%88 = load i32, i32* %12
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %22, i64 %89
%91 = load double, double* %90
%92 = load i32, i32* %12
%93 = sext i32 %92 to i64
%94 = getelementptr inbounds double, double* %22, i64 %93
%95 = load double, double* %94
%96 = load double, double* %15
%97 = call double @llvm.fmuladd.f64(double %91, double %95, double %96)
store double %97, double* %15
br label %98
98: ; preds = %87
%99 = load i32, i32* %12
%100 = add nsw i32 %99, 1
store i32 %100, i32* %12
br label %83
101: ; preds = %83
%102 = load double, double* %15
%103 = fmul double 0xC00921FB54442D18, %102
%104 = load i8*, i8** %9
ret double %103
}

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