int int_px_interp(int **t, double x_, double y_) {
  int x = x_;
  int y = y_;
  return (x_ - x) * (y_ - y) * t[x + 1][y + 1] +
         (x_ - x) * (1 - y_ + y) * t[x + 1][y] +
         (1 - x_ + x) * (y_ - y) * t[x][y + 1] +
         (1 - x_ + x) * (1 - y_ + y) * t[x][y];
}


x86 assembly (clang, -O0):
<int_px_interp>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
cvttsd2si -0x10(%rbp),%eax
mov    %eax,-0x1c(%rbp)
cvttsd2si -0x18(%rbp),%eax
mov    %eax,-0x20(%rbp)
movsd  -0x10(%rbp),%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  -0x18(%rbp),%xmm1
cvtsi2sdl -0x20(%rbp),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm2
movsd  -0x10(%rbp),%xmm1
cvtsi2sdl -0x1c(%rbp),%xmm3
subsd  %xmm3,%xmm1
movsd  0xe89(%rip),%xmm3
subsd  -0x18(%rbp),%xmm3
cvtsi2sdl -0x20(%rbp),%xmm4
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xe49(%rip),%xmm1
subsd  -0x10(%rbp),%xmm1
cvtsi2sdl -0x1c(%rbp),%xmm2
addsd  %xmm2,%xmm1
movsd  -0x18(%rbp),%xmm2
cvtsi2sdl -0x20(%rbp),%xmm3
subsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  0xdff(%rip),%xmm0
subsd  -0x10(%rbp),%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm2
addsd  %xmm2,%xmm0
movsd  0xde9(%rip),%xmm2
subsd  -0x18(%rbp),%xmm2
cvtsi2sdl -0x20(%rbp),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<int_px_interp>:
unpcklpd %xmm1,%xmm0
cvttpd2dq %xmm0,%xmm1
cvtdq2pd %xmm1,%xmm2
movapd 0xeec(%rip),%xmm3
subpd  %xmm0,%xmm3
subpd  %xmm2,%xmm0
movapd %xmm0,%xmm4
unpckhpd %xmm0,%xmm4
movapd %xmm0,%xmm5
mulsd  %xmm4,%xmm5
movd   %xmm1,%eax
cltq
mov    (%rdi,%rax,8),%rcx
mov    0x8(%rdi,%rax,8),%rax
pshufd $0x55,%xmm1,%xmm1
movd   %xmm1,%edx
movslq %edx,%rdx
xorps  %xmm1,%xmm1
cvtsi2sdl 0x4(%rax,%rdx,4),%xmm1
mulsd  %xmm5,%xmm1
addpd  %xmm2,%xmm3
movapd %xmm3,%xmm2
unpckhpd %xmm3,%xmm2
mulsd  %xmm2,%xmm0
xorps  %xmm5,%xmm5
cvtsi2sdl (%rax,%rdx,4),%xmm5
mulsd  %xmm0,%xmm5
addsd  %xmm1,%xmm5
mulsd  %xmm3,%xmm4
xorps  %xmm0,%xmm0
cvtsi2sdl 0x4(%rcx,%rdx,4),%xmm0
mulsd  %xmm4,%xmm0
addsd  %xmm5,%xmm0
mulsd  %xmm3,%xmm2
xorps  %xmm1,%xmm1
cvtsi2sdl (%rcx,%rdx,4),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
cvttsd2si %xmm1,%eax
ret

LLVM IR (-O0):
define i32 @int_px_interp(i32** %0, double %1, double %2) {
%4 = alloca i32**
%5 = alloca double
%6 = alloca double
%7 = alloca i32
%8 = alloca i32
store i32** %0, i32*** %4
store double %1, double* %5
store double %2, double* %6
%9 = load double, double* %5
%10 = fptosi double %9 to i32
store i32 %10, i32* %7
%11 = load double, double* %6
%12 = fptosi double %11 to i32
store i32 %12, i32* %8
%13 = load double, double* %5
%14 = load i32, i32* %7
%15 = sitofp i32 %14 to double
%16 = fsub double %13, %15
%17 = load double, double* %6
%18 = load i32, i32* %8
%19 = sitofp i32 %18 to double
%20 = fsub double %17, %19
%21 = fmul double %16, %20
%22 = load i32**, i32*** %4
%23 = load i32, i32* %7
%24 = add nsw i32 %23, 1
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds i32*, i32** %22, i64 %25
%27 = load i32*, i32** %26
%28 = load i32, i32* %8
%29 = add nsw i32 %28, 1
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds i32, i32* %27, i64 %30
%32 = load i32, i32* %31
%33 = sitofp i32 %32 to double
%34 = load double, double* %5
%35 = load i32, i32* %7
%36 = sitofp i32 %35 to double
%37 = fsub double %34, %36
%38 = load double, double* %6
%39 = fsub double 1.000000e+00, %38
%40 = load i32, i32* %8
%41 = sitofp i32 %40 to double
%42 = fadd double %39, %41
%43 = fmul double %37, %42
%44 = load i32**, i32*** %4
%45 = load i32, i32* %7
%46 = add nsw i32 %45, 1
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds i32*, i32** %44, i64 %47
%49 = load i32*, i32** %48
%50 = load i32, i32* %8
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds i32, i32* %49, i64 %51
%53 = load i32, i32* %52
%54 = sitofp i32 %53 to double
%55 = fmul double %43, %54
%56 = call double @llvm.fmuladd.f64(double %21, double %33, double %55)
%57 = load double, double* %5
%58 = fsub double 1.000000e+00, %57
%59 = load i32, i32* %7
%60 = sitofp i32 %59 to double
%61 = fadd double %58, %60
%62 = load double, double* %6
%63 = load i32, i32* %8
%64 = sitofp i32 %63 to double
%65 = fsub double %62, %64
%66 = fmul double %61, %65
%67 = load i32**, i32*** %4
%68 = load i32, i32* %7
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds i32*, i32** %67, i64 %69
%71 = load i32*, i32** %70
%72 = load i32, i32* %8
%73 = add nsw i32 %72, 1
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds i32, i32* %71, i64 %74
%76 = load i32, i32* %75
%77 = sitofp i32 %76 to double
%78 = call double @llvm.fmuladd.f64(double %66, double %77, double %56)
%79 = load double, double* %5
%80 = fsub double 1.000000e+00, %79
%81 = load i32, i32* %7
%82 = sitofp i32 %81 to double
%83 = fadd double %80, %82
%84 = load double, double* %6
%85 = fsub double 1.000000e+00, %84
%86 = load i32, i32* %8
%87 = sitofp i32 %86 to double
%88 = fadd double %85, %87
%89 = fmul double %83, %88
%90 = load i32**, i32*** %4
%91 = load i32, i32* %7
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds i32*, i32** %90, i64 %92
%94 = load i32*, i32** %93
%95 = load i32, i32* %8
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds i32, i32* %94, i64 %96
%98 = load i32, i32* %97
%99 = sitofp i32 %98 to double
%100 = call double @llvm.fmuladd.f64(double %89, double %99, double %78)
%101 = fptosi double %100 to i32
ret i32 %101
}

LLVM IR (-O2):
define i32 @int_px_interp(i32** nocapture %0, double %1, double %2) {
%4 = insertelement <2 x double> poison, double %1, i64 0
%5 = insertelement <2 x double> %4, double %2, i64 1
%6 = fptosi <2 x double> %5 to <2 x i32>
%7 = sitofp <2 x i32> %6 to <2 x double>
%8 = fsub <2 x double> %5, %7
%9 = extractelement <2 x double> %8, i64 0
%10 = extractelement <2 x double> %8, i64 1
%11 = fmul double %9, %10
%12 = extractelement <2 x i32> %6, i64 0
%13 = add nsw i32 %12, 1
%14 = sext i32 %13 to i64
%15 = getelementptr inbounds i32*, i32** %0, i64 %14
%16 = load i32*, i32** %15
%17 = extractelement <2 x i32> %6, i64 1
%18 = add nsw i32 %17, 1
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds i32, i32* %16, i64 %19
%21 = load i32, i32* %20
%22 = sitofp i32 %21 to double
%23 = fsub <2 x double> <double 1.000000e+00, double 1.000000e+00>, %5
%24 = fadd <2 x double> %23, %7
%25 = extractelement <2 x double> %24, i64 1
%26 = fmul double %9, %25
%27 = sext i32 %17 to i64
%28 = getelementptr inbounds i32, i32* %16, i64 %27
%29 = load i32, i32* %28
%30 = sitofp i32 %29 to double
%31 = fmul double %26, %30
%32 = tail call double @llvm.fmuladd.f64(double %11, double %22, double %31)
%33 = extractelement <2 x double> %24, i64 0
%34 = fmul double %33, %10
%35 = sext i32 %12 to i64
%36 = getelementptr inbounds i32*, i32** %0, i64 %35
%37 = load i32*, i32** %36
%38 = getelementptr inbounds i32, i32* %37, i64 %19
%39 = load i32, i32* %38
%40 = sitofp i32 %39 to double
%41 = tail call double @llvm.fmuladd.f64(double %34, double %40, double %32)
%42 = fmul double %33, %25
%43 = getelementptr inbounds i32, i32* %37, i64 %27
%44 = load i32, i32* %43
%45 = sitofp i32 %44 to double
%46 = tail call double @llvm.fmuladd.f64(double %42, double %45, double %41)
%47 = fptosi double %46 to i32
ret i32 %47
}