void copytranspose(float (*a)[4], float (*b)[4]) {
  int i, j;
  for (i = 0; i < 4; i++)
    for (j = 0; j < 4; j++)
      a[i][j] = b[j][i];
}


x86 assembly (clang, -O0):
<copytranspose>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
cmpl   $0x4,-0x14(%rbp)
jge    117f <copytranspose+0x7f>
movl   $0x0,-0x18(%rbp)
cmpl   $0x4,-0x18(%rbp)
jge    116c <copytranspose+0x6c>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x18(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1124 <copytranspose+0x24>
jmp    1171 <copytranspose+0x71>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1113 <copytranspose+0x13>
pop    %rbp
ret

x86 assembly (clang, -O2):
<copytranspose>:
movss  (%rsi),%xmm0
movss  %xmm0,(%rdi)
movss  0x10(%rsi),%xmm0
movss  %xmm0,0x4(%rdi)
movss  0x20(%rsi),%xmm0
movss  %xmm0,0x8(%rdi)
movss  0x30(%rsi),%xmm0
movss  %xmm0,0xc(%rdi)
movss  0x4(%rsi),%xmm0
movss  %xmm0,0x10(%rdi)
movss  0x14(%rsi),%xmm0
movss  %xmm0,0x14(%rdi)
movss  0x24(%rsi),%xmm0
movss  %xmm0,0x18(%rdi)
movss  0x34(%rsi),%xmm0
movss  %xmm0,0x1c(%rdi)
movss  0x8(%rsi),%xmm0
movss  %xmm0,0x20(%rdi)
movss  0x18(%rsi),%xmm0
movss  %xmm0,0x24(%rdi)
movss  0x28(%rsi),%xmm0
movss  %xmm0,0x28(%rdi)
movss  0x38(%rsi),%xmm0
movss  %xmm0,0x2c(%rdi)
movss  0xc(%rsi),%xmm0
movss  %xmm0,0x30(%rdi)
movss  0x1c(%rsi),%xmm0
movss  %xmm0,0x34(%rdi)
movss  0x2c(%rsi),%xmm0
movss  %xmm0,0x38(%rdi)
movss  0x3c(%rsi),%xmm0
movss  %xmm0,0x3c(%rdi)
ret

LLVM IR (-O0):
define void @copytranspose([4 x float]* %0, [4 x float]* %1) {
%3 = alloca [4 x float]*
%4 = alloca [4 x float]*
%5 = alloca i32
%6 = alloca i32
store [4 x float]* %0, [4 x float]** %3
store [4 x float]* %1, [4 x float]** %4
store i32 0, i32* %5
br label %7
7: ; preds = %34, %2
%8 = load i32, i32* %5
%9 = icmp slt i32 %8, 4
br i1 %9, label %10, label %37
10: ; preds = %7
store i32 0, i32* %6
br label %11
11: ; preds = %30, %10
%12 = load i32, i32* %6
%13 = icmp slt i32 %12, 4
br i1 %13, label %14, label %33
14: ; preds = %11
%15 = load [4 x float]*, [4 x float]** %4
%16 = load i32, i32* %6
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds [4 x float], [4 x float]* %15, i64 %17
%19 = load i32, i32* %5
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds [4 x float], [4 x float]* %18, i64 0, i64 %20
%22 = load float, float* %21
%23 = load [4 x float]*, [4 x float]** %3
%24 = load i32, i32* %5
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds [4 x float], [4 x float]* %23, i64 %25
%27 = load i32, i32* %6
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds [4 x float], [4 x float]* %26, i64 0, i64 %28
store float %22, float* %29
br label %30
30: ; preds = %14
%31 = load i32, i32* %6
%32 = add nsw i32 %31, 1
store i32 %32, i32* %6
br label %11
33: ; preds = %11
br label %34
34: ; preds = %33
%35 = load i32, i32* %5
%36 = add nsw i32 %35, 1
store i32 %36, i32* %5
br label %7
37: ; preds = %7
ret void
}

LLVM IR (-O2):
define void @copytranspose([4 x float]* nocapture writeonly %0, [4 x float]* nocapture %1) {
%3 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 0, i64 0
%4 = load float, float* %3
%5 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 0, i64 0
store float %4, float* %5
%6 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 1, i64 0
%7 = load float, float* %6
%8 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 0, i64 1
store float %7, float* %8
%9 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 2, i64 0
%10 = load float, float* %9
%11 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 0, i64 2
store float %10, float* %11
%12 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 3, i64 0
%13 = load float, float* %12
%14 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 0, i64 3
store float %13, float* %14
%15 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 0, i64 1
%16 = load float, float* %15
%17 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 1, i64 0
store float %16, float* %17
%18 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 1, i64 1
%19 = load float, float* %18
%20 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 1, i64 1
store float %19, float* %20
%21 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 2, i64 1
%22 = load float, float* %21
%23 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 1, i64 2
store float %22, float* %23
%24 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 3, i64 1
%25 = load float, float* %24
%26 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 1, i64 3
store float %25, float* %26
%27 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 0, i64 2
%28 = load float, float* %27
%29 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 2, i64 0
store float %28, float* %29
%30 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 1, i64 2
%31 = load float, float* %30
%32 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 2, i64 1
store float %31, float* %32
%33 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 2, i64 2
%34 = load float, float* %33
%35 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 2, i64 2
store float %34, float* %35
%36 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 3, i64 2
%37 = load float, float* %36
%38 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 2, i64 3
store float %37, float* %38
%39 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 0, i64 3
%40 = load float, float* %39
%41 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 3, i64 0
store float %40, float* %41
%42 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 1, i64 3
%43 = load float, float* %42
%44 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 3, i64 1
store float %43, float* %44
%45 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 2, i64 3
%46 = load float, float* %45
%47 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 3, i64 2
store float %46, float* %47
%48 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 3, i64 3
%49 = load float, float* %48
%50 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 3, i64 3
store float %49, float* %50
ret void
}