void opencl_nd_address(int linear_id, const unsigned int *dims,
                       unsigned int *id) {
  if (dims[2] == 1) {
    id[2] = 0;
    id[1] = linear_id / dims[0];
    id[0] = linear_id - (id[1] * dims[0]);
  } else {
    unsigned int items_per_slice = (dims[1] * dims[0]);
    id[2] = linear_id / items_per_slice;
    unsigned int tmp = linear_id - (id[2] * items_per_slice);
    id[1] = tmp / dims[0];
    id[0] = tmp - (id[1] * dims[0]);
  }
  return;
}


x86 assembly (clang, -O0):
<opencl_nd_address>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
cmpl   $0x1,0x8(%rax)
jne    115a <opencl_nd_address+0x5a>
mov    -0x18(%rbp),%rax
movl   $0x0,0x8(%rax)
mov    -0x4(%rbp),%eax
mov    -0x10(%rbp),%rcx
xor    %edx,%edx
divl   (%rcx)
mov    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,0x4(%rax)
mov    -0x4(%rbp),%ecx
mov    -0x18(%rbp),%rax
mov    0x4(%rax),%eax
mov    -0x10(%rbp),%rdx
imul   (%rdx),%eax
sub    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
jmp    11bc <opencl_nd_address+0xbc>
mov    -0x10(%rbp),%rax
mov    0x4(%rax),%eax
mov    -0x10(%rbp),%rcx
imul   (%rcx),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x4(%rbp),%eax
xor    %edx,%edx
divl   -0x1c(%rbp)
mov    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,0x8(%rax)
mov    -0x4(%rbp),%eax
mov    -0x18(%rbp),%rcx
mov    0x8(%rcx),%ecx
imul   -0x1c(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x10(%rbp),%rcx
xor    %edx,%edx
divl   (%rcx)
mov    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,0x4(%rax)
mov    -0x20(%rbp),%ecx
mov    -0x18(%rbp),%rax
mov    0x4(%rax),%eax
mov    -0x10(%rbp),%rdx
imul   (%rdx),%eax
sub    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<opencl_nd_address>:
mov    %rdx,%rcx
cmpl   $0x1,0x8(%rsi)
jne    1112 <opencl_nd_address+0x12>
movl   $0x0,0x8(%rcx)
jmp    1126 <opencl_nd_address+0x26>
mov    (%rsi),%r8d
imul   0x4(%rsi),%r8d
mov    %edi,%eax
xor    %edx,%edx
div    %r8d
mov    %eax,0x8(%rcx)
mov    %edx,%edi
mov    %edi,%eax
xor    %edx,%edx
divl   (%rsi)
mov    %eax,0x4(%rcx)
imul   (%rsi),%eax
sub    %eax,%edi
mov    %edi,(%rcx)
ret

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

LLVM IR (-O2):
define void @opencl_nd_address(i32 %0, i32* nocapture %1, i32* nocapture writeonly %2) {
%4 = getelementptr inbounds i32, i32* %1, i64 2
%5 = load i32, i32* %4
%6 = icmp eq i32 %5, 1
br i1 %6, label %7, label %9
7: ; preds = %3
%8 = getelementptr inbounds i32, i32* %2, i64 2
store i32 0, i32* %8
br label %18
9: ; preds = %3
%10 = getelementptr inbounds i32, i32* %1, i64 1
%11 = load i32, i32* %10
%12 = load i32, i32* %1
%13 = mul i32 %12, %11
%14 = udiv i32 %0, %13
%15 = getelementptr inbounds i32, i32* %2, i64 2
store i32 %14, i32* %15
%16 = mul i32 %14, %13
%17 = urem i32 %0, %13
br label %18
18: ; preds = %9, %7
%19 = phi i32 [ %17, %9 ], [ %0, %7 ]
%20 = load i32, i32* %1
%21 = udiv i32 %19, %20
%22 = getelementptr inbounds i32, i32* %2, i64 1
store i32 %21, i32* %22
%23 = load i32, i32* %1
%24 = mul i32 %23, %21
%25 = sub i32 %19, %24
store i32 %25, i32* %2
ret void
}