float wshape(float x, double *f, int len) {
  int i1;
  float x1, diff;
  x1 = (x + 1) * (float)((len - 1) / 2.);
  i1 = x1;
  diff = (float)(x1 - i1);
  return (*(f + i1) + (*(f + i1 + 1) - *(f + i1)) * diff);
}


x86 assembly (clang, -O0):
<wshape>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x4(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
movss  0xef0(%rip),%xmm0
addss  -0x4(%rbp),%xmm0
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  0xed1(%rip),%xmm2
divsd  %xmm2,%xmm1
cvtsd2ss %xmm1,%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
cvttss2si -0x1c(%rbp),%eax
mov    %eax,-0x18(%rbp)
movss  -0x1c(%rbp),%xmm0
cvtsi2ssl -0x18(%rbp),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  0x8(%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movss  -0x20(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<wshape>:
addss  0xef8(%rip),%xmm0
add    $0xffffffff,%esi
cvtsi2sd %esi,%xmm1
mulsd  0xef1(%rip),%xmm1
cvtsd2ss %xmm1,%xmm1
mulss  %xmm0,%xmm1
cvttss2si %xmm1,%eax
cvttps2dq %xmm1,%xmm0
cvtdq2ps %xmm0,%xmm0
subss  %xmm0,%xmm1
cltq
movsd  (%rdi,%rax,8),%xmm0
movsd  0x8(%rdi,%rax,8),%xmm2
subsd  %xmm0,%xmm2
cvtss2sd %xmm1,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
ret

LLVM IR (-O0):
define float @wshape(float %0, double* %1, i32 %2) {
%4 = alloca float
%5 = alloca double*
%6 = alloca i32
%7 = alloca i32
%8 = alloca float
%9 = alloca float
store float %0, float* %4
store double* %1, double** %5
store i32 %2, i32* %6
%10 = load float, float* %4
%11 = fadd float %10, 1.000000e+00
%12 = load i32, i32* %6
%13 = sub nsw i32 %12, 1
%14 = sitofp i32 %13 to double
%15 = fdiv double %14, 2.000000e+00
%16 = fptrunc double %15 to float
%17 = fmul float %11, %16
store float %17, float* %8
%18 = load float, float* %8
%19 = fptosi float %18 to i32
store i32 %19, i32* %7
%20 = load float, float* %8
%21 = load i32, i32* %7
%22 = sitofp i32 %21 to float
%23 = fsub float %20, %22
store float %23, float* %9
%24 = load double*, double** %5
%25 = load i32, i32* %7
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double, double* %24, i64 %26
%28 = load double, double* %27
%29 = load double*, double** %5
%30 = load i32, i32* %7
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds double, double* %29, i64 %31
%33 = getelementptr inbounds double, double* %32, i64 1
%34 = load double, double* %33
%35 = load double*, double** %5
%36 = load i32, i32* %7
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds double, double* %35, i64 %37
%39 = load double, double* %38
%40 = fsub double %34, %39
%41 = load float, float* %9
%42 = fpext float %41 to double
%43 = call double @llvm.fmuladd.f64(double %40, double %42, double %28)
%44 = fptrunc double %43 to float
ret float %44
}

LLVM IR (-O2):
define float @wshape(float %0, double* nocapture %1, i32 %2) {
%4 = fadd float %0, 1.000000e+00
%5 = add nsw i32 %2, -1
%6 = sitofp i32 %5 to double
%7 = fmul double %6, 5.000000e-01
%8 = fptrunc double %7 to float
%9 = fmul float %4, %8
%10 = fptosi float %9 to i32
%11 = sitofp i32 %10 to float
%12 = fsub float %9, %11
%13 = sext i32 %10 to i64
%14 = getelementptr inbounds double, double* %1, i64 %13
%15 = load double, double* %14
%16 = getelementptr inbounds double, double* %14, i64 1
%17 = load double, double* %16
%18 = fsub double %17, %15
%19 = fpext float %12 to double
%20 = tail call double @llvm.fmuladd.f64(double %18, double %19, double %15)
%21 = fptrunc double %20 to float
ret float %21
}