double _XcmsPolynomial(order, coeffs, x)
register int order;
double const *coeffs;
double x;
{
  auto double rtn_value;
  coeffs += order;
  rtn_value = *coeffs--;
  while (order-- > 0)
    rtn_value = *coeffs-- + (x * rtn_value);
  return (rtn_value);
}


x86 assembly (clang, -O0):
<_XcmsPolynomial>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    -0x4(%rbp),%ecx
mov    -0x10(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x10(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0xfffffffffffffff8,%rcx
mov    %rcx,-0x10(%rbp)
movsd  (%rax),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x4(%rbp)
cmp    $0x0,%eax
jle    1180 <_XcmsPolynomial+0x80>
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0xfffffffffffffff8,%rcx
mov    %rcx,-0x10(%rbp)
movsd  (%rax),%xmm1
movsd  -0x18(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    113d <_XcmsPolynomial+0x3d>
movsd  -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<_XcmsPolynomial>:
movapd %xmm0,%xmm1
movslq %edi,%rdx
movsd  (%rsi,%rdx,8),%xmm0
test   %edx,%edx
jle    117f <_XcmsPolynomial+0x7f>
lea    (%rsi,%rdx,8),%rax
lea    -0x1(%rdx),%ecx
and    $0x3,%edx
je     1136 <_XcmsPolynomial+0x36>
xor    %esi,%esi
xchg   %ax,%ax
mulsd  %xmm1,%xmm0
addsd  -0x8(%rax),%xmm0
add    $0xfffffffffffffff8,%rax
add    $0x1,%esi
cmp    %esi,%edx
jne    1120 <_XcmsPolynomial+0x20>
sub    %esi,%edi
cmp    $0x3,%ecx
jb     117f <_XcmsPolynomial+0x7f>
add    $0xfffffffffffffff8,%rax
add    $0xffffffff,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mulsd  %xmm1,%xmm0
addsd  (%rax),%xmm0
mulsd  %xmm1,%xmm0
addsd  -0x8(%rax),%xmm0
mulsd  %xmm1,%xmm0
addsd  -0x10(%rax),%xmm0
mulsd  %xmm1,%xmm0
addsd  -0x18(%rax),%xmm0
add    $0xffffffffffffffe0,%rax
add    $0xfffffffc,%edi
cmp    $0xfffffffe,%edi
jb     1150 <_XcmsPolynomial+0x50>
ret

LLVM IR (-O0):
define double @_XcmsPolynomial(i32 %0, double* %1, double %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca double
%7 = alloca double
store i32 %0, i32* %4
store double* %1, double** %5
store double %2, double* %6
%8 = load i32, i32* %4
%9 = load double*, double** %5
%10 = sext i32 %8 to i64
%11 = getelementptr inbounds double, double* %9, i64 %10
store double* %11, double** %5
%12 = load double*, double** %5
%13 = getelementptr inbounds double, double* %12, i32 -1
store double* %13, double** %5
%14 = load double, double* %12
store double %14, double* %7
br label %15
15: ; preds = %19, %3
%16 = load i32, i32* %4
%17 = add nsw i32 %16, -1
store i32 %17, i32* %4
%18 = icmp sgt i32 %16, 0
br i1 %18, label %19, label %26
19: ; preds = %15
%20 = load double*, double** %5
%21 = getelementptr inbounds double, double* %20, i32 -1
store double* %21, double** %5
%22 = load double, double* %20
%23 = load double, double* %6
%24 = load double, double* %7
%25 = call double @llvm.fmuladd.f64(double %23, double %24, double %22)
store double %25, double* %7
br label %15
26: ; preds = %15
%27 = load double, double* %7
ret double %27
}

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