void HTS_b2mc(const double *b, double *mc, int m, const double a) {
  double d, o;
  d = mc[m] = b[m];
  for (m--; m >= 0; m--) {
    o = b[m] + a * d;
    d = b[m];
    mc[m] = o;
  }
}


x86 assembly (clang, -O0):
<HTS_b2mc>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jl     119c <HTS_b2mc+0x9c>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x20(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    113c <HTS_b2mc+0x3c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<HTS_b2mc>:
movslq %edx,%rax
movsd  (%rdi,%rax,8),%xmm1
movsd  %xmm1,(%rsi,%rax,8)
test   %eax,%eax
jle    1193 <HTS_b2mc+0x93>
mov    %edx,%r8d
cmp    $0x4,%edx
jb     1162 <HTS_b2mc+0x62>
lea    -0x1(%r8),%r11
add    $0xffffffff,%edx
cmp    %r11d,%edx
jb     1162 <HTS_b2mc+0x62>
mov    %r11,%rax
shr    $0x20,%rax
jne    1162 <HTS_b2mc+0x62>
mov    %r8,%rax
neg    %rax
lea    (%rdx,%rax,1),%r10
add    $0x1,%r10
lea    (%rsi,%r10,8),%r9
lea    (%rdi,%rdx,8),%rax
add    $0x8,%rax
cmp    %rax,%r9
jae    1194 <HTS_b2mc+0x94>
lea    (%rsi,%rdx,8),%rax
add    $0x8,%rax
lea    (%rdi,%r10,8),%rdx
cmp    %rax,%rdx
jae    1194 <HTS_b2mc+0x94>
mov    %r8,%rax
add    $0x1,%rax
nopl   0x0(%rax)
lea    -0x2(%rax),%ecx
movsd  (%rdi,%rcx,8),%xmm2
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,(%rsi,%rcx,8)
add    $0xffffffffffffffff,%rax
movapd %xmm2,%xmm1
cmp    $0x1,%rax
ja     1170 <HTS_b2mc+0x70>
ret
mov    %r8d,%r9d
and    $0xfffffffc,%r9d
mov    %r8d,%eax
and    $0x3,%eax
movlhps %xmm1,%xmm4
movapd %xmm0,%xmm2
unpcklpd %xmm0,%xmm2
movapd %xmm2,%xmm3
shufpd $0x1,%xmm2,%xmm3
mov    %r9,%rdx
nopl   0x0(%rax,%rax,1)
mov    %r11d,%ecx
shl    $0x3,%rcx
movupd -0x18(%rdi,%rcx,1),%xmm1
movupd -0x8(%rdi,%rcx,1),%xmm5
movapd %xmm1,%xmm6
shufpd $0x1,%xmm1,%xmm6
movapd %xmm1,%xmm7
movsd  %xmm5,%xmm7
mulpd  %xmm2,%xmm7
addpd  %xmm6,%xmm7
movhlps %xmm5,%xmm4
mulpd  %xmm3,%xmm4
addpd  %xmm5,%xmm4
movupd %xmm4,-0x8(%rsi,%rcx,1)
shufpd $0x1,%xmm7,%xmm7
movupd %xmm7,-0x18(%rsi,%rcx,1)
add    $0xfffffffffffffffc,%r11
movapd %xmm6,%xmm4
add    $0xfffffffffffffffc,%rdx
jne    11c0 <HTS_b2mc+0xc0>
cmp    %r8,%r9
je     1193 <HTS_b2mc+0x93>
jmp    1165 <HTS_b2mc+0x65>

LLVM IR (-O0):
define void @HTS_b2mc(double* %0, double* %1, i32 %2, double %3) {
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca double
%9 = alloca double
%10 = alloca double
store double* %0, double** %5
store double* %1, double** %6
store i32 %2, i32* %7
store double %3, double* %8
%11 = load double*, double** %5
%12 = load i32, i32* %7
%13 = sext i32 %12 to i64
%14 = getelementptr inbounds double, double* %11, i64 %13
%15 = load double, double* %14
%16 = load double*, double** %6
%17 = load i32, i32* %7
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds double, double* %16, i64 %18
store double %15, double* %19
store double %15, double* %9
%20 = load i32, i32* %7
%21 = add nsw i32 %20, -1
store i32 %21, i32* %7
br label %22
22: ; preds = %44, %4
%23 = load i32, i32* %7
%24 = icmp sge i32 %23, 0
br i1 %24, label %25, label %47
25: ; preds = %22
%26 = load double*, double** %5
%27 = load i32, i32* %7
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds double, double* %26, i64 %28
%30 = load double, double* %29
%31 = load double, double* %8
%32 = load double, double* %9
%33 = call double @llvm.fmuladd.f64(double %31, double %32, double %30)
store double %33, double* %10
%34 = load double*, double** %5
%35 = load i32, i32* %7
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
%38 = load double, double* %37
store double %38, double* %9
%39 = load double, double* %10
%40 = load double*, double** %6
%41 = load i32, i32* %7
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %40, i64 %42
store double %39, double* %43
br label %44
44: ; preds = %25
%45 = load i32, i32* %7
%46 = add nsw i32 %45, -1
store i32 %46, i32* %7
br label %22
47: ; preds = %22
ret void
}

LLVM IR (-O2):
define void @HTS_b2mc(double* nocapture %0, double* nocapture writeonly %1, i32 %2, double %3) {
%5 = sext i32 %2 to i64
%6 = getelementptr inbounds double, double* %0, i64 %5
%7 = load double, double* %6
%8 = getelementptr inbounds double, double* %1, i64 %5
store double %7, double* %8
%9 = icmp sgt i32 %2, 0
br i1 %9, label %10, label %87
10: ; preds = %4
%11 = zext i32 %2 to i64
%12 = icmp ult i32 %2, 4
br i1 %12, label %74, label %13
13: ; preds = %10
%14 = add nsw i64 %11, -1
%15 = add i32 %2, -1
%16 = trunc i64 %14 to i32
%17 = icmp ult i32 %15, %16
%18 = icmp ugt i64 %14, 4294967295
%19 = or i1 %17, %18
br i1 %19, label %74, label %20
20: ; preds = %13
%21 = add i32 %2, -1
%22 = zext i32 %21 to i64
%23 = add nuw nsw i64 %22, 1
%24 = sub nsw i64 %23, %11
%25 = getelementptr double, double* %1, i64 %24
%26 = add nuw nsw i64 %22, 1
%27 = getelementptr double, double* %1, i64 %26
%28 = getelementptr double, double* %0, i64 %24
%29 = getelementptr double, double* %0, i64 %26
%30 = icmp ult double* %25, %29
%31 = icmp ult double* %28, %27
%32 = and i1 %30, %31
br i1 %32, label %74, label %33
33: ; preds = %20
%34 = and i64 %11, 4294967292
%35 = and i64 %11, 3
%36 = insertelement <2 x double> poison, double %7, i64 1
%37 = insertelement <2 x double> poison, double %3, i64 0
%38 = shufflevector <2 x double> %37, <2 x double> poison, <2 x i32> zeroinitializer
%39 = insertelement <2 x double> poison, double %3, i64 0
%40 = shufflevector <2 x double> %39, <2 x double> poison, <2 x i32> zeroinitializer
br label %41
41: ; preds = %41, %33
%42 = phi i64 [ 0, %33 ], [ %69, %41 ]
%43 = phi <2 x double> [ %36, %33 ], [ %56, %41 ]
%44 = xor i64 %42, -1
%45 = add i64 %44, %11
%46 = and i64 %45, 4294967295
%47 = getelementptr inbounds double, double* %0, i64 %46
%48 = getelementptr inbounds double, double* %47, i64 -1
%49 = bitcast double* %48 to <2 x double>*
%50 = load <2 x double>, <2 x double>* %49
%51 = shufflevector <2 x double> %50, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%52 = getelementptr inbounds double, double* %47, i64 -2
%53 = getelementptr inbounds double, double* %52, i64 -1
%54 = bitcast double* %53 to <2 x double>*
%55 = load <2 x double>, <2 x double>* %54
%56 = shufflevector <2 x double> %55, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%57 = shufflevector <2 x double> %43, <2 x double> %51, <2 x i32> <i32 1, i32 2>
%58 = shufflevector <2 x double> %51, <2 x double> %56, <2 x i32> <i32 1, i32 2>
%59 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %38, <2 x double> %57, <2 x double> %51)
%60 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %40, <2 x double> %58, <2 x double> %56)
%61 = getelementptr inbounds double, double* %1, i64 %46
%62 = shufflevector <2 x double> %59, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%63 = getelementptr inbounds double, double* %61, i64 -1
%64 = bitcast double* %63 to <2 x double>*
store <2 x double> %62, <2 x double>* %64
%65 = shufflevector <2 x double> %60, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%66 = getelementptr inbounds double, double* %61, i64 -2
%67 = getelementptr inbounds double, double* %66, i64 -1
%68 = bitcast double* %67 to <2 x double>*
store <2 x double> %65, <2 x double>* %68
%69 = add nuw i64 %42, 4
%70 = icmp eq i64 %69, %34
br i1 %70, label %71, label %41
71: ; preds = %41
%72 = icmp eq i64 %34, %11
%73 = extractelement <2 x double> %55, i64 0
br i1 %72, label %87, label %74
74: ; preds = %20, %13, %10, %71
%75 = phi i64 [ %35, %71 ], [ %11, %10 ], [ %11, %13 ], [ %11, %20 ]
%76 = phi double [ %73, %71 ], [ %7, %10 ], [ %7, %13 ], [ %7, %20 ]
br label %77
77: ; preds = %74, %77
%78 = phi i64 [ %80, %77 ], [ %75, %74 ]
%79 = phi double [ %83, %77 ], [ %76, %74 ]
%80 = add nsw i64 %78, -1
%81 = and i64 %80, 4294967295
%82 = getelementptr inbounds double, double* %0, i64 %81
%83 = load double, double* %82
%84 = tail call double @llvm.fmuladd.f64(double %3, double %79, double %83)
%85 = getelementptr inbounds double, double* %1, i64 %81
store double %84, double* %85
%86 = icmp ugt i64 %78, 1
br i1 %86, label %77, label %87
87: ; preds = %77, %71, %4
ret void
}