float cubic_interp(float p1, float p2, float p3, float p4, float amount) {
  float amount_sqrd = amount * amount;
  float amount_cubd = amount * amount * amount;
  float a1 = p4 - p3 - p2 + p1;
  float a2 = p1 - p2 - a1;
  float a3 = p3 - p1;
  float a4 = p2;
  return (a1 * amount_cubd) + (a2 * amount_sqrd) + (a3 * amount) + a4;
}


x86 assembly (clang, -O0):
<cubic_interp>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x4(%rbp)
movss  %xmm1,-0x8(%rbp)
movss  %xmm2,-0xc(%rbp)
movss  %xmm3,-0x10(%rbp)
movss  %xmm4,-0x14(%rbp)
movss  -0x14(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x14(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x10(%rbp),%xmm0
subss  -0xc(%rbp),%xmm0
subss  -0x8(%rbp),%xmm0
addss  -0x4(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x4(%rbp),%xmm0
subss  -0x8(%rbp),%xmm0
subss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0xc(%rbp),%xmm0
subss  -0x4(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x8(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x20(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm2
movss  -0x24(%rbp),%xmm0
mulss  -0x18(%rbp),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x28(%rbp),%xmm0
movss  -0x14(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
addss  -0x2c(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<cubic_interp>:
movaps %xmm4,%xmm5
mulss  %xmm4,%xmm5
subss  %xmm2,%xmm3
subss  %xmm1,%xmm3
addss  %xmm0,%xmm3
subss  %xmm0,%xmm2
subss  %xmm1,%xmm0
subss  %xmm3,%xmm0
mulss  %xmm5,%xmm0
mulss  %xmm4,%xmm5
mulss  %xmm3,%xmm5
addss  %xmm5,%xmm0
mulss  %xmm4,%xmm2
addss  %xmm2,%xmm0
addss  %xmm1,%xmm0
ret

LLVM IR (-O0):
define float @cubic_interp(float %0, float %1, float %2, float %3, float %4) {
%6 = alloca float
%7 = alloca float
%8 = alloca float
%9 = alloca float
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca float
%14 = alloca float
%15 = alloca float
%16 = alloca float
store float %0, float* %6
store float %1, float* %7
store float %2, float* %8
store float %3, float* %9
store float %4, float* %10
%17 = load float, float* %10
%18 = load float, float* %10
%19 = fmul float %17, %18
store float %19, float* %11
%20 = load float, float* %10
%21 = load float, float* %10
%22 = fmul float %20, %21
%23 = load float, float* %10
%24 = fmul float %22, %23
store float %24, float* %12
%25 = load float, float* %9
%26 = load float, float* %8
%27 = fsub float %25, %26
%28 = load float, float* %7
%29 = fsub float %27, %28
%30 = load float, float* %6
%31 = fadd float %29, %30
store float %31, float* %13
%32 = load float, float* %6
%33 = load float, float* %7
%34 = fsub float %32, %33
%35 = load float, float* %13
%36 = fsub float %34, %35
store float %36, float* %14
%37 = load float, float* %8
%38 = load float, float* %6
%39 = fsub float %37, %38
store float %39, float* %15
%40 = load float, float* %7
store float %40, float* %16
%41 = load float, float* %13
%42 = load float, float* %12
%43 = load float, float* %14
%44 = load float, float* %11
%45 = fmul float %43, %44
%46 = call float @llvm.fmuladd.f32(float %41, float %42, float %45)
%47 = load float, float* %15
%48 = load float, float* %10
%49 = call float @llvm.fmuladd.f32(float %47, float %48, float %46)
%50 = load float, float* %16
%51 = fadd float %49, %50
ret float %51
}

LLVM IR (-O2):
define float @cubic_interp(float %0, float %1, float %2, float %3, float %4) {
%6 = fmul float %4, %4
%7 = fmul float %6, %4
%8 = fsub float %3, %2
%9 = fsub float %8, %1
%10 = fadd float %9, %0
%11 = fsub float %0, %1
%12 = fsub float %11, %10
%13 = fsub float %2, %0
%14 = fmul float %6, %12
%15 = tail call float @llvm.fmuladd.f32(float %10, float %7, float %14)
%16 = tail call float @llvm.fmuladd.f32(float %13, float %4, float %15)
%17 = fadd float %16, %1
ret float %17
}