float hermite_kernel(float mu1, float a, float b, float c, float d) {
  float mu2, mu3, m0, m1, a0, a1, a2, a3;
  mu2 = mu1 * mu1;
  mu3 = mu2 * mu1;
  m0 = (c - a) * 0.5f;
  m1 = (d - b) * 0.5f;
  a0 = +2 * mu3 - 3 * mu2 + 1;
  a1 = mu3 - 2 * mu2 + mu1;
  a2 = mu3 - mu2;
  a3 = -2 * mu3 + 3 * mu2;
  return (a0 * b) + (a1 * m0) + (a2 * m1) + (a3 * c);
}


x86 assembly (clang, -O0):
<hermite_kernel>:
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  -0x4(%rbp),%xmm0
mulss  -0x4(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x18(%rbp),%xmm0
mulss  -0x4(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x10(%rbp),%xmm0
subss  -0x8(%rbp),%xmm0
movss  0xebb(%rip),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x14(%rbp),%xmm0
subss  -0xc(%rbp),%xmm0
movss  0xea0(%rip),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x1c(%rbp),%xmm0
movss  0xe82(%rip),%xmm1
mulss  -0x18(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
addss  %xmm0,%xmm0
addss  %xmm1,%xmm0
movss  0xe64(%rip),%xmm1
addss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x1c(%rbp),%xmm0
movss  -0x18(%rbp),%xmm1
addss  %xmm1,%xmm1
subss  %xmm1,%xmm0
addss  -0x4(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x1c(%rbp),%xmm0
subss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x1c(%rbp),%xmm1
movss  0xe1f(%rip),%xmm0
mulss  -0x18(%rbp),%xmm0
addss  %xmm1,%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x28(%rbp),%xmm0
movss  -0xc(%rbp),%xmm2
movss  -0x2c(%rbp),%xmm1
mulss  -0x20(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x30(%rbp),%xmm1
movss  -0x24(%rbp),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x34(%rbp),%xmm0
movss  -0x10(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<hermite_kernel>:
movaps %xmm0,%xmm8
mulss  %xmm0,%xmm8
movaps %xmm8,%xmm5
mulss  %xmm0,%xmm5
movss  0xeeb(%rip),%xmm7
mulss  %xmm8,%xmm7
movaps %xmm5,%xmm9
addss  %xmm5,%xmm9
movaps %xmm9,%xmm6
subss  %xmm7,%xmm6
subss  %xmm9,%xmm7
mulss  %xmm3,%xmm7
subss  %xmm1,%xmm3
movss  0xebc(%rip),%xmm1
mulss  %xmm1,%xmm3
subss  %xmm2,%xmm4
mulss  %xmm4,%xmm1
addss  0xeb0(%rip),%xmm6
movaps %xmm5,%xmm4
subss  %xmm8,%xmm5
addss  %xmm8,%xmm8
subss  %xmm8,%xmm4
addss  %xmm4,%xmm0
mulss  %xmm3,%xmm0
mulss  %xmm1,%xmm5
mulss  %xmm2,%xmm6
addss  %xmm6,%xmm0
addss  %xmm5,%xmm0
addss  %xmm7,%xmm0
ret

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

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