float spline_cubic_val(int n, float t[], float tval, float y[], float ypp[]) {
  float dt;
  float h;
  int i;
  int ival;
  float yval;
  ival = n - 2;
  for (i = 0; i < n - 1; i++) {
    if (tval < t[i + 1]) {
      ival = i;
      break;
    }
  }
  dt = tval - t[ival];
  h = t[ival + 1] - t[ival];
  yval = y[ival] +
         dt * ((y[ival + 1] - y[ival]) / h -
               (ypp[ival + 1] / 6.0E+00 + ypp[ival] / 3.0E+00) * h +
               dt * (0.5E+00 * ypp[ival] +
                     dt * ((ypp[ival + 1] - ypp[ival]) / (6.0E+00 * h))));
  return yval;
}


x86 assembly (clang, -O0):
<spline_cubic_val>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1177 <spline_cubic_val+0x77>
movss  -0x14(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss %xmm1,%xmm0
jbe    1164 <spline_cubic_val+0x64>
mov    -0x34(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    1177 <spline_cubic_val+0x77>
jmp    1169 <spline_cubic_val+0x69>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1128 <spline_cubic_val+0x28>
movss  -0x14(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm2
divss  -0x30(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm3
cvtss2sd %xmm3,%xmm3
movsd  0xdf6(%rip),%xmm4
divsd  %xmm4,%xmm3
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm4
cvtss2sd %xmm4,%xmm4
movsd  0xde1(%rip),%xmm5
divsd  %xmm5,%xmm4
addsd  %xmm4,%xmm3
movss  -0x30(%rbp),%xmm4
cvtss2sd %xmm4,%xmm4
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm4,%xmm3
addsd  %xmm2,%xmm3
movss  -0x2c(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm4
cvtss2sd %xmm4,%xmm4
movss  -0x2c(%rbp),%xmm5
cvtss2sd %xmm5,%xmm5
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm6
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm6
cvtss2sd %xmm6,%xmm6
movss  -0x30(%rbp),%xmm7
cvtss2sd %xmm7,%xmm8
movsd  0xd51(%rip),%xmm7
mulsd  %xmm8,%xmm7
divsd  %xmm7,%xmm6
mulsd  %xmm6,%xmm5
movsd  0xd34(%rip),%xmm6
mulsd  %xmm6,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x3c(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<spline_cubic_val>:
lea    -0x2(%rdi),%r8d
add    $0xffffffff,%edi
xor    %eax,%eax
test   %edi,%edi
cmovle %eax,%edi
xchg   %ax,%ax
cmp    %rax,%rdi
je     112a <spline_cubic_val+0x2a>
movss  0x4(%rsi,%rax,4),%xmm1
add    $0x1,%rax
ucomiss %xmm0,%xmm1
jbe    1110 <spline_cubic_val+0x10>
add    $0xffffffff,%eax
mov    %eax,%r8d
movslq %r8d,%rax
movss  (%rsi,%rax,4),%xmm1
movss  0x4(%rsi,%rax,4),%xmm3
subss  %xmm1,%xmm0
subss  %xmm1,%xmm3
movss  (%rdx,%rax,4),%xmm2
movss  0x4(%rdx,%rax,4),%xmm4
xorps  %xmm1,%xmm1
cvtss2sd %xmm2,%xmm1
cvtss2sd %xmm0,%xmm0
subss  %xmm2,%xmm4
divss  %xmm3,%xmm4
movsd  (%rcx,%rax,4),%xmm5
cvtps2pd %xmm5,%xmm2
cvtss2sd %xmm4,%xmm4
movapd %xmm2,%xmm6
divpd  0xe8a(%rip),%xmm6
movapd %xmm6,%xmm7
unpckhpd %xmm6,%xmm7
addsd  %xmm6,%xmm7
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm7
subsd  %xmm7,%xmm4
movaps %xmm5,%xmm6
shufps $0x55,%xmm5,%xmm6
subss  %xmm5,%xmm6
mulsd  0xe6f(%rip),%xmm3
xorps  %xmm5,%xmm5
cvtss2sd %xmm6,%xmm5
divsd  %xmm3,%xmm5
mulsd  0xe64(%rip),%xmm2
mulsd  %xmm0,%xmm5
addsd  %xmm5,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm4,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
xorps  %xmm0,%xmm0
cvtsd2ss %xmm2,%xmm0
ret

LLVM IR (-O0):
define float @spline_cubic_val(i32 %0, float* %1, float %2, float* %3, float* %4) {
%6 = alloca i32
%7 = alloca float*
%8 = alloca float
%9 = alloca float*
%10 = alloca float*
%11 = alloca float
%12 = alloca float
%13 = alloca i32
%14 = alloca i32
%15 = alloca float
store i32 %0, i32* %6
store float* %1, float** %7
store float %2, float* %8
store float* %3, float** %9
store float* %4, float** %10
%16 = load i32, i32* %6
%17 = sub nsw i32 %16, 2
store i32 %17, i32* %14
store i32 0, i32* %13
br label %18
18: ; preds = %35, %5
%19 = load i32, i32* %13
%20 = load i32, i32* %6
%21 = sub nsw i32 %20, 1
%22 = icmp slt i32 %19, %21
br i1 %22, label %23, label %38
23: ; preds = %18
%24 = load float, float* %8
%25 = load float*, float** %7
%26 = load i32, i32* %13
%27 = add nsw i32 %26, 1
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds float, float* %25, i64 %28
%30 = load float, float* %29
%31 = fcmp olt float %24, %30
br i1 %31, label %32, label %34
32: ; preds = %23
%33 = load i32, i32* %13
store i32 %33, i32* %14
br label %38
34: ; preds = %23
br label %35
35: ; preds = %34
%36 = load i32, i32* %13
%37 = add nsw i32 %36, 1
store i32 %37, i32* %13
br label %18
38: ; preds = %32, %18
%39 = load float, float* %8
%40 = load float*, float** %7
%41 = load i32, i32* %14
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds float, float* %40, i64 %42
%44 = load float, float* %43
%45 = fsub float %39, %44
store float %45, float* %11
%46 = load float*, float** %7
%47 = load i32, i32* %14
%48 = add nsw i32 %47, 1
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds float, float* %46, i64 %49
%51 = load float, float* %50
%52 = load float*, float** %7
%53 = load i32, i32* %14
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds float, float* %52, i64 %54
%56 = load float, float* %55
%57 = fsub float %51, %56
store float %57, float* %12
%58 = load float*, float** %9
%59 = load i32, i32* %14
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds float, float* %58, i64 %60
%62 = load float, float* %61
%63 = fpext float %62 to double
%64 = load float, float* %11
%65 = fpext float %64 to double
%66 = load float*, float** %9
%67 = load i32, i32* %14
%68 = add nsw i32 %67, 1
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds float, float* %66, i64 %69
%71 = load float, float* %70
%72 = load float*, float** %9
%73 = load i32, i32* %14
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds float, float* %72, i64 %74
%76 = load float, float* %75
%77 = fsub float %71, %76
%78 = load float, float* %12
%79 = fdiv float %77, %78
%80 = fpext float %79 to double
%81 = load float*, float** %10
%82 = load i32, i32* %14
%83 = add nsw i32 %82, 1
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float, float* %81, i64 %84
%86 = load float, float* %85
%87 = fpext float %86 to double
%88 = fdiv double %87, 6.000000e+00
%89 = load float*, float** %10
%90 = load i32, i32* %14
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds float, float* %89, i64 %91
%93 = load float, float* %92
%94 = fpext float %93 to double
%95 = fdiv double %94, 3.000000e+00
%96 = fadd double %88, %95
%97 = load float, float* %12
%98 = fpext float %97 to double
%99 = fneg double %96
%100 = call double @llvm.fmuladd.f64(double %99, double %98, double %80)
%101 = load float, float* %11
%102 = fpext float %101 to double
%103 = load float*, float** %10
%104 = load i32, i32* %14
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds float, float* %103, i64 %105
%107 = load float, float* %106
%108 = fpext float %107 to double
%109 = load float, float* %11
%110 = fpext float %109 to double
%111 = load float*, float** %10
%112 = load i32, i32* %14
%113 = add nsw i32 %112, 1
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds float, float* %111, i64 %114
%116 = load float, float* %115
%117 = load float*, float** %10
%118 = load i32, i32* %14
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds float, float* %117, i64 %119
%121 = load float, float* %120
%122 = fsub float %116, %121
%123 = fpext float %122 to double
%124 = load float, float* %12
%125 = fpext float %124 to double
%126 = fmul double 6.000000e+00, %125
%127 = fdiv double %123, %126
%128 = fmul double %110, %127
%129 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %108, double %128)
%130 = call double @llvm.fmuladd.f64(double %102, double %129, double %100)
%131 = call double @llvm.fmuladd.f64(double %65, double %130, double %63)
%132 = fptrunc double %131 to float
store float %132, float* %15
%133 = load float, float* %15
ret float %133
}

LLVM IR (-O2):
define float @spline_cubic_val(i32 %0, float* nocapture %1, float %2, float* nocapture %3, float* nocapture %4) {
%6 = add nsw i32 %0, -2
%7 = add i32 %0, -1
%8 = call i32 @llvm.smax.i32(i32 %7, i32 0)
%9 = zext i32 %8 to i64
br label %10
10: ; preds = %13, %5
%11 = phi i64 [ %14, %13 ], [ 0, %5 ]
%12 = icmp eq i64 %11, %9
br i1 %12, label %20, label %13
13: ; preds = %10
%14 = add nuw nsw i64 %11, 1
%15 = getelementptr inbounds float, float* %1, i64 %14
%16 = load float, float* %15
%17 = fcmp ogt float %16, %2
br i1 %17, label %18, label %10
18: ; preds = %13
%19 = trunc i64 %11 to i32
br label %20
20: ; preds = %10, %18
%21 = phi i32 [ %19, %18 ], [ %6, %10 ]
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds float, float* %1, i64 %22
%24 = load float, float* %23
%25 = fsub float %2, %24
%26 = add nsw i32 %21, 1
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds float, float* %1, i64 %27
%29 = load float, float* %28
%30 = fsub float %29, %24
%31 = getelementptr inbounds float, float* %3, i64 %22
%32 = load float, float* %31
%33 = fpext float %32 to double
%34 = fpext float %25 to double
%35 = getelementptr inbounds float, float* %3, i64 %27
%36 = load float, float* %35
%37 = fsub float %36, %32
%38 = fdiv float %37, %30
%39 = fpext float %38 to double
%40 = getelementptr inbounds float, float* %4, i64 %22
%41 = bitcast float* %40 to <2 x float>*
%42 = load <2 x float>, <2 x float>* %41
%43 = fpext <2 x float> %42 to <2 x double>
%44 = fdiv <2 x double> %43, <double 3.000000e+00, double 6.000000e+00>
%45 = shufflevector <2 x double> %44, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%46 = fadd <2 x double> %45, %44
%47 = extractelement <2 x double> %46, i64 0
%48 = fpext float %30 to double
%49 = fneg double %47
%50 = tail call double @llvm.fmuladd.f64(double %49, double %48, double %39)
%51 = shufflevector <2 x float> %42, <2 x float> poison, <2 x i32> <i32 1, i32 undef>
%52 = fsub <2 x float> %51, %42
%53 = extractelement <2 x float> %52, i64 0
%54 = fpext float %53 to double
%55 = fmul double %48, 6.000000e+00
%56 = fdiv double %54, %55
%57 = fmul double %56, %34
%58 = extractelement <2 x double> %43, i64 0
%59 = tail call double @llvm.fmuladd.f64(double %58, double 5.000000e-01, double %57)
%60 = tail call double @llvm.fmuladd.f64(double %34, double %59, double %50)
%61 = tail call double @llvm.fmuladd.f64(double %34, double %60, double %33)
%62 = fptrunc double %61 to float
ret float %62
}