void unpermute_indices(int *i, const int *pi, const int ifa) {
  const int axis_permut[6][4] = {{0, 2, 1, 0}, {1, 2, 0, 1}, {2, 0, 1, 1},
                                 {2, 1, 0, 0}, {0, 1, 2, 1}, {1, 0, 2, 0}};
  i[axis_permut[ifa][0]] = pi[0];
  i[axis_permut[ifa][1]] = pi[1];
  i[axis_permut[ifa][2]] = pi[2];
}


x86 assembly (clang, -O0):
<unpermute_indices>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
lea    -0x80(%rbp),%rdi
lea    0xecf(%rip),%rsi
mov    $0x60,%edx
call   1030 <memcpy@plt>
mov    -0x10(%rbp),%rax
mov    (%rax),%edx
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rsi
lea    -0x80(%rbp),%rcx
shl    $0x4,%rsi
add    %rsi,%rcx
movslq (%rcx),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    0x4(%rax),%edx
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rsi
lea    -0x80(%rbp),%rcx
shl    $0x4,%rsi
add    %rsi,%rcx
movslq 0x4(%rcx),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%edx
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rsi
lea    -0x80(%rbp),%rcx
shl    $0x4,%rsi
add    %rsi,%rcx
movslq 0x8(%rcx),%rcx
mov    %edx,(%rax,%rcx,4)
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<unpermute_indices>:
mov    (%rsi),%eax
movslq %edx,%rcx
shl    $0x4,%rcx
lea    0xef0(%rip),%r8
movslq (%rcx,%r8,1),%rdx
mov    %eax,(%rdi,%rdx,4)
mov    0x4(%rsi),%eax
movslq 0x4(%rcx,%r8,1),%rdx
mov    %eax,(%rdi,%rdx,4)
mov    0x8(%rsi),%eax
movslq 0x8(%rcx,%r8,1),%rcx
mov    %eax,(%rdi,%rcx,4)
ret

LLVM IR (-O0):
@__const.unpermute_indices.axis_permut = [6 x [4 x i32]] [[4 x i32] [i32 0, i32 2, i32 1, i32 0], [4 x i32] [i32 1, i32 2, i32 0, i32 1], [4 x i32] [i32 2, i32 0, i32 1, i32 1], [4 x i32] [i32 2, i32 1, i32 0, i32 0], [4 x i32] [i32 0, i32 1, i32 2, i32 1], [4 x i32] [i32 1, i32 0, i32 2, i32 0]]
define void @unpermute_indices(i32* %0, i32* %1, i32 %2) {
%4 = alloca i32*
%5 = alloca i32*
%6 = alloca i32
%7 = alloca [6 x [4 x i32]]
store i32* %0, i32** %4
store i32* %1, i32** %5
store i32 %2, i32* %6
%8 = bitcast [6 x [4 x i32]]* %7 to i8*
%9 = load i32*, i32** %5
%10 = getelementptr inbounds i32, i32* %9, i64 0
%11 = load i32, i32* %10
%12 = load i32*, i32** %4
%13 = load i32, i32* %6
%14 = sext i32 %13 to i64
%15 = getelementptr inbounds [6 x [4 x i32]], [6 x [4 x i32]]* %7, i64 0, i64 %14
%16 = getelementptr inbounds [4 x i32], [4 x i32]* %15, i64 0, i64 0
%17 = load i32, i32* %16
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds i32, i32* %12, i64 %18
store i32 %11, i32* %19
%20 = load i32*, i32** %5
%21 = getelementptr inbounds i32, i32* %20, i64 1
%22 = load i32, i32* %21
%23 = load i32*, i32** %4
%24 = load i32, i32* %6
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds [6 x [4 x i32]], [6 x [4 x i32]]* %7, i64 0, i64 %25
%27 = getelementptr inbounds [4 x i32], [4 x i32]* %26, i64 0, i64 1
%28 = load i32, i32* %27
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds i32, i32* %23, i64 %29
store i32 %22, i32* %30
%31 = load i32*, i32** %5
%32 = getelementptr inbounds i32, i32* %31, i64 2
%33 = load i32, i32* %32
%34 = load i32*, i32** %4
%35 = load i32, i32* %6
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds [6 x [4 x i32]], [6 x [4 x i32]]* %7, i64 0, i64 %36
%38 = getelementptr inbounds [4 x i32], [4 x i32]* %37, i64 0, i64 2
%39 = load i32, i32* %38
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds i32, i32* %34, i64 %40
store i32 %33, i32* %41
ret void
}

LLVM IR (-O2):
@__const.unpermute_indices.axis_permut = [6 x [4 x i32]] [[4 x i32] [i32 0, i32 2, i32 1, i32 0], [4 x i32] [i32 1, i32 2, i32 0, i32 1], [4 x i32] [i32 2, i32 0, i32 1, i32 1], [4 x i32] [i32 2, i32 1, i32 0, i32 0], [4 x i32] [i32 0, i32 1, i32 2, i32 1], [4 x i32] [i32 1, i32 0, i32 2, i32 0]]
define void @unpermute_indices(i32* nocapture writeonly %0, i32* nocapture %1, i32 %2) {
%4 = load i32, i32* %1
%5 = sext i32 %2 to i64
%6 = getelementptr inbounds [6 x [4 x i32]], [6 x [4 x i32]]* @__const.unpermute_indices.axis_permut, i64 0, i64 %5, i64 0
%7 = load i32, i32* %6
%8 = sext i32 %7 to i64
%9 = getelementptr inbounds i32, i32* %0, i64 %8
store i32 %4, i32* %9
%10 = getelementptr inbounds i32, i32* %1, i64 1
%11 = load i32, i32* %10
%12 = getelementptr inbounds [6 x [4 x i32]], [6 x [4 x i32]]* @__const.unpermute_indices.axis_permut, i64 0, i64 %5, i64 1
%13 = load i32, i32* %12
%14 = sext i32 %13 to i64
%15 = getelementptr inbounds i32, i32* %0, i64 %14
store i32 %11, i32* %15
%16 = getelementptr inbounds i32, i32* %1, i64 2
%17 = load i32, i32* %16
%18 = getelementptr inbounds [6 x [4 x i32]], [6 x [4 x i32]]* @__const.unpermute_indices.axis_permut, i64 0, i64 %5, i64 2
%19 = load i32, i32* %18
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds i32, i32* %0, i64 %20
store i32 %17, i32* %21
ret void
}