user: What is the c source code of the assembly code below?

```asm
<Matrix_Multiply>:
  endbr64
  push   %r15
  mov    %rdi,%r10
  mov    %edx,%r11d
  movzwl %r8w,%edi
  push   %r14
  movzwl %dx,%eax
  and    $0x3,%r11d
  push   %r13
  push   %r12
  mov    %edx,%r12d
  push   %rbp
  shr    $0x2,%r12w
  lea    0x0(,%rax,4),%rbp
  push   %rbx
  mov    %rcx,%rbx
  lea    0x0(,%rdi,4),%rcx
  mov    %r9,-0x10(%rsp)
  mov    %si,-0x14(%rsp)
  cmp    $0x1,%r8w
  jne    L147
  lea    -0x1(%r12),%eax
  lea    (%rdi,%rdi,2),%r8
  shl    $0x4,%rdi
  xor    %r9d,%r9d
  movzwl %ax,%eax
  shl    $0x2,%r8
  pxor   %xmm2,%xmm2
  add    $0x1,%rax
  mov    %rax,%rsi
  imul   %rdi,%rax
  shl    $0x4,%rsi
  add    %r10,%rsi
  lea    (%rbx,%rax,1),%r15
  mov    %r15,-0x20(%rsp)
  mov    -0x10(%rsp),%r15
  nopw   0x0(%rax,%rax,1)
L121:
  movzwl %r9w,%eax
  lea    (%r15,%rax,4),%r14
  test   %r12w,%r12w
  je     L123
  mov    %rbx,%rdx
  mov    %r10,%rax
  movaps %xmm2,%xmm0
  mov    %rsi,%r13
  xchg   %ax,%ax
L130:
  movss  (%rax),%xmm1
  mulss  (%rdx),%xmm1
  add    $0x10,%rax
  addss  %xmm1,%xmm0
  movss  (%rdx,%rcx,1),%xmm1
  mulss  -0xc(%rax),%xmm1
  addss  %xmm1,%xmm0
  movss  -0x8(%rax),%xmm1
  mulss  (%rdx,%rcx,2),%xmm1
  addss  %xmm0,%xmm1
  movss  (%rdx,%r8,1),%xmm0
  mulss  -0x4(%rax),%xmm0
  add    %rdi,%rdx
  addss  %xmm1,%xmm0
  cmp    %rsi,%rax
  jne    L130
  mov    -0x20(%rsp),%rax
L95:
  test   %r11w,%r11w
  je     L17
  movss  (%rax),%xmm1
  mulss  0x0(%r13),%xmm1
  add    %rcx,%rax
  addss  %xmm1,%xmm0
  cmp    $0x1,%r11w
  je     L17
  movss  0x4(%r13),%xmm1
  mulss  (%rax),%xmm1
  add    %rcx,%rax
  addss  %xmm1,%xmm0
  cmp    $0x2,%r11w
  je     L17
  movss  (%rax),%xmm1
  mulss  0x8(%r13),%xmm1
  addss  %xmm1,%xmm0
L17:
  movss  %xmm0,(%r14)
  add    $0x1,%r9d
  add    %rbp,%r10
  add    %rbp,%rsi
  cmp    %r9w,-0x14(%rsp)
  jne    L121
  pop    %rbx
  xor    %eax,%eax
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L123:
  mov    %rbx,%rax
  mov    %r10,%r13
  movaps %xmm2,%xmm0
  jmp    L95
L147:
  lea    -0x1(%r12),%eax
  lea    (%rdi,%rdi,2),%r9
  shl    $0x4,%rdi
  mov    %r8d,%r15d
  movzwl %ax,%eax
  mov    %rbp,-0x8(%rsp)
  shl    $0x2,%r9
  mov    %rdi,%rsi
  add    $0x1,%rax
  pxor   %xmm2,%xmm2
  mov    %rax,%r8
  imul   %rdi,%rax
  shl    $0x4,%r8
  add    %r10,%r8
  mov    %rax,-0x20(%rsp)
  xor    %eax,%eax
  mov    %ax,-0x12(%rsp)
L18:
  mov    -0x10(%rsp),%rdi
  movzwl -0x12(%rsp),%eax
  xor    %ebp,%ebp
  lea    (%rdi,%rax,4),%r13
  mov    %rbx,%rdi
  nopl   0x0(%rax,%rax,1)
L1:
  test   %r12w,%r12w
  je     L87
  mov    %rdi,%rdx
  mov    %r10,%rax
  movaps %xmm2,%xmm0
  mov    %r8,%r14
  xchg   %ax,%ax
L126:
  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    %rsi,%rdx
  addss  %xmm1,%xmm0
  cmp    %r8,%rax
  jne    L126
  add    -0x20(%rsp),%rdi
L19:
  test   %r11w,%r11w
  je     L70
  movss  (%r14),%xmm1
  mulss  (%rdi),%xmm1
  add    %rcx,%rdi
  addss  %xmm1,%xmm0
  cmp    $0x1,%r11w
  je     L70
  movss  0x4(%r14),%xmm1
  mulss  (%rdi),%xmm1
  add    %rcx,%rdi
  addss  %xmm1,%xmm0
  cmp    $0x2,%r11w
  je     L70
  movss  0x8(%r14),%xmm1
  mulss  (%rdi),%xmm1
  addss  %xmm1,%xmm0
L70:
  add    $0x1,%ebp
  add    $0x4,%r13
  movss  %xmm0,-0x4(%r13)
  movzwl %bp,%eax
  lea    (%rbx,%rax,4),%rdi
  cmp    %bp,%r15w
  jne    L1
  mov    -0x8(%rsp),%rax
  add    %r15w,-0x12(%rsp)
  add    %rax,%r10
  add    %rax,%r8
  subw   $0x1,-0x14(%rsp)
  jne    L18
  pop    %rbx
  xor    %eax,%eax
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L87:
  mov    %r10,%r14
  movaps %xmm2,%xmm0
  jmp    L19
```

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);
}
```