float getsample(float *fx, int w, int h, int pd, int i, int j, int l) {
  if (i < 0 || i >= w || j < 0 || j >= h || l < 0 || l >= pd)
    return 0;
  float(*x)[w][pd] = (void *)fx;
  return x[j][i][l];
}


x86 assembly (clang, -O0):
<getsample>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%eax
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
cmpl   $0x0,-0x20(%rbp)
jl     115e <getsample+0x5e>
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    115e <getsample+0x5e>
cmpl   $0x0,-0x24(%rbp)
jl     115e <getsample+0x5e>
mov    -0x24(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    115e <getsample+0x5e>
cmpl   $0x0,0x10(%rbp)
jl     115e <getsample+0x5e>
mov    0x10(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     116b <getsample+0x6b>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    11b1 <getsample+0xb1>
mov    -0x14(%rbp),%eax
mov    %eax,%esi
mov    -0x1c(%rbp),%eax
mov    %eax,%edx
mov    -0x10(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   %rdx,%rsi
imul   %rsi,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq 0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<getsample>:
xorps  %xmm0,%xmm0
test   %r8d,%r8d
js     1148 <getsample+0x48>
cmp    %esi,%r8d
jge    1148 <getsample+0x48>
test   %r9d,%r9d
js     1148 <getsample+0x48>
mov    0x8(%rsp),%eax
cmp    %ecx,%eax
jge    1148 <getsample+0x48>
cmp    %edx,%r9d
jge    1148 <getsample+0x48>
test   %eax,%eax
js     1148 <getsample+0x48>
mov    %esi,%edx
mov    %ecx,%ecx
mov    %r9d,%esi
imul   %rcx,%rdx
imul   %rsi,%rdx
lea    (%rdi,%rdx,4),%rdx
mov    %r8d,%esi
imul   %rcx,%rsi
lea    (%rdx,%rsi,4),%rcx
mov    %eax,%eax
movss  (%rcx,%rax,4),%xmm0
ret

LLVM IR (-O0):
define float @getsample(float* %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6) {
%8 = alloca float
%9 = alloca float*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca float*
store float* %0, float** %9
store i32 %1, i32* %10
store i32 %2, i32* %11
store i32 %3, i32* %12
store i32 %4, i32* %13
store i32 %5, i32* %14
store i32 %6, i32* %15
%17 = load i32, i32* %13
%18 = icmp slt i32 %17, 0
br i1 %18, label %37, label %19
19: ; preds = %7
%20 = load i32, i32* %13
%21 = load i32, i32* %10
%22 = icmp sge i32 %20, %21
br i1 %22, label %37, label %23
23: ; preds = %19
%24 = load i32, i32* %14
%25 = icmp slt i32 %24, 0
br i1 %25, label %37, label %26
26: ; preds = %23
%27 = load i32, i32* %14
%28 = load i32, i32* %11
%29 = icmp sge i32 %27, %28
br i1 %29, label %37, label %30
30: ; preds = %26
%31 = load i32, i32* %15
%32 = icmp slt i32 %31, 0
br i1 %32, label %37, label %33
33: ; preds = %30
%34 = load i32, i32* %15
%35 = load i32, i32* %12
%36 = icmp sge i32 %34, %35
br i1 %36, label %37, label %38
37: ; preds = %33, %30, %26, %23, %19, %7
store float 0.000000e+00, float* %8
br label %60
38: ; preds = %33
%39 = load i32, i32* %10
%40 = zext i32 %39 to i64
%41 = load i32, i32* %12
%42 = zext i32 %41 to i64
%43 = load float*, float** %9
%44 = bitcast float* %43 to i8*
%45 = bitcast i8* %44 to float*
store float* %45, float** %16
%46 = load float*, float** %16
%47 = load i32, i32* %14
%48 = sext i32 %47 to i64
%49 = mul nuw i64 %40, %42
%50 = mul nsw i64 %48, %49
%51 = getelementptr inbounds float, float* %46, i64 %50
%52 = load i32, i32* %13
%53 = sext i32 %52 to i64
%54 = mul nsw i64 %53, %42
%55 = getelementptr inbounds float, float* %51, i64 %54
%56 = load i32, i32* %15
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds float, float* %55, i64 %57
%59 = load float, float* %58
store float %59, float* %8
br label %60
60: ; preds = %38, %37
%61 = load float, float* %8
ret float %61
}

LLVM IR (-O2):
define float @getsample(float* nocapture %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6) {
%8 = icmp slt i32 %4, 0
br i1 %8, label %32, label %9
9: ; preds = %7
%10 = icmp sge i32 %4, %1
%11 = icmp slt i32 %5, 0
%12 = or i1 %10, %11
br i1 %12, label %32, label %13
13: ; preds = %9
%14 = icmp slt i32 %5, %2
%15 = icmp sgt i32 %6, -1
%16 = and i1 %14, %15
%17 = icmp slt i32 %6, %3
%18 = and i1 %17, %16
br i1 %18, label %19, label %32
19: ; preds = %13
%20 = zext i32 %1 to i64
%21 = zext i32 %3 to i64
%22 = zext i32 %5 to i64
%23 = mul nuw i64 %21, %20
%24 = mul nsw i64 %23, %22
%25 = getelementptr inbounds float, float* %0, i64 %24
%26 = zext i32 %4 to i64
%27 = mul nuw nsw i64 %26, %21
%28 = getelementptr inbounds float, float* %25, i64 %27
%29 = zext i32 %6 to i64
%30 = getelementptr inbounds float, float* %28, i64 %29
%31 = load float, float* %30
br label %32
32: ; preds = %7, %9, %13, %19
%33 = phi float [ %31, %19 ], [ 0.000000e+00, %13 ], [ 0.000000e+00, %9 ], [ 0.000000e+00, %7 ]
ret float %33
}