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

```asm
<kernel_2mm>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %edi,-0x14(%rbp)
  mov    %esi,-0x18(%rbp)
  mov    %edx,-0x1c(%rbp)
  mov    %ecx,-0x20(%rbp)
  movsd  %xmm0,-0x28(%rbp)
  movsd  %xmm1,-0x30(%rbp)
  mov    %r8,-0x38(%rbp)
  mov    %r9,-0x40(%rbp)
  movl   $0x0,-0xc(%rbp)
  jmp    L6
L95:
  movl   $0x0,-0x8(%rbp)
  jmp    L13
L82:
  mov    -0xc(%rbp),%eax
  cltq
  imul   $0x1c20,%rax,%rdx
  mov    -0x38(%rbp),%rax
  add    %rax,%rdx
  mov    -0x8(%rbp),%eax
  cltq
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  movl   $0x0,-0x4(%rbp)
  jmp    L45
L80:
  mov    -0xc(%rbp),%eax
  cltq
  imul   $0x1c20,%rax,%rdx
  mov    -0x38(%rbp),%rax
  add    %rax,%rdx
  mov    -0x8(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm1
  mov    -0xc(%rbp),%eax
  cltq
  imul   $0x2260,%rax,%rdx
  mov    -0x40(%rbp),%rax
  add    %rax,%rdx
  mov    -0x4(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0x28(%rbp),%xmm2
  mov    -0x4(%rbp),%eax
  cltq
  imul   $0x1c20,%rax,%rdx
  mov    0x10(%rbp),%rax
  add    %rax,%rdx
  mov    -0x8(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm0
  mulsd  %xmm2,%xmm0
  mov    -0xc(%rbp),%eax
  cltq
  imul   $0x1c20,%rax,%rdx
  mov    -0x38(%rbp),%rax
  add    %rax,%rdx
  addsd  %xmm1,%xmm0
  mov    -0x8(%rbp),%eax
  cltq
  movsd  %xmm0,(%rdx,%rax,8)
  addl   $0x1,-0x4(%rbp)
L45:
  mov    -0x4(%rbp),%eax
  cmp    -0x1c(%rbp),%eax
  jl     L80
  addl   $0x1,-0x8(%rbp)
L13:
  mov    -0x8(%rbp),%eax
  cmp    -0x18(%rbp),%eax
  jl     L82
  addl   $0x1,-0xc(%rbp)
L6:
  mov    -0xc(%rbp),%eax
  cmp    -0x14(%rbp),%eax
  jl     L95
  movl   $0x0,-0xc(%rbp)
  jmp    L114
L55:
  movl   $0x0,-0x8(%rbp)
  jmp    L79
L16:
  mov    -0xc(%rbp),%eax
  cltq
  imul   $0x2580,%rax,%rdx
  mov    0x20(%rbp),%rax
  add    %rax,%rdx
  mov    -0x8(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm0
  mov    -0xc(%rbp),%eax
  cltq
  imul   $0x2580,%rax,%rdx
  mov    0x20(%rbp),%rax
  add    %rax,%rdx
  mulsd  -0x30(%rbp),%xmm0
  mov    -0x8(%rbp),%eax
  cltq
  movsd  %xmm0,(%rdx,%rax,8)
  movl   $0x0,-0x4(%rbp)
  jmp    L36
L44:
  mov    -0xc(%rbp),%eax
  cltq
  imul   $0x2580,%rax,%rdx
  mov    0x20(%rbp),%rax
  add    %rax,%rdx
  mov    -0x8(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm1
  mov    -0xc(%rbp),%eax
  cltq
  imul   $0x1c20,%rax,%rdx
  mov    -0x38(%rbp),%rax
  add    %rax,%rdx
  mov    -0x4(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm2
  mov    -0x4(%rbp),%eax
  cltq
  imul   $0x2580,%rax,%rdx
  mov    0x18(%rbp),%rax
  add    %rax,%rdx
  mov    -0x8(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm0
  mulsd  %xmm2,%xmm0
  mov    -0xc(%rbp),%eax
  cltq
  imul   $0x2580,%rax,%rdx
  mov    0x20(%rbp),%rax
  add    %rax,%rdx
  addsd  %xmm1,%xmm0
  mov    -0x8(%rbp),%eax
  cltq
  movsd  %xmm0,(%rdx,%rax,8)
  addl   $0x1,-0x4(%rbp)
L36:
  mov    -0x4(%rbp),%eax
  cmp    -0x18(%rbp),%eax
  jl     L44
  addl   $0x1,-0x8(%rbp)
L79:
  mov    -0x8(%rbp),%eax
  cmp    -0x20(%rbp),%eax
  jl     L16
  addl   $0x1,-0xc(%rbp)
L114:
  mov    -0xc(%rbp),%eax
  cmp    -0x14(%rbp),%eax
  jl     L55
  nop
  nop
  pop    %rbp
  ret
```

assistant: ```c
void kernel_2mm(int ni, int nj, int nk, int nl, double alpha, double beta, double tmp[800 + 0][900 + 0], double A[800 + 0][1100 + 0], double B[1100 + 0][900 + 0], double C[900 + 0][1200 + 0], double D[800 + 0][1200 + 0]) {
  int i, j, k;
#pragma scop
  for (i = 0; i < ni; i++)
    for (j = 0; j < nj; j++) {
      tmp[i][j] = 0.0;
      for (k = 0; k < nk; ++k)
        tmp[i][j] += alpha * A[i][k] * B[k][j];
    }
  for (i = 0; i < ni; i++)
    for (j = 0; j < nl; j++) {
      D[i][j] *= beta;
      for (k = 0; k < nj; ++k)
        D[i][j] += tmp[i][k] * C[k][j];
    }
#pragma endscop
}
```