float cubic(float t, float pl, float p0, float p1, float ph) {
  float m0 = (p1 - pl) / 2;
  float m1 = (ph - p0) / 2;
  return (1 + 2 * t) * (1 - t) * (1 - t) * p0 + t * (1 - t) * (1 - t) * m0 +
         t * t * (3 - 2 * t) * p1 + t * t * (t - 1) * m1;
}


x86 assembly (clang, -O0):
<cubic>:
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  -0x10(%rbp),%xmm0
subss  -0x8(%rbp),%xmm0
movss  0xed9(%rip),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x14(%rbp),%xmm0
subss  -0xc(%rbp),%xmm0
movss  0xebe(%rip),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x4(%rbp),%xmm0
addss  %xmm0,%xmm0
movss  0xe9c(%rip),%xmm1
addss  %xmm1,%xmm0
movss  0xe90(%rip),%xmm1
subss  -0x4(%rbp),%xmm1
mulss  %xmm1,%xmm0
movss  0xe7f(%rip),%xmm1
subss  -0x4(%rbp),%xmm1
mulss  %xmm1,%xmm0
movss  -0xc(%rbp),%xmm2
movss  -0x4(%rbp),%xmm1
movss  0xe64(%rip),%xmm3
subss  -0x4(%rbp),%xmm3
mulss  %xmm3,%xmm1
movss  0xe53(%rip),%xmm3
subss  -0x4(%rbp),%xmm3
mulss  %xmm3,%xmm1
mulss  -0x18(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x4(%rbp),%xmm1
mulss  -0x4(%rbp),%xmm1
movss  -0x4(%rbp),%xmm3
addss  %xmm3,%xmm3
movss  0xe26(%rip),%xmm2
subss  %xmm3,%xmm2
mulss  %xmm2,%xmm1
movss  -0x10(%rbp),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x4(%rbp),%xmm0
mulss  -0x4(%rbp),%xmm0
movss  -0x4(%rbp),%xmm2
movss  0xdf6(%rip),%xmm3
subss  %xmm3,%xmm2
mulss  %xmm2,%xmm0
movss  -0x1c(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<cubic>:
movaps %xmm0,%xmm5
addss  %xmm0,%xmm5
movaps %xmm0,%xmm8
mulss  %xmm0,%xmm8
movss  0xef0(%rip),%xmm7
subss  %xmm5,%xmm7
mulss  %xmm8,%xmm7
mulss  %xmm3,%xmm7
subss  %xmm1,%xmm3
movss  0xecf(%rip),%xmm1
mulss  %xmm1,%xmm3
subss  %xmm2,%xmm4
mulss  %xmm4,%xmm1
movss  0xebf(%rip),%xmm6
addss  %xmm6,%xmm5
subss  %xmm0,%xmm6
mulss  %xmm6,%xmm5
mulss  %xmm6,%xmm5
movaps %xmm6,%xmm4
mulss  %xmm0,%xmm4
mulss  %xmm6,%xmm4
mulss  %xmm3,%xmm4
mulss  %xmm2,%xmm5
addss  %xmm4,%xmm5
addss  0xe98(%rip),%xmm0
addss  %xmm5,%xmm7
mulss  %xmm8,%xmm0
mulss  %xmm1,%xmm0
addss  %xmm7,%xmm0
ret

LLVM IR (-O0):
define float @cubic(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
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
%13 = load float, float* %9
%14 = load float, float* %7
%15 = fsub float %13, %14
%16 = fdiv float %15, 2.000000e+00
store float %16, float* %11
%17 = load float, float* %10
%18 = load float, float* %8
%19 = fsub float %17, %18
%20 = fdiv float %19, 2.000000e+00
store float %20, float* %12
%21 = load float, float* %6
%22 = call float @llvm.fmuladd.f32(float 2.000000e+00, float %21, float 1.000000e+00)
%23 = load float, float* %6
%24 = fsub float 1.000000e+00, %23
%25 = fmul float %22, %24
%26 = load float, float* %6
%27 = fsub float 1.000000e+00, %26
%28 = fmul float %25, %27
%29 = load float, float* %8
%30 = load float, float* %6
%31 = load float, float* %6
%32 = fsub float 1.000000e+00, %31
%33 = fmul float %30, %32
%34 = load float, float* %6
%35 = fsub float 1.000000e+00, %34
%36 = fmul float %33, %35
%37 = load float, float* %11
%38 = fmul float %36, %37
%39 = call float @llvm.fmuladd.f32(float %28, float %29, float %38)
%40 = load float, float* %6
%41 = load float, float* %6
%42 = fmul float %40, %41
%43 = load float, float* %6
%44 = call float @llvm.fmuladd.f32(float -2.000000e+00, float %43, float 3.000000e+00)
%45 = fmul float %42, %44
%46 = load float, float* %9
%47 = call float @llvm.fmuladd.f32(float %45, float %46, float %39)
%48 = load float, float* %6
%49 = load float, float* %6
%50 = fmul float %48, %49
%51 = load float, float* %6
%52 = fsub float %51, 1.000000e+00
%53 = fmul float %50, %52
%54 = load float, float* %12
%55 = call float @llvm.fmuladd.f32(float %53, float %54, float %47)
ret float %55
}

LLVM IR (-O2):
define float @cubic(float %0, float %1, float %2, float %3, float %4) {
%6 = fsub float %3, %1
%7 = fmul float %6, 5.000000e-01
%8 = fsub float %4, %2
%9 = fmul float %8, 5.000000e-01
%10 = tail call float @llvm.fmuladd.f32(float %0, float 2.000000e+00, float 1.000000e+00)
%11 = fsub float 1.000000e+00, %0
%12 = fmul float %10, %11
%13 = fmul float %11, %12
%14 = fmul float %11, %0
%15 = fmul float %11, %14
%16 = fmul float %15, %7
%17 = tail call float @llvm.fmuladd.f32(float %13, float %2, float %16)
%18 = fmul float %0, %0
%19 = tail call float @llvm.fmuladd.f32(float %0, float -2.000000e+00, float 3.000000e+00)
%20 = fmul float %18, %19
%21 = tail call float @llvm.fmuladd.f32(float %20, float %3, float %17)
%22 = fadd float %0, -1.000000e+00
%23 = fmul float %18, %22
%24 = tail call float @llvm.fmuladd.f32(float %23, float %9, float %21)
ret float %24
}