user: 这段汇编代码用C语言表示是什么样的?

```asm
<matmat>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %rsi,-0x50(%rbp)
  mov    %rdx,-0x58(%rbp)
  movsd  %xmm0,-0x60(%rbp)
  movsd  %xmm1,-0x68(%rbp)
  mov    %ecx,-0x6c(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movb   $0x74,-0x32(%rbp)
  movb   $0x74,-0x31(%rbp)
  mov    -0x6c(%rbp),%eax
  mov    %eax,-0x30(%rbp)
  mov    -0x6c(%rbp),%eax
  mov    %eax,-0x2c(%rbp)
  mov    -0x6c(%rbp),%eax
  mov    %eax,-0x28(%rbp)
  mov    -0x6c(%rbp),%eax
  mov    %eax,-0x24(%rbp)
  mov    -0x6c(%rbp),%eax
  mov    %eax,-0x20(%rbp)
  mov    -0x6c(%rbp),%eax
  mov    %eax,-0x1c(%rbp)
  mov    -0x6c(%rbp),%eax
  cltq
  mov    %eax,%edx
  mov    -0x6c(%rbp),%eax
  cltq
  imul   %edx,%eax
  shl    $0x3,%eax
  mov    %eax,%edi
  call   <malloc@plt>
  mov    %rax,-0x10(%rbp)
  movl   $0x0,-0x18(%rbp)
  jmp    L164
L30:
  movl   $0x0,-0x14(%rbp)
  jmp    L20
L123:
  mov    -0x6c(%rbp),%eax
  imul   -0x18(%rbp),%eax
  mov    %eax,%edx
  mov    -0x14(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rax,%rdx
  mov    -0x6c(%rbp),%eax
  imul   -0x14(%rbp),%eax
  mov    %eax,%ecx
  mov    -0x18(%rbp),%eax
  add    %ecx,%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x10(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x14(%rbp)
L20:
  mov    -0x14(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L123
  addl   $0x1,-0x18(%rbp)
L164:
  mov    -0x18(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L30
  movl   $0x0,-0x18(%rbp)
  jmp    L40
L3:
  movl   $0x0,-0x14(%rbp)
  jmp    L80
L127:
  mov    -0x6c(%rbp),%eax
  imul   -0x14(%rbp),%eax
  mov    %eax,%edx
  mov    -0x18(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rax,%rdx
  mov    -0x6c(%rbp),%eax
  imul   -0x14(%rbp),%eax
  mov    %eax,%ecx
  mov    -0x18(%rbp),%eax
  add    %ecx,%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x58(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x14(%rbp)
L80:
  mov    -0x14(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L127
  addl   $0x1,-0x18(%rbp)
L40:
  mov    -0x18(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L3
  lea    -0x60(%rbp),%r9
  lea    -0x28(%rbp),%r8
  lea    -0x2c(%rbp),%rcx
  lea    -0x30(%rbp),%rdx
  lea    -0x31(%rbp),%rsi
  lea    -0x32(%rbp),%rax
  sub    $0x8,%rsp
  lea    -0x1c(%rbp),%rdi
  push   %rdi
  push   -0x58(%rbp)
  lea    -0x68(%rbp),%rdi
  push   %rdi
  lea    -0x20(%rbp),%rdi
  push   %rdi
  push   -0x50(%rbp)
  lea    -0x24(%rbp),%rdi
  push   %rdi
  push   -0x48(%rbp)
  mov    %rax,%rdi
  call   <dgemm_@plt>
  add    $0x40,%rsp
  movl   $0x0,-0x18(%rbp)
  jmp    L169
L149:
  movl   $0x0,-0x14(%rbp)
  jmp    L119
L153:
  mov    -0x6c(%rbp),%eax
  imul   -0x18(%rbp),%eax
  mov    %eax,%edx
  mov    -0x14(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rax,%rdx
  mov    -0x6c(%rbp),%eax
  imul   -0x14(%rbp),%eax
  mov    %eax,%ecx
  mov    -0x18(%rbp),%eax
  add    %ecx,%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x10(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x14(%rbp)
L119:
  mov    -0x14(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L153
  addl   $0x1,-0x18(%rbp)
L169:
  mov    -0x18(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L149
  movl   $0x0,-0x18(%rbp)
  jmp    L37
L96:
  movl   $0x0,-0x14(%rbp)
  jmp    L160
L116:
  mov    -0x6c(%rbp),%eax
  imul   -0x14(%rbp),%eax
  mov    %eax,%edx
  mov    -0x18(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rax,%rdx
  mov    -0x6c(%rbp),%eax
  imul   -0x14(%rbp),%eax
  mov    %eax,%ecx
  mov    -0x18(%rbp),%eax
  add    %ecx,%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x58(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x14(%rbp)
L160:
  mov    -0x14(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L116
  addl   $0x1,-0x18(%rbp)
L37:
  mov    -0x18(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  jl     L96
  mov    -0x10(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L150
  call   <__stack_chk_fail@plt>
L150:
  leave
  ret
```

assistant: ```c
void matmat(double *A, double *B, double *C, double alpha, double beta, int nA) {
  int i, j;
  char TRANSA = 't';
  char TRANSB = 't';
  int m = nA;
  int n = nA;
  int k = nA;
  int LDA = nA;
  int LDB = nA;
  int LDC = nA;
  double *work = (double *)malloc(sizeof(double) * nA * nA);
  for (i = 0; i < nA; i++) {
    for (j = 0; j < nA; j++)
      work[i + nA * j] = C[j + nA * i];
  }
  for (i = 0; i < nA; i++) {
    for (j = 0; j < nA; j++)
      C[i + nA * j] = work[i + nA * j];
  }
  dgemm_(&TRANSA, &TRANSB, &m, &n, &k, &alpha, A, &LDA, B, &LDB, &beta, C, &LDC);
  for (i = 0; i < nA; i++) {
    for (j = 0; j < nA; j++)
      work[i + nA * j] = C[j + nA * i];
  }
  for (i = 0; i < nA; i++) {
    for (j = 0; j < nA; j++)
      C[i + nA * j] = work[i + nA * j];
  }
  free(work);
  return;
}
```