double astro_polynomial(double t, int termc, const double *termv) {
  double T = 1.0;
  double v = 0.0;
  int i;
  for (i = 0; i < termc; i++) {
    v += *termv++ * T;
    T *= t;
  }
  return v;
}


x86 assembly (clang, -O0):
<astro_polynomial>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
movsd  0xee8(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    117f <astro_polynomial+0x7f>
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x18(%rbp)
movsd  (%rax),%xmm0
movsd  -0x20(%rbp),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x8(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    112c <astro_polynomial+0x2c>
movsd  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<astro_polynomial>:
test   %edi,%edi
jle    111a <astro_polynomial+0x1a>
lea    -0x1(%rdi),%eax
cmp    $0x3,%eax
jae    111e <astro_polynomial+0x1e>
movsd  0xeec(%rip),%xmm2
xorpd  %xmm1,%xmm1
jmp    117c <astro_polynomial+0x7c>
xorps  %xmm0,%xmm0
ret
mov    %edi,%eax
and    $0xfffffffc,%eax
movsd  0xed5(%rip),%xmm2
xorpd  %xmm1,%xmm1
nop
movsd  (%rsi),%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm1,%xmm3
mulsd  %xmm0,%xmm2
movsd  0x8(%rsi),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm3,%xmm1
mulsd  %xmm0,%xmm2
movsd  0x10(%rsi),%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm1,%xmm3
mulsd  %xmm0,%xmm2
movsd  0x18(%rsi),%xmm1
mulsd  %xmm2,%xmm1
add    $0x20,%rsi
addsd  %xmm3,%xmm1
mulsd  %xmm0,%xmm2
add    $0xfffffffc,%eax
jne    1130 <astro_polynomial+0x30>
test   $0x3,%dil
je     11a9 <astro_polynomial+0xa9>
and    $0x3,%edi
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movsd  (%rsi,%rax,8),%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm3,%xmm1
mulsd  %xmm0,%xmm2
add    $0x1,%rax
cmp    %eax,%edi
jne    1190 <astro_polynomial+0x90>
movapd %xmm1,%xmm0
ret

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

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