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


x86 assembly (clang, -O0):
<lmafirt>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x8(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x1,-0x24(%rbp)
jle    1195 <lmafirt+0x95>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x8(%rbp),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1151 <lmafirt+0x51>
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
addsd  0x8(%rax),%xmm0
movsd  %xmm0,0x8(%rax)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11e6 <lmafirt+0xe6>
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11ad <lmafirt+0xad>
movsd  -0x30(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<lmafirt>:
push   %rbx
sub    $0x10,%rsp
movsd  (%rdx),%xmm6
movslq %esi,%rax
movsd  (%rdi,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rax,8)
mov    %esi,%r8d
cmp    $0x3,%eax
jl     11e0 <lmafirt+0xd0>
xor    %r9d,%r9d
mov    %r8,%rcx
sub    $0x3,%rcx
mov    $0x0,%r10d
cmovae %rcx,%r10
add    $0x1,%r10
mov    %r8,%rax
cmp    $0x4,%r10
jb     11b5 <lmafirt+0xa5>
sub    $0x3,%rax
mov    $0x0,%ecx
cmovae %rax,%rcx
mov    %r8d,%eax
sub    $0x3,%eax
cmovb  %r9d,%eax
lea    -0x1(%rsi),%r9d
cmp    %eax,%r9d
jb     11b2 <lmafirt+0xa2>
shr    $0x20,%rcx
jne    11b2 <lmafirt+0xa2>
cmp    $0x3,%r8
mov    $0x3,%ecx
cmovb  %r8,%rcx
add    %r9,%rcx
sub    %r8,%rcx
lea    (%rdx,%rcx,8),%r11
lea    (%rdi,%r9,8),%rax
add    $0x8,%rax
cmp    %rax,%r11
jae    121a <lmafirt+0x10a>
lea    (%rdx,%r9,8),%rax
add    $0x8,%rax
lea    (%rdi,%rcx,8),%rcx
cmp    %rax,%rcx
jae    121a <lmafirt+0x10a>
mov    %r8,%rax
add    $0x1,%rax
nopl   0x0(%rax)
lea    -0x2(%rax),%ecx
movsd  (%rdi,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdx,%rcx,8),%xmm1
movsd  %xmm1,(%rdx,%rcx,8)
add    $0xffffffffffffffff,%rax
cmp    $0x3,%rax
ja     11c0 <lmafirt+0xb0>
xorpd  %xmm0,%xmm0
addsd  0x8(%rdx),%xmm0
movsd  %xmm0,0x8(%rdx)
test   %esi,%esi
jle    1211 <lmafirt+0x101>
lea    0x8(%rdx),%rsi
shl    $0x3,%r8
mov    %rdx,%rdi
mov    %r8,%rdx
movsd  %xmm6,0x8(%rsp)
call   1030 <memmove@plt>
movsd  0x8(%rsp),%xmm6
movaps %xmm6,%xmm0
add    $0x10,%rsp
pop    %rbx
ret
mov    %r10,%r9
and    $0xfffffffffffffffc,%r9
mov    %r8,%rax
sub    %r9,%rax
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x1(%r8),%rcx
shufpd $0x1,%xmm1,%xmm1
mov    %r9,%r11
nopl   0x0(%rax,%rax,1)
mov    %ecx,%ebx
shl    $0x3,%rbx
movupd -0x18(%rdi,%rbx,1),%xmm2
movupd -0x8(%rdi,%rbx,1),%xmm3
movupd -0x18(%rdx,%rbx,1),%xmm4
movupd -0x8(%rdx,%rbx,1),%xmm5
mulpd  %xmm1,%xmm3
addpd  %xmm5,%xmm3
movupd %xmm3,-0x8(%rdx,%rbx,1)
mulpd  %xmm1,%xmm2
addpd  %xmm4,%xmm2
movupd %xmm2,-0x18(%rdx,%rbx,1)
add    $0xfffffffffffffffc,%rcx
add    $0xfffffffffffffffc,%r11
jne    1240 <lmafirt+0x130>
cmp    %r9,%r10
je     11e0 <lmafirt+0xd0>
jmp    11b5 <lmafirt+0xa5>

LLVM IR (-O0):
define double @lmafirt(double %0, double* %1, i32 %2, double* %3) {
%5 = alloca double
%6 = alloca double*
%7 = alloca i32
%8 = alloca double*
%9 = alloca i32
%10 = alloca double
store double %0, double* %5
store double* %1, double** %6
store i32 %2, i32* %7
store double* %3, double** %8
store double 0.000000e+00, double* %10
%11 = load double*, double** %8
%12 = getelementptr inbounds double, double* %11, i64 0
%13 = load double, double* %12
store double %13, double* %10
%14 = load double*, double** %6
%15 = load i32, i32* %7
%16 = sext i32 %15 to i64
%17 = getelementptr inbounds double, double* %14, i64 %16
%18 = load double, double* %17
%19 = load double, double* %5
%20 = fmul double %18, %19
%21 = load double*, double** %8
%22 = load i32, i32* %7
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %21, i64 %23
store double %20, double* %24
%25 = load i32, i32* %7
%26 = sub nsw i32 %25, 1
store i32 %26, i32* %9
br label %27
27: ; preds = %43, %4
%28 = load i32, i32* %9
%29 = icmp sgt i32 %28, 1
br i1 %29, label %30, label %46
30: ; preds = %27
%31 = load double*, double** %6
%32 = load i32, i32* %9
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double, double* %31, i64 %33
%35 = load double, double* %34
%36 = load double, double* %5
%37 = load double*, double** %8
%38 = load i32, i32* %9
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double, double* %37, i64 %39
%41 = load double, double* %40
%42 = call double @llvm.fmuladd.f64(double %35, double %36, double %41)
store double %42, double* %40
br label %43
43: ; preds = %30
%44 = load i32, i32* %9
%45 = add nsw i32 %44, -1
store i32 %45, i32* %9
br label %27
46: ; preds = %27
%47 = load double*, double** %8
%48 = getelementptr inbounds double, double* %47, i64 1
%49 = load double, double* %48
%50 = fadd double %49, 0.000000e+00
store double %50, double* %48
store i32 0, i32* %9
br label %51
51: ; preds = %66, %46
%52 = load i32, i32* %9
%53 = load i32, i32* %7
%54 = icmp slt i32 %52, %53
br i1 %54, label %55, label %69
55: ; preds = %51
%56 = load double*, double** %8
%57 = load i32, i32* %9
%58 = add nsw i32 %57, 1
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %56, i64 %59
%61 = load double, double* %60
%62 = load double*, double** %8
%63 = load i32, i32* %9
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds double, double* %62, i64 %64
store double %61, double* %65
br label %66
66: ; preds = %55
%67 = load i32, i32* %9
%68 = add nsw i32 %67, 1
store i32 %68, i32* %9
br label %51
69: ; preds = %51
%70 = load double, double* %10
ret double %70
}

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