#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 */
void transpose(double *table, size_t nRow, size_t nCol) {
  size_t i;
  for (i = 1; i < nRow * nCol - 1; i++) {
    size_t x = nRow * (i % nCol) + i / nCol;
    if (x <= i) {
      continue;
    }
    while (x > i) {
      x = nRow * (x % nCol) + x / nCol;
    }
    if (x < i) {
      continue;
    }
    {
      double tmp = table[i];
      size_t s = i;
      x = nRow * (i % nCol) + i / nCol;
      while (x != i) {
        table[s] = table[x];
        s = x;
        x = nRow * (x % nCol) + x / nCol;
      }
      table[s] = tmp;
    }
  }
}


x86 assembly (clang, -O0):
<transpose>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movq   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
imul   -0x18(%rbp),%rcx
sub    $0x1,%rcx
cmp    %rcx,%rax
jae    12ca <transpose+0x1ca>
mov    -0x10(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x20(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x18(%rbp)
mov    -0x48(%rbp),%rax
imul   %rdx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x20(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x18(%rbp)
mov    %rax,%rcx
mov    -0x40(%rbp),%rax
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x20(%rbp),%rax
ja     117f <transpose+0x7f>
jmp    12b9 <transpose+0x1b9>
jmp    1184 <transpose+0x84>
mov    -0x28(%rbp),%rax
cmp    -0x20(%rbp),%rax
jbe    11d1 <transpose+0xd1>
mov    -0x10(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x28(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x18(%rbp)
mov    -0x58(%rbp),%rax
imul   %rdx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x28(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x18(%rbp)
mov    %rax,%rcx
mov    -0x50(%rbp),%rax
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
jmp    1184 <transpose+0x84>
mov    -0x28(%rbp),%rax
cmp    -0x20(%rbp),%rax
jae    11e4 <transpose+0xe4>
jmp    12b9 <transpose+0x1b9>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x68(%rbp)
mov    -0x20(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x18(%rbp)
mov    -0x68(%rbp),%rax
imul   %rdx,%rax
mov    %rax,-0x60(%rbp)
mov    -0x20(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x18(%rbp)
mov    %rax,%rcx
mov    -0x60(%rbp),%rax
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x20(%rbp),%rax
je     12a7 <transpose+0x1a7>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x78(%rbp)
mov    -0x28(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x18(%rbp)
mov    -0x78(%rbp),%rax
imul   %rdx,%rax
mov    %rax,-0x70(%rbp)
mov    -0x28(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x18(%rbp)
mov    %rax,%rcx
mov    -0x70(%rbp),%rax
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
jmp    1238 <transpose+0x138>
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    1118 <transpose+0x18>
pop    %rbp
ret

x86 assembly (clang, -O2):
<transpose>:
mov    %rdx,%r11
imul   %rsi,%r11
add    $0xffffffffffffffff,%r11
cmp    $0x2,%r11
jb     11f3 <transpose+0xf3>
mov    %rdx,%r9
mov    $0x1,%r8d
jmp    1133 <transpose+0x33>
movsd  %xmm0,(%rdi,%r10,8)
add    $0x1,%r8
cmp    %r11,%r8
je     11f3 <transpose+0xf3>
mov    %r8,%rax
or     %r9,%rax
shr    $0x20,%rax
je     1150 <transpose+0x50>
mov    %r8,%rax
xor    %edx,%edx
div    %r9
mov    %rdx,%rcx
jmp    115a <transpose+0x5a>
nopl   0x0(%rax)
mov    %r8d,%eax
xor    %edx,%edx
div    %r9d
mov    %edx,%ecx
imul   %rsi,%rcx
add    %rax,%rcx
mov    %rcx,%rax
cmp    %r8,%rcx
ja     1187 <transpose+0x87>
jmp    1126 <transpose+0x26>
nopl   0x0(%rax,%rax,1)
xor    %edx,%edx
div    %r9
mov    %rax,%r10
mov    %rdx,%rax
imul   %rsi,%rax
add    %r10,%rax
cmp    %r8,%rax
jbe    11a0 <transpose+0xa0>
mov    %rax,%rdx
or     %r9,%rdx
shr    $0x20,%rdx
jne    1170 <transpose+0x70>
xor    %edx,%edx
div    %r9d
mov    %eax,%r10d
mov    %edx,%eax
jmp    117b <transpose+0x7b>
nop
jb     1126 <transpose+0x26>
movsd  (%rdi,%r8,8),%xmm0
mov    %r8,%rax
jmp    11ce <transpose+0xce>
nopl   (%rax)
mov    %r10,%rax
xor    %edx,%edx
div    %r9
mov    %rdx,%rcx
imul   %rsi,%rcx
add    %rax,%rcx
mov    %r10,%rax
cmp    %r8,%rcx
je     1120 <transpose+0x20>
mov    %rcx,%r10
movsd  (%rdi,%rcx,8),%xmm1
movsd  %xmm1,(%rdi,%rax,8)
mov    %rcx,%rax
or     %r9,%rax
shr    $0x20,%rax
jne    11b0 <transpose+0xb0>
mov    %r10d,%eax
xor    %edx,%edx
div    %r9d
mov    %edx,%ecx
jmp    11bb <transpose+0xbb>
ret

LLVM IR (-O0):
define void @transpose(double* %0, i64 %1, i64 %2) {
%4 = alloca double*
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca i64
%9 = alloca double
%10 = alloca i64
store double* %0, double** %4
store i64 %1, i64* %5
store i64 %2, i64* %6
store i64 1, i64* %7
br label %11
11: ; preds = %94, %3
%12 = load i64, i64* %7
%13 = load i64, i64* %5
%14 = load i64, i64* %6
%15 = mul i64 %13, %14
%16 = sub i64 %15, 1
%17 = icmp ult i64 %12, %16
br i1 %17, label %18, label %97
18: ; preds = %11
%19 = load i64, i64* %5
%20 = load i64, i64* %7
%21 = load i64, i64* %6
%22 = urem i64 %20, %21
%23 = mul i64 %19, %22
%24 = load i64, i64* %7
%25 = load i64, i64* %6
%26 = udiv i64 %24, %25
%27 = add i64 %23, %26
store i64 %27, i64* %8
%28 = load i64, i64* %8
%29 = load i64, i64* %7
%30 = icmp ule i64 %28, %29
br i1 %30, label %31, label %32
31: ; preds = %18
br label %94
32: ; preds = %18
br label %33
33: ; preds = %37, %32
%34 = load i64, i64* %8
%35 = load i64, i64* %7
%36 = icmp ugt i64 %34, %35
br i1 %36, label %37, label %47
37: ; preds = %33
%38 = load i64, i64* %5
%39 = load i64, i64* %8
%40 = load i64, i64* %6
%41 = urem i64 %39, %40
%42 = mul i64 %38, %41
%43 = load i64, i64* %8
%44 = load i64, i64* %6
%45 = udiv i64 %43, %44
%46 = add i64 %42, %45
store i64 %46, i64* %8
br label %33
47: ; preds = %33
%48 = load i64, i64* %8
%49 = load i64, i64* %7
%50 = icmp ult i64 %48, %49
br i1 %50, label %51, label %52
51: ; preds = %47
br label %94
52: ; preds = %47
%53 = load double*, double** %4
%54 = load i64, i64* %7
%55 = getelementptr inbounds double, double* %53, i64 %54
%56 = load double, double* %55
store double %56, double* %9
%57 = load i64, i64* %7
store i64 %57, i64* %10
%58 = load i64, i64* %5
%59 = load i64, i64* %7
%60 = load i64, i64* %6
%61 = urem i64 %59, %60
%62 = mul i64 %58, %61
%63 = load i64, i64* %7
%64 = load i64, i64* %6
%65 = udiv i64 %63, %64
%66 = add i64 %62, %65
store i64 %66, i64* %8
br label %67
67: ; preds = %71, %52
%68 = load i64, i64* %8
%69 = load i64, i64* %7
%70 = icmp ne i64 %68, %69
br i1 %70, label %71, label %89
71: ; preds = %67
%72 = load double*, double** %4
%73 = load i64, i64* %8
%74 = getelementptr inbounds double, double* %72, i64 %73
%75 = load double, double* %74
%76 = load double*, double** %4
%77 = load i64, i64* %10
%78 = getelementptr inbounds double, double* %76, i64 %77
store double %75, double* %78
%79 = load i64, i64* %8
store i64 %79, i64* %10
%80 = load i64, i64* %5
%81 = load i64, i64* %8
%82 = load i64, i64* %6
%83 = urem i64 %81, %82
%84 = mul i64 %80, %83
%85 = load i64, i64* %8
%86 = load i64, i64* %6
%87 = udiv i64 %85, %86
%88 = add i64 %84, %87
store i64 %88, i64* %8
br label %67
89: ; preds = %67
%90 = load double, double* %9
%91 = load double*, double** %4
%92 = load i64, i64* %10
%93 = getelementptr inbounds double, double* %91, i64 %92
store double %90, double* %93
br label %94
94: ; preds = %89, %51, %31
%95 = load i64, i64* %7
%96 = add i64 %95, 1
store i64 %96, i64* %7
br label %11
97: ; preds = %11
ret void
}

LLVM IR (-O2):
define void @transpose(double* nocapture %0, i64 %1, i64 %2) {
%4 = mul i64 %2, %1
%5 = add i64 %4, -1
%6 = icmp ugt i64 %5, 1
br i1 %6, label %7, label %42
7: ; preds = %3, %39
%8 = phi i64 [ %40, %39 ], [ 1, %3 ]
%9 = urem i64 %8, %2
%10 = mul i64 %9, %1
%11 = udiv i64 %8, %2
%12 = add i64 %10, %11
%13 = icmp ugt i64 %12, %8
br i1 %13, label %14, label %39
14: ; preds = %7, %14
%15 = phi i64 [ %19, %14 ], [ %12, %7 ]
%16 = urem i64 %15, %2
%17 = mul i64 %16, %1
%18 = udiv i64 %15, %2
%19 = add i64 %17, %18
%20 = icmp ugt i64 %19, %8
br i1 %20, label %14, label %21
21: ; preds = %14
%22 = icmp ult i64 %19, %8
br i1 %22, label %39, label %23
23: ; preds = %21
%24 = getelementptr inbounds double, double* %0, i64 %8
%25 = load double, double* %24
br label %26
26: ; preds = %23, %26
%27 = phi i64 [ %28, %26 ], [ %8, %23 ]
%28 = phi i64 [ %35, %26 ], [ %12, %23 ]
%29 = getelementptr inbounds double, double* %0, i64 %28
%30 = load double, double* %29
%31 = getelementptr inbounds double, double* %0, i64 %27
store double %30, double* %31
%32 = urem i64 %28, %2
%33 = mul i64 %32, %1
%34 = udiv i64 %28, %2
%35 = add i64 %33, %34
%36 = icmp eq i64 %35, %8
br i1 %36, label %37, label %26
37: ; preds = %26
%38 = getelementptr inbounds double, double* %0, i64 %28
store double %25, double* %38
br label %39
39: ; preds = %21, %7, %37
%40 = add nuw i64 %8, 1
%41 = icmp eq i64 %40, %5
br i1 %41, label %42, label %7
42: ; preds = %39, %3
ret void
}