user: Demonstrate the C code that generates the assembly output below:

```asm
<matrix_multiply>:
  endbr64
  push   %r15
  mov    %r9d,%r15d
  push   %r14
  mov    %r8d,%r14d
  push   %r13
  mov    %rdi,%r13
  push   %r12
  mov    %rdx,%r12
  mov    %ecx,%edx
  push   %rbp
  imul   %r8d,%edx
  push   %rbx
  mov    %ecx,%ebx
  shl    $0x3,%edx
  sub    $0x48,%rsp
  mov    %rsi,0x28(%rsp)
  xor    %esi,%esi
  mov    %ecx,0x10(%rsp)
  call   <memset@plt>
  mov    %r14d,%eax
  or     %r15d,%eax
  or     %ebx,%eax
  test   $0x1f,%al
  jne    L11
  test   %r15d,%r15d
  je     L85
  mov    0x10(%rsp),%eax
  test   %eax,%eax
  je     L85
  test   %r14d,%r14d
  je     L85
  mov    %r14d,%eax
  movl   $0x0,0x3c(%rsp)
  shl    $0x5,%eax
  movl   $0x0,0x38(%rsp)
  mov    %eax,0x30(%rsp)
  mov    %r15d,%eax
  shl    $0x5,%eax
  mov    %r12,0x8(%rsp)
  mov    %eax,0x34(%rsp)
  mov    %r13,(%rsp)
  mov    %r15d,%r13d
  mov    %r14d,%r15d
L87:
  mov    0x38(%rsp),%eax
  movl   $0x0,0x18(%rsp)
  movl   $0x0,0x1c(%rsp)
  mov    0x3c(%rsp),%r12d
  mov    %rax,0x20(%rsp)
  movl   $0x0,0x14(%rsp)
L122:
  mov    0x28(%rsp),%rbx
  mov    0x1c(%rsp),%eax
  xor    %r14d,%r14d
  add    0x20(%rsp),%rax
  lea    (%rbx,%rax,8),%rbp
  mov    %r13d,%eax
  mov    0x18(%rsp),%ebx
  mov    %r14d,%r13d
  mov    %r15d,%r14d
  mov    %eax,%r15d
  nopl   (%rax)
L63:
  mov    0x8(%rsp),%rcx
  mov    %r13d,%eax
  mov    %r14d,%r9d
  mov    %r15d,%r8d
  lea    (%r12,%rax,1),%rdx
  add    %rbx,%rax
  mov    %rbp,%rsi
  add    $0x20,%r13d
  lea    (%rcx,%rdx,8),%rdx
  mov    (%rsp),%rcx
  lea    (%rcx,%rax,8),%rdi
  mov    %r14d,%ecx
  call   <matrix_multiply_block@plt>
  cmp    %r13d,%r14d
  ja     L63
  mov    0x34(%rsp),%esi
  addl   $0x20,0x14(%rsp)
  mov    %r15d,%r13d
  mov    %r14d,%r15d
  add    %esi,0x1c(%rsp)
  mov    0x14(%rsp),%eax
  mov    0x30(%rsp),%esi
  add    %esi,0x18(%rsp)
  cmp    %eax,0x10(%rsp)
  ja     L122
  addl   $0x20,0x38(%rsp)
  mov    0x38(%rsp),%eax
  add    %esi,0x3c(%rsp)
  cmp    %eax,%r13d
  ja     L87
L85:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L11:
  test   %r15d,%r15d
  je     L85
  mov    0x10(%rsp),%edx
  test   %edx,%edx
  je     L85
  movl   $0x0,0x8(%rsp)
  xor    %r8d,%r8d
  lea    (%r14,%r14,1),%ecx
  test   %r14d,%r14d
  je     L85
  mov    %r8d,(%rsp)
  mov    0x28(%rsp),%rbp
  mov    %r8d,%r11d
  nopl   0x0(%rax,%rax,1)
L90:
  mov    0x8(%rsp),%edi
  mov    %r14d,%r8d
  xor    %r10d,%r10d
  nopw   0x0(%rax,%rax,1)
L130:
  mov    %r8d,%ebx
  mov    %r11d,%eax
  sub    %r14d,%ebx
  lea    0x0(%rbp,%rax,8),%r9
  mov    %ebx,%eax
L71:
  mov    %eax,%edx
  lea    (%rax,%rdi,1),%esi
  add    $0x1,%eax
  movsd  (%r12,%rsi,8),%xmm0
  lea    0x0(%r13,%rdx,8),%rdx
  mulsd  (%r9),%xmm0
  addsd  (%rdx),%xmm0
  movsd  %xmm0,(%rdx)
  cmp    %eax,%r8d
  jne    L71
  add    $0x1,%r10d
  add    %r15d,%r11d
  lea    (%rcx,%rbx,1),%r8d
  sub    %r14d,%edi
  cmp    %r10d,0x10(%rsp)
  jne    L130
  addl   $0x1,(%rsp)
  mov    (%rsp),%eax
  add    %r14d,0x8(%rsp)
  cmp    %eax,%r15d
  je     L85
  mov    %eax,%r11d
  jmp    L90
```

assistant: ```c
void matrix_multiply(double C[], const double A[], const double B[], const unsigned int m, const unsigned int n, const unsigned int q) {
  unsigned int i, j, k;
  memset(C, 0, m * n * sizeof(double));
  if (m % 32 || n % 32 || q % 32) {
    for (k = 0; k < q; ++k) {
      for (i = 0; i < m; ++i) {
        for (j = 0; j < n; ++j) {
          C[i * n + j] += A[i * q + k] * B[k * n + j];
        }
      }
    }
    return;
  }
  for (k = 0; k < q; k += 32) {
    const double *B_k = B + k * n;
    for (i = 0; i < m; i += 32) {
      const double *A_i_k = A + i * q + k;
      double *C_i = C + i * n;
      for (j = 0; j < n; j += 32) {
        matrix_multiply_block(C_i + j, A_i_k, B_k + j, n, q, n);
      }
    }
  }
}
```