void agbnp3_cspline_interpolate_soa(float *kv, float *xh, float dx, int m,
                                    float *yp, float *y, float *y2p, float *y2,
                                    float *f, float *fp) {
  int i;
  float b, a, a2, b2, kf, dp1, dp2, yk, ypk, y2k, y2pk, dxinv;
  dp1 = dx / 6.0f;
  dp2 = dx * dp1;
  dxinv = 1.f / dx;
#pragma vector aligned
#pragma ivdep
  for (i = 0; i < m; i++) {
    kf = kv[i];
    yk = y[i];
    ypk = yp[i];
    y2k = y2[i];
    y2pk = y2p[i];
    a = (kf + 1.0f) - xh[i];
    a2 = a * a;
    b = xh[i] - kf;
    b2 = b * b;
    f[i] = a * yk + b * ypk + ((a2 * a - a) * y2k + (b2 * b - b) * y2pk) * dp2;
    fp[i] = (ypk - yk) * dxinv -
            ((3.0f * a2 - 1.f) * y2k + (3.0f * b2 - 1.f) * y2pk) * dp1;
  }
}


x86 assembly (clang, -O0):
<agbnp3_cspline_interpolate_soa>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movss  -0x14(%rbp),%xmm0
movss  0xecb(%rip),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x14(%rbp),%xmm0
mulss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  0xea7(%rip),%xmm0
divss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x64(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1344 <agbnp3_cspline_interpolate_soa+0x244>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x54(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x58(%rbp)
mov    0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x5c(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x60(%rbp)
movss  0xe28(%rip),%xmm0
addss  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x3c(%rbp),%xmm0
mulss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm0
mulss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x3c(%rbp),%xmm1
movss  -0x54(%rbp),%xmm2
movss  -0x38(%rbp),%xmm0
mulss  -0x58(%rbp),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x40(%rbp),%xmm0
movss  -0x3c(%rbp),%xmm3
movss  -0x3c(%rbp),%xmm2
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  %xmm3,%xmm0
addss  %xmm2,%xmm0
movss  -0x5c(%rbp),%xmm3
movss  -0x44(%rbp),%xmm2
movss  -0x38(%rbp),%xmm5
movss  -0x38(%rbp),%xmm4
movd   %xmm4,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm4
mulss  %xmm5,%xmm2
addss  %xmm4,%xmm2
mulss  -0x60(%rbp),%xmm2
mulss  %xmm3,%xmm0
addss  %xmm2,%xmm0
movss  -0x50(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x58(%rbp),%xmm0
subss  -0x54(%rbp),%xmm0
movss  -0x64(%rbp),%xmm2
movss  -0x40(%rbp),%xmm1
movss  0xd38(%rip),%xmm3
mulss  %xmm3,%xmm1
movss  0xd30(%rip),%xmm3
addss  %xmm3,%xmm1
movss  -0x5c(%rbp),%xmm4
movss  -0x44(%rbp),%xmm3
movss  0xd16(%rip),%xmm5
mulss  %xmm5,%xmm3
movss  0xd0e(%rip),%xmm5
addss  %xmm5,%xmm3
mulss  -0x60(%rbp),%xmm3
mulss  %xmm4,%xmm1
addss  %xmm3,%xmm1
mulss  -0x4c(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    116a <agbnp3_cspline_interpolate_soa+0x6a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<agbnp3_cspline_interpolate_soa>:
push   %rbx
test   %edx,%edx
jle    1228 <agbnp3_cspline_interpolate_soa+0x128>
mov    0x20(%rsp),%r10
mov    0x18(%rsp),%r11
movaps %xmm0,%xmm11
divss  0xee0(%rip),%xmm11
mov    0x10(%rsp),%rax
movaps %xmm11,%xmm9
mulss  %xmm0,%xmm9
movss  0xecd(%rip),%xmm8
movaps %xmm8,%xmm13
divss  %xmm0,%xmm13
mov    %edx,%edx
xor    %ebx,%ebx
movss  0xebb(%rip),%xmm10
movss  0xeb6(%rip),%xmm12
cs nopw 0x0(%rax,%rax,1)
movss  (%rdi,%rbx,4),%xmm4
movss  (%r8,%rbx,4),%xmm1
movss  (%rcx,%rbx,4),%xmm7
movss  (%rax,%rbx,4),%xmm14
movss  (%r9,%rbx,4),%xmm15
movaps %xmm4,%xmm0
addss  %xmm8,%xmm0
movss  (%rsi,%rbx,4),%xmm5
subss  %xmm5,%xmm0
movaps %xmm0,%xmm2
mulss  %xmm0,%xmm2
subss  %xmm4,%xmm5
movaps %xmm5,%xmm4
mulss  %xmm5,%xmm4
movaps %xmm7,%xmm6
mulss  %xmm5,%xmm6
movaps %xmm2,%xmm3
mulss  %xmm0,%xmm3
subss  %xmm0,%xmm3
mulss  %xmm1,%xmm0
addss  %xmm6,%xmm0
movaps %xmm4,%xmm6
mulss  %xmm5,%xmm6
subss  %xmm5,%xmm6
mulss  %xmm15,%xmm6
mulss  %xmm14,%xmm3
addss  %xmm6,%xmm3
mulss  %xmm9,%xmm3
addss  %xmm0,%xmm3
movss  %xmm3,(%r11,%rbx,4)
subss  %xmm1,%xmm7
mulss  %xmm13,%xmm7
mulss  %xmm10,%xmm2
addss  %xmm12,%xmm2
mulss  %xmm10,%xmm4
addss  %xmm12,%xmm4
mulss  %xmm15,%xmm4
mulss  %xmm14,%xmm2
addss  %xmm4,%xmm2
mulss  %xmm11,%xmm2
subss  %xmm2,%xmm7
movss  %xmm7,(%r10,%rbx,4)
add    $0x1,%rbx
cmp    %rbx,%rdx
jne    1160 <agbnp3_cspline_interpolate_soa+0x60>
pop    %rbx
ret

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

LLVM IR (-O2):
define void @agbnp3_cspline_interpolate_soa(float* nocapture %0, float* nocapture %1, float %2, i32 %3, float* nocapture %4, float* nocapture %5, float* nocapture %6, float* nocapture %7, float* nocapture writeonly %8, float* nocapture writeonly %9) {
%11 = fdiv float %2, 6.000000e+00
%12 = fmul float %11, %2
%13 = fdiv float 1.000000e+00, %2
%14 = icmp sgt i32 %3, 0
br i1 %14, label %15, label %57
15: ; preds = %10
%16 = zext i32 %3 to i64
br label %17
17: ; preds = %15, %17
%18 = phi i64 [ 0, %15 ], [ %55, %17 ]
%19 = getelementptr inbounds float, float* %0, i64 %18
%20 = load float, float* %19
%21 = getelementptr inbounds float, float* %5, i64 %18
%22 = load float, float* %21
%23 = getelementptr inbounds float, float* %4, i64 %18
%24 = load float, float* %23
%25 = getelementptr inbounds float, float* %7, i64 %18
%26 = load float, float* %25
%27 = getelementptr inbounds float, float* %6, i64 %18
%28 = load float, float* %27
%29 = fadd float %20, 1.000000e+00
%30 = getelementptr inbounds float, float* %1, i64 %18
%31 = load float, float* %30
%32 = fsub float %29, %31
%33 = fmul float %32, %32
%34 = fsub float %31, %20
%35 = fmul float %34, %34
%36 = fmul float %24, %34
%37 = tail call float @llvm.fmuladd.f32(float %32, float %22, float %36)
%38 = fneg float %32
%39 = tail call float @llvm.fmuladd.f32(float %33, float %32, float %38)
%40 = fneg float %34
%41 = tail call float @llvm.fmuladd.f32(float %35, float %34, float %40)
%42 = fmul float %28, %41
%43 = tail call float @llvm.fmuladd.f32(float %39, float %26, float %42)
%44 = tail call float @llvm.fmuladd.f32(float %43, float %12, float %37)
%45 = getelementptr inbounds float, float* %8, i64 %18
store float %44, float* %45
%46 = fsub float %24, %22
%47 = tail call float @llvm.fmuladd.f32(float %33, float 3.000000e+00, float -1.000000e+00)
%48 = tail call float @llvm.fmuladd.f32(float %35, float 3.000000e+00, float -1.000000e+00)
%49 = fmul float %28, %48
%50 = tail call float @llvm.fmuladd.f32(float %47, float %26, float %49)
%51 = fneg float %50
%52 = fmul float %11, %51
%53 = tail call float @llvm.fmuladd.f32(float %46, float %13, float %52)
%54 = getelementptr inbounds float, float* %9, i64 %18
store float %53, float* %54
%55 = add nuw nsw i64 %18, 1
%56 = icmp eq i64 %55, %16
br i1 %56, label %57, label %17
57: ; preds = %17, %10
ret void
}