#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int floor(double);
double BGBMID_SplinePolate(short *samp, int len, double t) {
  double f, g, u, v;
  double a, b, c, d;
  double p;
  int i;
  i = floor(t);
  if (i < 1) {
    if (i < 0)
      return (samp[0]);
    f = t - i;
    g = 1.0 - f;
    a = samp[0];
    b = samp[1];
    p = (a * g) + (b * f);
    return (p);
  }
  if (i >= (len - 2)) {
    if (i >= (len - 1))
      return (samp[len - 1]);
    f = t - i;
    g = 1.0 - f;
    a = samp[i + 0];
    b = samp[i + 1];
    p = (a * g) + (b * f);
    return (p);
  }
  f = t - i;
  g = 1.0 - f;
  a = samp[i - 1];
  b = samp[i + 0];
  c = samp[i + 1];
  d = samp[i + 2];
  u = b + f * (b - a);
  v = c + g * (c - d);
  p = (u * g) + (v * f);
  return (p);
}


x86 assembly (clang, -O0):
<BGBMID_SplinePolate>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
call   1030 <floor@plt>
mov    %eax,-0x6c(%rbp)
cmpl   $0x1,-0x6c(%rbp)
jge    11d4 <BGBMID_SplinePolate+0xc4>
cmpl   $0x0,-0x6c(%rbp)
jge    115a <BGBMID_SplinePolate+0x4a>
mov    -0x10(%rbp),%rax
movswl (%rax),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    138a <BGBMID_SplinePolate+0x27a>
movsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x6c(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm1
movsd  0xe86(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movswl (%rax),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movswl 0x2(%rax),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    138a <BGBMID_SplinePolate+0x27a>
mov    -0x6c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jl     1299 <BGBMID_SplinePolate+0x189>
mov    -0x6c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jl     1215 <BGBMID_SplinePolate+0x105>
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0xffffffff,%ecx
movslq %ecx,%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    138a <BGBMID_SplinePolate+0x27a>
movsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x6c(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm1
movsd  0xdcb(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movswl 0x2(%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    138a <BGBMID_SplinePolate+0x27a>
movsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x6c(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm1
movsd  0xd47(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
add    $0xffffffff,%ecx
movslq %ecx,%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movswl 0x2(%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movswl 0x4(%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x50(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
subsd  -0x48(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x58(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
subsd  -0x60(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<BGBMID_SplinePolate>:
push   %r14
push   %rbx
push   %rax
mov    %esi,%ebx
mov    %rdi,%r14
movsd  %xmm0,(%rsp)
call   1030 <floor@plt>
test   %eax,%eax
jle    11ae <BGBMID_SplinePolate+0x9e>
lea    -0x2(%rbx),%ecx
cmp    %ecx,%eax
jge    11d4 <BGBMID_SplinePolate+0xc4>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  (%rsp),%xmm6
subsd  %xmm0,%xmm6
movsd  0xeb2(%rip),%xmm1
lea    -0x1(%rax),%ecx
movswl (%r14,%rcx,2),%ecx
cvtsi2sd %ecx,%xmm2
mov    %eax,%eax
movswl (%r14,%rax,2),%ecx
cvtsi2sd %ecx,%xmm3
subsd  %xmm6,%xmm1
movswl 0x2(%r14,%rax,2),%ecx
cvtsi2sd %ecx,%xmm4
movswl 0x4(%r14,%rax,2),%eax
cvtsi2sd %eax,%xmm5
movapd %xmm3,%xmm0
subsd  %xmm2,%xmm0
mulsd  %xmm6,%xmm0
addsd  %xmm3,%xmm0
movapd %xmm4,%xmm2
subsd  %xmm5,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm4,%xmm2
mulsd  %xmm6,%xmm2
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
jmp    1232 <BGBMID_SplinePolate+0x122>
js     121d <BGBMID_SplinePolate+0x10d>
movsd  0xe48(%rip),%xmm0
movsd  (%rsp),%xmm2
subsd  %xmm2,%xmm0
movswl (%r14),%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm0,%xmm1
movswl 0x2(%r14),%eax
jmp    120c <BGBMID_SplinePolate+0xfc>
add    $0xffffffff,%ebx
cmp    %ebx,%eax
jge    1223 <BGBMID_SplinePolate+0x113>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  (%rsp),%xmm2
subsd  %xmm0,%xmm2
movsd  0xe0d(%rip),%xmm0
subsd  %xmm2,%xmm0
mov    %eax,%eax
movswl (%r14,%rax,2),%ecx
cvtsi2sd %ecx,%xmm1
mulsd  %xmm0,%xmm1
movswl 0x2(%r14,%rax,2),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
jmp    1232 <BGBMID_SplinePolate+0x122>
movswl (%r14),%eax
jmp    122b <BGBMID_SplinePolate+0x11b>
movslq %ebx,%rax
movswl (%r14,%rax,2),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
add    $0x8,%rsp
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
define double @BGBMID_SplinePolate(i16* nocapture %0, i32 %1, double %2) {
%4 = tail call i32 @floor(double %2)
%5 = icmp slt i32 %4, 1
br i1 %5, label %6, label %20
6: ; preds = %3
%7 = icmp slt i32 %4, 0
br i1 %7, label %8, label %11
8: ; preds = %6
%9 = load i16, i16* %0
%10 = sitofp i16 %9 to double
br label %75
11: ; preds = %6
%12 = fsub double 1.000000e+00, %2
%13 = load i16, i16* %0
%14 = sitofp i16 %13 to double
%15 = getelementptr inbounds i16, i16* %0, i64 1
%16 = load i16, i16* %15
%17 = sitofp i16 %16 to double
%18 = fmul double %17, %2
%19 = tail call double @llvm.fmuladd.f64(double %14, double %12, double %18)
br label %75
20: ; preds = %3
%21 = add nsw i32 %1, -2
%22 = icmp slt i32 %4, %21
br i1 %22, label %46, label %23
23: ; preds = %20
%24 = add nsw i32 %1, -1
%25 = icmp slt i32 %4, %24
br i1 %25, label %31, label %26
26: ; preds = %23
%27 = sext i32 %24 to i64
%28 = getelementptr inbounds i16, i16* %0, i64 %27
%29 = load i16, i16* %28
%30 = sitofp i16 %29 to double
br label %75
31: ; preds = %23
%32 = sitofp i32 %4 to double
%33 = fsub double %2, %32
%34 = fsub double 1.000000e+00, %33
%35 = zext i32 %4 to i64
%36 = getelementptr inbounds i16, i16* %0, i64 %35
%37 = load i16, i16* %36
%38 = sitofp i16 %37 to double
%39 = add nuw nsw i32 %4, 1
%40 = zext i32 %39 to i64
%41 = getelementptr inbounds i16, i16* %0, i64 %40
%42 = load i16, i16* %41
%43 = sitofp i16 %42 to double
%44 = fmul double %33, %43
%45 = tail call double @llvm.fmuladd.f64(double %38, double %34, double %44)
br label %75
46: ; preds = %20
%47 = sitofp i32 %4 to double
%48 = fsub double %2, %47
%49 = fsub double 1.000000e+00, %48
%50 = add nsw i32 %4, -1
%51 = zext i32 %50 to i64
%52 = getelementptr inbounds i16, i16* %0, i64 %51
%53 = load i16, i16* %52
%54 = sitofp i16 %53 to double
%55 = zext i32 %4 to i64
%56 = getelementptr inbounds i16, i16* %0, i64 %55
%57 = load i16, i16* %56
%58 = sitofp i16 %57 to double
%59 = add nuw nsw i32 %4, 1
%60 = zext i32 %59 to i64
%61 = getelementptr inbounds i16, i16* %0, i64 %60
%62 = load i16, i16* %61
%63 = sitofp i16 %62 to double
%64 = add nuw nsw i32 %4, 2
%65 = zext i32 %64 to i64
%66 = getelementptr inbounds i16, i16* %0, i64 %65
%67 = load i16, i16* %66
%68 = sitofp i16 %67 to double
%69 = fsub double %58, %54
%70 = tail call double @llvm.fmuladd.f64(double %48, double %69, double %58)
%71 = fsub double %63, %68
%72 = tail call double @llvm.fmuladd.f64(double %49, double %71, double %63)
%73 = fmul double %48, %72
%74 = tail call double @llvm.fmuladd.f64(double %70, double %49, double %73)
br label %75
75: ; preds = %46, %31, %26, %11, %8
%76 = phi double [ %10, %8 ], [ %19, %11 ], [ %30, %26 ], [ %45, %31 ], [ %74, %46 ]
ret double %76
}