user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<Matrix_Multiply>:
  endbr64
  push   %r15
  movzwl %dx,%eax
  pxor   %xmm2,%xmm2
  push   %r14
  shl    $0x2,%rax
  mov    %edx,%r14d
  push   %r13
  mov    %edx,%r13d
  shr    $0x2,%r14w
  push   %r12
  and    $0x3,%r13d
  lea    -0x1(%r14),%r15d
  mov    %r8d,%r12d
  push   %rbp
  movzwl %r8w,%r8d
  movzwl %r15w,%r15d
  mov    %rcx,%rbp
  push   %rbx
  add    $0x1,%r15
  lea    0x0(,%r8,4),%rcx
  mov    %rdi,%rbx
  mov    %r15,%rdi
  shl    $0x4,%rdi
  mov    %rax,-0x18(%rsp)
  lea    -0x1(%r13),%eax
  add    %rbx,%rdi
  movzwl %ax,%eax
  mov    %r9,-0x8(%rsp)
  lea    (%r8,%r8,2),%r9
  shl    $0x4,%r8
  mov    %si,-0xa(%rsp)
  lea    0x4(,%rax,4),%rax
  shl    $0x2,%r9
  mov    %rax,-0x20(%rsp)
  xor    %eax,%eax
  imul   %r8,%r15
  mov    %ax,-0xc(%rsp)
  nopl   (%rax)
L5:
  mov    -0x8(%rsp),%rsi
  movzwl -0xc(%rsp),%eax
  xor    %r10d,%r10d
  lea    (%rsi,%rax,4),%r11
  mov    %rbp,%rsi
  nopl   0x0(%rax)
L42:
  test   %r14w,%r14w
  je     L63
  mov    %rsi,%rdx
  mov    %rbx,%rax
  movaps %xmm2,%xmm0
  nopl   0x0(%rax,%rax,1)
L59:
  movss  (%rdx),%xmm1
  mulss  (%rax),%xmm1
  add    $0x10,%rax
  addss  %xmm1,%xmm0
  movss  -0xc(%rax),%xmm1
  mulss  (%rdx,%rcx,1),%xmm1
  addss  %xmm1,%xmm0
  movss  (%rdx,%rcx,2),%xmm1
  mulss  -0x8(%rax),%xmm1
  addss  %xmm0,%xmm1
  movss  (%rdx,%r9,1),%xmm0
  mulss  -0x4(%rax),%xmm0
  add    %r8,%rdx
  addss  %xmm1,%xmm0
  cmp    %rdi,%rax
  jne    L59
  add    %r15,%rsi
L30:
  test   %r13w,%r13w
  je     L60
  mov    -0x20(%rsp),%rdx
  add    %rax,%rdx
L25:
  movss  (%rax),%xmm1
  mulss  (%rsi),%xmm1
  add    $0x4,%rax
  add    %rcx,%rsi
  addss  %xmm1,%xmm0
  cmp    %rdx,%rax
  jne    L25
L60:
  add    $0x1,%r10d
  add    $0x4,%r11
  movss  %xmm0,-0x4(%r11)
  movzwl %r10w,%eax
  lea    0x0(%rbp,%rax,4),%rsi
  cmp    %r10w,%r12w
  jne    L42
  mov    -0x18(%rsp),%rax
  add    %r12w,-0xc(%rsp)
  add    %rax,%rbx
  add    %rax,%rdi
  subw   $0x1,-0xa(%rsp)
  jne    L5
  pop    %rbx
  xor    %eax,%eax
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L63:
  mov    %rbx,%rax
  movaps %xmm2,%xmm0
  jmp    L30
```

assistant: ```c
int Matrix_Multiply(float *pSrcA, unsigned short numRowsA, unsigned short numColsA, float *pSrcB, unsigned short numColsB, float *pDst) {
  float *pIn1 = pSrcA;
  float *pIn2 = pSrcB;
  float *pInA = pSrcA;
  float *pOut = pDst;
  float *px;
  float sum;
  float in1, in2, in3, in4;
  unsigned short col, i = 0u, j, row = numRowsA, colCnt;
  int status;
  do {
    px = pOut + i;
    col = numColsB;
    pIn2 = pSrcB;
    j = 0u;
    do {
      sum = 0.0f;
      pIn1 = pInA;
      colCnt = numColsA >> 2u;
      while (colCnt > 0u) {
        in3 = *pIn2;
        pIn2 += numColsB;
        in1 = pIn1[0];
        in2 = pIn1[1];
        sum += in1 * in3;
        in4 = *pIn2;
        pIn2 += numColsB;
        sum += in2 * in4;
        in3 = *pIn2;
        pIn2 += numColsB;
        in1 = pIn1[2];
        in2 = pIn1[3];
        sum += in1 * in3;
        in4 = *pIn2;
        pIn2 += numColsB;
        sum += in2 * in4;
        pIn1 += 4u;
        colCnt--;
      }
      colCnt = numColsA & 0x3u;
      while (colCnt > 0u) {
        sum += *pIn1++ * (*pIn2);
        pIn2 += numColsB;
        colCnt--;
      }
      *px++ = sum;
      j++;
      pIn2 = pSrcB + j;
      col--;
    } while (col > 0u);
    i = i + numColsB;
    pInA = pInA + numColsA;
    row--;
  } while (row > 0u);
  status = 0;
  return (status);
}
```