void solution_update_upwind15(int n, double *curr, double const *last,
                              double alpha, double dx, double dt, double bc_0,
                              double bc_1) {
  double const f2 = 1.0 / 24;
  double const f1 = 1.0 / 6;
  double const f0 = 1.0 / 4;
  double const k = alpha * alpha * dt / (dx * dx);
  double const k2 = k * k;
  int i;
  curr[0] = bc_0;
  curr[1] = last[1] + k * (last[0] - 2 * last[1] + last[2]);
  curr[n - 2] = last[n - 2] + k * (last[n - 3] - 2 * last[n - 2] + last[n - 1]);
  curr[n - 1] = bc_1;
  for (i = 2; i < n - 2; i++)
    curr[i] = f2 * (12 * k2 - 2 * k) * last[i - 2] +
              f2 * (12 * k2 - 2 * k) * last[i + 2] -
              f1 * (12 * k2 - 8 * k) * last[i - 1] -
              f1 * (12 * k2 - 8 * k) * last[i + 1] +
              f0 * (12 * k2 - 10 * k + 4) * last[i];
}


x86 assembly (clang, -O0):
<solution_update_upwind15>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
movsd  %xmm2,-0x30(%rbp)
movsd  %xmm3,-0x38(%rbp)
movsd  %xmm4,-0x40(%rbp)
movsd  0xee0(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0xecb(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xeb6(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x20(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
mulsd  -0x28(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm1
movsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm3
addsd  %xmm3,%xmm3
subsd  %xmm3,%xmm2
mov    -0x18(%rbp),%rax
addsd  0x10(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x8(%rax)
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x3,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm3
addsd  %xmm3,%xmm3
subsd  %xmm3,%xmm2
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x2,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jge    148f <solution_update_upwind15+0x38f>
movsd  -0x68(%rbp),%xmm0
movsd  0xdbf(%rip),%xmm1
mulsd  -0x60(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xd83(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xd63(%rip),%xmm1
mulsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
mov    -0x6c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x68(%rbp),%xmm3
movsd  0xd68(%rip),%xmm0
mulsd  -0x60(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xd2c(%rip),%xmm4
mulsd  %xmm4,%xmm3
addsd  %xmm0,%xmm3
movsd  0xd0c(%rip),%xmm0
mulsd  %xmm3,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x6c(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm2
movsd  0xd01(%rip),%xmm0
mulsd  -0x60(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xccd(%rip),%xmm3
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
movsd  0xca5(%rip),%xmm0
mulsd  %xmm2,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x6c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x68(%rbp),%xmm2
movsd  0xc8b(%rip),%xmm1
mulsd  -0x60(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xc57(%rip),%xmm3
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  0xc2f(%rip),%xmm1
mulsd  %xmm2,%xmm1
mov    -0x18(%rbp),%rax
mov    -0x6c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm2
movsd  0xc0d(%rip),%xmm0
mulsd  -0x60(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xbe1(%rip),%xmm3
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
movsd  0xbc9(%rip),%xmm0
addsd  %xmm0,%xmm2
movsd  0xba5(%rip),%xmm0
mulsd  %xmm2,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    125b <solution_update_upwind15+0x15b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<solution_update_upwind15>:
mulsd  %xmm0,%xmm0
mulsd  %xmm2,%xmm0
movsd  %xmm3,(%rsi)
movsd  (%rdx),%xmm2
movsd  0x8(%rdx),%xmm3
movapd %xmm3,%xmm5
addsd  %xmm3,%xmm5
subsd  %xmm5,%xmm2
mulsd  %xmm1,%xmm1
addsd  0x10(%rdx),%xmm2
divsd  %xmm1,%xmm0
mulsd  %xmm0,%xmm2
addsd  %xmm3,%xmm2
movsd  %xmm2,0x8(%rsi)
movslq %edi,%rax
movsd  -0x18(%rdx,%rax,8),%xmm1
movsd  -0x10(%rdx,%rax,8),%xmm2
movapd %xmm2,%xmm3
addsd  %xmm2,%xmm3
subsd  %xmm3,%xmm1
addsd  -0x8(%rdx,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,-0x10(%rsi,%rax,8)
movsd  %xmm4,-0x8(%rsi,%rax,8)
cmp    $0x5,%edi
jl     12f2 <solution_update_upwind15+0x1f2>
movapd %xmm0,%xmm2
mulsd  %xmm0,%xmm2
add    $0xfffffffe,%edi
movapd %xmm0,%xmm3
mulsd  0xe70(%rip),%xmm2
addsd  %xmm0,%xmm3
movsd  0xe74(%rip),%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
mulsd  0xe74(%rip),%xmm0
addsd  %xmm2,%xmm0
subsd  %xmm3,%xmm2
mulsd  0xe4c(%rip),%xmm2
mulsd  0xe54(%rip),%xmm1
addsd  0xe5c(%rip),%xmm0
mulsd  0xe5c(%rip),%xmm0
mov    %edi,%eax
lea    -0x2(%rax),%r8
mov    $0x2,%ecx
cmp    $0x2,%r8
jb     12a0 <solution_update_upwind15+0x1a0>
lea    0x10(%rsi),%r9
lea    (%rdx,%rax,8),%rdi
add    $0x10,%rdi
cmp    %rdi,%r9
jae    1207 <solution_update_upwind15+0x107>
lea    (%rsi,%rax,8),%rdi
cmp    %rdx,%rdi
ja     12a0 <solution_update_upwind15+0x1a0>
mov    %r8,%r9
and    $0xfffffffffffffffe,%r9
lea    0x2(%r9),%rcx
movapd %xmm2,%xmm9
unpcklpd %xmm2,%xmm9
movapd %xmm1,%xmm10
unpcklpd %xmm1,%xmm10
movapd %xmm0,%xmm8
unpcklpd %xmm0,%xmm8
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd 0x20(%rdx,%rdi,8),%xmm6
mulpd  %xmm9,%xmm6
movupd (%rdx,%rdi,8),%xmm7
movupd 0x8(%rdx,%rdi,8),%xmm5
movupd 0x10(%rdx,%rdi,8),%xmm3
mulpd  %xmm8,%xmm3
movupd 0x18(%rdx,%rdi,8),%xmm4
mulpd  %xmm9,%xmm7
addpd  %xmm6,%xmm7
mulpd  %xmm10,%xmm5
addpd  %xmm7,%xmm5
mulpd  %xmm10,%xmm4
addpd  %xmm5,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm4,0x10(%rsi,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%r9
jne    1240 <solution_update_upwind15+0x140>
cmp    %r9,%r8
je     12f2 <solution_update_upwind15+0x1f2>
nopw   0x0(%rax,%rax,1)
movsd  0x10(%rdx,%rcx,8),%xmm3
mulsd  %xmm2,%xmm3
movsd  -0x10(%rdx,%rcx,8),%xmm4
mulsd  %xmm2,%xmm4
addsd  %xmm3,%xmm4
movsd  -0x8(%rdx,%rcx,8),%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm4,%xmm3
movsd  0x8(%rdx,%rcx,8),%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm3,%xmm4
movsd  (%rdx,%rcx,8),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm4,%xmm3
movsd  %xmm3,(%rsi,%rcx,8)
lea    0x1(%rcx),%rdi
mov    %rdi,%rcx
cmp    %rdi,%rax
jne    12a0 <solution_update_upwind15+0x1a0>
ret

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

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