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

```asm
<Matrix4Double>:
  endbr64
  push   %r15
  push   %r14
  mov    %edx,%r14d
  push   %r13
  mov    %esi,%r13d
  mov    $0x8,%esi
  push   %r12
  push   %rbp
  mov    %edi,%ebp
  push   %rbx
  mov    %ebp,%r15d
  mov    %ecx,%ebx
  imul   %r13d,%r15d
  sub    $0x8,%rsp
  call   <calloc@plt>
  mov    %r15d,%edi
  mov    $0x8,%esi
  mov    %rax,%r12
  imul   %r14d,%r15d
  call   <calloc@plt>
  mov    $0x8,%esi
  mov    %rax,(%r12)
  mov    %r15d,%edi
  call   <calloc@plt>
  mov    (%r12),%rdx
  mov    %r15d,%edi
  mov    $0x8,%esi
  imul   %ebx,%edi
  mov    %rax,(%rdx)
  call   <calloc@plt>
  mov    %rax,%r8
  mov    (%r12),%rax
  mov    (%rax),%rax
  mov    %r8,(%rax)
  cmp    $0x1,%ebp
  jle    L25
  mov    %r13d,%eax
  movslq %r13d,%r9
  lea    -0x1(%rbp),%edx
  imul   %r14d,%eax
  shl    $0x3,%r9
  lea    (%r12,%rdx,8),%r10
  movslq %eax,%r8
  imul   %ebx,%eax
  shl    $0x3,%r8
  movslq %eax,%rdi
  mov    %r12,%rax
  shl    $0x3,%rdi
  nopl   0x0(%rax,%rax,1)
L73:
  mov    %rax,%rdx
  add    $0x8,%rax
  mov    (%rdx),%rsi
  add    %r9,%rsi
  mov    %rsi,(%rax)
  mov    (%rdx),%rcx
  mov    (%rcx),%rdx
  add    %r8,%rdx
  mov    %rdx,(%rsi)
  mov    (%rcx),%rcx
  mov    (%rcx),%rsi
  add    %rdi,%rsi
  mov    %rsi,(%rdx)
  cmp    %r10,%rax
  jne    L73
L22:
  mov    %ebx,%edi
  movslq %r14d,%r8
  imul   %r14d,%edi
  shl    $0x3,%r8
  movslq %edi,%rdi
  shl    $0x3,%rdi
  cmp    $0x1,%r13d
  jle    L104
  lea    -0x1(%r13),%r11d
  mov    %r12,%r10
  xor    %r9d,%r9d
  shl    $0x3,%r11
  nopl   0x0(%rax)
L91:
  mov    (%r10),%rax
  lea    (%r11,%rax,1),%rsi
  nopw   0x0(%rax,%rax,1)
L48:
  mov    (%rax),%rcx
  add    $0x8,%rax
  mov    (%rcx),%r15
  lea    (%rcx,%r8,1),%rdx
  mov    %rdx,(%rax)
  add    %rdi,%r15
  mov    %r15,(%rdx)
  cmp    %rax,%rsi
  jne    L48
  add    $0x1,%r9d
  add    $0x8,%r10
  cmp    %r9d,%ebp
  jg     L91
  movslq %ebx,%rsi
  shl    $0x3,%rsi
L88:
  cmp    $0x1,%r14d
  jle    L28
  mov    %r13d,%r9d
  lea    -0x2(%r14),%eax
  mov    %r12,%r11
  xor    %ebx,%ebx
  shl    $0x3,%r9
  lea    0x10(,%rax,8),%r10
  xchg   %ax,%ax
L19:
  mov    (%r11),%r8
  xor    %edi,%edi
  nopl   (%rax)
L106:
  mov    (%r8,%rdi,1),%rcx
  mov    (%rcx),%rdx
  lea    0x8(%rcx),%rax
  add    %r10,%rcx
  cs nopw 0x0(%rax,%rax,1)
L49:
  add    %rsi,%rdx
  add    $0x8,%rax
  mov    %rdx,-0x8(%rax)
  cmp    %rax,%rcx
  jne    L49
  add    $0x8,%rdi
  cmp    %rdi,%r9
  jne    L106
  add    $0x1,%ebx
  add    $0x8,%r11
  cmp    %ebx,%ebp
  jg     L19
L28:
  add    $0x8,%rsp
  mov    %r12,%rax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L25:
  jne    L28
  jmp    L22
  nopl   (%rax)
L104:
  movslq %ebx,%rsi
  shl    $0x3,%rsi
  cmp    $0x1,%r13d
  jne    L28
  jmp    L88
```

assistant: ```c
double ****Matrix4Double(int M, int N, int P, int Q) {
  double ****out;
  int i, j, k;
  out = (double ****)calloc(M, sizeof(double ***));
  out[0] = (double ***)calloc(M * N, sizeof(double **));
  out[0][0] = (double **)calloc(M * N * P, sizeof(double *));
  out[0][0][0] = (double *)calloc(M * N * P * Q, sizeof(double));
  for (i = 1; i < M; i++) {
    out[i] = out[i - 1] + N;
    out[i][0] = out[i - 1][0] + N * P;
    out[i][0][0] = out[i - 1][0][0] + N * P * Q;
  }
  for (i = 0; i < M; i++)
    for (j = 1; j < N; j++) {
      out[i][j] = out[i][j - 1] + P;
      out[i][j][0] = out[i][j - 1][0] + P * Q;
    }
  for (i = 0; i < M; i++)
    for (j = 0; j < N; j++)
      for (k = 1; k < P; k++)
        out[i][j][k] = out[i][j][k - 1] + Q;
  return out;
}
```