void ref_transpose_column(long double **in, long height, long double *out,
                          long in_offset, long cmp) {
  long i, j;
  for (i = 0; i < height; i++) {
    for (j = 0; j < cmp; j++) {
      out[i * cmp + j] = in[i][in_offset + j];
    }
  }
}


x86 assembly (clang, -O0):
<ref_transpose_column>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    119f <ref_transpose_column+0x9f>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    1189 <ref_transpose_column+0x89>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rsi
mov    -0x38(%rbp),%rdx
add    %rdx,%rsi
shl    $0x4,%rsi
fldt   (%rax,%rsi,1)
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%rsi
imul   %rsi,%rcx
add    %rdx,%rcx
shl    $0x4,%rcx
fstpt  (%rax,%rcx,1)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1136 <ref_transpose_column+0x36>
jmp    118e <ref_transpose_column+0x8e>
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1120 <ref_transpose_column+0x20>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ref_transpose_column>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %rsi,%rsi
jle    11a6 <ref_transpose_column+0xa6>
mov    %r8,%r15
and    $0xfffffffffffffffe,%r15
mov    %r8,%r9
shl    $0x4,%r9
xor    %r10d,%r10d
mov    %rcx,%r11
shl    $0x4,%r11
mov    %rdx,%r12
jmp    113c <ref_transpose_column+0x3c>
nop
add    $0x1,%r10
add    %r9,%r12
cmp    %rsi,%r10
je     11a6 <ref_transpose_column+0xa6>
test   %r8,%r8
jle    1130 <ref_transpose_column+0x30>
mov    (%rdi,%r10,8),%r14
cmp    $0x1,%r8
jne    1150 <ref_transpose_column+0x50>
xor    %eax,%eax
jmp    1180 <ref_transpose_column+0x80>
nop
lea    (%r14,%r11,1),%r13
mov    $0x10,%ebx
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
fldt   -0x10(%r13,%rbx,1)
fstpt  -0x10(%r12,%rbx,1)
fldt   0x0(%r13,%rbx,1)
fstpt  (%r12,%rbx,1)
add    $0x2,%rax
add    $0x20,%rbx
cmp    %rax,%r15
jne    1160 <ref_transpose_column+0x60>
test   $0x1,%r8b
je     1130 <ref_transpose_column+0x30>
mov    %r10,%r13
imul   %r8,%r13
lea    (%rax,%rcx,1),%rbx
shl    $0x4,%rbx
fldt   (%r14,%rbx,1)
add    %rax,%r13
shl    $0x4,%r13
fstpt  (%rdx,%r13,1)
jmp    1130 <ref_transpose_column+0x30>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

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

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