user: Translate the below assembly code into C programming language:

```asm
<mulMat_park>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rsi,-0x20(%rsp)
  mov    %rdx,-0x18(%rsp)
  mov    %rcx,-0x10(%rsp)
  test   %edi,%edi
  jle    L21
  mov    %edi,%eax
  movl   $0x0,-0x50(%rsp)
  mov    %edi,%r8d
  movq   $0x0,-0x48(%rsp)
  shl    $0x5,%eax
  mov    %eax,-0x8(%rsp)
  cltq
  mov    %rax,-0x30(%rsp)
  movslq %edi,%rax
  lea    0x0(,%rax,8),%r14
  lea    -0x1(%rdi),%eax
  and    $0xffffffe0,%eax
  mov    %r14,%r15
  add    $0x3f,%eax
  mov    %eax,-0x4(%rsp)
L43:
  mov    -0x48(%rsp),%rax
  mov    -0x10(%rsp),%rbx
  movq   $0x0,-0x58(%rsp)
  movl   $0x0,-0x4c(%rsp)
  shl    $0x3,%rax
  add    %rax,%rbx
  add    -0x20(%rsp),%rax
  mov    %rax,-0x40(%rsp)
  mov    -0x50(%rsp),%eax
  mov    %rbx,-0x28(%rsp)
  add    $0x1f,%eax
  mov    %eax,-0x5c(%rsp)
L37:
  mov    -0x58(%rsp),%rax
  mov    %eax,%r13d
  cmp    %eax,%r8d
  jle    L104
  mov    %rax,%rbx
  mov    -0x40(%rsp),%rax
  mov    $0x1f,%r9d
  mov    $0x1,%r11d
  lea    (%rax,%rbx,8),%rax
  mov    -0x18(%rsp),%rbx
  mov    %rax,-0x38(%rsp)
  movslq -0x4c(%rsp),%rax
  lea    (%rbx,%rax,8),%r14
  lea    0x1f(%r13),%ebx
L38:
  mov    -0x38(%rsp),%r12
  mov    -0x28(%rsp),%rcx
  mov    -0x50(%rsp),%ebp
  xchg   %ax,%ax
L55:
  mov    %r14,%rdi
  mov    %r12,%rsi
  mov    %r13d,%r10d
  nopl   0x0(%rax)
L105:
  mov    %r11,%rax
  jmp    L19
  nopl   (%rax)
L39:
  cmp    %edx,%r8d
  jle    L68
L19:
  movsd  (%rsi),%xmm0
  mulsd  -0x8(%rdi,%rax,8),%xmm0
  mov    %eax,%edx
  addsd  -0x8(%rcx,%rax,8),%xmm0
  movsd  %xmm0,-0x8(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %edx,%r9d
  jge    L39
L68:
  add    $0x1,%r10d
  add    $0x8,%rsi
  add    %r15,%rdi
  cmp    %ebx,%r10d
  jg     L91
  cmp    %r10d,%r8d
  jg     L105
L91:
  add    $0x1,%ebp
  add    %r15,%rcx
  add    %r15,%r12
  cmp    -0x5c(%rsp),%ebp
  jg     L57
  cmp    %ebp,%r8d
  jg     L55
L57:
  add    $0x20,%r11
  add    $0x20,%r9d
  cmp    %r9d,-0x4(%rsp)
  jne    L38
L104:
  addq   $0x20,-0x58(%rsp)
  mov    -0x58(%rsp),%rax
  mov    -0x8(%rsp),%ebx
  add    %ebx,-0x4c(%rsp)
  cmp    %eax,%r8d
  jg     L37
  addl   $0x20,-0x50(%rsp)
  mov    -0x50(%rsp),%eax
  mov    -0x30(%rsp),%rsi
  add    %rsi,-0x48(%rsp)
  cmp    %eax,%r8d
  jg     L43
L21:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void mulMat_park(const int matsize, const double *mata, const double *matb, double *matc) {
  int i, j, k, ii, jj, kk;
  int ixmatsize;
  int kxmatsize;
  for (ii = 0; ii < matsize; ii += 32) {
    for (kk = 0; kk < matsize; kk += 32) {
      for (jj = 0; jj < matsize; jj += 32) {
        ixmatsize = ii * matsize;
        for (i = ii; i < ii + 32 && i < matsize; i++) {
          kxmatsize = kk * matsize;
          for (k = kk; k < kk + 32 && k < matsize; k++) {
            for (j = jj; j < jj + 32 && j < matsize; j++) {
              matc[ixmatsize + j] += mata[ixmatsize + k] * matb[kxmatsize + j];
            }
            kxmatsize += matsize;
          }
          ixmatsize += matsize;
        }
      }
    }
  }
  return;
}
```