double cubic(double dx, int jm1, int j, int jp1, int jp2) {
  return (double)((((-jm1 + 3 * j - 3 * jp1 + jp2) * dx +
                    (2 * jm1 - 5 * j + 4 * jp1 - jp2)) *
                       dx +
                   (-jm1 + jp1)) *
                      dx +
                  (j + j)) /
         2.0;
}


x86 assembly (clang, -O0):
<cubic>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
xor    %eax,%eax
sub    -0xc(%rbp),%eax
imul   $0x3,-0x10(%rbp),%ecx
add    %ecx,%eax
imul   $0x3,-0x14(%rbp),%ecx
sub    %ecx,%eax
add    -0x18(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  -0x8(%rbp),%xmm2
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
imul   $0x5,-0x10(%rbp),%ecx
sub    %ecx,%eax
mov    -0x14(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
sub    -0x18(%rbp),%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x8(%rbp),%xmm2
xor    %eax,%eax
sub    -0xc(%rbp),%eax
add    -0x14(%rbp),%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x8(%rbp),%xmm2
mov    -0x10(%rbp),%eax
add    -0x10(%rbp),%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xe73(%rip),%xmm1
divsd  %xmm1,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<cubic>:
lea    (%rdi,%rdi,1),%r8d
lea    (%rsi,%rsi,4),%eax
sub    %eax,%r8d
mov    %esi,%eax
sub    %edx,%eax
lea    (%rax,%rax,2),%r9d
lea    (%r8,%rdx,4),%eax
sub    %ecx,%eax
sub    %edi,%ecx
add    %r9d,%ecx
cvtsi2sd %ecx,%xmm1
cvtsi2sd %eax,%xmm2
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
sub    %edi,%edx
xorps  %xmm2,%xmm2
cvtsi2sd %edx,%xmm2
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
add    %esi,%esi
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
mulsd  0xea9(%rip),%xmm0
ret

LLVM IR (-O0):
define double @cubic(double %0, i32 %1, i32 %2, i32 %3, i32 %4) {
%6 = alloca double
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store double %0, double* %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store i32 %4, i32* %10
%11 = load i32, i32* %7
%12 = sub nsw i32 0, %11
%13 = load i32, i32* %8
%14 = mul nsw i32 3, %13
%15 = add nsw i32 %12, %14
%16 = load i32, i32* %9
%17 = mul nsw i32 3, %16
%18 = sub nsw i32 %15, %17
%19 = load i32, i32* %10
%20 = add nsw i32 %18, %19
%21 = sitofp i32 %20 to double
%22 = load double, double* %6
%23 = load i32, i32* %7
%24 = mul nsw i32 2, %23
%25 = load i32, i32* %8
%26 = mul nsw i32 5, %25
%27 = sub nsw i32 %24, %26
%28 = load i32, i32* %9
%29 = mul nsw i32 4, %28
%30 = add nsw i32 %27, %29
%31 = load i32, i32* %10
%32 = sub nsw i32 %30, %31
%33 = sitofp i32 %32 to double
%34 = call double @llvm.fmuladd.f64(double %21, double %22, double %33)
%35 = load double, double* %6
%36 = load i32, i32* %7
%37 = sub nsw i32 0, %36
%38 = load i32, i32* %9
%39 = add nsw i32 %37, %38
%40 = sitofp i32 %39 to double
%41 = call double @llvm.fmuladd.f64(double %34, double %35, double %40)
%42 = load double, double* %6
%43 = load i32, i32* %8
%44 = load i32, i32* %8
%45 = add nsw i32 %43, %44
%46 = sitofp i32 %45 to double
%47 = call double @llvm.fmuladd.f64(double %41, double %42, double %46)
%48 = fdiv double %47, 2.000000e+00
ret double %48
}

LLVM IR (-O2):
define double @cubic(double %0, i32 %1, i32 %2, i32 %3, i32 %4) {
%6 = sub i32 %2, %3
%7 = mul i32 %6, 3
%8 = sub i32 %4, %1
%9 = add i32 %8, %7
%10 = sitofp i32 %9 to double
%11 = shl nsw i32 %1, 1
%12 = mul i32 %2, -5
%13 = add i32 %12, %11
%14 = shl nsw i32 %3, 2
%15 = add nsw i32 %13, %14
%16 = sub i32 %15, %4
%17 = sitofp i32 %16 to double
%18 = tail call double @llvm.fmuladd.f64(double %10, double %0, double %17)
%19 = sub i32 %3, %1
%20 = sitofp i32 %19 to double
%21 = tail call double @llvm.fmuladd.f64(double %18, double %0, double %20)
%22 = shl nsw i32 %2, 1
%23 = sitofp i32 %22 to double
%24 = tail call double @llvm.fmuladd.f64(double %21, double %0, double %23)
%25 = fmul double %24, 5.000000e-01
ret double %25
}