void Matrix_Zero(float *A, unsigned short numRows, unsigned short numCols) {
  float *pIn = A;
  unsigned int numSamples = numRows * numCols;
  unsigned int blkCnt = numSamples >> 2u;
  while (blkCnt > 0u) {
    (*pIn++) = 0.0f;
    (*pIn++) = 0.0f;
    (*pIn++) = 0.0f;
    (*pIn++) = 0.0f;
    blkCnt--;
  }
  blkCnt = numSamples & 0x03u;
  while (blkCnt > 0u) {
    (*pIn++) = 0.0f;
    blkCnt--;
  }
}


x86 assembly (clang, -O0):
<Matrix_Zero>:
push   %rbp
mov    %rsp,%rbp
mov    %dx,%ax
mov    %si,%cx
mov    %rdi,-0x8(%rbp)
mov    %cx,-0xa(%rbp)
mov    %ax,-0xc(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-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    11a5 <Matrix_Zero+0xa5>
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    1135 <Matrix_Zero+0x35>
mov    -0x1c(%rbp),%eax
and    $0x3,%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
jbe    11dc <Matrix_Zero+0xdc>
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    11ae <Matrix_Zero+0xae>
pop    %rbp
ret

x86 assembly (clang, -O2):
<Matrix_Zero>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
mov    %edx,%r12d
mov    %esi,%r15d
mov    %rdi,%r14
mov    %edx,%ebp
imul   %esi,%ebp
cmp    $0x4,%ebp
jb     114c <Matrix_Zero+0x3c>
mov    %ebp,%ebx
shr    $0x2,%ebx
add    $0xffffffff,%ebx
shl    $0x4,%rbx
lea    0x10(%rbx),%rdx
mov    %r14,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
add    %rbx,%r14
add    $0x10,%r14
test   $0x3,%bpl
je     116b <Matrix_Zero+0x5b>
imul   %r15d,%r12d
shl    $0x2,%r12d
and    $0xc,%r12d
mov    %r14,%rdi
xor    %esi,%esi
mov    %r12,%rdx
call   1030 <memset@plt>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @Matrix_Zero(float* %0, i16 zeroext %1, i16 zeroext %2) {
%4 = alloca float*
%5 = alloca i16
%6 = alloca i16
%7 = alloca float*
%8 = alloca i32
%9 = alloca i32
store float* %0, float** %4
store i16 %1, i16* %5
store i16 %2, i16* %6
%10 = load float*, float** %4
store float* %10, float** %7
%11 = load i16, i16* %5
%12 = zext i16 %11 to i32
%13 = load i16, i16* %6
%14 = zext i16 %13 to i32
%15 = mul nsw i32 %12, %14
store i32 %15, i32* %8
%16 = load i32, i32* %8
%17 = lshr i32 %16, 2
store i32 %17, i32* %9
br label %18
18: ; preds = %21, %3
%19 = load i32, i32* %9
%20 = icmp ugt i32 %19, 0
br i1 %20, label %21, label %32
21: ; preds = %18
%22 = load float*, float** %7
%23 = getelementptr inbounds float, float* %22, i32 1
store float* %23, float** %7
store float 0.000000e+00, float* %22
%24 = load float*, float** %7
%25 = getelementptr inbounds float, float* %24, i32 1
store float* %25, float** %7
store float 0.000000e+00, float* %24
%26 = load float*, float** %7
%27 = getelementptr inbounds float, float* %26, i32 1
store float* %27, float** %7
store float 0.000000e+00, float* %26
%28 = load float*, float** %7
%29 = getelementptr inbounds float, float* %28, i32 1
store float* %29, float** %7
store float 0.000000e+00, float* %28
%30 = load i32, i32* %9
%31 = add i32 %30, -1
store i32 %31, i32* %9
br label %18
32: ; preds = %18
%33 = load i32, i32* %8
%34 = and i32 %33, 3
store i32 %34, i32* %9
br label %35
35: ; preds = %38, %32
%36 = load i32, i32* %9
%37 = icmp ugt i32 %36, 0
br i1 %37, label %38, label %43
38: ; preds = %35
%39 = load float*, float** %7
%40 = getelementptr inbounds float, float* %39, i32 1
store float* %40, float** %7
store float 0.000000e+00, float* %39
%41 = load i32, i32* %9
%42 = add i32 %41, -1
store i32 %42, i32* %9
br label %35
43: ; preds = %35
ret void
}

LLVM IR (-O2):
define void @Matrix_Zero(float* nocapture writeonly %0, i16 zeroext %1, i16 zeroext %2) {
%4 = bitcast float* %0 to i8*
%5 = zext i16 %1 to i32
%6 = zext i16 %2 to i32
%7 = mul nuw nsw i32 %6, %5
%8 = icmp ult i32 %7, 4
br i1 %8, label %19, label %9
9: ; preds = %3
%10 = lshr i32 %7, 2
%11 = add nsw i32 %10, -1
%12 = zext i32 %11 to i64
%13 = shl nuw nsw i64 %12, 4
%14 = add nuw nsw i64 %13, 16
%15 = shl nuw nsw i64 %12, 2
%16 = add nuw nsw i64 %15, 4
%17 = getelementptr float, float* %0, i64 %16
%18 = bitcast float* %17 to i8*
br label %19
19: ; preds = %9, %3
%20 = phi i8* [ %18, %9 ], [ %4, %3 ]
%21 = and i32 %7, 3
%22 = icmp eq i32 %21, 0
br i1 %22, label %29, label %23
23: ; preds = %19
%24 = zext i16 %2 to i64
%25 = zext i16 %1 to i64
%26 = mul nuw nsw i64 %24, %25
%27 = shl nuw nsw i64 %26, 2
%28 = and i64 %27, 12
br label %29
29: ; preds = %23, %19
ret void
}