void cubic_coef(float frac, float interp[4]) {
  interp[0] = -0.16667f * frac + 0.16667f * frac * frac * frac;
  interp[1] = frac + 0.5f * frac * frac - 0.5f * frac * frac * frac;
  interp[3] =
      -0.33333f * frac + 0.5f * frac * frac - 0.16667f * frac * frac * frac;
  interp[2] = 1. - interp[0] - interp[1] - interp[3];
}


x86 assembly (clang, -O0):
<cubic_coef>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x4(%rbp)
mov    %rdi,-0x10(%rbp)
movss  -0x4(%rbp),%xmm0
movss  0xeee(%rip),%xmm1
mulss  -0x4(%rbp),%xmm1
mulss  -0x4(%rbp),%xmm1
mulss  -0x4(%rbp),%xmm1
movss  0xee3(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x4(%rbp),%xmm0
movss  0xec2(%rip),%xmm1
mulss  -0x4(%rbp),%xmm1
movss  -0x4(%rbp),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  0xea8(%rip),%xmm0
mulss  -0x4(%rbp),%xmm0
mulss  -0x4(%rbp),%xmm0
movss  -0x4(%rbp),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x4(%rax)
movss  -0x4(%rbp),%xmm1
movss  0xe6e(%rip),%xmm0
mulss  -0x4(%rbp),%xmm0
mulss  -0x4(%rbp),%xmm0
movss  0xe58(%rip),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  0xe44(%rip),%xmm0
mulss  -0x4(%rbp),%xmm0
mulss  -0x4(%rbp),%xmm0
movss  -0x4(%rbp),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0xc(%rax)
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdfb(%rip),%xmm0
subsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm1
cvtss2sd %xmm1,%xmm1
subsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  0xc(%rax),%xmm1
cvtss2sd %xmm1,%xmm1
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x8(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<cubic_coef>:
movss  0xef8(%rip),%xmm1
mulss  %xmm0,%xmm1
mulss  %xmm0,%xmm1
mulss  %xmm0,%xmm1
movss  0xee8(%rip),%xmm2
mulss  %xmm0,%xmm2
addss  %xmm1,%xmm2
movss  %xmm2,(%rdi)
movss  0xed8(%rip),%xmm3
mulss  %xmm0,%xmm3
mulss  %xmm0,%xmm3
movaps %xmm3,%xmm4
addss  %xmm0,%xmm4
movaps %xmm3,%xmm5
mulss  %xmm0,%xmm5
subss  %xmm5,%xmm4
movss  %xmm4,0x4(%rdi)
mulss  0xeb5(%rip),%xmm0
addss  %xmm3,%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,0xc(%rdi)
xorps  %xmm1,%xmm1
cvtss2sd %xmm2,%xmm1
movsd  0xe9d(%rip),%xmm2
subsd  %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtss2sd %xmm4,%xmm1
subsd  %xmm1,%xmm2
cvtss2sd %xmm0,%xmm0
subsd  %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtsd2ss %xmm2,%xmm0
movss  %xmm0,0x8(%rdi)
ret

LLVM IR (-O0):
define void @cubic_coef(float %0, float* %1) {
%3 = alloca float
%4 = alloca float*
store float %0, float* %3
store float* %1, float** %4
%5 = load float, float* %3
%6 = load float, float* %3
%7 = fmul float 0x3FC5557140000000, %6
%8 = load float, float* %3
%9 = fmul float %7, %8
%10 = load float, float* %3
%11 = fmul float %9, %10
%12 = call float @llvm.fmuladd.f32(float 0xBFC5557140000000, float %5, float %11)
%13 = load float*, float** %4
%14 = getelementptr inbounds float, float* %13, i64 0
store float %12, float* %14
%15 = load float, float* %3
%16 = load float, float* %3
%17 = fmul float 5.000000e-01, %16
%18 = load float, float* %3
%19 = call float @llvm.fmuladd.f32(float %17, float %18, float %15)
%20 = load float, float* %3
%21 = fmul float 5.000000e-01, %20
%22 = load float, float* %3
%23 = fmul float %21, %22
%24 = load float, float* %3
%25 = fneg float %23
%26 = call float @llvm.fmuladd.f32(float %25, float %24, float %19)
%27 = load float*, float** %4
%28 = getelementptr inbounds float, float* %27, i64 1
store float %26, float* %28
%29 = load float, float* %3
%30 = load float, float* %3
%31 = fmul float 5.000000e-01, %30
%32 = load float, float* %3
%33 = fmul float %31, %32
%34 = call float @llvm.fmuladd.f32(float 0xBFD5554760000000, float %29, float %33)
%35 = load float, float* %3
%36 = fmul float 0x3FC5557140000000, %35
%37 = load float, float* %3
%38 = fmul float %36, %37
%39 = load float, float* %3
%40 = fneg float %38
%41 = call float @llvm.fmuladd.f32(float %40, float %39, float %34)
%42 = load float*, float** %4
%43 = getelementptr inbounds float, float* %42, i64 3
store float %41, float* %43
%44 = load float*, float** %4
%45 = getelementptr inbounds float, float* %44, i64 0
%46 = load float, float* %45
%47 = fpext float %46 to double
%48 = fsub double 1.000000e+00, %47
%49 = load float*, float** %4
%50 = getelementptr inbounds float, float* %49, i64 1
%51 = load float, float* %50
%52 = fpext float %51 to double
%53 = fsub double %48, %52
%54 = load float*, float** %4
%55 = getelementptr inbounds float, float* %54, i64 3
%56 = load float, float* %55
%57 = fpext float %56 to double
%58 = fsub double %53, %57
%59 = fptrunc double %58 to float
%60 = load float*, float** %4
%61 = getelementptr inbounds float, float* %60, i64 2
store float %59, float* %61
ret void
}

LLVM IR (-O2):
define void @cubic_coef(float %0, float* nocapture %1) {
%3 = fmul float %0, 0x3FC5557140000000
%4 = fmul float %3, %0
%5 = fmul float %4, %0
%6 = tail call float @llvm.fmuladd.f32(float %0, float 0xBFC5557140000000, float %5)
store float %6, float* %1
%7 = fmul float %0, 5.000000e-01
%8 = tail call float @llvm.fmuladd.f32(float %7, float %0, float %0)
%9 = fmul float %7, %0
%10 = fneg float %9
%11 = tail call float @llvm.fmuladd.f32(float %10, float %0, float %8)
%12 = getelementptr inbounds float, float* %1, i64 1
store float %11, float* %12
%13 = tail call float @llvm.fmuladd.f32(float %0, float 0xBFD5554760000000, float %9)
%14 = fneg float %4
%15 = tail call float @llvm.fmuladd.f32(float %14, float %0, float %13)
%16 = getelementptr inbounds float, float* %1, i64 3
store float %15, float* %16
%17 = fpext float %6 to double
%18 = fsub double 1.000000e+00, %17
%19 = fpext float %11 to double
%20 = fsub double %18, %19
%21 = fpext float %15 to double
%22 = fsub double %20, %21
%23 = fptrunc double %22 to float
%24 = getelementptr inbounds float, float* %1, i64 2
store float %23, float* %24
ret void
}