#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 */
float fmod(float, double);
float trilinear_interp(float *c_n, int c_n_x, float *v_n, int v_n_x,
                       float *pos_radar, int p_r_x, float radar_bin_height) {
  float **data;
  int i;
  float v[4];
  float interp = 0;
  float diff_i = 0;
  float diff_j = 0;
  float pos_i = 0.0;
  float pos_j = 0.0;
  int j = 0;
  for (j = 0; j < c_n_x / 2; j++) {
    v[j] = v_n[2 * j + 1] - (v_n[2 * j + 1] - v_n[2 * j]) /
                                (c_n[2 * j] - c_n[2 * j + 1]) *
                                (radar_bin_height - c_n[2 * j + 1]);
  }
  pos_i = fmod(pos_radar[1], 1.0);
  pos_j = fmod(pos_radar[0], 1.0);
  diff_i = 1.0 - pos_i;
  diff_j = 1.0 - pos_j;
  interp = diff_i * diff_j * v[0] + pos_i * v[2] * diff_j +
           diff_i * v[1] * pos_j + pos_i * pos_j * v[3];
  return interp;
}


x86 assembly (clang, -O0):
<trilinear_interp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movss  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x5c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
movl   $0x0,-0x68(%rbp)
movl   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%eax
mov    %eax,-0x6c(%rbp)
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x6c(%rbp),%eax
cmp    %ecx,%eax
jge    1235 <trilinear_interp+0x125>
mov    -0x18(%rbp),%rax
mov    -0x68(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x68(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x68(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm2
divss  %xmm2,%xmm0
movss  -0x30(%rbp),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movslq -0x68(%rbp),%rax
movss  %xmm0,-0x50(%rbp,%rax,4)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    1169 <trilinear_interp+0x59>
mov    -0x28(%rbp),%rax
movss  0x4(%rax),%xmm0
movsd  0xdba(%rip),%xmm1
call   1030 <fmod@plt>
movss  %xmm0,-0x60(%rbp)
mov    -0x28(%rbp),%rax
movss  (%rax),%xmm0
movsd  0xda0(%rip),%xmm1
call   1030 <fmod@plt>
movss  %xmm0,-0x64(%rbp)
movss  -0x60(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd85(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x64(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd67(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x5c(%rbp)
movss  -0x58(%rbp),%xmm0
mulss  -0x5c(%rbp),%xmm0
movss  -0x50(%rbp),%xmm2
movss  -0x60(%rbp),%xmm1
mulss  -0x48(%rbp),%xmm1
mulss  -0x5c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x58(%rbp),%xmm1
mulss  -0x4c(%rbp),%xmm1
movss  -0x64(%rbp),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x60(%rbp),%xmm0
mulss  -0x64(%rbp),%xmm0
movss  -0x44(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x54(%rbp),%xmm0
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<trilinear_interp>:
push   %r14
push   %rbx
sub    $0x38,%rsp
mov    %r8,%r14
cmp    $0x2,%esi
jl     121f <trilinear_interp+0x10f>
mov    %esi,%eax
shr    %eax
cmp    $0x8,%esi
jae    1133 <trilinear_interp+0x23>
xor    %ecx,%ecx
jmp    11b0 <trilinear_interp+0xa0>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
and    $0xfffffff8,%esi
add    %rsi,%rsi
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movups (%rdx,%rbx,2),%xmm2
movups 0x10(%rdx,%rbx,2),%xmm3
movaps %xmm2,%xmm4
shufps $0x88,%xmm3,%xmm4
shufps $0xdd,%xmm3,%xmm2
movaps %xmm2,%xmm3
subps  %xmm4,%xmm3
movups (%rdi,%rbx,2),%xmm4
movups 0x10(%rdi,%rbx,2),%xmm5
movaps %xmm4,%xmm6
shufps $0x88,%xmm5,%xmm6
shufps $0xdd,%xmm5,%xmm4
subps  %xmm4,%xmm6
divps  %xmm6,%xmm3
movaps %xmm1,%xmm5
subps  %xmm4,%xmm5
mulps  %xmm3,%xmm5
subps  %xmm5,%xmm2
movaps %xmm2,0x20(%rsp,%rbx,1)
add    $0x10,%rbx
cmp    %rbx,%rsi
jne    1150 <trilinear_interp+0x40>
cmp    %rax,%rcx
je     11ef <trilinear_interp+0xdf>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  0x4(%rdx,%rcx,8),%xmm1
movaps %xmm1,%xmm2
subss  (%rdx,%rcx,8),%xmm2
movss  (%rdi,%rcx,8),%xmm3
movss  0x4(%rdi,%rcx,8),%xmm4
subss  %xmm4,%xmm3
divss  %xmm3,%xmm2
movaps %xmm0,%xmm3
subss  %xmm4,%xmm3
mulss  %xmm2,%xmm3
subss  %xmm3,%xmm1
movss  %xmm1,0x20(%rsp,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    11b0 <trilinear_interp+0xa0>
movss  0x20(%rsp),%xmm0
movss  %xmm0,0xc(%rsp)
movss  0x24(%rsp),%xmm0
movss  %xmm0,0x10(%rsp)
movss  0x28(%rsp),%xmm0
movss  %xmm0,0x14(%rsp)
movss  0x2c(%rsp),%xmm0
movss  %xmm0,0x18(%rsp)
movss  0x4(%r14),%xmm0
movsd  0xdd3(%rip),%xmm1
call   1030 <fmod@plt>
movss  %xmm0,0x1c(%rsp)
movss  (%r14),%xmm0
movsd  0xdbb(%rip),%xmm1
call   1030 <fmod@plt>
movss  0xdb6(%rip),%xmm1
movaps %xmm1,%xmm2
movss  0x1c(%rsp),%xmm5
subss  %xmm5,%xmm2
subss  %xmm0,%xmm1
movaps %xmm2,%xmm3
mulss  %xmm1,%xmm3
mulss  0xc(%rsp),%xmm3
movss  0x14(%rsp),%xmm4
mulss  %xmm5,%xmm4
mulss  %xmm1,%xmm4
addss  %xmm3,%xmm4
mulss  0x10(%rsp),%xmm2
mulss  %xmm0,%xmm2
addss  %xmm4,%xmm2
mulss  %xmm5,%xmm0
mulss  0x18(%rsp),%xmm0
addss  %xmm2,%xmm0
add    $0x38,%rsp
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
define float @trilinear_interp(float* nocapture %0, i32 %1, float* nocapture %2, i32 %3, float* nocapture %4, i32 %5, float %6) {
%8 = alloca [4 x float]
%9 = bitcast [4 x float]* %8 to i8*
%10 = icmp sgt i32 %1, 1
br i1 %10, label %11, label %78
11: ; preds = %7
%12 = lshr i32 %1, 1
%13 = zext i32 %12 to i64
%14 = icmp ult i32 %1, 8
br i1 %14, label %46, label %15
15: ; preds = %11
%16 = and i64 %13, 2147483644
%17 = insertelement <4 x float> poison, float %6, i64 0
%18 = shufflevector <4 x float> %17, <4 x float> poison, <4 x i32> zeroinitializer
%19 = getelementptr inbounds float, float* %2, i64 -1
br label %20
20: ; preds = %20, %15
%21 = phi i64 [ 0, %15 ], [ %42, %20 ]
%22 = shl nuw nsw i64 %21, 1
%23 = or i64 %22, 1
%24 = getelementptr inbounds float, float* %19, i64 %23
%25 = bitcast float* %24 to <8 x float>*
%26 = load <8 x float>, <8 x float>* %25
%27 = shufflevector <8 x float> %26, <8 x float> poison, <4 x i32> <i32 0, i32 2, i32 4, i32 6>
%28 = shufflevector <8 x float> %26, <8 x float> poison, <4 x i32> <i32 1, i32 3, i32 5, i32 7>
%29 = fsub <4 x float> %28, %27
%30 = getelementptr inbounds float, float* %0, i64 %22
%31 = bitcast float* %30 to <8 x float>*
%32 = load <8 x float>, <8 x float>* %31
%33 = shufflevector <8 x float> %32, <8 x float> poison, <4 x i32> <i32 0, i32 2, i32 4, i32 6>
%34 = shufflevector <8 x float> %32, <8 x float> poison, <4 x i32> <i32 1, i32 3, i32 5, i32 7>
%35 = fsub <4 x float> %33, %34
%36 = fsub <4 x float> %18, %34
%37 = fneg <4 x float> %29
%38 = fdiv <4 x float> %37, %35
%39 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %38, <4 x float> %36, <4 x float> %28)
%40 = getelementptr inbounds [4 x float], [4 x float]* %8, i64 0, i64 %21
%41 = bitcast float* %40 to <4 x float>*
store <4 x float> %39, <4 x float>* %41
%42 = add nuw i64 %21, 4
%43 = icmp eq i64 %42, %16
br i1 %43, label %44, label %20
44: ; preds = %20
%45 = icmp eq i64 %16, %13
br i1 %45, label %69, label %46
46: ; preds = %11, %44
%47 = phi i64 [ 0, %11 ], [ %16, %44 ]
br label %48
48: ; preds = %46, %48
%49 = phi i64 [ %67, %48 ], [ %47, %46 ]
%50 = shl nuw nsw i64 %49, 1
%51 = or i64 %50, 1
%52 = getelementptr inbounds float, float* %2, i64 %51
%53 = load float, float* %52
%54 = getelementptr inbounds float, float* %2, i64 %50
%55 = load float, float* %54
%56 = fsub float %53, %55
%57 = getelementptr inbounds float, float* %0, i64 %50
%58 = load float, float* %57
%59 = getelementptr inbounds float, float* %0, i64 %51
%60 = load float, float* %59
%61 = fsub float %58, %60
%62 = fsub float %6, %60
%63 = fneg float %56
%64 = fdiv float %63, %61
%65 = tail call float @llvm.fmuladd.f32(float %64, float %62, float %53)
%66 = getelementptr inbounds [4 x float], [4 x float]* %8, i64 0, i64 %49
store float %65, float* %66
%67 = add nuw nsw i64 %49, 1
%68 = icmp eq i64 %67, %13
br i1 %68, label %69, label %48
69: ; preds = %48, %44
%70 = getelementptr inbounds [4 x float], [4 x float]* %8, i64 0, i64 0
%71 = load float, float* %70
%72 = getelementptr inbounds [4 x float], [4 x float]* %8, i64 0, i64 2
%73 = load float, float* %72
%74 = getelementptr inbounds [4 x float], [4 x float]* %8, i64 0, i64 1
%75 = load float, float* %74
%76 = getelementptr inbounds [4 x float], [4 x float]* %8, i64 0, i64 3
%77 = load float, float* %76
br label %78
78: ; preds = %69, %7
%79 = phi float [ %77, %69 ], [ undef, %7 ]
%80 = phi float [ %75, %69 ], [ undef, %7 ]
%81 = phi float [ %73, %69 ], [ undef, %7 ]
%82 = phi float [ %71, %69 ], [ undef, %7 ]
%83 = getelementptr inbounds float, float* %4, i64 1
%84 = load float, float* %83
%85 = tail call float @fmod(float %84, double 1.000000e+00)
%86 = load float, float* %4
%87 = tail call float @fmod(float %86, double 1.000000e+00)
%88 = fsub float 1.000000e+00, %85
%89 = fsub float 1.000000e+00, %87
%90 = fmul float %88, %89
%91 = fmul float %85, %81
%92 = fmul float %89, %91
%93 = tail call float @llvm.fmuladd.f32(float %90, float %82, float %92)
%94 = fmul float %88, %80
%95 = tail call float @llvm.fmuladd.f32(float %94, float %87, float %93)
%96 = fmul float %85, %87
%97 = tail call float @llvm.fmuladd.f32(float %96, float %79, float %95)
ret float %97
}