#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ array_swap_d(double *, double *, unsigned int const,
                                unsigned int const);
int /*<<< orphan*/ array_swap_ui(unsigned int *, unsigned int *, int, int);
int /*<<< orphan*/ free(unsigned int *);
unsigned int *malloc(unsigned int const);
void array_rearrange_d(const unsigned int NRows, const unsigned int NCols,
                       const unsigned int *Ordering, const char layout,
                       double *A) {
  unsigned int i, RowInd, RowSub, ReOrder, *RowTracker;
  RowTracker = malloc(NRows * sizeof *RowTracker);
  for (i = 0; i < NRows; i++)
    RowTracker[i] = i;
  for (RowInd = 0; RowInd < NRows; RowInd++) {
    ReOrder = Ordering[RowInd];
    for (RowSub = ReOrder; RowTracker[RowSub] != ReOrder;
         RowSub = RowTracker[RowSub])
      ;
    if (RowInd != RowSub) {
      if (layout == 'C')
        array_swap_d(&A[RowInd], &A[RowSub], NCols, NRows);
      else
        array_swap_d(&A[RowInd * NCols], &A[RowSub * NCols], NCols, 1);
      array_swap_ui(&RowTracker[RowInd], &RowTracker[RowSub], 1, 1);
    }
  }
  free(RowTracker);
}


x86 assembly (clang, -O0):
<array_rearrange_d>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %cl,%al
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %al,-0x11(%rbp)
mov    %r8,-0x20(%rbp)
mov    -0x4(%rbp),%eax
shl    $0x2,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jae    119b <array_rearrange_d+0x5b>
mov    -0x24(%rbp),%edx
mov    -0x38(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    %edx,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1174 <array_rearrange_d+0x34>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jae    12a0 <array_rearrange_d+0x160>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    (%rax,%rcx,4),%eax
cmp    -0x30(%rbp),%eax
je     11eb <array_rearrange_d+0xab>
jmp    11d9 <array_rearrange_d+0x99>
mov    -0x38(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
jmp    11c1 <array_rearrange_d+0x81>
mov    -0x28(%rbp),%eax
cmp    -0x2c(%rbp),%eax
je     128d <array_rearrange_d+0x14d>
movsbl -0x11(%rbp),%eax
cmp    $0x43,%eax
jne    1230 <array_rearrange_d+0xf0>
mov    -0x20(%rbp),%rdi
mov    -0x28(%rbp),%eax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x20(%rbp),%rsi
mov    -0x2c(%rbp),%eax
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x8(%rbp),%edx
mov    -0x4(%rbp),%ecx
call   1040 <array_swap_d@plt>
jmp    1265 <array_rearrange_d+0x125>
mov    -0x20(%rbp),%rdi
mov    -0x28(%rbp),%eax
imul   -0x8(%rbp),%eax
mov    %eax,%eax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x20(%rbp),%rsi
mov    -0x2c(%rbp),%eax
imul   -0x8(%rbp),%eax
mov    %eax,%eax
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x8(%rbp),%edx
mov    $0x1,%ecx
call   1040 <array_swap_d@plt>
mov    -0x38(%rbp),%rdi
mov    -0x28(%rbp),%eax
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x38(%rbp),%rsi
mov    -0x2c(%rbp),%eax
shl    $0x2,%rax
add    %rax,%rsi
mov    $0x1,%ecx
mov    %ecx,%edx
call   1050 <array_swap_ui@plt>
jmp    1292 <array_rearrange_d+0x152>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11a2 <array_rearrange_d+0x62>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<array_rearrange_d>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r8,0x10(%rsp)
mov    %ecx,0xc(%rsp)
mov    %rdx,%r13
mov    %esi,%r14d
mov    %edi,%ebp
lea    0x0(,%rbp,4),%edi
call   1060 <malloc@plt>
mov    %rax,%rbx
mov    %rbp,(%rsp)
test   %ebp,%ebp
je     137e <array_rearrange_d+0x23e>
mov    %r14d,%ebp
mov    (%rsp),%rax
mov    %eax,%r14d
cmp    $0x8,%eax
jae    1190 <array_rearrange_d+0x50>
xor    %eax,%eax
jmp    12e0 <array_rearrange_d+0x1a0>
mov    %r14d,%eax
and    $0xfffffff8,%eax
lea    -0x8(%rax),%rsi
mov    %rsi,%rdx
shr    $0x3,%rdx
add    $0x1,%rdx
mov    %edx,%ecx
and    $0x3,%ecx
cmp    $0x18,%rsi
jae    11bf <array_rearrange_d+0x7f>
movdqa 0xe48(%rip),%xmm0
xor    %esi,%esi
jmp    1288 <array_rearrange_d+0x148>
and    $0xfffffffffffffffc,%rdx
movdqa 0xe35(%rip),%xmm0
xor    %esi,%esi
movdqa 0xe3a(%rip),%xmm8
movdqa 0xe41(%rip),%xmm9
movdqa 0xe48(%rip),%xmm10
movdqa 0xe50(%rip),%xmm4
movdqa 0xe58(%rip),%xmm5
movdqa 0xe60(%rip),%xmm6
movdqa 0xe68(%rip),%xmm7
movdqa 0xe70(%rip),%xmm1
movdqa %xmm0,%xmm2
paddd  %xmm8,%xmm2
movdqu %xmm0,(%rbx,%rsi,4)
movdqu %xmm2,0x10(%rbx,%rsi,4)
movdqa %xmm0,%xmm2
paddd  %xmm9,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm10,%xmm3
movdqu %xmm2,0x20(%rbx,%rsi,4)
movdqu %xmm3,0x30(%rbx,%rsi,4)
movdqa %xmm0,%xmm2
paddd  %xmm4,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm5,%xmm3
movdqu %xmm2,0x40(%rbx,%rsi,4)
movdqu %xmm3,0x50(%rbx,%rsi,4)
movdqa %xmm0,%xmm2
paddd  %xmm6,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm7,%xmm3
movdqu %xmm2,0x60(%rbx,%rsi,4)
movdqu %xmm3,0x70(%rbx,%rsi,4)
add    $0x20,%rsi
paddd  %xmm1,%xmm0
add    $0xfffffffffffffffc,%rdx
jne    1210 <array_rearrange_d+0xd0>
test   %rcx,%rcx
je     12d0 <array_rearrange_d+0x190>
lea    (%rbx,%rsi,4),%rdx
add    $0x10,%rdx
shl    $0x5,%rcx
xor    %esi,%esi
movdqa 0xd6d(%rip),%xmm1
movdqa 0xd75(%rip),%xmm2
nopl   0x0(%rax,%rax,1)
movdqa %xmm0,%xmm3
paddd  %xmm1,%xmm3
movdqu %xmm0,-0x10(%rdx,%rsi,1)
movdqu %xmm3,(%rdx,%rsi,1)
paddd  %xmm2,%xmm0
add    $0x20,%rsi
cmp    %rsi,%rcx
jne    12b0 <array_rearrange_d+0x170>
cmp    %r14,%rax
je     12ec <array_rearrange_d+0x1ac>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %eax,(%rbx,%rax,4)
add    $0x1,%rax
cmp    %rax,%r14
jne    12e0 <array_rearrange_d+0x1a0>
cmpl   $0x0,(%rsp)
je     137e <array_rearrange_d+0x23e>
xor    %r12d,%r12d
jmp    1342 <array_rearrange_d+0x202>
nopl   0x0(%rax,%rax,1)
mov    %r12d,%ecx
imul   %ebp,%ecx
mov    0x10(%rsp),%rdx
lea    (%rdx,%rcx,8),%rdi
imul   %ebp,%eax
lea    (%rdx,%rax,8),%rsi
mov    %ebp,%edx
mov    $0x1,%ecx
call   1040 <array_swap_d@plt>
lea    (%rbx,%r15,4),%rsi
lea    (%rbx,%r12,4),%rdi
mov    $0x1,%edx
mov    $0x1,%ecx
call   1050 <array_swap_ui@plt>
add    $0x1,%r12
cmp    %r14,%r12
je     137e <array_rearrange_d+0x23e>
mov    0x0(%r13,%r12,4),%ecx
mov    %ecx,%edx
nopl   0x0(%rax)
mov    %edx,%eax
mov    %edx,%r15d
mov    (%rbx,%r15,4),%edx
cmp    %ecx,%edx
jne    1350 <array_rearrange_d+0x210>
cmp    %r15,%r12
je     1339 <array_rearrange_d+0x1f9>
cmpb   $0x43,0xc(%rsp)
jne    1300 <array_rearrange_d+0x1c0>
mov    0x10(%rsp),%rax
lea    (%rax,%r12,8),%rdi
lea    (%rax,%r15,8),%rsi
mov    %ebp,%edx
mov    (%rsp),%rcx
jmp    131d <array_rearrange_d+0x1dd>
mov    %rbx,%rdi
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @array_rearrange_d(i32 %0, i32 %1, i32* %2, i8 signext %3, double* %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32*
%9 = alloca i8
%10 = alloca double*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32*
store i32 %0, i32* %6
store i32 %1, i32* %7
store i32* %2, i32** %8
store i8 %3, i8* %9
store double* %4, double** %10
%16 = load i32, i32* %6
%17 = zext i32 %16 to i64
%18 = mul i64 %17, 4
%19 = trunc i64 %18 to i32
%20 = call i32* @malloc(i32 %19)
store i32* %20, i32** %15
store i32 0, i32* %11
br label %21
21: ; preds = %31, %5
%22 = load i32, i32* %11
%23 = load i32, i32* %6
%24 = icmp ult i32 %22, %23
br i1 %24, label %25, label %34
25: ; preds = %21
%26 = load i32, i32* %11
%27 = load i32*, i32** %15
%28 = load i32, i32* %11
%29 = zext i32 %28 to i64
%30 = getelementptr inbounds i32, i32* %27, i64 %29
store i32 %26, i32* %30
br label %31
31: ; preds = %25
%32 = load i32, i32* %11
%33 = add i32 %32, 1
store i32 %33, i32* %11
br label %21
34: ; preds = %21
store i32 0, i32* %12
br label %35
35: ; preds = %107, %34
%36 = load i32, i32* %12
%37 = load i32, i32* %6
%38 = icmp ult i32 %36, %37
br i1 %38, label %39, label %110
39: ; preds = %35
%40 = load i32*, i32** %8
%41 = load i32, i32* %12
%42 = zext i32 %41 to i64
%43 = getelementptr inbounds i32, i32* %40, i64 %42
%44 = load i32, i32* %43
store i32 %44, i32* %14
%45 = load i32, i32* %14
store i32 %45, i32* %13
br label %46
46: ; preds = %55, %39
%47 = load i32*, i32** %15
%48 = load i32, i32* %13
%49 = zext i32 %48 to i64
%50 = getelementptr inbounds i32, i32* %47, i64 %49
%51 = load i32, i32* %50
%52 = load i32, i32* %14
%53 = icmp ne i32 %51, %52
br i1 %53, label %54, label %61
54: ; preds = %46
br label %55
55: ; preds = %54
%56 = load i32*, i32** %15
%57 = load i32, i32* %13
%58 = zext i32 %57 to i64
%59 = getelementptr inbounds i32, i32* %56, i64 %58
%60 = load i32, i32* %59
store i32 %60, i32* %13
br label %46
61: ; preds = %46
%62 = load i32, i32* %12
%63 = load i32, i32* %13
%64 = icmp ne i32 %62, %63
br i1 %64, label %65, label %106
65: ; preds = %61
%66 = load i8, i8* %9
%67 = sext i8 %66 to i32
%68 = icmp eq i32 %67, 67
br i1 %68, label %69, label %81
69: ; preds = %65
%70 = load double*, double** %10
%71 = load i32, i32* %12
%72 = zext i32 %71 to i64
%73 = getelementptr inbounds double, double* %70, i64 %72
%74 = load double*, double** %10
%75 = load i32, i32* %13
%76 = zext i32 %75 to i64
%77 = getelementptr inbounds double, double* %74, i64 %76
%78 = load i32, i32* %7
%79 = load i32, i32* %6
%80 = call i32 @array_swap_d(double* %73, double* %77, i32 %78, i32 %79)
br label %96
81: ; preds = %65
%82 = load double*, double** %10
%83 = load i32, i32* %12
%84 = load i32, i32* %7
%85 = mul i32 %83, %84
%86 = zext i32 %85 to i64
%87 = getelementptr inbounds double, double* %82, i64 %86
%88 = load double*, double** %10
%89 = load i32, i32* %13
%90 = load i32, i32* %7
%91 = mul i32 %89, %90
%92 = zext i32 %91 to i64
%93 = getelementptr inbounds double, double* %88, i64 %92
%94 = load i32, i32* %7
%95 = call i32 @array_swap_d(double* %87, double* %93, i32 %94, i32 1)
br label %96
96: ; preds = %81, %69
%97 = load i32*, i32** %15
%98 = load i32, i32* %12
%99 = zext i32 %98 to i64
%100 = getelementptr inbounds i32, i32* %97, i64 %99
%101 = load i32*, i32** %15
%102 = load i32, i32* %13
%103 = zext i32 %102 to i64
%104 = getelementptr inbounds i32, i32* %101, i64 %103
%105 = call i32 @array_swap_ui(i32* %100, i32* %104, i32 1, i32 1)
br label %106
106: ; preds = %96, %61
br label %107
107: ; preds = %106
%108 = load i32, i32* %12
%109 = add i32 %108, 1
store i32 %109, i32* %12
br label %35
110: ; preds = %35
%111 = load i32*, i32** %15
%112 = call i32 @free(i32* %111)
ret void
}

LLVM IR (-O2):
define void @array_rearrange_d(i32 %0, i32 %1, i32* nocapture %2, i8 signext %3, double* %4) {
%6 = shl i32 %0, 2
%7 = tail call i32* @malloc(i32 %6)
%8 = icmp eq i32 %0, 0
br i1 %8, label %120, label %9
9: ; preds = %5
%10 = zext i32 %0 to i64
%11 = icmp ult i32 %0, 8
br i1 %11, label %74, label %12
12: ; preds = %9
%13 = and i64 %10, 4294967288
%14 = add nsw i64 %13, -8
%15 = lshr exact i64 %14, 3
%16 = add nuw nsw i64 %15, 1
%17 = and i64 %16, 3
%18 = icmp ult i64 %14, 24
br i1 %18, label %55, label %19
19: ; preds = %12
%20 = and i64 %16, 4611686018427387900
br label %21
21: ; preds = %21, %19
%22 = phi i64 [ 0, %19 ], [ %51, %21 ]
%23 = phi <4 x i32> [ <i32 0, i32 1, i32 2, i32 3>, %19 ], [ %52, %21 ]
%24 = phi i64 [ 0, %19 ], [ %53, %21 ]
%25 = add <4 x i32> %23, <i32 4, i32 4, i32 4, i32 4>
%26 = getelementptr inbounds i32, i32* %7, i64 %22
%27 = bitcast i32* %26 to <4 x i32>*
store <4 x i32> %23, <4 x i32>* %27
%28 = getelementptr inbounds i32, i32* %26, i64 4
%29 = bitcast i32* %28 to <4 x i32>*
store <4 x i32> %25, <4 x i32>* %29
%30 = or i64 %22, 8
%31 = add <4 x i32> %23, <i32 8, i32 8, i32 8, i32 8>
%32 = add <4 x i32> %23, <i32 12, i32 12, i32 12, i32 12>
%33 = getelementptr inbounds i32, i32* %7, i64 %30
%34 = bitcast i32* %33 to <4 x i32>*
store <4 x i32> %31, <4 x i32>* %34
%35 = getelementptr inbounds i32, i32* %33, i64 4
%36 = bitcast i32* %35 to <4 x i32>*
store <4 x i32> %32, <4 x i32>* %36
%37 = or i64 %22, 16
%38 = add <4 x i32> %23, <i32 16, i32 16, i32 16, i32 16>
%39 = add <4 x i32> %23, <i32 20, i32 20, i32 20, i32 20>
%40 = getelementptr inbounds i32, i32* %7, i64 %37
%41 = bitcast i32* %40 to <4 x i32>*
store <4 x i32> %38, <4 x i32>* %41
%42 = getelementptr inbounds i32, i32* %40, i64 4
%43 = bitcast i32* %42 to <4 x i32>*
store <4 x i32> %39, <4 x i32>* %43
%44 = or i64 %22, 24
%45 = add <4 x i32> %23, <i32 24, i32 24, i32 24, i32 24>
%46 = add <4 x i32> %23, <i32 28, i32 28, i32 28, i32 28>
%47 = getelementptr inbounds i32, i32* %7, i64 %44
%48 = bitcast i32* %47 to <4 x i32>*
store <4 x i32> %45, <4 x i32>* %48
%49 = getelementptr inbounds i32, i32* %47, i64 4
%50 = bitcast i32* %49 to <4 x i32>*
store <4 x i32> %46, <4 x i32>* %50
%51 = add nuw i64 %22, 32
%52 = add <4 x i32> %23, <i32 32, i32 32, i32 32, i32 32>
%53 = add i64 %24, 4
%54 = icmp eq i64 %53, %20
br i1 %54, label %55, label %21
55: ; preds = %21, %12
%56 = phi i64 [ 0, %12 ], [ %51, %21 ]
%57 = phi <4 x i32> [ <i32 0, i32 1, i32 2, i32 3>, %12 ], [ %52, %21 ]
%58 = icmp eq i64 %17, 0
br i1 %58, label %72, label %59
59: ; preds = %55, %59
%60 = phi i64 [ %68, %59 ], [ %56, %55 ]
%61 = phi <4 x i32> [ %69, %59 ], [ %57, %55 ]
%62 = phi i64 [ %70, %59 ], [ 0, %55 ]
%63 = add <4 x i32> %61, <i32 4, i32 4, i32 4, i32 4>
%64 = getelementptr inbounds i32, i32* %7, i64 %60
%65 = bitcast i32* %64 to <4 x i32>*
store <4 x i32> %61, <4 x i32>* %65
%66 = getelementptr inbounds i32, i32* %64, i64 4
%67 = bitcast i32* %66 to <4 x i32>*
store <4 x i32> %63, <4 x i32>* %67
%68 = add nuw i64 %60, 8
%69 = add <4 x i32> %61, <i32 8, i32 8, i32 8, i32 8>
%70 = add i64 %62, 1
%71 = icmp eq i64 %70, %17
br i1 %71, label %72, label %59
72: ; preds = %59, %55
%73 = icmp eq i64 %13, %10
br i1 %73, label %76, label %74
74: ; preds = %9, %72
%75 = phi i64 [ 0, %9 ], [ %13, %72 ]
br label %80
76: ; preds = %80, %72
br i1 %8, label %120, label %77
77: ; preds = %76
%78 = icmp eq i8 %3, 67
%79 = zext i32 %0 to i64
br label %86
80: ; preds = %74, %80
%81 = phi i64 [ %84, %80 ], [ %75, %74 ]
%82 = getelementptr inbounds i32, i32* %7, i64 %81
%83 = trunc i64 %81 to i32
store i32 %83, i32* %82
%84 = add nuw nsw i64 %81, 1
%85 = icmp eq i64 %84, %10
br i1 %85, label %76, label %80
86: ; preds = %77, %117
%87 = phi i64 [ 0, %77 ], [ %118, %117 ]
%88 = getelementptr inbounds i32, i32* %2, i64 %87
%89 = load i32, i32* %88
br label %90
90: ; preds = %90, %86
%91 = phi i32 [ %89, %86 ], [ %94, %90 ]
%92 = zext i32 %91 to i64
%93 = getelementptr inbounds i32, i32* %7, i64 %92
%94 = load i32, i32* %93
%95 = icmp eq i32 %94, %89
br i1 %95, label %96, label %90
96: ; preds = %90
%97 = zext i32 %91 to i64
%98 = getelementptr inbounds i32, i32* %7, i64 %97
%99 = icmp eq i64 %87, %97
br i1 %99, label %117, label %100
100: ; preds = %96
br i1 %78, label %101, label %105
101: ; preds = %100
%102 = getelementptr inbounds double, double* %4, i64 %87
%103 = getelementptr inbounds double, double* %4, i64 %97
%104 = tail call i32 @array_swap_d(double* %102, double* %103, i32 %1, i32 %0)
br label %114
105: ; preds = %100
%106 = trunc i64 %87 to i32
%107 = mul i32 %106, %1
%108 = zext i32 %107 to i64
%109 = getelementptr inbounds double, double* %4, i64 %108
%110 = mul i32 %91, %1
%111 = zext i32 %110 to i64
%112 = getelementptr inbounds double, double* %4, i64 %111
%113 = tail call i32 @array_swap_d(double* %109, double* %112, i32 %1, i32 1)
br label %114
114: ; preds = %105, %101
%115 = getelementptr inbounds i32, i32* %7, i64 %87
%116 = tail call i32 @array_swap_ui(i32* %115, i32* %98, i32 1, i32 1)
br label %117
117: ; preds = %96, %114
%118 = add nuw nsw i64 %87, 1
%119 = icmp eq i64 %118, %79
br i1 %119, label %120, label %86
120: ; preds = %117, %5, %76
%121 = tail call i32 @free(i32* %7)
ret void
}