int bilinear_interpolate(the_image, x, y, rows, cols)
double x, y;
long cols, rows;
short **the_image;
{
  double fraction_x, fraction_y, one_minus_x, one_minus_y, tmp_double;
  int ceil_x, ceil_y, floor_x, floor_y;
  short p1, p2, p3, result = 150;
  if (x < 0.0 || x >= (double)(cols - 1) || y < 0.0 || y >= (double)(rows - 1))
    return (result);
  tmp_double = floor(x);
  floor_x = tmp_double;
  tmp_double = floor(y);
  floor_y = tmp_double;
  tmp_double = ceil(x);
  ceil_x = tmp_double;
  tmp_double = ceil(y);
  ceil_y = tmp_double;
  fraction_x = x - floor(x);
  fraction_y = y - floor(y);
  one_minus_x = 1.0 - fraction_x;
  one_minus_y = 1.0 - fraction_y;
  tmp_double = one_minus_x * (double)(the_image[floor_y][floor_x]) +
               fraction_x * (double)(the_image[floor_y][ceil_x]);
  p1 = tmp_double;
  tmp_double = one_minus_x * (double)(the_image[ceil_y][floor_x]) +
               fraction_x * (double)(the_image[ceil_y][ceil_x]);
  p2 = tmp_double;
  tmp_double = one_minus_y * (double)(p1) + fraction_y * (double)(p2);
  p3 = tmp_double;
  return (p3);
}


x86 assembly (clang, -O0):
<bilinear_interpolate>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    %rsi,-0x28(%rbp)
mov    %rdx,-0x30(%rbp)
movw   $0x96,-0x70(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x18(%rbp),%xmm0
ja     119b <bilinear_interpolate+0x7b>
movsd  -0x18(%rbp),%xmm0
mov    -0x30(%rbp),%rax
sub    $0x1,%rax
cvtsi2sd %rax,%xmm1
ucomisd %xmm1,%xmm0
jae    119b <bilinear_interpolate+0x7b>
xorps  %xmm0,%xmm0
ucomisd -0x20(%rbp),%xmm0
ja     119b <bilinear_interpolate+0x7b>
movsd  -0x20(%rbp),%xmm0
mov    -0x28(%rbp),%rax
sub    $0x1,%rax
cvtsi2sd %rax,%xmm1
ucomisd %xmm1,%xmm0
jb     11a7 <bilinear_interpolate+0x87>
movswl -0x70(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1339 <bilinear_interpolate+0x219>
movsd  -0x18(%rbp),%xmm0
call   1040 <floor@plt>
movsd  %xmm0,-0x58(%rbp)
cvttsd2si -0x58(%rbp),%eax
mov    %eax,-0x64(%rbp)
movsd  -0x20(%rbp),%xmm0
call   1040 <floor@plt>
movsd  %xmm0,-0x58(%rbp)
cvttsd2si -0x58(%rbp),%eax
mov    %eax,-0x68(%rbp)
movsd  -0x18(%rbp),%xmm0
call   1030 <ceil@plt>
movsd  %xmm0,-0x58(%rbp)
cvttsd2si -0x58(%rbp),%eax
mov    %eax,-0x5c(%rbp)
movsd  -0x20(%rbp),%xmm0
call   1030 <ceil@plt>
movsd  %xmm0,-0x58(%rbp)
cvttsd2si -0x58(%rbp),%eax
mov    %eax,-0x60(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
call   1040 <floor@plt>
movaps %xmm0,%xmm1
movsd  -0x80(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
call   1040 <floor@plt>
movaps %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm2
movsd  0xdb0(%rip),%xmm0
movaps %xmm0,%xmm1
subsd  %xmm2,%xmm1
movsd  %xmm1,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x68(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x64(%rbp),%rcx
movswl (%rax,%rcx,2),%ecx
cvtsi2sd %ecx,%xmm2
movsd  -0x38(%rbp),%xmm1
movslq -0x5c(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
cvttsd2si -0x58(%rbp),%eax
mov    %ax,-0x6a(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x60(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x64(%rbp),%rcx
movswl (%rax,%rcx,2),%ecx
cvtsi2sd %ecx,%xmm2
movsd  -0x38(%rbp),%xmm1
movslq -0x5c(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
cvttsd2si -0x58(%rbp),%eax
mov    %ax,-0x6c(%rbp)
movsd  -0x50(%rbp),%xmm0
movswl -0x6a(%rbp),%eax
cvtsi2sd %eax,%xmm2
movsd  -0x40(%rbp),%xmm1
movswl -0x6c(%rbp),%eax
cvtsi2sd %eax,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
cvttsd2si %xmm0,%eax
mov    %ax,-0x6e(%rbp)
movswl -0x6e(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bilinear_interpolate>:
mov    $0x96,%eax
xorpd  %xmm2,%xmm2
ucomisd %xmm0,%xmm2
ja     127f <bilinear_interpolate+0x15f>
ucomisd %xmm1,%xmm2
ja     127f <bilinear_interpolate+0x15f>
add    $0xffffffffffffffff,%rdx
xorps  %xmm2,%xmm2
cvtsi2sd %rdx,%xmm2
ucomisd %xmm2,%xmm0
jae    127f <bilinear_interpolate+0x15f>
add    $0xffffffffffffffff,%rsi
xorps  %xmm2,%xmm2
cvtsi2sd %rsi,%xmm2
ucomisd %xmm2,%xmm1
jae    127f <bilinear_interpolate+0x15f>
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x28,%rsp
mov    %rdi,%rbx
movsd  %xmm1,0x8(%rsp)
movsd  %xmm0,0x10(%rsp)
call   1040 <floor@plt>
movsd  %xmm0,0x20(%rsp)
cvttsd2si %xmm0,%r14d
movsd  0x8(%rsp),%xmm0
call   1040 <floor@plt>
movsd  %xmm0,0x18(%rsp)
cvttsd2si %xmm0,%r15d
movsd  0x10(%rsp),%xmm0
call   1030 <ceil@plt>
cvttsd2si %xmm0,%ebp
movsd  0x8(%rsp),%xmm0
call   1030 <ceil@plt>
cvttsd2si %xmm0,%eax
movsd  0x10(%rsp),%xmm2
subsd  0x20(%rsp),%xmm2
movsd  0x8(%rsp),%xmm4
subsd  0x18(%rsp),%xmm4
movsd  0xe1a(%rip),%xmm0
movapd %xmm0,%xmm1
subsd  %xmm2,%xmm1
movapd %xmm2,%xmm5
subsd  %xmm4,%xmm0
movslq %r15d,%rcx
mov    (%rbx,%rcx,8),%rcx
movslq %r14d,%rdx
movswl (%rcx,%rdx,2),%esi
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
movslq %ebp,%rsi
movswl (%rcx,%rsi,2),%ecx
cvtsi2sd %ecx,%xmm3
mulsd  %xmm5,%xmm3
mulsd  %xmm1,%xmm2
addsd  %xmm3,%xmm2
cvttsd2si %xmm2,%ecx
cltq
mov    (%rbx,%rax,8),%rax
movswl (%rax,%rdx,2),%edx
xorps  %xmm2,%xmm2
cvtsi2sd %edx,%xmm2
movswl (%rax,%rsi,2),%eax
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
mulsd  %xmm5,%xmm3
mulsd  %xmm1,%xmm2
addsd  %xmm3,%xmm2
cvttsd2si %xmm2,%eax
movswl %cx,%ecx
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
mulsd  %xmm0,%xmm1
cwtl
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  %xmm4,%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
cwtl
add    $0x28,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

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