int boscili(float amp, float si, double *farray, int len, float *phs,
            float *array, int alen) {
  register int i, j, k;
  float frac, temp;
  float *fp = array;
  temp = *phs;
  for (j = alen; j > 0; --j) {
    i = temp;
    k = (i + 1) % len;
    frac = temp - i;
    temp += si;
    while (temp >= len)
      temp -= len;
    *fp++ = ((*(farray + i) + (*(farray + k) - *(farray + i)) * frac) * amp);
  }
  *phs = temp;
  return (j);
}


x86 assembly (clang, -O0):
<boscili>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x4(%rbp)
movss  %xmm1,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8d,-0x2c(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jle    1212 <boscili+0x112>
cvttss2si -0x40(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
cltd
idivl  -0x14(%rbp)
mov    %edx,-0x38(%rbp)
movss  -0x40(%rbp),%xmm0
cvtsi2ssl -0x30(%rbp),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x8(%rbp),%xmm0
addss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x40(%rbp),%xmm0
cvtsi2ssl -0x14(%rbp),%xmm1
ucomiss %xmm1,%xmm0
jb     11a8 <boscili+0xa8>
cvtsi2ssl -0x14(%rbp),%xmm1
movss  -0x40(%rbp),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
jmp    117d <boscili+0x7d>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movss  -0x3c(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movss  -0x4(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x48(%rbp)
movss  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    113c <boscili+0x3c>
movss  -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<boscili>:
mov    %rdx,%r9
movss  (%rdx),%xmm4
test   %r8d,%r8d
jle    1190 <boscili+0x90>
cvtsi2ss %esi,%xmm2
cvtss2sd %xmm0,%xmm3
movaps %xmm4,%xmm0
jmp    116b <boscili+0x6b>
nopl   (%rax)
cvttps2dq %xmm4,%xmm5
cvtdq2ps %xmm5,%xmm5
subss  %xmm5,%xmm4
movslq %r10d,%rax
movsd  (%rdi,%rax,8),%xmm5
movslq %edx,%rax
movsd  (%rdi,%rax,8),%xmm6
subsd  %xmm5,%xmm6
cvtss2sd %xmm4,%xmm4
mulsd  %xmm6,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm3,%xmm4
cvtsd2ss %xmm4,%xmm4
movss  %xmm4,(%rcx)
add    $0x4,%rcx
lea    -0x1(%r8),%eax
movaps %xmm0,%xmm4
cmp    $0x1,%r8d
mov    %eax,%r8d
jle    118b <boscili+0x8b>
cvttss2si %xmm4,%r10d
lea    0x1(%r10),%eax
cltd
idiv   %esi
addss  %xmm1,%xmm0
ucomiss %xmm2,%xmm0
jb     1120 <boscili+0x20>
subss  %xmm2,%xmm0
ucomiss %xmm2,%xmm0
jae    1180 <boscili+0x80>
jmp    1120 <boscili+0x20>
xor    %r8d,%r8d
jmp    1193 <boscili+0x93>
movaps %xmm4,%xmm0
movss  %xmm0,(%r9)
mov    %r8d,%eax
ret

LLVM IR (-O0):
define i32 @boscili(float %0, float %1, double* %2, i32 %3, float* %4, float* %5, i32 %6) {
%8 = alloca float
%9 = alloca float
%10 = alloca double*
%11 = alloca i32
%12 = alloca float*
%13 = alloca float*
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca float
%19 = alloca float
%20 = alloca float*
store float %0, float* %8
store float %1, float* %9
store double* %2, double** %10
store i32 %3, i32* %11
store float* %4, float** %12
store float* %5, float** %13
store i32 %6, i32* %14
%21 = load float*, float** %13
store float* %21, float** %20
%22 = load float*, float** %12
%23 = load float, float* %22
store float %23, float* %19
%24 = load i32, i32* %14
store i32 %24, i32* %16
br label %25
25: ; preds = %78, %7
%26 = load i32, i32* %16
%27 = icmp sgt i32 %26, 0
br i1 %27, label %28, label %81
28: ; preds = %25
%29 = load float, float* %19
%30 = fptosi float %29 to i32
store i32 %30, i32* %15
%31 = load i32, i32* %15
%32 = add nsw i32 %31, 1
%33 = load i32, i32* %11
%34 = srem i32 %32, %33
store i32 %34, i32* %17
%35 = load float, float* %19
%36 = load i32, i32* %15
%37 = sitofp i32 %36 to float
%38 = fsub float %35, %37
store float %38, float* %18
%39 = load float, float* %9
%40 = load float, float* %19
%41 = fadd float %40, %39
store float %41, float* %19
br label %42
42: ; preds = %47, %28
%43 = load float, float* %19
%44 = load i32, i32* %11
%45 = sitofp i32 %44 to float
%46 = fcmp oge float %43, %45
br i1 %46, label %47, label %52
47: ; preds = %42
%48 = load i32, i32* %11
%49 = sitofp i32 %48 to float
%50 = load float, float* %19
%51 = fsub float %50, %49
store float %51, float* %19
br label %42
52: ; preds = %42
%53 = load double*, double** %10
%54 = load i32, i32* %15
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %53, i64 %55
%57 = load double, double* %56
%58 = load double*, double** %10
%59 = load i32, i32* %17
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %58, i64 %60
%62 = load double, double* %61
%63 = load double*, double** %10
%64 = load i32, i32* %15
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
%67 = load double, double* %66
%68 = fsub double %62, %67
%69 = load float, float* %18
%70 = fpext float %69 to double
%71 = call double @llvm.fmuladd.f64(double %68, double %70, double %57)
%72 = load float, float* %8
%73 = fpext float %72 to double
%74 = fmul double %71, %73
%75 = fptrunc double %74 to float
%76 = load float*, float** %20
%77 = getelementptr inbounds float, float* %76, i32 1
store float* %77, float** %20
store float %75, float* %76
br label %78
78: ; preds = %52
%79 = load i32, i32* %16
%80 = add nsw i32 %79, -1
store i32 %80, i32* %16
br label %25
81: ; preds = %25
%82 = load float, float* %19
%83 = load float*, float** %12
store float %82, float* %83
%84 = load i32, i32* %16
ret i32 %84
}

LLVM IR (-O2):
define i32 @boscili(float %0, float %1, double* nocapture %2, i32 %3, float* nocapture %4, float* nocapture writeonly %5, i32 %6) {
%8 = load float, float* %4
%9 = icmp sgt i32 %6, 0
br i1 %9, label %10, label %44
10: ; preds = %7
%11 = sitofp i32 %3 to float
%12 = fpext float %0 to double
br label %13
13: ; preds = %10, %28
%14 = phi float* [ %5, %10 ], [ %41, %28 ]
%15 = phi float [ %8, %10 ], [ %29, %28 ]
%16 = phi i32 [ %6, %10 ], [ %42, %28 ]
%17 = fptosi float %15 to i32
%18 = add nsw i32 %17, 1
%19 = srem i32 %18, %3
%20 = sitofp i32 %17 to float
%21 = fsub float %15, %20
%22 = fadd float %15, %1
%23 = fcmp ult float %22, %11
br i1 %23, label %28, label %24
24: ; preds = %13, %24
%25 = phi float [ %26, %24 ], [ %22, %13 ]
%26 = fsub float %25, %11
%27 = fcmp ult float %26, %11
br i1 %27, label %28, label %24
28: ; preds = %24, %13
%29 = phi float [ %22, %13 ], [ %26, %24 ]
%30 = sext i32 %17 to i64
%31 = getelementptr inbounds double, double* %2, i64 %30
%32 = load double, double* %31
%33 = sext i32 %19 to i64
%34 = getelementptr inbounds double, double* %2, i64 %33
%35 = load double, double* %34
%36 = fsub double %35, %32
%37 = fpext float %21 to double
%38 = tail call double @llvm.fmuladd.f64(double %36, double %37, double %32)
%39 = fmul double %38, %12
%40 = fptrunc double %39 to float
%41 = getelementptr inbounds float, float* %14, i64 1
store float %40, float* %14
%42 = add nsw i32 %16, -1
%43 = icmp sgt i32 %16, 1
br i1 %43, label %13, label %44
44: ; preds = %28, %7
%45 = phi i32 [ %6, %7 ], [ 0, %28 ]
%46 = phi float [ %8, %7 ], [ %29, %28 ]
store float %46, float* %4
ret i32 %45
}