void transpose_matrix02(float **m, long nrow, long ncol, float **transposed) {
  long i, j;
  float m0, m1;
  float *pm;
  for (i = 0; i < nrow; i++) {
    pm = m[i];
    for (j = 0; j < ncol; j++) {
      m0 = *pm++;
      m1 = *pm++;
      transposed[j++][i] = m0;
      transposed[j][i] = m1;
    }
  }
}


x86 assembly (clang, -O0):
<transpose_matrix02>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    11e6 <transpose_matrix02+0xe6>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x40(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    11d0 <transpose_matrix02+0xd0>
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x40(%rbp)
movss  (%rax),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x40(%rbp)
movss  (%rax),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x34(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    %rcx,%rdx
add    $0x1,%rdx
mov    %rdx,-0x30(%rbp)
mov    (%rax,%rcx,8),%rax
mov    -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1142 <transpose_matrix02+0x42>
jmp    11d5 <transpose_matrix02+0xd5>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    111c <transpose_matrix02+0x1c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<transpose_matrix02>:
push   %r15
push   %r14
push   %rbx
test   %rsi,%rsi
jle    11e6 <transpose_matrix02+0xe6>
lea    -0x1(%rdx),%r8
mov    %r8,%r10
shr    %r10
add    $0x1,%r10
mov    %r10,%r9
and    $0xfffffffffffffffe,%r9
xor    %r11d,%r11d
jmp    113d <transpose_matrix02+0x3d>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r11
cmp    %rsi,%r11
je     11e6 <transpose_matrix02+0xe6>
test   %rdx,%rdx
jle    1130 <transpose_matrix02+0x30>
mov    (%rdi,%r11,8),%r14
cmp    $0x2,%r8
jae    1150 <transpose_matrix02+0x50>
xor    %eax,%eax
jmp    11b0 <transpose_matrix02+0xb0>
mov    %r9,%r15
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%r14),%xmm0
movss  0x4(%r14),%xmm1
mov    (%rcx,%rax,8),%rbx
movss  %xmm0,(%rbx,%r11,4)
mov    0x8(%rcx,%rax,8),%rbx
movss  %xmm1,(%rbx,%r11,4)
movss  0x8(%r14),%xmm0
movss  0xc(%r14),%xmm1
add    $0x10,%r14
mov    0x10(%rcx,%rax,8),%rbx
movss  %xmm0,(%rbx,%r11,4)
mov    0x18(%rcx,%rax,8),%rbx
movss  %xmm1,(%rbx,%r11,4)
add    $0x4,%rax
add    $0xfffffffffffffffe,%r15
jne    1160 <transpose_matrix02+0x60>
test   $0x1,%r10b
je     1130 <transpose_matrix02+0x30>
movss  (%r14),%xmm0
movss  0x4(%r14),%xmm1
mov    (%rcx,%rax,8),%rbx
shl    $0x3,%rax
movss  %xmm0,(%rbx,%r11,4)
or     $0x8,%rax
mov    (%rcx,%rax,1),%rax
movss  %xmm1,(%rax,%r11,4)
jmp    1130 <transpose_matrix02+0x30>
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @transpose_matrix02(float** %0, i64 %1, i64 %2, float** %3) {
%5 = alloca float**
%6 = alloca i64
%7 = alloca i64
%8 = alloca float**
%9 = alloca i64
%10 = alloca i64
%11 = alloca float
%12 = alloca float
%13 = alloca float*
store float** %0, float*** %5
store i64 %1, i64* %6
store i64 %2, i64* %7
store float** %3, float*** %8
store i64 0, i64* %9
br label %14
14: ; preds = %53, %4
%15 = load i64, i64* %9
%16 = load i64, i64* %6
%17 = icmp slt i64 %15, %16
br i1 %17, label %18, label %56
18: ; preds = %14
%19 = load float**, float*** %5
%20 = load i64, i64* %9
%21 = getelementptr inbounds float*, float** %19, i64 %20
%22 = load float*, float** %21
store float* %22, float** %13
store i64 0, i64* %10
br label %23
23: ; preds = %49, %18
%24 = load i64, i64* %10
%25 = load i64, i64* %7
%26 = icmp slt i64 %24, %25
br i1 %26, label %27, label %52
27: ; preds = %23
%28 = load float*, float** %13
%29 = getelementptr inbounds float, float* %28, i32 1
store float* %29, float** %13
%30 = load float, float* %28
store float %30, float* %11
%31 = load float*, float** %13
%32 = getelementptr inbounds float, float* %31, i32 1
store float* %32, float** %13
%33 = load float, float* %31
store float %33, float* %12
%34 = load float, float* %11
%35 = load float**, float*** %8
%36 = load i64, i64* %10
%37 = add nsw i64 %36, 1
store i64 %37, i64* %10
%38 = getelementptr inbounds float*, float** %35, i64 %36
%39 = load float*, float** %38
%40 = load i64, i64* %9
%41 = getelementptr inbounds float, float* %39, i64 %40
store float %34, float* %41
%42 = load float, float* %12
%43 = load float**, float*** %8
%44 = load i64, i64* %10
%45 = getelementptr inbounds float*, float** %43, i64 %44
%46 = load float*, float** %45
%47 = load i64, i64* %9
%48 = getelementptr inbounds float, float* %46, i64 %47
store float %42, float* %48
br label %49
49: ; preds = %27
%50 = load i64, i64* %10
%51 = add nsw i64 %50, 1
store i64 %51, i64* %10
br label %23
52: ; preds = %23
br label %53
53: ; preds = %52
%54 = load i64, i64* %9
%55 = add nsw i64 %54, 1
store i64 %55, i64* %9
br label %14
56: ; preds = %14
ret void
}

LLVM IR (-O2):
define void @transpose_matrix02(float** nocapture %0, i64 %1, i64 %2, float** nocapture %3) {
%5 = icmp sgt i64 %1, 0
br i1 %5, label %6, label %67
6: ; preds = %4
%7 = icmp sgt i64 %2, 0
%8 = add i64 %2, -1
%9 = lshr i64 %8, 1
%10 = add nuw i64 %9, 1
%11 = and i64 %10, 1
%12 = icmp ult i64 %8, 2
%13 = and i64 %10, -2
%14 = icmp eq i64 %11, 0
br label %15
15: ; preds = %6, %64
%16 = phi i64 [ 0, %6 ], [ %65, %64 ]
br i1 %7, label %17, label %64
17: ; preds = %15
%18 = getelementptr inbounds float*, float** %0, i64 %16
%19 = load float*, float** %18
br i1 %12, label %50, label %20
20: ; preds = %17, %20
%21 = phi float* [ %38, %20 ], [ %19, %17 ]
%22 = phi i64 [ %47, %20 ], [ 0, %17 ]
%23 = phi i64 [ %48, %20 ], [ 0, %17 ]
%24 = getelementptr inbounds float, float* %21, i64 1
%25 = load float, float* %21
%26 = getelementptr inbounds float, float* %21, i64 2
%27 = load float, float* %24
%28 = or i64 %22, 1
%29 = getelementptr inbounds float*, float** %3, i64 %22
%30 = load float*, float** %29
%31 = getelementptr inbounds float, float* %30, i64 %16
store float %25, float* %31
%32 = getelementptr inbounds float*, float** %3, i64 %28
%33 = load float*, float** %32
%34 = getelementptr inbounds float, float* %33, i64 %16
store float %27, float* %34
%35 = or i64 %22, 2
%36 = getelementptr inbounds float, float* %21, i64 3
%37 = load float, float* %26
%38 = getelementptr inbounds float, float* %21, i64 4
%39 = load float, float* %36
%40 = or i64 %22, 3
%41 = getelementptr inbounds float*, float** %3, i64 %35
%42 = load float*, float** %41
%43 = getelementptr inbounds float, float* %42, i64 %16
store float %37, float* %43
%44 = getelementptr inbounds float*, float** %3, i64 %40
%45 = load float*, float** %44
%46 = getelementptr inbounds float, float* %45, i64 %16
store float %39, float* %46
%47 = add nuw nsw i64 %22, 4
%48 = add nuw nsw i64 %23, 2
%49 = icmp eq i64 %48, %13
br i1 %49, label %50, label %20
50: ; preds = %20, %17
%51 = phi float* [ %19, %17 ], [ %38, %20 ]
%52 = phi i64 [ 0, %17 ], [ %47, %20 ]
br i1 %14, label %64, label %53
53: ; preds = %50
%54 = getelementptr inbounds float, float* %51, i64 1
%55 = load float, float* %51
%56 = load float, float* %54
%57 = or i64 %52, 1
%58 = getelementptr inbounds float*, float** %3, i64 %52
%59 = load float*, float** %58
%60 = getelementptr inbounds float, float* %59, i64 %16
store float %55, float* %60
%61 = getelementptr inbounds float*, float** %3, i64 %57
%62 = load float*, float** %61
%63 = getelementptr inbounds float, float* %62, i64 %16
store float %56, float* %63
br label %64
64: ; preds = %53, %50, %15
%65 = add nuw nsw i64 %16, 1
%66 = icmp eq i64 %65, %1
br i1 %66, label %67, label %15
67: ; preds = %64, %4
ret void
}