double imglsadff1(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];
  }
  y = g * 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):
<imglsadff1>:
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 <imglsadff1+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 <imglsadff1+0x71>
movsd  -0x28(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  -0x8(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jle    1258 <imglsadff1+0x158>
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    1221 <imglsadff1+0x121>
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  -0x40(%rbp),%xmm0
pop    %rbp
ret

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

LLVM IR (-O0):
define double @imglsadff1(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 = call double @llvm.fmuladd.f64(double %69, double %70, double %71)
store double %72, double* %14
%73 = load i32, i32* %9
store i32 %73, i32* %13
br label %74
74: ; preds = %88, %68
%75 = load i32, i32* %13
%76 = icmp sgt i32 %75, 0
br i1 %76, label %77, label %91
77: ; preds = %74
%78 = load double*, double** %12
%79 = load i32, i32* %13
%80 = sub nsw i32 %79, 1
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %78, i64 %81
%83 = load double, double* %82
%84 = load double*, double** %12
%85 = load i32, i32* %13
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
store double %83, double* %87
br label %88
88: ; preds = %77
%89 = load i32, i32* %13
%90 = add nsw i32 %89, -1
store i32 %90, i32* %13
br label %74
91: ; preds = %74
%92 = load double, double* %10
%93 = load double*, double** %12
%94 = getelementptr inbounds double, double* %93, i64 0
%95 = load double, double* %94
%96 = load double, double* %15
%97 = load double, double* %7
%98 = fmul double %96, %97
%99 = call double @llvm.fmuladd.f64(double %92, double %95, double %98)
%100 = load double*, double** %12
%101 = getelementptr inbounds double, double* %100, i64 0
store double %99, double* %101
%102 = load double, double* %14
ret double %102
}

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