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


x86 assembly (clang, -O0):
<imglsadff>:
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)
movsd  -0x20(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xeba(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%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    11fa <imglsadff+0xfa>
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    116c <imglsadff+0x6c>
movsd  -0x8(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jle    1246 <imglsadff+0x146>
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    120f <imglsadff+0x10f>
movsd  -0x20(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  -0x8(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x38(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<imglsadff>:
push   %r14
push   %rbx
mov    %rdx,%r8
movsd  (%rdx),%xmm8
movsd  0x8(%rdi),%xmm2
mulsd  %xmm8,%xmm2
mov    %esi,%ecx
cmp    $0x2,%esi
jl     11f5 <imglsadff+0xf5>
movsd  0x8(%r8),%xmm4
lea    -0x1(%rcx),%r9
cmp    $0x2,%esi
jne    113e <imglsadff+0x3e>
mov    $0x1,%edx
movapd %xmm8,%xmm5
jmp    11cc <imglsadff+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),%xmm3
movsd  0x10(%r8,%rax,8),%xmm4
movapd %xmm3,%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  %xmm2,%xmm7
lea    0x2(%rax),%rdx
mulsd  %xmm1,%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,0x8(%r8,%rax,8)
movsd  0x10(%rdi,%rax,8),%xmm2
mulsd  %xmm5,%xmm2
addsd  %xmm7,%xmm2
lea    (%r10,%rax,1),%rbx
add    $0x2,%rbx
movapd %xmm4,%xmm6
mov    %rdx,%rax
cmp    $0x1,%rbx
jne    1160 <imglsadff+0x60>
test   $0x1,%r9b
je     11f5 <imglsadff+0xf5>
movsd  0x8(%r8,%rdx,8),%xmm6
subsd  %xmm5,%xmm6
mulsd  %xmm1,%xmm6
addsd  %xmm4,%xmm6
movsd  %xmm6,(%r8,%rdx,8)
mulsd  0x8(%rdi,%rdx,8),%xmm6
addsd  %xmm6,%xmm2
test   %esi,%esi
jle    122a <imglsadff+0x12a>
cmp    $0x4,%esi
jae    125d <imglsadff+0x15d>
add    $0x1,%rcx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    -0x2(%rcx),%eax
movsd  (%r8,%rax,8),%xmm3
movsd  %xmm3,-0x8(%r8,%rcx,8)
add    $0xffffffffffffffff,%rcx
cmp    $0x1,%rcx
ja     1210 <imglsadff+0x110>
addsd  %xmm0,%xmm2
movapd %xmm1,%xmm3
mulsd  %xmm1,%xmm3
movsd  0xdc2(%rip),%xmm4
subsd  %xmm3,%xmm4
mulsd  %xmm0,%xmm4
mulsd  %xmm1,%xmm8
addsd  %xmm4,%xmm8
movsd  %xmm8,(%r8)
movapd %xmm2,%xmm0
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     11fe <imglsadff+0xfe>
mov    %r14,%rdi
shr    $0x20,%rdi
jne    11fe <imglsadff+0xfe>
lea    (%r8,%rcx,8),%rbx
mov    %rbx,%rdi
sub    %rax,%rdi
cmp    %rbx,%rdi
ja     11fe <imglsadff+0xfe>
test   %dl,%dl
jne    11fe <imglsadff+0xfe>
lea    (%r8,%rsi,8),%rbx
mov    %rbx,%rdi
sub    %rax,%rdi
cmp    %rbx,%rdi
ja     11fe <imglsadff+0xfe>
test   %dl,%dl
jne    11fe <imglsadff+0xfe>
lea    0x8(%r8),%rax
lea    (%r8,%rsi,8),%rdx
add    $0x8,%rdx
cmp    %rdx,%rax
jae    12e7 <imglsadff+0x1e7>
lea    (%r8,%rcx,8),%rax
add    $0x8,%rax
add    $0x1,%rsi
sub    %rcx,%rsi
lea    (%r8,%rsi,8),%rdx
cmp    %rax,%rdx
jb     11fe <imglsadff+0xfe>
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     13af <imglsadff+0x2af>
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),%xmm4
movups -0x8(%r8,%rdx,8),%xmm5
movups %xmm5,(%rax,%rdi,8)
movups %xmm4,-0x10(%rax,%rdi,8)
lea    (%r14,%rdi,1),%edx
movupd -0x18(%r8,%rdx,8),%xmm4
movupd -0x8(%r8,%rdx,8),%xmm5
movupd %xmm5,-0x20(%rax,%rdi,8)
movupd %xmm4,-0x30(%rax,%rdi,8)
add    $0xfffffffffffffff8,%rdi
add    $0xfffffffffffffffe,%rsi
jne    1330 <imglsadff+0x230>
neg    %rdi
test   $0x1,%r11b
je     139e <imglsadff+0x29e>
mov    %rcx,%rax
sub    %rdi,%rax
add    %eax,%r9d
movupd -0x18(%r8,%r9,8),%xmm4
movupd -0x8(%r8,%r9,8),%xmm5
movupd %xmm5,-0x8(%r8,%rax,8)
movupd %xmm4,-0x18(%r8,%rax,8)
cmp    %rcx,%r10
je     122a <imglsadff+0x12a>
and    $0x3,%ecx
jmp    11fe <imglsadff+0xfe>
xor    %edi,%edi
test   $0x1,%r11b
jne    1379 <imglsadff+0x279>
jmp    139e <imglsadff+0x29e>

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

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