#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 */
scalar_t__ floor(double);
float LinInterp1D(float *arr, double index, long nelem) {
  long intpart;
  double fracpart;
  intpart = (long)floor(index);
  fracpart = index - intpart;
  if (intpart < 0) {
    return (arr[0]);
  } else if (intpart >= nelem - 1) {
    return (arr[nelem - 1]);
  } else {
    return (((1 - fracpart) * arr[intpart] + fracpart * arr[intpart + 1]) /
            2.0);
  }
}


x86 assembly (clang, -O0):
<LinInterp1D>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
movsd  -0x18(%rbp),%xmm0
call   1030 <floor@plt>
mov    %rax,-0x28(%rbp)
movsd  -0x18(%rbp),%xmm0
cvtsi2sdq -0x28(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
cmpq   $0x0,-0x28(%rbp)
jge    1164 <LinInterp1D+0x54>
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    11ea <LinInterp1D+0xda>
mov    -0x28(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    $0x1,%rcx
cmp    %rcx,%rax
jl     1194 <LinInterp1D+0x84>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    $0x1,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    11ea <LinInterp1D+0xda>
movsd  0xe6c(%rip),%xmm0
subsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  -0x30(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xe23(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<LinInterp1D>:
push   %r14
push   %rbx
push   %rax
mov    %rsi,%rbx
mov    %rdi,%r14
movsd  %xmm0,(%rsp)
call   1030 <floor@plt>
test   %rax,%rax
js     1181 <LinInterp1D+0x71>
lea    -0x1(%rbx),%rcx
cmp    %rcx,%rax
jge    1188 <LinInterp1D+0x78>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  (%rsp),%xmm3
subsd  %xmm0,%xmm3
movsd  0xeb5(%rip),%xmm0
subsd  %xmm3,%xmm0
movss  (%r14,%rax,4),%xmm1
movss  0x4(%r14,%rax,4),%xmm2
cvtss2sd %xmm1,%xmm1
mulsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm1,%xmm0
mulsd  0xe8d(%rip),%xmm0
cvtsd2ss %xmm0,%xmm0
jmp    118f <LinInterp1D+0x7f>
movss  (%r14),%xmm0
jmp    118f <LinInterp1D+0x7f>
movss  -0x4(%r14,%rbx,4),%xmm0
add    $0x8,%rsp
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define float @LinInterp1D(float* %0, double %1, i64 %2) {
%4 = alloca float
%5 = alloca float*
%6 = alloca double
%7 = alloca i64
%8 = alloca i64
%9 = alloca double
store float* %0, float** %5
store double %1, double* %6
store i64 %2, i64* %7
%10 = load double, double* %6
%11 = call i64 @floor(double %10)
store i64 %11, i64* %8
%12 = load double, double* %6
%13 = load i64, i64* %8
%14 = sitofp i64 %13 to double
%15 = fsub double %12, %14
store double %15, double* %9
%16 = load i64, i64* %8
%17 = icmp slt i64 %16, 0
br i1 %17, label %18, label %22
18: ; preds = %3
%19 = load float*, float** %5
%20 = getelementptr inbounds float, float* %19, i64 0
%21 = load float, float* %20
store float %21, float* %4
br label %52
22: ; preds = %3
%23 = load i64, i64* %8
%24 = load i64, i64* %7
%25 = sub nsw i64 %24, 1
%26 = icmp sge i64 %23, %25
br i1 %26, label %27, label %33
27: ; preds = %22
%28 = load float*, float** %5
%29 = load i64, i64* %7
%30 = sub nsw i64 %29, 1
%31 = getelementptr inbounds float, float* %28, i64 %30
%32 = load float, float* %31
store float %32, float* %4
br label %52
33: ; preds = %22
%34 = load double, double* %9
%35 = fsub double 1.000000e+00, %34
%36 = load float*, float** %5
%37 = load i64, i64* %8
%38 = getelementptr inbounds float, float* %36, i64 %37
%39 = load float, float* %38
%40 = fpext float %39 to double
%41 = load double, double* %9
%42 = load float*, float** %5
%43 = load i64, i64* %8
%44 = add nsw i64 %43, 1
%45 = getelementptr inbounds float, float* %42, i64 %44
%46 = load float, float* %45
%47 = fpext float %46 to double
%48 = fmul double %41, %47
%49 = call double @llvm.fmuladd.f64(double %35, double %40, double %48)
%50 = fdiv double %49, 2.000000e+00
%51 = fptrunc double %50 to float
store float %51, float* %4
br label %52
52: ; preds = %33, %27, %18
%53 = load float, float* %4
ret float %53
}

LLVM IR (-O2):
define float @LinInterp1D(float* nocapture %0, double %1, i64 %2) {
%4 = tail call i64 @floor(double %1)
%5 = sitofp i64 %4 to double
%6 = fsub double %1, %5
%7 = icmp slt i64 %4, 0
br i1 %7, label %8, label %10
8: ; preds = %3
%9 = load float, float* %0
br label %29
10: ; preds = %3
%11 = add nsw i64 %2, -1
%12 = icmp slt i64 %4, %11
br i1 %12, label %16, label %13
13: ; preds = %10
%14 = getelementptr inbounds float, float* %0, i64 %11
%15 = load float, float* %14
br label %29
16: ; preds = %10
%17 = fsub double 1.000000e+00, %6
%18 = getelementptr inbounds float, float* %0, i64 %4
%19 = load float, float* %18
%20 = fpext float %19 to double
%21 = add nuw nsw i64 %4, 1
%22 = getelementptr inbounds float, float* %0, i64 %21
%23 = load float, float* %22
%24 = fpext float %23 to double
%25 = fmul double %6, %24
%26 = tail call double @llvm.fmuladd.f64(double %17, double %20, double %25)
%27 = fmul double %26, 5.000000e-01
%28 = fptrunc double %27 to float
br label %29
29: ; preds = %16, %13, %8
%30 = phi float [ %9, %8 ], [ %15, %13 ], [ %28, %16 ]
ret float %30
}