unsigned int part_nd_address_rec(unsigned int stride, int cur, int dim,
                                 unsigned int addr, const unsigned int *size,
                                 unsigned int *pos) {
  if (cur == dim)
    return addr;
  else {
    unsigned int rem =
        part_nd_address_rec(size[cur] * stride, cur + 1, dim, addr, size, pos);
    pos[cur] = rem / stride;
    return rem % stride;
  }
}


x86 assembly (clang, -O0):
<part_nd_address_rec>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x10(%rbp),%eax
jne    1143 <part_nd_address_rec+0x33>
mov    -0x14(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    118e <part_nd_address_rec+0x7e>
mov    -0x20(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,4),%edi
imul   -0x8(%rbp),%edi
mov    -0xc(%rbp),%esi
add    $0x1,%esi
mov    -0x10(%rbp),%edx
mov    -0x14(%rbp),%ecx
mov    -0x20(%rbp),%r8
mov    -0x28(%rbp),%r9
call   1030 <part_nd_address_rec@plt>
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
xor    %edx,%edx
divl   -0x8(%rbp)
mov    %eax,%edx
mov    -0x28(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
xor    %edx,%edx
divl   -0x8(%rbp)
mov    %edx,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<part_nd_address_rec>:
cmp    %edx,%esi
je     113d <part_nd_address_rec+0x2d>
push   %rbp
push   %r14
push   %rbx
mov    %r9,%rbx
mov    %edi,%ebp
movslq %esi,%r14
mov    (%r8,%r14,4),%edi
imul   %ebp,%edi
add    $0x1,%esi
call   1030 <part_nd_address_rec@plt>
xor    %edx,%edx
div    %ebp
mov    %edx,%ecx
mov    %eax,(%rbx,%r14,4)
pop    %rbx
pop    %r14
pop    %rbp
mov    %ecx,%eax
ret

LLVM IR (-O0):
define i32 @part_nd_address_rec(i32 %0, i32 %1, i32 %2, i32 %3, i32* %4, i32* %5) {
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32*
%13 = alloca i32*
%14 = alloca i32
store i32 %0, i32* %8
store i32 %1, i32* %9
store i32 %2, i32* %10
store i32 %3, i32* %11
store i32* %4, i32** %12
store i32* %5, i32** %13
%15 = load i32, i32* %9
%16 = load i32, i32* %10
%17 = icmp eq i32 %15, %16
br i1 %17, label %18, label %20
18: ; preds = %6
%19 = load i32, i32* %11
store i32 %19, i32* %7
br label %45
20: ; preds = %6
%21 = load i32*, i32** %12
%22 = load i32, i32* %9
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds i32, i32* %21, i64 %23
%25 = load i32, i32* %24
%26 = load i32, i32* %8
%27 = mul i32 %25, %26
%28 = load i32, i32* %9
%29 = add nsw i32 %28, 1
%30 = load i32, i32* %10
%31 = load i32, i32* %11
%32 = load i32*, i32** %12
%33 = load i32*, i32** %13
%34 = call i32 @part_nd_address_rec(i32 %27, i32 %29, i32 %30, i32 %31, i32* %32, i32* %33)
store i32 %34, i32* %14
%35 = load i32, i32* %14
%36 = load i32, i32* %8
%37 = udiv i32 %35, %36
%38 = load i32*, i32** %13
%39 = load i32, i32* %9
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds i32, i32* %38, i64 %40
store i32 %37, i32* %41
%42 = load i32, i32* %14
%43 = load i32, i32* %8
%44 = urem i32 %42, %43
store i32 %44, i32* %7
br label %45
45: ; preds = %20, %18
%46 = load i32, i32* %7
ret i32 %46
}

LLVM IR (-O2):
define i32 @part_nd_address_rec(i32 %0, i32 %1, i32 %2, i32 %3, i32* nocapture %4, i32* nocapture writeonly %5) {
%7 = icmp eq i32 %1, %2
br i1 %7, label %8, label %10
8: ; preds = %6, %10
%9 = phi i32 [ %19, %10 ], [ %3, %6 ]
ret i32 %9
10: ; preds = %6
%11 = sext i32 %1 to i64
%12 = getelementptr inbounds i32, i32* %4, i64 %11
%13 = load i32, i32* %12
%14 = mul i32 %13, %0
%15 = add nsw i32 %1, 1
%16 = tail call i32 @part_nd_address_rec(i32 %14, i32 %15, i32 %2, i32 %3, i32* %4, i32* %5)
%17 = udiv i32 %16, %0
%18 = getelementptr inbounds i32, i32* %5, i64 %11
store i32 %17, i32* %18
%19 = urem i32 %16, %0
br label %8
}