void Matrix_Copy(float *pSrc, unsigned short numRows, unsigned short numCols,
                 float *pDst) {
  unsigned int numSamples;
  unsigned int blkCnt;
  float in1, in2, in3, in4;
  numSamples = (unsigned int)numRows * numCols;
  blkCnt = numSamples >> 2u;
  while (blkCnt > 0u) {
    in1 = *pSrc++;
    in2 = *pSrc++;
    in3 = *pSrc++;
    in4 = *pSrc++;
    *pDst++ = in1;
    *pDst++ = in2;
    *pDst++ = in3;
    *pDst++ = in4;
    blkCnt--;
  }
  blkCnt = numSamples & 0x3u;
  while (blkCnt > 0u) {
    *pDst++ = *pSrc++;
    blkCnt--;
  }
}


x86 assembly (clang, -O0):
<Matrix_Copy>:
push   %rbp
mov    %rsp,%rbp
mov    %dx,%ax
mov    %si,%dx
mov    %rdi,-0x8(%rbp)
mov    %dx,-0xa(%rbp)
mov    %ax,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
movzwl -0xa(%rbp),%eax
movzwl -0xc(%rbp),%ecx
imul   %ecx,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
shr    $0x2,%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
jbe    1209 <Matrix_Copy+0x109>
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x8(%rbp)
movss  (%rax),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x8(%rbp)
movss  (%rax),%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x8(%rbp)
movss  (%rax),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x8(%rbp)
movss  (%rax),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x24(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
movss  %xmm0,(%rax)
movss  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
movss  %xmm0,(%rax)
movss  -0x2c(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
movss  %xmm0,(%rax)
movss  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
movss  %xmm0,(%rax)
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    1131 <Matrix_Copy+0x31>
mov    -0x1c(%rbp),%eax
and    $0x3,%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
jbe    1250 <Matrix_Copy+0x150>
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x8(%rbp)
movss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
movss  %xmm0,(%rax)
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    1212 <Matrix_Copy+0x112>
pop    %rbp
ret

x86 assembly (clang, -O2):
<Matrix_Copy>:
imul   %edx,%esi
cmp    $0x4,%esi
jb     118b <Matrix_Copy+0x8b>
mov    %esi,%eax
shr    $0x2,%eax
lea    -0x1(%rax),%r8d
mov    %eax,%r9d
and    $0x7,%r9d
je     1138 <Matrix_Copy+0x38>
xor    %edx,%edx
movups (%rdi),%xmm0
add    $0x10,%rdi
movups %xmm0,(%rcx)
add    $0x10,%rcx
add    $0x1,%edx
cmp    %edx,%r9d
jne    1120 <Matrix_Copy+0x20>
sub    %edx,%eax
cmp    $0x7,%r8d
jb     118b <Matrix_Copy+0x8b>
xchg   %ax,%ax
movups (%rdi),%xmm0
movups %xmm0,(%rcx)
movups 0x10(%rdi),%xmm0
movups %xmm0,0x10(%rcx)
movups 0x20(%rdi),%xmm0
movups %xmm0,0x20(%rcx)
movups 0x30(%rdi),%xmm0
movups %xmm0,0x30(%rcx)
movups 0x40(%rdi),%xmm0
movups %xmm0,0x40(%rcx)
movups 0x50(%rdi),%xmm0
movups %xmm0,0x50(%rcx)
movups 0x60(%rdi),%xmm0
movups %xmm0,0x60(%rcx)
movups 0x70(%rdi),%xmm0
sub    $0xffffffffffffff80,%rdi
movups %xmm0,0x70(%rcx)
sub    $0xffffffffffffff80,%rcx
add    $0xfffffff8,%eax
jne    1140 <Matrix_Copy+0x40>
and    $0x3,%esi
je     11b6 <Matrix_Copy+0xb6>
movss  (%rdi),%xmm0
movss  %xmm0,(%rcx)
cmp    $0x1,%esi
je     11b6 <Matrix_Copy+0xb6>
movss  0x4(%rdi),%xmm0
movss  %xmm0,0x4(%rcx)
cmp    $0x2,%esi
je     11b6 <Matrix_Copy+0xb6>
movss  0x8(%rdi),%xmm0
movss  %xmm0,0x8(%rcx)
ret

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

LLVM IR (-O2):
define void @Matrix_Copy(float* nocapture %0, i16 zeroext %1, i16 zeroext %2, float* nocapture writeonly %3) {
%5 = zext i16 %1 to i32
%6 = zext i16 %2 to i32
%7 = mul nuw i32 %6, %5
%8 = icmp ult i32 %7, 4
br i1 %8, label %80, label %9
9: ; preds = %4
%10 = lshr i32 %7, 2
%11 = add nsw i32 %10, -1
%12 = and i32 %10, 7
%13 = icmp eq i32 %12, 0
br i1 %13, label %27, label %14
14: ; preds = %9, %14
%15 = phi float* [ %19, %14 ], [ %0, %9 ]
%16 = phi float* [ %22, %14 ], [ %3, %9 ]
%17 = phi i32 [ %24, %14 ], [ %10, %9 ]
%18 = phi i32 [ %25, %14 ], [ 0, %9 ]
%19 = getelementptr inbounds float, float* %15, i64 4
%20 = bitcast float* %15 to <4 x float>*
%21 = load <4 x float>, <4 x float>* %20
%22 = getelementptr inbounds float, float* %16, i64 4
%23 = bitcast float* %16 to <4 x float>*
store <4 x float> %21, <4 x float>* %23
%24 = add i32 %17, -1
%25 = add i32 %18, 1
%26 = icmp eq i32 %25, %12
br i1 %26, label %27, label %14
27: ; preds = %14, %9
%28 = phi float* [ undef, %9 ], [ %19, %14 ]
%29 = phi float* [ undef, %9 ], [ %22, %14 ]
%30 = phi float* [ %0, %9 ], [ %19, %14 ]
%31 = phi float* [ %3, %9 ], [ %22, %14 ]
%32 = phi i32 [ %10, %9 ], [ %24, %14 ]
%33 = icmp ult i32 %11, 7
br i1 %33, label %80, label %34
34: ; preds = %27, %34
%35 = phi float* [ %73, %34 ], [ %30, %27 ]
%36 = phi float* [ %76, %34 ], [ %31, %27 ]
%37 = phi i32 [ %78, %34 ], [ %32, %27 ]
%38 = getelementptr inbounds float, float* %35, i64 4
%39 = bitcast float* %35 to <4 x float>*
%40 = load <4 x float>, <4 x float>* %39
%41 = getelementptr inbounds float, float* %36, i64 4
%42 = bitcast float* %36 to <4 x float>*
store <4 x float> %40, <4 x float>* %42
%43 = getelementptr inbounds float, float* %35, i64 8
%44 = bitcast float* %38 to <4 x float>*
%45 = load <4 x float>, <4 x float>* %44
%46 = getelementptr inbounds float, float* %36, i64 8
%47 = bitcast float* %41 to <4 x float>*
store <4 x float> %45, <4 x float>* %47
%48 = getelementptr inbounds float, float* %35, i64 12
%49 = bitcast float* %43 to <4 x float>*
%50 = load <4 x float>, <4 x float>* %49
%51 = getelementptr inbounds float, float* %36, i64 12
%52 = bitcast float* %46 to <4 x float>*
store <4 x float> %50, <4 x float>* %52
%53 = getelementptr inbounds float, float* %35, i64 16
%54 = bitcast float* %48 to <4 x float>*
%55 = load <4 x float>, <4 x float>* %54
%56 = getelementptr inbounds float, float* %36, i64 16
%57 = bitcast float* %51 to <4 x float>*
store <4 x float> %55, <4 x float>* %57
%58 = getelementptr inbounds float, float* %35, i64 20
%59 = bitcast float* %53 to <4 x float>*
%60 = load <4 x float>, <4 x float>* %59
%61 = getelementptr inbounds float, float* %36, i64 20
%62 = bitcast float* %56 to <4 x float>*
store <4 x float> %60, <4 x float>* %62
%63 = getelementptr inbounds float, float* %35, i64 24
%64 = bitcast float* %58 to <4 x float>*
%65 = load <4 x float>, <4 x float>* %64
%66 = getelementptr inbounds float, float* %36, i64 24
%67 = bitcast float* %61 to <4 x float>*
store <4 x float> %65, <4 x float>* %67
%68 = getelementptr inbounds float, float* %35, i64 28
%69 = bitcast float* %63 to <4 x float>*
%70 = load <4 x float>, <4 x float>* %69
%71 = getelementptr inbounds float, float* %36, i64 28
%72 = bitcast float* %66 to <4 x float>*
store <4 x float> %70, <4 x float>* %72
%73 = getelementptr inbounds float, float* %35, i64 32
%74 = bitcast float* %68 to <4 x float>*
%75 = load <4 x float>, <4 x float>* %74
%76 = getelementptr inbounds float, float* %36, i64 32
%77 = bitcast float* %71 to <4 x float>*
store <4 x float> %75, <4 x float>* %77
%78 = add i32 %37, -8
%79 = icmp eq i32 %78, 0
br i1 %79, label %80, label %34
80: ; preds = %27, %34, %4
%81 = phi float* [ %3, %4 ], [ %29, %27 ], [ %76, %34 ]
%82 = phi float* [ %0, %4 ], [ %28, %27 ], [ %73, %34 ]
%83 = and i32 %7, 3
%84 = icmp eq i32 %83, 0
br i1 %84, label %97, label %85
85: ; preds = %80
%86 = load float, float* %82
store float %86, float* %81
%87 = icmp eq i32 %83, 1
br i1 %87, label %97, label %88
88: ; preds = %85
%89 = getelementptr inbounds float, float* %81, i64 1
%90 = getelementptr inbounds float, float* %82, i64 1
%91 = load float, float* %90
store float %91, float* %89
%92 = icmp eq i32 %83, 2
br i1 %92, label %97, label %93
93: ; preds = %88
%94 = getelementptr inbounds float, float* %81, i64 2
%95 = getelementptr inbounds float, float* %82, i64 2
%96 = load float, float* %95
store float %96, float* %94
br label %97
97: ; preds = %85, %88, %93, %80
ret void
}