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


x86 assembly (clang, -O0):
<transpose_matrix01>:
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    11a9 <transpose_matrix01+0xa9>
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    1193 <transpose_matrix01+0x93>
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x40(%rbp)
movss  (%rax),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x34(%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_matrix01+0x42>
jmp    1198 <transpose_matrix01+0x98>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    111c <transpose_matrix01+0x1c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<transpose_matrix01>:
push   %r15
push   %r14
push   %rbx
test   %rsi,%rsi
jle    11ce <transpose_matrix01+0xce>
lea    -0x1(%rdx),%r8
mov    %edx,%r10d
and    $0x3,%r10d
mov    %rdx,%r9
and    $0xfffffffffffffffc,%r9
xor    %r11d,%r11d
jmp    113d <transpose_matrix01+0x3d>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r11
cmp    %rsi,%r11
je     11ce <transpose_matrix01+0xce>
test   %rdx,%rdx
jle    1130 <transpose_matrix01+0x30>
mov    (%rdi,%r11,8),%r15
xor    %ebx,%ebx
cmp    $0x3,%r8
jb     119f <transpose_matrix01+0x9f>
xchg   %ax,%ax
movss  (%r15),%xmm0
mov    (%rcx,%rbx,8),%rax
movss  %xmm0,(%rax,%r11,4)
movss  0x4(%r15),%xmm0
mov    0x8(%rcx,%rbx,8),%rax
movss  %xmm0,(%rax,%r11,4)
movss  0x8(%r15),%xmm0
mov    0x10(%rcx,%rbx,8),%rax
movss  %xmm0,(%rax,%r11,4)
movss  0xc(%r15),%xmm0
add    $0x10,%r15
mov    0x18(%rcx,%rbx,8),%rax
movss  %xmm0,(%rax,%r11,4)
add    $0x4,%rbx
cmp    %rbx,%r9
jne    1150 <transpose_matrix01+0x50>
test   %r10,%r10
je     1130 <transpose_matrix01+0x30>
lea    (%rcx,%rbx,8),%r14
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movss  (%r15,%rbx,4),%xmm0
mov    (%r14,%rbx,8),%rax
movss  %xmm0,(%rax,%r11,4)
add    $0x1,%rbx
cmp    %rbx,%r10
jne    11b0 <transpose_matrix01+0xb0>
jmp    1130 <transpose_matrix01+0x30>
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @transpose_matrix01(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*
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 %13
13: ; preds = %41, %4
%14 = load i64, i64* %9
%15 = load i64, i64* %6
%16 = icmp slt i64 %14, %15
br i1 %16, label %17, label %44
17: ; preds = %13
%18 = load float**, float*** %5
%19 = load i64, i64* %9
%20 = getelementptr inbounds float*, float** %18, i64 %19
%21 = load float*, float** %20
store float* %21, float** %12
store i64 0, i64* %10
br label %22
22: ; preds = %37, %17
%23 = load i64, i64* %10
%24 = load i64, i64* %7
%25 = icmp slt i64 %23, %24
br i1 %25, label %26, label %40
26: ; preds = %22
%27 = load float*, float** %12
%28 = getelementptr inbounds float, float* %27, i32 1
store float* %28, float** %12
%29 = load float, float* %27
store float %29, float* %11
%30 = load float, float* %11
%31 = load float**, float*** %8
%32 = load i64, i64* %10
%33 = getelementptr inbounds float*, float** %31, i64 %32
%34 = load float*, float** %33
%35 = load i64, i64* %9
%36 = getelementptr inbounds float, float* %34, i64 %35
store float %30, float* %36
br label %37
37: ; preds = %26
%38 = load i64, i64* %10
%39 = add nsw i64 %38, 1
store i64 %39, i64* %10
br label %22
40: ; preds = %22
br label %41
41: ; preds = %40
%42 = load i64, i64* %9
%43 = add nsw i64 %42, 1
store i64 %43, i64* %9
br label %13
44: ; preds = %13
ret void
}

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