user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<matprod_mat_mat>:
  endbr64
  push   %r15
  mov    %rsi,%r15
  push   %r14
  mov    %r8d,%r14d
  push   %r13
  push   %r12
  mov    %ecx,%r12d
  push   %rbp
  push   %rbx
  movslq %r9d,%rbx
  sub    $0x48,%rsp
  mov    %rdi,0x18(%rsp)
  mov    $0x1,%edi
  mov    %rdx,0x10(%rsp)
  call   <assert@plt>
  mov    $0x1,%edi
  call   <assert@plt>
  mov    $0x1,%edi
  call   <assert@plt>
  cmp    $0x1,%r14d
  jle    L7
  cmp    $0x2,%r12d
  jle    L222
  cmp    $0x1,%ebx
  jle    L188
  cmp    $0x3c0,%r12d
  jg     L156
  cmp    $0x20,%r14d
  jle    L146
L156:
  mov    %ebx,0x24(%rsp)
  cmp    $0x3c0,%r12d
  jg     L55
L82:
  mov    0x24(%rsp),%eax
  mov    %ebx,%r10d
  mov    %r14d,%r8d
  mov    %r12d,%ebx
  add    %eax,%eax
  mov    %eax,0x38(%rsp)
  movslq %r14d,%rax
  mov    %rax,0x28(%rsp)
  movslq %r12d,%rax
  mov    %rax,0x30(%rsp)
  jmp    L148
  nop
L111:
  cmp    %r10d,0x38(%rsp)
  jle    L173
  mov    %r10d,%ebp
  shr    $0x1f,%ebp
  add    %r10d,%ebp
  sar    %ebp
L173:
  add    $0x7,%ebp
  and    $0xfffffff8,%ebp
  cmp    %r10d,%ebp
  cmovg  %r10d,%ebp
  sub    %ebp,%r10d
  cmp    $0x3c0,%ebx
  jg     L198
  mov    0x10(%rsp),%r14
  mov    0x18(%rsp),%r12
  mov    %ebx,%r13d
L136:
  mov    %r10d,0x20(%rsp)
  mov    %ebx,%ecx
  mov    %r15,%rsi
  mov    %ebp,%r9d
  push   %rbx
  mov    %r14,%rdx
  mov    %r12,%rdi
  push   %r13
  mov    %r8d,0x1c(%rsp)
  call   <matprod_mat_mat_sub_xrows@plt>
  pop    %rcx
  pop    %rsi
  mov    0x20(%rsp),%r10d
  test   %r10d,%r10d
  je     L89
  mov    0x28(%rsp),%rax
  movslq %ebp,%rbp
  mov    0xc(%rsp),%r8d
  imul   %rbp,%rax
  imul   0x30(%rsp),%rbp
  lea    (%r15,%rax,8),%r15
  mov    0x10(%rsp),%rax
  lea    (%rax,%rbp,8),%rax
  mov    %rax,0x10(%rsp)
L148:
  mov    0x24(%rsp),%ebp
  cmp    %ebp,%r10d
  jg     L111
  cmp    $0x3c0,%ebx
  jle    L236
  movl   $0x0,0x20(%rsp)
L166:
  cmp    $0x77f,%ebx
  jle    L224
  mov    0x10(%rsp),%r14
  mov    0x18(%rsp),%r12
  mov    %ebx,%r13d
  mov    %r8d,%ebp
  mov    %r10d,%r9d
  nopl   0x0(%rax)
L13:
  push   %rbx
  mov    %r14,%rdx
  mov    %r12,%rdi
  mov    %ebp,%r8d
  push   $0x3c0
  mov    %ebx,%ecx
  mov    %r15,%rsi
  sub    $0x3c0,%r13d
  add    $0x1e00,%r12
  add    $0x1e00,%r14
  mov    %r9d,0x1c(%rsp)
  call   <matprod_mat_mat_sub_xrows@plt>
  cmp    $0x77f,%r13d
  pop    %r9
  pop    %r10
  mov    0xc(%rsp),%r9d
  jg     L13
  mov    %ebp,%r8d
  mov    %r9d,%r10d
  cmp    $0x3c0,%r13d
  jne    L38
  mov    0x20(%rsp),%r10d
  mov    %r9d,%ebp
  jmp    L136
  nopl   0x0(%rax)
L222:
  je     L209
  cmp    $0x1,%r12d
  je     L59
L89:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nop
L224:
  mov    0x10(%rsp),%r14
  mov    0x18(%rsp),%r12
  mov    %ebx,%r13d
L38:
  lea    0x3(%r13),%eax
  push   %rbx
  mov    %r12,%rdi
  mov    %r10d,%r9d
  sar    %eax
  mov    %r14,%rdx
  mov    %ebx,%ecx
  mov    %r15,%rsi
  and    $0xfffffffc,%eax
  push   %rax
  mov    %eax,%ebp
  sub    %ebp,%r13d
  mov    %r10d,0x4c(%rsp)
  mov    %r8d,0x1c(%rsp)
  call   <matprod_mat_mat_sub_xrows@plt>
  movslq %ebp,%rax
  pop    %rdi
  pop    %r8
  shl    $0x3,%rax
  mov    0x3c(%rsp),%ebp
  mov    0x20(%rsp),%r10d
  mov    0xc(%rsp),%r8d
  add    %rax,%r12
  add    %rax,%r14
  jmp    L136
  nopl   (%rax)
L236:
  push   %rbx
  mov    %ebx,%ecx
  mov    %r15,%rsi
  mov    %r10d,%r9d
  push   %rbx
  mov    0x20(%rsp),%rdx
  mov    0x28(%rsp),%rdi
  call   <matprod_mat_mat_sub_xrows@plt>
  pop    %rax
  pop    %rdx
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L7:
  je     L139
  movslq %r12d,%rsi
  mov    0x10(%rsp),%rdi
  add    $0x48,%rsp
  imul   %rbx,%rsi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    10c0 <set_to_zeros@plt>
  nopl   0x0(%rax)
L55:
  pxor   %xmm1,%xmm1
  movsd  D1(%rip),%xmm0
  cvtsi2sd %r14d,%xmm1
  addsd  D2(%rip),%xmm1
  divsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%eax
  add    $0x1,%eax
  mov    %eax,0x24(%rsp)
  jmp    L82
  nopl   (%rax)
L139:
  mov    0x10(%rsp),%rdx
  mov    0x18(%rsp),%rdi
  mov    %ebx,%r8d
  mov    %r12d,%ecx
  add    $0x48,%rsp
  mov    %r15,%rsi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1110 <matprod_outer@plt>
  xchg   %ax,%ax
L146:
  push   $0x0
  mov    %ebx,%r9d
  mov    %r14d,%r8d
  mov    %r12d,%ecx
  push   %r14
  mov    %r15,%rsi
  push   %r12
  push   %r12
  mov    0x30(%rsp),%rdx
  mov    0x38(%rsp),%rdi
  call   <matprod_mat_mat_sub_xrowscols@plt>
  add    $0x20,%rsp
  jmp    L89
L209:
  mov    0x10(%rsp),%rdx
  mov    0x18(%rsp),%rdi
  mov    %ebx,%r8d
  mov    %r14d,%ecx
  add    $0x48,%rsp
  mov    %r15,%rsi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1130 <matprod_mat_mat_n2@plt>
  nopw   0x0(%rax,%rax,1)
L188:
  jne    L89
  mov    0x10(%rsp),%rdx
  mov    0x18(%rsp),%rdi
  mov    %r14d,%r8d
  mov    %r12d,%ecx
  add    $0x48,%rsp
  mov    %r15,%rsi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    10d0 <matprod_mat_vec@plt>
L59:
  mov    0x10(%rsp),%rdx
  mov    0x18(%rsp),%rdi
  mov    %ebx,%r8d
  mov    %r14d,%ecx
  add    $0x48,%rsp
  mov    %r15,%rsi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1100 <matprod_vec_mat@plt>
L198:
  mov    %r10d,0x20(%rsp)
  mov    %ebp,%r10d
  jmp    L166
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D1
type: byte[8]
name: extract_symbolic_data
arguments: data_tag: D2
type: byte[8]

tool: {"D1": [0x00,0x00,0x00,0x00,0x80,0x31,0xf7,0x40]}

tool: {"D2": [0x00,0x00,0x00,0x00,0x00,0x00,0x8e,0x40]}

assistant: ```c
void matprod_mat_mat(double *restrict x, double *restrict y, double *restrict z, int n, int k, int m) {
  assert(((uintptr_t)(x) & (1 - 1)) == 0);
  assert(((uintptr_t)(y) & (1 - 1)) == 0);
  assert(((uintptr_t)(z) & (1 - 1)) == 0);
  x = (x);
  y = (y);
  z = (z);
  if (k <= 1) {
    if (k == 1) {
      matprod_outer(x, y, z, n, m);
    } else {
      set_to_zeros(z, (size_t)n * m);
    }
    return;
  }
  if (n <= 2) {
    if (n == 2) {
      matprod_mat_mat_n2(x, y, z, k, m);
    } else if (n == 1) {
      matprod_vec_mat(x, y, z, k, m);
    }
    return;
  }
  if (m <= 1) {
    if (m == 1) {
      matprod_mat_vec(x, y, z, n, k);
    }
    return;
  }
  if (n <= (1024 - 64) && k <= 32) {
    matprod_mat_mat_sub_xrowscols(x, y, z, n, k, m, n, n, k, 0);
    return;
  }
  int cachable = n <= (1024 - 64) ? m : 1 + (int)(95000 / ((1024 - 64) + (double)k));
  int mm = m;
  for (;;) {
    double *xx = x;
    double *zz = z;
    int xrows = n;
    int m1 = mm;
    if (m1 > cachable) {
      m1 = mm < 2 * cachable ? mm / 2 : cachable;
      m1 = (m1 + 7) & ~7;
      if (m1 > mm)
        m1 = mm;
    }
    if (xrows > (1024 - 64)) {
      while (xrows >= 2 * (1024 - 64)) {
        matprod_mat_mat_sub_xrows(xx, y, zz, n, k, m1, (1024 - 64), n);
        xx += (1024 - 64);
        zz += (1024 - 64);
        xrows -= (1024 - 64);
      }
      if (xrows > (1024 - 64)) {
        int nr = (((xrows + 3) >> 1) & ~3);
        matprod_mat_mat_sub_xrows(xx, y, zz, n, k, m1, nr, n);
        xx += nr;
        zz += nr;
        xrows -= nr;
      }
    }
    matprod_mat_mat_sub_xrows(xx, y, zz, n, k, m1, xrows, n);
    mm -= m1;
    if (mm == 0)
      break;
    y += (size_t)m1 * k;
    z += (size_t)m1 * n;
  }
}
```