void extract(float (*b)[4], float (*a)[4][4], int k) {
  b[0][0] = a[0][0][k];
  b[0][1] = a[1][0][k];
  b[0][2] = a[2][0][k];
  b[0][3] = a[3][0][k];
  b[1][0] = a[0][1][k];
  b[1][1] = a[1][1][k];
  b[1][2] = a[2][1][k];
  b[1][3] = a[3][1][k];
  b[2][0] = a[0][2][k];
  b[2][1] = a[1][2][k];
  b[2][2] = a[2][2][k];
  b[2][3] = a[3][2][k];
  b[3][0] = a[0][3][k];
  b[3][1] = a[1][3][k];
  b[3][2] = a[2][3][k];
  b[3][3] = a[3][3][k];
}


x86 assembly (clang, -O0):
<extract>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0x40(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x4(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0x80(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x8(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0xc0(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0xc(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0x10(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x10(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0x50(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x14(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0x90(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x18(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0xd0(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x1c(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0x20(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x20(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0x60(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x24(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0xa0(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x28(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0xe0(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x2c(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0x30(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x30(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0x70(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x34(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0xb0(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x38(%rax)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  0xf0(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x3c(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<extract>:
movslq %edx,%rax
movss  (%rsi,%rax,4),%xmm0
movss  %xmm0,(%rdi)
movss  0x40(%rsi,%rax,4),%xmm0
movss  %xmm0,0x4(%rdi)
movss  0x80(%rsi,%rax,4),%xmm0
movss  %xmm0,0x8(%rdi)
movss  0xc0(%rsi,%rax,4),%xmm0
movss  %xmm0,0xc(%rdi)
movss  0x10(%rsi,%rax,4),%xmm0
movss  %xmm0,0x10(%rdi)
movss  0x50(%rsi,%rax,4),%xmm0
movss  %xmm0,0x14(%rdi)
movss  0x90(%rsi,%rax,4),%xmm0
movss  %xmm0,0x18(%rdi)
movss  0xd0(%rsi,%rax,4),%xmm0
movss  %xmm0,0x1c(%rdi)
movss  0x20(%rsi,%rax,4),%xmm0
movss  %xmm0,0x20(%rdi)
movss  0x60(%rsi,%rax,4),%xmm0
movss  %xmm0,0x24(%rdi)
movss  0xa0(%rsi,%rax,4),%xmm0
movss  %xmm0,0x28(%rdi)
movss  0xe0(%rsi,%rax,4),%xmm0
movss  %xmm0,0x2c(%rdi)
movss  0x30(%rsi,%rax,4),%xmm0
movss  %xmm0,0x30(%rdi)
movss  0x70(%rsi,%rax,4),%xmm0
movss  %xmm0,0x34(%rdi)
movss  0xb0(%rsi,%rax,4),%xmm0
movss  %xmm0,0x38(%rdi)
movss  0xf0(%rsi,%rax,4),%xmm0
movss  %xmm0,0x3c(%rdi)
ret

LLVM IR (-O0):
define void @extract([4 x float]* %0, [4 x [4 x float]]* %1, i32 %2) {
%4 = alloca [4 x float]*
%5 = alloca [4 x [4 x float]]*
%6 = alloca i32
store [4 x float]* %0, [4 x float]** %4
store [4 x [4 x float]]* %1, [4 x [4 x float]]** %5
store i32 %2, i32* %6
%7 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%8 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %7, i64 0
%9 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %8, i64 0, i64 0
%10 = load i32, i32* %6
%11 = sext i32 %10 to i64
%12 = getelementptr inbounds [4 x float], [4 x float]* %9, i64 0, i64 %11
%13 = load float, float* %12
%14 = load [4 x float]*, [4 x float]** %4
%15 = getelementptr inbounds [4 x float], [4 x float]* %14, i64 0
%16 = getelementptr inbounds [4 x float], [4 x float]* %15, i64 0, i64 0
store float %13, float* %16
%17 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%18 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %17, i64 1
%19 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %18, i64 0, i64 0
%20 = load i32, i32* %6
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds [4 x float], [4 x float]* %19, i64 0, i64 %21
%23 = load float, float* %22
%24 = load [4 x float]*, [4 x float]** %4
%25 = getelementptr inbounds [4 x float], [4 x float]* %24, i64 0
%26 = getelementptr inbounds [4 x float], [4 x float]* %25, i64 0, i64 1
store float %23, float* %26
%27 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%28 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %27, i64 2
%29 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %28, i64 0, i64 0
%30 = load i32, i32* %6
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds [4 x float], [4 x float]* %29, i64 0, i64 %31
%33 = load float, float* %32
%34 = load [4 x float]*, [4 x float]** %4
%35 = getelementptr inbounds [4 x float], [4 x float]* %34, i64 0
%36 = getelementptr inbounds [4 x float], [4 x float]* %35, i64 0, i64 2
store float %33, float* %36
%37 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%38 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %37, i64 3
%39 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %38, i64 0, i64 0
%40 = load i32, i32* %6
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds [4 x float], [4 x float]* %39, i64 0, i64 %41
%43 = load float, float* %42
%44 = load [4 x float]*, [4 x float]** %4
%45 = getelementptr inbounds [4 x float], [4 x float]* %44, i64 0
%46 = getelementptr inbounds [4 x float], [4 x float]* %45, i64 0, i64 3
store float %43, float* %46
%47 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%48 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %47, i64 0
%49 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %48, i64 0, i64 1
%50 = load i32, i32* %6
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds [4 x float], [4 x float]* %49, i64 0, i64 %51
%53 = load float, float* %52
%54 = load [4 x float]*, [4 x float]** %4
%55 = getelementptr inbounds [4 x float], [4 x float]* %54, i64 1
%56 = getelementptr inbounds [4 x float], [4 x float]* %55, i64 0, i64 0
store float %53, float* %56
%57 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%58 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %57, i64 1
%59 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %58, i64 0, i64 1
%60 = load i32, i32* %6
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds [4 x float], [4 x float]* %59, i64 0, i64 %61
%63 = load float, float* %62
%64 = load [4 x float]*, [4 x float]** %4
%65 = getelementptr inbounds [4 x float], [4 x float]* %64, i64 1
%66 = getelementptr inbounds [4 x float], [4 x float]* %65, i64 0, i64 1
store float %63, float* %66
%67 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%68 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %67, i64 2
%69 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %68, i64 0, i64 1
%70 = load i32, i32* %6
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds [4 x float], [4 x float]* %69, i64 0, i64 %71
%73 = load float, float* %72
%74 = load [4 x float]*, [4 x float]** %4
%75 = getelementptr inbounds [4 x float], [4 x float]* %74, i64 1
%76 = getelementptr inbounds [4 x float], [4 x float]* %75, i64 0, i64 2
store float %73, float* %76
%77 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%78 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %77, i64 3
%79 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %78, i64 0, i64 1
%80 = load i32, i32* %6
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds [4 x float], [4 x float]* %79, i64 0, i64 %81
%83 = load float, float* %82
%84 = load [4 x float]*, [4 x float]** %4
%85 = getelementptr inbounds [4 x float], [4 x float]* %84, i64 1
%86 = getelementptr inbounds [4 x float], [4 x float]* %85, i64 0, i64 3
store float %83, float* %86
%87 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%88 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %87, i64 0
%89 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %88, i64 0, i64 2
%90 = load i32, i32* %6
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds [4 x float], [4 x float]* %89, i64 0, i64 %91
%93 = load float, float* %92
%94 = load [4 x float]*, [4 x float]** %4
%95 = getelementptr inbounds [4 x float], [4 x float]* %94, i64 2
%96 = getelementptr inbounds [4 x float], [4 x float]* %95, i64 0, i64 0
store float %93, float* %96
%97 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%98 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %97, i64 1
%99 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %98, i64 0, i64 2
%100 = load i32, i32* %6
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds [4 x float], [4 x float]* %99, i64 0, i64 %101
%103 = load float, float* %102
%104 = load [4 x float]*, [4 x float]** %4
%105 = getelementptr inbounds [4 x float], [4 x float]* %104, i64 2
%106 = getelementptr inbounds [4 x float], [4 x float]* %105, i64 0, i64 1
store float %103, float* %106
%107 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%108 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %107, i64 2
%109 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %108, i64 0, i64 2
%110 = load i32, i32* %6
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds [4 x float], [4 x float]* %109, i64 0, i64 %111
%113 = load float, float* %112
%114 = load [4 x float]*, [4 x float]** %4
%115 = getelementptr inbounds [4 x float], [4 x float]* %114, i64 2
%116 = getelementptr inbounds [4 x float], [4 x float]* %115, i64 0, i64 2
store float %113, float* %116
%117 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%118 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %117, i64 3
%119 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %118, i64 0, i64 2
%120 = load i32, i32* %6
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds [4 x float], [4 x float]* %119, i64 0, i64 %121
%123 = load float, float* %122
%124 = load [4 x float]*, [4 x float]** %4
%125 = getelementptr inbounds [4 x float], [4 x float]* %124, i64 2
%126 = getelementptr inbounds [4 x float], [4 x float]* %125, i64 0, i64 3
store float %123, float* %126
%127 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%128 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %127, i64 0
%129 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %128, i64 0, i64 3
%130 = load i32, i32* %6
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds [4 x float], [4 x float]* %129, i64 0, i64 %131
%133 = load float, float* %132
%134 = load [4 x float]*, [4 x float]** %4
%135 = getelementptr inbounds [4 x float], [4 x float]* %134, i64 3
%136 = getelementptr inbounds [4 x float], [4 x float]* %135, i64 0, i64 0
store float %133, float* %136
%137 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%138 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %137, i64 1
%139 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %138, i64 0, i64 3
%140 = load i32, i32* %6
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds [4 x float], [4 x float]* %139, i64 0, i64 %141
%143 = load float, float* %142
%144 = load [4 x float]*, [4 x float]** %4
%145 = getelementptr inbounds [4 x float], [4 x float]* %144, i64 3
%146 = getelementptr inbounds [4 x float], [4 x float]* %145, i64 0, i64 1
store float %143, float* %146
%147 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%148 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %147, i64 2
%149 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %148, i64 0, i64 3
%150 = load i32, i32* %6
%151 = sext i32 %150 to i64
%152 = getelementptr inbounds [4 x float], [4 x float]* %149, i64 0, i64 %151
%153 = load float, float* %152
%154 = load [4 x float]*, [4 x float]** %4
%155 = getelementptr inbounds [4 x float], [4 x float]* %154, i64 3
%156 = getelementptr inbounds [4 x float], [4 x float]* %155, i64 0, i64 2
store float %153, float* %156
%157 = load [4 x [4 x float]]*, [4 x [4 x float]]** %5
%158 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %157, i64 3
%159 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %158, i64 0, i64 3
%160 = load i32, i32* %6
%161 = sext i32 %160 to i64
%162 = getelementptr inbounds [4 x float], [4 x float]* %159, i64 0, i64 %161
%163 = load float, float* %162
%164 = load [4 x float]*, [4 x float]** %4
%165 = getelementptr inbounds [4 x float], [4 x float]* %164, i64 3
%166 = getelementptr inbounds [4 x float], [4 x float]* %165, i64 0, i64 3
store float %163, float* %166
ret void
}

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