float interpolateHermite(float *data, double address, long arraysize) {
  long pos = (long)address;
  float posFract = (float)(address - (double)pos);
  long posMinus1 = (pos == 0) ? arraysize - 1 : pos - 1;
  long posPlus1 = (pos + 1) % arraysize;
  long posPlus2 = (pos + 2) % arraysize;
  float a = ((3.0f * (data[pos] - data[posPlus1])) - data[posMinus1] +
             data[posPlus2]) *
            0.5f;
  float b = (2.0f * data[posPlus1]) + data[posMinus1] - (2.5f * data[pos]) -
            (data[posPlus2] * 0.5f);
  float c = (data[posPlus1] - data[posMinus1]) * 0.5f;
  return ((((a * posFract) + b) * posFract + c) * posFract) + data[pos];
}


x86 assembly (clang, -O0):
<interpolateHermite>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
cvttsd2si -0x10(%rbp),%rax
mov    %rax,-0x20(%rbp)
movsd  -0x10(%rbp),%xmm0
cvtsi2sdq -0x20(%rbp),%xmm1
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
cmpq   $0x0,-0x20(%rbp)
jne    114f <interpolateHermite+0x4f>
mov    -0x18(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x58(%rbp)
jmp    115b <interpolateHermite+0x5b>
mov    -0x20(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x58(%rbp)
mov    -0x58(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
cqto
idivq  -0x18(%rbp)
mov    %rdx,-0x38(%rbp)
mov    -0x20(%rbp),%rax
add    $0x2,%rax
cqto
idivq  -0x18(%rbp)
mov    %rdx,-0x40(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
movss  0xe45(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  0xe20(%rip),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
addss  %xmm0,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movss  0xde4(%rip),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  0xdb6(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  0xd87(%rip),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x44(%rbp),%xmm0
movss  -0x24(%rbp),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x24(%rbp),%xmm2
movss  -0x4c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x24(%rbp),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<interpolateHermite>:
cvttsd2si %xmm0,%rcx
cvtsi2sd %rcx,%xmm1
test   %rcx,%rcx
mov    %rcx,%r9
cmove  %rsi,%r9
subsd  %xmm1,%xmm0
lea    0x1(%rcx),%rax
mov    %rax,%rdx
or     %rsi,%rdx
shr    $0x20,%rdx
je     1132 <interpolateHermite+0x32>
cqto
idiv   %rsi
mov    %rdx,%r8
jmp    1139 <interpolateHermite+0x39>
xor    %edx,%edx
div    %esi
mov    %edx,%r8d
cvtsd2ss %xmm0,%xmm8
lea    0x2(%rcx),%rax
mov    %rax,%rdx
or     %rsi,%rdx
shr    $0x20,%rdx
je     1155 <interpolateHermite+0x55>
cqto
idiv   %rsi
jmp    1159 <interpolateHermite+0x59>
xor    %edx,%edx
div    %esi
movss  (%rdi,%rcx,4),%xmm2
movss  (%rdi,%r8,4),%xmm3
movaps %xmm2,%xmm0
subss  %xmm3,%xmm0
movss  -0x4(%rdi,%r9,4),%xmm4
mulss  0xe86(%rip),%xmm0
subss  %xmm4,%xmm0
movss  (%rdi,%rdx,4),%xmm5
addss  %xmm5,%xmm0
movss  0xe75(%rip),%xmm6
mulss  %xmm6,%xmm0
movaps %xmm3,%xmm7
addss  %xmm3,%xmm7
addss  %xmm4,%xmm7
movss  0xe62(%rip),%xmm1
mulss  %xmm2,%xmm1
addss  %xmm7,%xmm1
mulss  %xmm6,%xmm5
subss  %xmm5,%xmm1
subss  %xmm4,%xmm3
mulss  %xmm6,%xmm3
mulss  %xmm8,%xmm0
addss  %xmm1,%xmm0
mulss  %xmm8,%xmm0
addss  %xmm3,%xmm0
mulss  %xmm8,%xmm0
addss  %xmm2,%xmm0
ret

LLVM IR (-O0):
define float @interpolateHermite(float* %0, double %1, i64 %2) {
%4 = alloca float*
%5 = alloca double
%6 = alloca i64
%7 = alloca i64
%8 = alloca float
%9 = alloca i64
%10 = alloca i64
%11 = alloca i64
%12 = alloca float
%13 = alloca float
%14 = alloca float
store float* %0, float** %4
store double %1, double* %5
store i64 %2, i64* %6
%15 = load double, double* %5
%16 = fptosi double %15 to i64
store i64 %16, i64* %7
%17 = load double, double* %5
%18 = load i64, i64* %7
%19 = sitofp i64 %18 to double
%20 = fsub double %17, %19
%21 = fptrunc double %20 to float
store float %21, float* %8
%22 = load i64, i64* %7
%23 = icmp eq i64 %22, 0
br i1 %23, label %24, label %27
24: ; preds = %3
%25 = load i64, i64* %6
%26 = sub nsw i64 %25, 1
br label %30
27: ; preds = %3
%28 = load i64, i64* %7
%29 = sub nsw i64 %28, 1
br label %30
30: ; preds = %27, %24
%31 = phi i64 [ %26, %24 ], [ %29, %27 ]
store i64 %31, i64* %9
%32 = load i64, i64* %7
%33 = add nsw i64 %32, 1
%34 = load i64, i64* %6
%35 = srem i64 %33, %34
store i64 %35, i64* %10
%36 = load i64, i64* %7
%37 = add nsw i64 %36, 2
%38 = load i64, i64* %6
%39 = srem i64 %37, %38
store i64 %39, i64* %11
%40 = load float*, float** %4
%41 = load i64, i64* %7
%42 = getelementptr inbounds float, float* %40, i64 %41
%43 = load float, float* %42
%44 = load float*, float** %4
%45 = load i64, i64* %10
%46 = getelementptr inbounds float, float* %44, i64 %45
%47 = load float, float* %46
%48 = fsub float %43, %47
%49 = load float*, float** %4
%50 = load i64, i64* %9
%51 = getelementptr inbounds float, float* %49, i64 %50
%52 = load float, float* %51
%53 = fneg float %52
%54 = call float @llvm.fmuladd.f32(float 3.000000e+00, float %48, float %53)
%55 = load float*, float** %4
%56 = load i64, i64* %11
%57 = getelementptr inbounds float, float* %55, i64 %56
%58 = load float, float* %57
%59 = fadd float %54, %58
%60 = fmul float %59, 5.000000e-01
store float %60, float* %12
%61 = load float*, float** %4
%62 = load i64, i64* %10
%63 = getelementptr inbounds float, float* %61, i64 %62
%64 = load float, float* %63
%65 = load float*, float** %4
%66 = load i64, i64* %9
%67 = getelementptr inbounds float, float* %65, i64 %66
%68 = load float, float* %67
%69 = call float @llvm.fmuladd.f32(float 2.000000e+00, float %64, float %68)
%70 = load float*, float** %4
%71 = load i64, i64* %7
%72 = getelementptr inbounds float, float* %70, i64 %71
%73 = load float, float* %72
%74 = call float @llvm.fmuladd.f32(float -2.500000e+00, float %73, float %69)
%75 = load float*, float** %4
%76 = load i64, i64* %11
%77 = getelementptr inbounds float, float* %75, i64 %76
%78 = load float, float* %77
%79 = fneg float %78
%80 = call float @llvm.fmuladd.f32(float %79, float 5.000000e-01, float %74)
store float %80, float* %13
%81 = load float*, float** %4
%82 = load i64, i64* %10
%83 = getelementptr inbounds float, float* %81, i64 %82
%84 = load float, float* %83
%85 = load float*, float** %4
%86 = load i64, i64* %9
%87 = getelementptr inbounds float, float* %85, i64 %86
%88 = load float, float* %87
%89 = fsub float %84, %88
%90 = fmul float %89, 5.000000e-01
store float %90, float* %14
%91 = load float, float* %12
%92 = load float, float* %8
%93 = load float, float* %13
%94 = call float @llvm.fmuladd.f32(float %91, float %92, float %93)
%95 = load float, float* %8
%96 = load float, float* %14
%97 = call float @llvm.fmuladd.f32(float %94, float %95, float %96)
%98 = load float, float* %8
%99 = load float*, float** %4
%100 = load i64, i64* %7
%101 = getelementptr inbounds float, float* %99, i64 %100
%102 = load float, float* %101
%103 = call float @llvm.fmuladd.f32(float %97, float %98, float %102)
ret float %103
}

LLVM IR (-O2):
define float @interpolateHermite(float* nocapture %0, double %1, i64 %2) {
%4 = fptosi double %1 to i64
%5 = sitofp i64 %4 to double
%6 = fsub double %1, %5
%7 = fptrunc double %6 to float
%8 = icmp eq i64 %4, 0
%9 = select i1 %8, i64 %2, i64 %4
%10 = add nsw i64 %9, -1
%11 = add nsw i64 %4, 1
%12 = srem i64 %11, %2
%13 = add nsw i64 %4, 2
%14 = srem i64 %13, %2
%15 = getelementptr inbounds float, float* %0, i64 %4
%16 = load float, float* %15
%17 = getelementptr inbounds float, float* %0, i64 %12
%18 = load float, float* %17
%19 = fsub float %16, %18
%20 = getelementptr inbounds float, float* %0, i64 %10
%21 = load float, float* %20
%22 = fneg float %21
%23 = tail call float @llvm.fmuladd.f32(float %19, float 3.000000e+00, float %22)
%24 = getelementptr inbounds float, float* %0, i64 %14
%25 = load float, float* %24
%26 = fadd float %23, %25
%27 = fmul float %26, 5.000000e-01
%28 = tail call float @llvm.fmuladd.f32(float %18, float 2.000000e+00, float %21)
%29 = tail call float @llvm.fmuladd.f32(float %16, float -2.500000e+00, float %28)
%30 = fneg float %25
%31 = tail call float @llvm.fmuladd.f32(float %30, float 5.000000e-01, float %29)
%32 = fsub float %18, %21
%33 = fmul float %32, 5.000000e-01
%34 = tail call float @llvm.fmuladd.f32(float %27, float %7, float %31)
%35 = tail call float @llvm.fmuladd.f32(float %34, float %7, float %33)
%36 = tail call float @llvm.fmuladd.f32(float %35, float %7, float %16)
ret float %36
}