double HTS_mglsadff(double x, double *b, const int m, const double a,
                    double *d) {
  int i;
  double y;
  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 -= y;
  for (i = m; i > 0; i--)
    d[i] = d[i - 1];
  d[0] = a * d[0] + (1 - a * a) * x;
  return x;
}


x86 assembly (clang, -O0):
<HTS_mglsadff>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11c4 <HTS_mglsadff+0xc4>
movsd  -0x20(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1136 <HTS_mglsadff+0x36>
movsd  -0x38(%rbp),%xmm1
movsd  -0x8(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jle    1214 <HTS_mglsadff+0x114>
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    11dd <HTS_mglsadff+0xdd>
movsd  -0x20(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x20(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm3
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm3,%xmm1
movsd  0xdb2(%rip),%xmm3
addsd  %xmm3,%xmm1
mulsd  -0x8(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<HTS_mglsadff>:
push   %r14
push   %rbx
mov    %rdx,%r8
movsd  (%rdx),%xmm8
movsd  0x8(%rdi),%xmm3
mulsd  %xmm8,%xmm3
mov    %esi,%ecx
cmp    $0x2,%esi
jl     11f5 <HTS_mglsadff+0xf5>
movsd  0x8(%r8),%xmm4
lea    -0x1(%rcx),%r9
cmp    $0x2,%esi
jne    113e <HTS_mglsadff+0x3e>
mov    $0x1,%edx
movapd %xmm8,%xmm5
jmp    11cc <HTS_mglsadff+0xcc>
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
neg    %r10
mov    $0x1,%eax
movapd %xmm4,%xmm6
movapd %xmm8,%xmm5
cs nopw 0x0(%rax,%rax,1)
movsd  0x8(%r8,%rax,8),%xmm2
movsd  0x10(%r8,%rax,8),%xmm4
movapd %xmm2,%xmm7
subsd  %xmm5,%xmm7
mulsd  %xmm1,%xmm7
addsd  %xmm6,%xmm7
movsd  %xmm7,(%r8,%rax,8)
movapd %xmm4,%xmm5
subsd  %xmm7,%xmm5
mulsd  0x8(%rdi,%rax,8),%xmm7
addsd  %xmm3,%xmm7
lea    0x2(%rax),%rdx
mulsd  %xmm1,%xmm5
addsd  %xmm2,%xmm5
movsd  %xmm5,0x8(%r8,%rax,8)
movsd  0x10(%rdi,%rax,8),%xmm3
mulsd  %xmm5,%xmm3
addsd  %xmm7,%xmm3
lea    (%r10,%rax,1),%rbx
add    $0x2,%rbx
movapd %xmm4,%xmm6
mov    %rdx,%rax
cmp    $0x1,%rbx
jne    1160 <HTS_mglsadff+0x60>
test   $0x1,%r9b
je     11f5 <HTS_mglsadff+0xf5>
movsd  0x8(%r8,%rdx,8),%xmm2
subsd  %xmm5,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm4,%xmm2
movsd  %xmm2,(%r8,%rdx,8)
mulsd  0x8(%rdi,%rdx,8),%xmm2
addsd  %xmm2,%xmm3
subsd  %xmm3,%xmm0
test   %esi,%esi
jle    122a <HTS_mglsadff+0x12a>
cmp    $0x4,%esi
jae    1255 <HTS_mglsadff+0x155>
add    $0x1,%rcx
cs nopw 0x0(%rax,%rax,1)
lea    -0x2(%rcx),%eax
movsd  (%r8,%rax,8),%xmm2
movsd  %xmm2,-0x8(%r8,%rcx,8)
add    $0xffffffffffffffff,%rcx
cmp    $0x1,%rcx
ja     1210 <HTS_mglsadff+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     1202 <HTS_mglsadff+0x102>
mov    %r14,%rdi
shr    $0x20,%rdi
jne    1202 <HTS_mglsadff+0x102>
lea    (%r8,%rcx,8),%rbx
mov    %rbx,%rdi
sub    %rax,%rdi
cmp    %rbx,%rdi
ja     1202 <HTS_mglsadff+0x102>
test   %dl,%dl
jne    1202 <HTS_mglsadff+0x102>
lea    (%r8,%rsi,8),%rbx
mov    %rbx,%rdi
sub    %rax,%rdi
cmp    %rbx,%rdi
ja     1202 <HTS_mglsadff+0x102>
test   %dl,%dl
jne    1202 <HTS_mglsadff+0x102>
lea    0x8(%r8),%rax
lea    (%r8,%rsi,8),%rdx
add    $0x8,%rdx
cmp    %rdx,%rax
jae    12db <HTS_mglsadff+0x1db>
lea    (%r8,%rcx,8),%rax
add    $0x8,%rax
add    $0x1,%rsi
sub    %rcx,%rsi
lea    (%r8,%rsi,8),%rdx
cmp    %rax,%rdx
jb     1202 <HTS_mglsadff+0x102>
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 <HTS_mglsadff+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
nopl   0x0(%rax)
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 <HTS_mglsadff+0x220>
neg    %rdi
test   $0x1,%r11b
je     138e <HTS_mglsadff+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 <HTS_mglsadff+0x12a>
and    $0x3,%ecx
jmp    1202 <HTS_mglsadff+0x102>
xor    %edi,%edi
test   $0x1,%r11b
jne    1369 <HTS_mglsadff+0x269>
jmp    138e <HTS_mglsadff+0x28e>

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

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