int get_grid_point(const int grid_double[3], const int mesh[3]) {
  int i, grid[3];
  for (i = 0; i < 3; i++) {
    if (grid_double[i] % 2 == 0) {
      grid[i] = grid_double[i] / 2;
    } else {
      grid[i] = (grid_double[i] - 1) / 2;
    }
  }
  return grid[2] * mesh[0] * mesh[1] + grid[1] * mesh[0] + grid[0];
}


x86 assembly (clang, -O0):
<get_grid_point>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
cmpl   $0x3,-0x14(%rbp)
jge    118e <get_grid_point+0x8e>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    115b <get_grid_point+0x5b>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
movslq -0x14(%rbp),%rax
mov    %ecx,-0x20(%rbp,%rax,4)
jmp    117b <get_grid_point+0x7b>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
sub    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
movslq -0x14(%rbp),%rax
mov    %ecx,-0x20(%rbp,%rax,4)
jmp    1180 <get_grid_point+0x80>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1113 <get_grid_point+0x13>
mov    -0x18(%rbp),%eax
mov    -0x10(%rbp),%rcx
imul   (%rcx),%eax
mov    -0x10(%rbp),%rcx
imul   0x4(%rcx),%eax
mov    -0x1c(%rbp),%ecx
mov    -0x10(%rbp),%rdx
imul   (%rdx),%ecx
add    %ecx,%eax
add    -0x20(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_grid_point>:
mov    (%rdi),%ecx
movq   0x4(%rdi),%xmm0
psrad  $0x1,%xmm0
sar    %ecx
pshufd $0x55,%xmm0,%xmm1
movd   %xmm1,%edx
imul   0x4(%rsi),%edx
movd   %xmm0,%eax
add    %edx,%eax
imul   (%rsi),%eax
add    %ecx,%eax
ret

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

LLVM IR (-O2):
define i32 @get_grid_point(i32* nocapture %0, i32* nocapture %1) {
%3 = load i32, i32* %0
%4 = and i32 %3, -2
%5 = getelementptr inbounds i32, i32* %0, i64 1
%6 = bitcast i32* %5 to <2 x i32>*
%7 = load <2 x i32>, <2 x i32>* %6
%8 = and <2 x i32> %7, <i32 -2, i32 -2>
%9 = sdiv <2 x i32> %8, <i32 2, i32 2>
%10 = sdiv i32 %4, 2
%11 = load i32, i32* %1
%12 = getelementptr inbounds i32, i32* %1, i64 1
%13 = load i32, i32* %12
%14 = extractelement <2 x i32> %9, i64 1
%15 = mul i32 %13, %14
%16 = extractelement <2 x i32> %9, i64 0
%17 = add i32 %16, %15
%18 = mul i32 %17, %11
%19 = add i32 %18, %10
ret i32 %19
}