#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);
double nlerp2(double *, double *);
double func_lerpf(int w, int h, float *f, double x, double y) {
  int x_f = (int)floor(x);
  int y_f = (int)floor(y);
  double fsub[4] = {f[y_f * w + x_f], f[y_f * w + x_f + 1],
                    f[(y_f + 1) * w + x_f], f[(y_f + 1) * w + x_f + 1]};
  double xin[2] = {x - x_f, y - y_f};
  return nlerp2(xin, fsub);
}


x86 assembly (clang, -O0):
<func_lerpf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
movsd  -0x18(%rbp),%xmm0
call   1040 <floor@plt>
mov    %eax,-0x24(%rbp)
movsd  -0x20(%rbp),%xmm0
call   1040 <floor@plt>
mov    %eax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
imul   -0x4(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
imul   -0x4(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x18(%rbp),%xmm0
cvtsi2sdl -0x24(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x28(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
lea    -0x60(%rbp),%rdi
lea    -0x50(%rbp),%rsi
call   1030 <nlerp2@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<func_lerpf>:
push   %rbp
push   %r14
push   %rbx
sub    $0x40,%rsp
movsd  %xmm1,(%rsp)
movsd  %xmm0,0x8(%rsp)
mov    %rdx,%r14
mov    %edi,%ebp
call   1040 <floor@plt>
mov    %rax,%rbx
movsd  (%rsp),%xmm0
call   1040 <floor@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    %eax,%ecx
imul   %ebp,%ecx
add    %ebx,%ecx
movslq %ecx,%rcx
movsd  (%r14,%rcx,4),%xmm1
cvtps2pd %xmm1,%xmm1
add    $0x1,%eax
movaps %xmm1,0x20(%rsp)
imul   %ebp,%eax
add    %ebx,%eax
cltq
movsd  (%r14,%rax,4),%xmm1
cvtps2pd %xmm1,%xmm1
movaps %xmm1,0x30(%rsp)
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
movsd  0x8(%rsp),%xmm2
subsd  %xmm1,%xmm2
movsd  %xmm2,0x10(%rsp)
movsd  (%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x18(%rsp)
lea    0x10(%rsp),%rdi
lea    0x20(%rsp),%rsi
call   1030 <nlerp2@plt>
add    $0x40,%rsp
pop    %rbx
pop    %r14
pop    %rbp
ret

LLVM IR (-O0):
define double @func_lerpf(i32 %0, i32 %1, float* %2, double %3, double %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca float*
%9 = alloca double
%10 = alloca double
%11 = alloca i32
%12 = alloca i32
%13 = alloca [4 x double]
%14 = alloca [2 x double]
store i32 %0, i32* %6
store i32 %1, i32* %7
store float* %2, float** %8
store double %3, double* %9
store double %4, double* %10
%15 = load double, double* %9
%16 = call i64 @floor(double %15)
%17 = trunc i64 %16 to i32
store i32 %17, i32* %11
%18 = load double, double* %10
%19 = call i64 @floor(double %18)
%20 = trunc i64 %19 to i32
store i32 %20, i32* %12
%21 = getelementptr inbounds [4 x double], [4 x double]* %13, i64 0, i64 0
%22 = load float*, float** %8
%23 = load i32, i32* %12
%24 = load i32, i32* %6
%25 = mul nsw i32 %23, %24
%26 = load i32, i32* %11
%27 = add nsw i32 %25, %26
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds float, float* %22, i64 %28
%30 = load float, float* %29
%31 = fpext float %30 to double
store double %31, double* %21
%32 = getelementptr inbounds double, double* %21, i64 1
%33 = load float*, float** %8
%34 = load i32, i32* %12
%35 = load i32, i32* %6
%36 = mul nsw i32 %34, %35
%37 = load i32, i32* %11
%38 = add nsw i32 %36, %37
%39 = add nsw i32 %38, 1
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds float, float* %33, i64 %40
%42 = load float, float* %41
%43 = fpext float %42 to double
store double %43, double* %32
%44 = getelementptr inbounds double, double* %32, i64 1
%45 = load float*, float** %8
%46 = load i32, i32* %12
%47 = add nsw i32 %46, 1
%48 = load i32, i32* %6
%49 = mul nsw i32 %47, %48
%50 = load i32, i32* %11
%51 = add nsw i32 %49, %50
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds float, float* %45, i64 %52
%54 = load float, float* %53
%55 = fpext float %54 to double
store double %55, double* %44
%56 = getelementptr inbounds double, double* %44, i64 1
%57 = load float*, float** %8
%58 = load i32, i32* %12
%59 = add nsw i32 %58, 1
%60 = load i32, i32* %6
%61 = mul nsw i32 %59, %60
%62 = load i32, i32* %11
%63 = add nsw i32 %61, %62
%64 = add nsw i32 %63, 1
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds float, float* %57, i64 %65
%67 = load float, float* %66
%68 = fpext float %67 to double
store double %68, double* %56
%69 = getelementptr inbounds [2 x double], [2 x double]* %14, i64 0, i64 0
%70 = load double, double* %9
%71 = load i32, i32* %11
%72 = sitofp i32 %71 to double
%73 = fsub double %70, %72
store double %73, double* %69
%74 = getelementptr inbounds double, double* %69, i64 1
%75 = load double, double* %10
%76 = load i32, i32* %12
%77 = sitofp i32 %76 to double
%78 = fsub double %75, %77
store double %78, double* %74
%79 = getelementptr inbounds [2 x double], [2 x double]* %14, i64 0, i64 0
%80 = getelementptr inbounds [4 x double], [4 x double]* %13, i64 0, i64 0
%81 = call double @nlerp2(double* %79, double* %80)
ret double %81
}

LLVM IR (-O2):
define double @func_lerpf(i32 %0, i32 %1, float* nocapture %2, double %3, double %4) {
%6 = alloca [4 x double]
%7 = alloca [2 x double]
%8 = tail call i64 @floor(double %3)
%9 = trunc i64 %8 to i32
%10 = tail call i64 @floor(double %4)
%11 = trunc i64 %10 to i32
%12 = bitcast [4 x double]* %6 to i8*
%13 = getelementptr inbounds [4 x double], [4 x double]* %6, i64 0, i64 0
%14 = mul nsw i32 %11, %0
%15 = add nsw i32 %14, %9
%16 = sext i32 %15 to i64
%17 = getelementptr inbounds float, float* %2, i64 %16
%18 = load float, float* %17
%19 = add nsw i32 %15, 1
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds float, float* %2, i64 %20
%22 = load float, float* %21
%23 = insertelement <2 x float> poison, float %18, i64 0
%24 = insertelement <2 x float> %23, float %22, i64 1
%25 = fpext <2 x float> %24 to <2 x double>
%26 = bitcast [4 x double]* %6 to <2 x double>*
store <2 x double> %25, <2 x double>* %26
%27 = getelementptr inbounds [4 x double], [4 x double]* %6, i64 0, i64 2
%28 = add nsw i32 %11, 1
%29 = mul nsw i32 %28, %0
%30 = add nsw i32 %29, %9
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds float, float* %2, i64 %31
%33 = load float, float* %32
%34 = add nsw i32 %30, 1
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds float, float* %2, i64 %35
%37 = load float, float* %36
%38 = insertelement <2 x float> poison, float %33, i64 0
%39 = insertelement <2 x float> %38, float %37, i64 1
%40 = fpext <2 x float> %39 to <2 x double>
%41 = bitcast double* %27 to <2 x double>*
store <2 x double> %40, <2 x double>* %41
%42 = bitcast [2 x double]* %7 to i8*
%43 = getelementptr inbounds [2 x double], [2 x double]* %7, i64 0, i64 0
%44 = sitofp i32 %9 to double
%45 = fsub double %3, %44
store double %45, double* %43
%46 = getelementptr inbounds [2 x double], [2 x double]* %7, i64 0, i64 1
%47 = sitofp i32 %11 to double
%48 = fsub double %4, %47
store double %48, double* %46
%49 = call double @nlerp2(double* %43, double* %13)
ret double %49
}