void grid_point_to_address_double(int address_double[3], const int grid_point,
                                  const int mesh[3], const int is_shift[3]) {
  int i;
  int address[3];
  address[2] = grid_point / (mesh[0] * mesh[1]);
  address[1] = (grid_point - address[2] * mesh[0] * mesh[1]) / mesh[0];
  address[0] = grid_point % mesh[0];
  for (i = 0; i < 3; i++) {
    address_double[i] = address[i] * 2 + is_shift[i];
  }
}


x86 assembly (clang, -O0):
<grid_point_to_address_double>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0xc(%rbp),%eax
mov    -0x18(%rbp),%rcx
mov    (%rcx),%ecx
mov    -0x18(%rbp),%rdx
imul   0x4(%rdx),%ecx
cltd
idiv   %ecx
mov    %eax,-0x28(%rbp)
mov    -0xc(%rbp),%eax
mov    -0x28(%rbp),%ecx
mov    -0x18(%rbp),%rdx
imul   (%rdx),%ecx
mov    -0x18(%rbp),%rdx
imul   0x4(%rdx),%ecx
sub    %ecx,%eax
mov    -0x18(%rbp),%rcx
cltd
idivl  (%rcx)
mov    %eax,-0x2c(%rbp)
mov    -0xc(%rbp),%eax
mov    -0x18(%rbp),%rcx
cltd
idivl  (%rcx)
mov    %edx,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
cmpl   $0x3,-0x24(%rbp)
jge    1198 <grid_point_to_address_double+0x98>
movslq -0x24(%rbp),%rax
mov    -0x30(%rbp,%rax,4),%edx
shl    $0x1,%edx
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
add    (%rax,%rcx,4),%edx
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    115f <grid_point_to_address_double+0x5f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<grid_point_to_address_double>:
mov    %esi,%r8d
mov    (%rdx),%r10d
mov    0x4(%rdx),%esi
imul   %r10d,%esi
mov    %r8d,%eax
cltd
idiv   %esi
mov    %eax,%r9d
mov    %edx,%eax
cltd
idiv   %r10d
mov    %eax,%esi
mov    %r8d,%eax
cltd
idiv   %r10d
add    %edx,%edx
add    (%rcx),%edx
mov    %edx,(%rdi)
add    %esi,%esi
add    0x4(%rcx),%esi
mov    %esi,0x4(%rdi)
add    %r9d,%r9d
add    0x8(%rcx),%r9d
mov    %r9d,0x8(%rdi)
ret

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

LLVM IR (-O2):
define void @grid_point_to_address_double(i32* nocapture writeonly %0, i32 %1, i32* nocapture %2, i32* nocapture %3) {
%5 = load i32, i32* %2
%6 = getelementptr inbounds i32, i32* %2, i64 1
%7 = load i32, i32* %6
%8 = mul i32 %7, %5
%9 = sdiv i32 %1, %8
%10 = mul i32 %8, %9
%11 = srem i32 %1, %8
%12 = sdiv i32 %11, %5
%13 = srem i32 %1, %5
%14 = shl nsw i32 %13, 1
%15 = load i32, i32* %3
%16 = add nsw i32 %14, %15
store i32 %16, i32* %0
%17 = shl nsw i32 %12, 1
%18 = getelementptr inbounds i32, i32* %3, i64 1
%19 = load i32, i32* %18
%20 = add nsw i32 %17, %19
%21 = getelementptr inbounds i32, i32* %0, i64 1
store i32 %20, i32* %21
%22 = shl nsw i32 %9, 1
%23 = getelementptr inbounds i32, i32* %3, i64 2
%24 = load i32, i32* %23
%25 = add nsw i32 %22, %24
%26 = getelementptr inbounds i32, i32* %0, i64 2
store i32 %25, i32* %26
ret void
}