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

```asm
<MyGemm>:
  endbr64
  push   %r15
  mov    %esi,%eax
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %edi,0x30(%rsp)
  mov    %edx,0x34(%rsp)
  mov    %rcx,0x50(%rsp)
  mov    %r9,0x58(%rsp)
  test   %esi,%esi
  jle    L56
  test   %edi,%edi
  jle    L56
  mov    %edx,%esi
  test   %edx,%edx
  jle    L56
  mov    0xb0(%rsp),%ecx
  lea    -0x4(%rax),%edx
  mov    %eax,0x38(%rsp)
  sub    $0x1,%eax
  and    $0xfffffffc,%eax
  movl   $0x0,0x40(%rsp)
  mov    %r8d,%r13d
  shl    $0x2,%ecx
  sub    %eax,%edx
  lea    0x0(,%r8,4),%eax
  movl   $0x0,0x3c(%rsp)
  mov    %ecx,0x44(%rsp)
  mov    0xa0(%rsp),%ecx
  cltq
  shl    $0x3,%rax
  mov    %edx,0x4c(%rsp)
  lea    -0x1(%rdi),%edx
  shl    $0x2,%ecx
  and    $0xfffffffc,%edx
  mov    %rax,0x10(%rsp)
  lea    -0x1(%rsi),%eax
  mov    %ecx,0x48(%rsp)
  lea    -0x4(%rsi),%ecx
  lea    0x4(%rdx),%edi
  and    $0xfffffffc,%eax
  mov    %ecx,%edx
  mov    %edi,0x2c(%rsp)
  sub    %eax,%edx
  mov    %edx,0xc(%rsp)
  nopl   0x0(%rax,%rax,1)
L31:
  mov    0x38(%rsp),%eax
  mov    $0x4,%r14d
  mov    0xa8(%rsp),%rcx
  movl   $0x0,0x28(%rsp)
  cmp    %r14d,%eax
  cmovle %eax,%r14d
  mov    0x50(%rsp),%rax
  mov    %rax,0x18(%rsp)
  movslq 0x3c(%rsp),%rax
  lea    (%rcx,%rax,8),%rax
  mov    0x58(%rsp),%rcx
  mov    %rax,(%rsp)
  movslq 0x40(%rsp),%rax
  lea    (%rcx,%rax,8),%rax
  mov    %rax,0x20(%rsp)
  nopl   0x0(%rax,%rax,1)
L32:
  mov    $0x4,%eax
  mov    0x30(%rsp),%r12d
  sub    0x28(%rsp),%r12d
  mov    0x20(%rsp),%r15
  mov    0x18(%rsp),%rbp
  cmp    %eax,%r12d
  cmovg  %eax,%r12d
  mov    0x34(%rsp),%ebx
  mov    %r13d,%eax
  mov    %rbp,%r13
  mov    %r15,%rbp
  mov    %eax,%r15d
  nopl   0x0(%rax,%rax,1)
L16:
  mov    $0x4,%edx
  mov    %rbp,%r9
  mov    %r13,%rcx
  mov    %r15d,%r8d
  cmp    %edx,%ebx
  mov    %r14d,%esi
  mov    %r12d,%edi
  cmovle %ebx,%edx
  sub    $0x8,%rsp
  sub    $0x4,%ebx
  add    $0x20,%rbp
  mov    0xb8(%rsp),%eax
  push   %rax
  push   0x10(%rsp)
  mov    0xb8(%rsp),%eax
  push   %rax
  xor    %eax,%eax
  call   <Gemm_P_Ger@plt>
  add    0x30(%rsp),%r13
  add    $0x20,%rsp
  cmp    0xc(%rsp),%ebx
  jne    L16
  addl   $0x4,0x28(%rsp)
  mov    %r15d,%r13d
  mov    0x28(%rsp),%eax
  addq   $0x20,0x18(%rsp)
  addq   $0x20,(%rsp)
  cmp    0x2c(%rsp),%eax
  jne    L32
  subl   $0x4,0x38(%rsp)
  mov    0x44(%rsp),%esi
  mov    0x48(%rsp),%edi
  add    %esi,0x3c(%rsp)
  add    %edi,0x40(%rsp)
  mov    0x38(%rsp),%eax
  cmp    0x4c(%rsp),%eax
  jne    L31
L56:
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void MyGemm(int m, int n, int k, double *A, int ldA, double *B, int ldB, double *C, int ldC) {
  for (int j = 0; j < n; j += 4) {
    int jb = ((n - j) < (4) ? (n - j) : (4));
    for (int i = 0; i < m; i += 4) {
      int ib = ((m - i) < (4) ? (m - i) : (4));
      for (int p = 0; p < k; p += 4) {
        int pb = ((k - p) < (4) ? (k - p) : (4));
        Gemm_P_Ger(ib, jb, pb, &A[(p)*ldA + (i)], ldA, &B[(j)*ldB + (p)], ldB, &C[(j)*ldC + (i)], ldC);
      }
    }
  }
}
```