double mglsadff1(double x, double *b, const int m, const double a,
                 const double g, double *d) {
  int i;
  double y, aa;
  aa = 1 - a * a;
  y = d[0] * b[1];
  for (i = 1; i < m; i++) {
    d[i] += a * (d[i + 1] - d[i - 1]);
    y += d[i] * b[i + 1];
  }
  x -= g * y;
  for (i = m; i > 0; i--)
    d[i] = d[i - 1];
  d[0] = a * d[0] + aa * x;
  return (x);
}


x86 assembly (clang, -O0):
<mglsadff1>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
movsd  %xmm1,-0x20(%rbp)
movsd  %xmm2,-0x28(%rbp)
mov    %rdx,-0x30(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xeb5(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11ff <mglsadff1+0xff>
movsd  -0x20(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm2
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1171 <mglsadff1+0x71>
movsd  -0x28(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  -0x8(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jle    126f <mglsadff1+0x16f>
mov    -0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    1238 <mglsadff1+0x138>
movsd  -0x20(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  -0x8(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<mglsadff1>:
push   %r14
push   %rbx
mov    %rdx,%r8
movsd  (%rdx),%xmm8
movsd  0x8(%rdi),%xmm4
mulsd  %xmm8,%xmm4
mov    %esi,%ecx
cmp    $0x2,%esi
jl     11f8 <mglsadff1+0xf8>
movsd  0x8(%r8),%xmm9
lea    -0x1(%rcx),%r9
cmp    $0x2,%esi
jne    113e <mglsadff1+0x3e>
mov    $0x1,%edx
movapd %xmm8,%xmm6
jmp    11ce <mglsadff1+0xce>
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
neg    %r10
mov    $0x1,%eax
movapd %xmm9,%xmm7
movapd %xmm8,%xmm6
nopw   0x0(%rax,%rax,1)
movsd  0x8(%r8,%rax,8),%xmm5
movsd  0x10(%r8,%rax,8),%xmm9
movapd %xmm5,%xmm3
subsd  %xmm6,%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm7,%xmm3
movsd  %xmm3,(%r8,%rax,8)
movapd %xmm9,%xmm6
subsd  %xmm3,%xmm6
mulsd  0x8(%rdi,%rax,8),%xmm3
addsd  %xmm4,%xmm3
lea    0x2(%rax),%rdx
mulsd  %xmm1,%xmm6
addsd  %xmm5,%xmm6
movsd  %xmm6,0x8(%r8,%rax,8)
movsd  0x10(%rdi,%rax,8),%xmm4
mulsd  %xmm6,%xmm4
addsd  %xmm3,%xmm4
lea    (%r10,%rax,1),%rbx
add    $0x2,%rbx
movapd %xmm9,%xmm7
mov    %rdx,%rax
cmp    $0x1,%rbx
jne    1160 <mglsadff1+0x60>
test   $0x1,%r9b
je     11f8 <mglsadff1+0xf8>
movsd  0x8(%r8,%rdx,8),%xmm3
subsd  %xmm6,%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm9,%xmm3
movsd  %xmm3,(%r8,%rdx,8)
mulsd  0x8(%rdi,%rdx,8),%xmm3
addsd  %xmm3,%xmm4
mulsd  %xmm4,%xmm2
subsd  %xmm2,%xmm0
test   %esi,%esi
jle    122a <mglsadff1+0x12a>
cmp    $0x4,%esi
jae    1255 <mglsadff1+0x155>
add    $0x1,%rcx
nopl   (%rax)
lea    -0x2(%rcx),%eax
movsd  (%r8,%rax,8),%xmm2
movsd  %xmm2,-0x8(%r8,%rcx,8)
add    $0xffffffffffffffff,%rcx
cmp    $0x1,%rcx
ja     1210 <mglsadff1+0x110>
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
movsd  0xdc6(%rip),%xmm3
subsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm3
mulsd  %xmm1,%xmm8
addsd  %xmm3,%xmm8
movsd  %xmm8,(%r8)
pop    %rbx
pop    %r14
ret
lea    -0x1(%rcx),%r14
add    $0xffffffff,%esi
mov    $0x8,%edx
mov    %r14,%rax
mul    %rdx
seto   %dl
cmp    %r14d,%esi
jb     1209 <mglsadff1+0x109>
mov    %r14,%rdi
shr    $0x20,%rdi
jne    1209 <mglsadff1+0x109>
lea    (%r8,%rcx,8),%rbx
mov    %rbx,%rdi
sub    %rax,%rdi
cmp    %rbx,%rdi
ja     1209 <mglsadff1+0x109>
test   %dl,%dl
jne    1209 <mglsadff1+0x109>
lea    (%r8,%rsi,8),%rbx
mov    %rbx,%rdi
sub    %rax,%rdi
cmp    %rbx,%rdi
ja     1209 <mglsadff1+0x109>
test   %dl,%dl
jne    1209 <mglsadff1+0x109>
lea    0x8(%r8),%rax
lea    (%r8,%rsi,8),%rdx
add    $0x8,%rdx
cmp    %rdx,%rax
jae    12d7 <mglsadff1+0x1d7>
lea    (%r8,%rcx,8),%rax
add    $0x8,%rax
add    $0x1,%rsi
sub    %rcx,%rsi
lea    (%r8,%rsi,8),%rdx
cmp    %rax,%rdx
jb     1209 <mglsadff1+0x109>
mov    $0xffffffff,%r9d
mov    %ecx,%r10d
and    $0xfffffffc,%r10d
lea    -0x4(%r10),%rax
mov    %rax,%r11
shr    $0x2,%r11
add    $0x1,%r11
test   %rax,%rax
je     139f <mglsadff1+0x29f>
mov    %r11,%rsi
and    $0xfffffffffffffffe,%rsi
lea    (%r8,%r14,8),%rax
lea    (%rcx,%r9,1),%r14
add    $0xfffffffffffffffc,%r14
lea    (%rcx,%r9,1),%rbx
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nop
lea    (%rbx,%rdi,1),%edx
movups -0x18(%r8,%rdx,8),%xmm2
movups -0x8(%r8,%rdx,8),%xmm3
movups %xmm3,(%rax,%rdi,8)
movups %xmm2,-0x10(%rax,%rdi,8)
lea    (%r14,%rdi,1),%edx
movupd -0x18(%r8,%rdx,8),%xmm2
movupd -0x8(%r8,%rdx,8),%xmm3
movupd %xmm3,-0x20(%rax,%rdi,8)
movupd %xmm2,-0x30(%rax,%rdi,8)
add    $0xfffffffffffffff8,%rdi
add    $0xfffffffffffffffe,%rsi
jne    1320 <mglsadff1+0x220>
neg    %rdi
test   $0x1,%r11b
je     138e <mglsadff1+0x28e>
mov    %rcx,%rax
sub    %rdi,%rax
add    %eax,%r9d
movupd -0x18(%r8,%r9,8),%xmm2
movupd -0x8(%r8,%r9,8),%xmm3
movupd %xmm3,-0x8(%r8,%rax,8)
movupd %xmm2,-0x18(%r8,%rax,8)
cmp    %rcx,%r10
je     122a <mglsadff1+0x12a>
and    $0x3,%ecx
jmp    1209 <mglsadff1+0x109>
xor    %edi,%edi
test   $0x1,%r11b
jne    1369 <mglsadff1+0x269>
jmp    138e <mglsadff1+0x28e>

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

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