double intpol(diffn, coeffs, k)
double *diffn;
double *coeffs;
int k;
{
  double s;
  int i;
  s = 0.0;
  coeffs += k;
  diffn += k;
  for (i = 0; i < k; i++) {
    s += (*--coeffs) * (*--diffn);
  }
  return (s);
}


x86 assembly (clang, -O0):
<intpol>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%ecx
mov    -0x10(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x10(%rbp)
mov    -0x14(%rbp),%ecx
mov    -0x8(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    119c <intpol+0x9c>
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0xfffffffffffffff8,%rcx
mov    %rcx,-0x10(%rbp)
movsd  -0x8(%rax),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0xfffffffffffffff8,%rcx
mov    %rcx,-0x8(%rbp)
movsd  -0x8(%rax),%xmm2
movsd  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1148 <intpol+0x48>
movsd  -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<intpol>:
test   %edx,%edx
jle    112b <intpol+0x2b>
movslq %edx,%rcx
lea    -0x1(%rdx),%r8d
mov    %edx,%eax
and    $0x3,%eax
cmp    $0x3,%r8d
jae    112f <intpol+0x2f>
lea    (%rdi,%rcx,8),%rdx
lea    (%rsi,%rcx,8),%rcx
xorpd  %xmm0,%xmm0
test   %eax,%eax
jne    1199 <intpol+0x99>
jmp    11b8 <intpol+0xb8>
xorps  %xmm0,%xmm0
ret
and    $0xfffffffc,%edx
xorpd  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
movsd  -0x8(%rsi,%rcx,8),%xmm1
mulsd  -0x8(%rdi,%rcx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  -0x10(%rsi,%rcx,8),%xmm2
mulsd  -0x10(%rdi,%rcx,8),%xmm2
addsd  %xmm1,%xmm2
movsd  -0x18(%rsi,%rcx,8),%xmm1
mulsd  -0x18(%rdi,%rcx,8),%xmm1
movsd  -0x20(%rsi,%rcx,8),%xmm0
addsd  %xmm2,%xmm1
mulsd  -0x20(%rdi,%rcx,8),%xmm0
addsd  %xmm1,%xmm0
add    $0xffffffffffffffe0,%rdi
add    $0xffffffffffffffe0,%rsi
add    $0xfffffffc,%edx
jne    1140 <intpol+0x40>
lea    (%rdi,%rcx,8),%rdx
lea    (%rsi,%rcx,8),%rcx
test   %eax,%eax
je     11b8 <intpol+0xb8>
neg    %eax
xor    %esi,%esi
nopl   (%rax)
movsd  -0x8(%rcx,%rsi,8),%xmm1
mulsd  -0x8(%rdx,%rsi,8),%xmm1
addsd  %xmm1,%xmm0
add    $0xffffffffffffffff,%rsi
cmp    %esi,%eax
jne    11a0 <intpol+0xa0>
ret

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

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