void l_coeff(double *CC, unsigned int order, double CENTER, double WIDTH) {
  CC[0] = 1;
  unsigned int size = (order + 1) * (order + 1);
  for (unsigned int i = 1; i < size; ++i)
    CC[i] = 0;
  if (order == 0)
    return;
  CC[order] = -CENTER / WIDTH;
  CC[order + 1] = 1 / WIDTH;
  double *prev = CC;
  double *current = CC + order + 1;
  for (int n = 1; n < order; ++n) {
    double *next = current + order + 1;
    next[0] = n * prev[0] / (n + 1);
    for (int i = 0; i < n; ++i)
      next[i + 1] = ((2 * n + 1) * current[i] + prev[i + 1]) / (n + 1);
    prev = current;
    current = next;
  }
}


x86 assembly (clang, -O0):
<l_coeff>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movsd  0xedf(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    %eax,-0x24(%rbp)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jae    1167 <l_coeff+0x67>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113e <l_coeff+0x3e>
cmpl   $0x0,-0xc(%rbp)
jne    1176 <l_coeff+0x76>
jmp    12be <l_coeff+0x1be>
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
movsd  %xmm0,(%rax,%rcx,8)
movsd  0xe55(%rip),%xmm0
divsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
shl    $0x3,%rcx
add    %rcx,%rax
add    $0x8,%rax
mov    %rax,-0x38(%rbp)
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    12be <l_coeff+0x1be>
mov    -0x38(%rbp),%rax
mov    -0xc(%rbp),%ecx
shl    $0x3,%rcx
add    %rcx,%rax
add    $0x8,%rax
mov    %rax,-0x48(%rbp)
cvtsi2sdl -0x3c(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x48(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    12a0 <l_coeff+0x1a0>
mov    -0x3c(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm0
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x30(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x48(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1232 <l_coeff+0x132>
mov    -0x38(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x48(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11e6 <l_coeff+0xe6>
pop    %rbp
ret

x86 assembly (clang, -O2):
<l_coeff>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %esi,%r14d
mov    %rdi,%rbx
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rdi)
lea    0x1(%r14),%ebp
mov    %ebp,%eax
imul   %eax,%eax
cmp    $0x2,%eax
jb     116c <l_coeff+0x5c>
lea    0x8(%rbx),%rdi
add    $0xfffffffe,%eax
lea    0x8(,%rax,8),%rdx
xor    %esi,%esi
movsd  %xmm1,0x8(%rsp)
movaps %xmm0,0x10(%rsp)
call   1030 <memset@plt>
movapd 0x10(%rsp),%xmm0
movsd  0x8(%rsp),%xmm1
test   %r14d,%r14d
je     11a3 <l_coeff+0x93>
xorpd  0xe87(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    %r14d,%eax
mov    %rax,0x30(%rsp)
movsd  %xmm0,(%rbx,%rax,8)
movsd  0xe7e(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    %ebp,%eax
movsd  %xmm0,(%rbx,%rax,8)
cmp    $0x1,%r14d
jne    11b2 <l_coeff+0xa2>
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    0x30(%rsp),%rcx
lea    (%rbx,%rcx,8),%r8
add    $0x8,%r8
mov    %rcx,%rax
shl    $0x4,%rax
lea    (%rax,%rbx,1),%r14
add    $0x18,%r14
lea    0x8(,%rcx,8),%r12
lea    0x20(%rax,%rbx,1),%rdx
mov    %rdx,0x58(%rsp)
lea    0x10(,%rcx,8),%rdx
mov    %rdx,0x50(%rsp)
lea    0x10(%rbx,%rcx,8),%r9
lea    0x8(%rbx),%rcx
mov    %rcx,0x40(%rsp)
lea    0x10(%rbx),%rcx
mov    %rcx,0x38(%rsp)
lea    (%rax,%rbx,1),%rdi
add    $0x28,%rdi
mov    $0x1,%r15d
mov    $0xffffffffffffffff,%r10
xor    %r13d,%r13d
mov    %r9,0x48(%rsp)
mov    %r14,0x60(%rsp)
jmp    1250 <l_coeff+0x140>
nopl   0x0(%rax,%rax,1)
add    %r12,%rdi
add    %r12,%rbx
add    $0x1,%r10
add    %r12,%r14
add    %r12,%r9
mov    0x8(%rsp),%r8
cmp    0x30(%rsp),%r15
je     11a3 <l_coeff+0x93>
mov    %r15,%rbp
mov    %r13,%rsi
add    $0x1,%r13
xorps  %xmm1,%xmm1
cvtsi2sd %ebp,%xmm1
lea    (%r8,%r12,1),%rax
mov    %rax,0x8(%rsp)
mulsd  (%rbx),%xmm1
add    $0x1,%r15
xorps  %xmm0,%xmm0
cvtsi2sd %r15d,%xmm0
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r8,%r12,1)
lea    0x1(%rbp,%rbp,1),%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
cmp    $0x2,%r13
jae    1330 <l_coeff+0x220>
xor    %r11d,%r11d
mov    %r11,%rax
test   $0x1,%bpl
je     12c9 <l_coeff+0x1b9>
mov    %r11,%rax
or     $0x1,%rax
movsd  (%r8,%r11,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rbx,%rax,8),%xmm2
divsd  %xmm0,%xmm2
mov    0x8(%rsp),%rcx
movsd  %xmm2,(%rcx,%rax,8)
cmp    %r11,%rsi
je     1230 <l_coeff+0x120>
sub    %rax,%rbp
lea    (%r14,%rax,8),%rcx
lea    (%rbx,%rax,8),%rdx
lea    (%r9,%rax,8),%rsi
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  -0x8(%rsi,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  0x8(%rdx,%rax,8),%xmm2
divsd  %xmm0,%xmm2
movsd  %xmm2,(%rcx,%rax,8)
movsd  (%rsi,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  0x10(%rdx,%rax,8),%xmm2
divsd  %xmm0,%xmm2
movsd  %xmm2,0x8(%rcx,%rax,8)
add    $0x2,%rax
cmp    %rax,%rbp
jne    12f0 <l_coeff+0x1e0>
jmp    1230 <l_coeff+0x120>
mov    %r15,0x10(%rsp)
mov    %r13,0x28(%rsp)
mov    %r12,%rcx
imul   %rsi,%rcx
mov    0x60(%rsp),%rax
lea    (%rax,%rcx,1),%r13
mov    0x50(%rsp),%rdx
mov    %rsi,0x20(%rsp)
imul   %rsi,%rdx
mov    0x58(%rsp),%rsi
add    %rdx,%rsi
mov    0x38(%rsp),%rax
lea    (%rax,%rdx,1),%r11
add    0x40(%rsp),%rcx
add    0x48(%rsp),%rdx
cmp    %rdx,%r13
setb   %r15b
cmp    %r8,%rsi
seta   %dl
cmp    %r11,%r13
setb   %al
cmp    %rsi,%rcx
setb   %cl
test   %dl,%r15b
jne    13d7 <l_coeff+0x2c7>
mov    $0x0,%r11d
and    %cl,%al
mov    0x28(%rsp),%r13
mov    0x10(%rsp),%r15
mov    0x20(%rsp),%rsi
jne    129c <l_coeff+0x18c>
lea    -0x1(%rsi),%rcx
mov    %rcx,%rax
shr    %rax
add    $0x1,%rax
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
cmp    $0x2,%rcx
jae    13ee <l_coeff+0x2de>
xor    %r15d,%r15d
jmp    144a <l_coeff+0x33a>
xor    %r11d,%r11d
mov    0x28(%rsp),%r13
mov    0x10(%rsp),%r15
mov    0x20(%rsp),%rsi
jmp    129c <l_coeff+0x18c>
mov    %r10,%rcx
shr    %rcx
add    $0x1,%rcx
and    $0xfffffffffffffffe,%rcx
xor    %r15d,%r15d
nop
movupd (%r8,%r15,8),%xmm4
movupd 0x8(%rbx,%r15,8),%xmm5
mulpd  %xmm2,%xmm4
addpd  %xmm5,%xmm4
divpd  %xmm3,%xmm4
movupd %xmm4,-0x10(%rdi,%r15,8)
movupd 0x10(%r8,%r15,8),%xmm4
movupd 0x18(%rbx,%r15,8),%xmm5
mulpd  %xmm2,%xmm4
addpd  %xmm5,%xmm4
divpd  %xmm3,%xmm4
movupd %xmm4,(%rdi,%r15,8)
add    $0x4,%r15
add    $0xfffffffffffffffe,%rcx
jne    1400 <l_coeff+0x2f0>
mov    %r13,%r11
and    $0xfffffffffffffffe,%r11
test   $0x1,%al
je     1480 <l_coeff+0x370>
movupd (%r8,%r15,8),%xmm4
shl    $0x3,%r15
or     $0x8,%r15
movupd (%rbx,%r15,1),%xmm5
mulpd  %xmm4,%xmm2
addpd  %xmm5,%xmm2
divpd  %xmm3,%xmm2
mov    0x8(%rsp),%rax
movupd %xmm2,(%rax,%r15,1)
cmp    %r11,%r13
mov    0x10(%rsp),%r15
je     1230 <l_coeff+0x120>
jmp    129c <l_coeff+0x18c>

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

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