user: 下面的汇编指令对应的C源代码是什么:

```asm
<dmatvec>:
  endbr64
  push   %r15
  push   %r14
  mov    %r9,%r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x60,%rsp
  mov    %edi,0x48(%rsp)
  mov    %esi,-0x44(%rsp)
  mov    %edx,0x2c(%rsp)
  mov    %rcx,0x20(%rsp)
  mov    %r8,0x30(%rsp)
  cmp    $0x7,%edx
  jle    L109
  movslq %edi,%rax
  mov    %edx,%r11d
  mov    %rcx,%r15
  mov    %r8,%rbp
  lea    0x0(,%rax,8),%rbx
  lea    (%rax,%rax,4),%rdx
  mov    %rax,%rdi
  lea    0x0(,%rdx,8),%r10
  mov    %rbx,%rdx
  shl    $0x6,%rdi
  mov    %rbx,-0x58(%rsp)
  sub    %rax,%rdx
  mov    %rdi,-0x30(%rsp)
  lea    (%rcx,%rdx,8),%r12
  lea    -0x8(%r11),%edx
  mov    %r10,-0x10(%rsp)
  lea    0x8(%rcx,%rbx,1),%rcx
  shr    $0x3,%edx
  mov    %rcx,-0x40(%rsp)
  mov    %edx,%ecx
  mov    %edx,0x4c(%rsp)
  mov    %rcx,%rdx
  shl    $0x6,%rdx
  mov    %rcx,0x50(%rsp)
  lea    0x40(%r8,%rdx,1),%rcx
  lea    (%rax,%rax,2),%rdx
  mov    %rcx,-0x28(%rsp)
  mov    %rax,%rcx
  lea    0x0(,%rdx,8),%r11
  mov    %rdx,%rdi
  shl    $0x4,%rcx
  mov    %esi,%edx
  shl    $0x4,%rdi
  mov    %r11,-0x8(%rsp)
  mov    %rcx,-0x18(%rsp)
  mov    %rax,%rcx
  shr    %edx
  shl    $0x5,%rcx
  sub    $0x1,%edx
  mov    %rdi,(%rsp)
  mov    %rcx,-0x20(%rsp)
  mov    %rbx,%rcx
  add    $0x1,%rdx
  add    %r11,%rcx
  shl    $0x4,%rdx
  mov    %rcx,0x10(%rsp)
  mov    %rbx,%rcx
  add    %r10,%rbx
  mov    %rbx,0x18(%rsp)
  mov    %esi,%ebx
  add    %rdi,%rcx
  and    $0xfffffffe,%ebx
  mov    %rdx,0x8(%rsp)
  movslq %ebx,%rdx
  mov    %rcx,0x38(%rsp)
  add    %rdx,%rax
  lea    0x0(,%rdx,8),%r13
  mov    %ebx,0x28(%rsp)
  shl    $0x3,%rax
  lea    (%r9,%r13,1),%rbx
  mov    %rax,-0x50(%rsp)
  lea    -0x1(%rsi),%eax
  mov    %rax,0x40(%rsp)
  nopl   0x0(%rax,%rax,1)
L54:
  mov    -0x18(%rsp),%rax
  mov    -0x58(%rsp),%r8
  mov    -0x10(%rsp),%r9
  mov    -0x44(%rsp),%r10d
  lea    (%r15,%rax,1),%rdx
  mov    -0x8(%rsp),%rax
  lea    (%r8,%r15,1),%rcx
  movsd  0x0(%rbp),%xmm0
  movsd  0x8(%rbp),%xmm3
  movsd  0x10(%rbp),%xmm10
  lea    (%r15,%r9,1),%rsi
  lea    (%r15,%rax,1),%rdi
  mov    -0x20(%rsp),%rax
  movsd  0x18(%rbp),%xmm11
  movsd  0x20(%rbp),%xmm12
  movsd  0x28(%rbp),%xmm13
  lea    (%r15,%rax,1),%r11
  mov    (%rsp),%rax
  movsd  0x30(%rbp),%xmm14
  movsd  0x38(%rbp),%xmm15
  add    %r15,%rax
  test   %r10d,%r10d
  jle    L16
  mov    -0x40(%rsp),%r10
  add    %r9,%r10
  mov    %r10,%r9
  sub    %r8,%r9
  mov    %r9,%r8
  sub    -0x58(%rsp),%r8
  mov    %r9,-0x78(%rsp)
  mov    %r8,%r9
  sub    -0x58(%rsp),%r9
  cmpl   $0x1,-0x44(%rsp)
  mov    %r8,-0x68(%rsp)
  mov    %r9,-0x38(%rsp)
  setne  %r9b
  cmp    -0x40(%rsp),%r14
  setne  %r8b
  and    %r8d,%r9d
  lea    0x8(%r12),%r8
  cmp    %r8,%r14
  setne  %r8b
  and    %r9d,%r8d
  cmp    %r10,%r14
  mov    -0x38(%rsp),%r10
  setne  %r9b
  and    %r8d,%r9d
  cmp    -0x78(%rsp),%r14
  setne  %r8b
  and    %r9d,%r8d
  cmp    -0x68(%rsp),%r14
  setne  %r9b
  and    %r9d,%r8d
  cmp    %r10,%r14
  setne  %r9b
  sub    -0x58(%rsp),%r10
  and    %r8d,%r9d
  cmp    %r10,%r14
  setne  %r8b
  test   %r8b,%r9b
  je     L277
  lea    0x8(%r15),%r8
  cmp    %r8,%r14
  je     L277
  movapd %xmm14,%xmm1
  mov    0x10(%rsp),%r10
  mov    0x18(%rsp),%r9
  mov    %rax,-0x38(%rsp)
  unpcklpd %xmm1,%xmm1
  movapd %xmm0,%xmm9
  movapd %xmm3,%xmm8
  mov    0x8(%rsp),%rax
  movaps %xmm1,-0x78(%rsp)
  movapd %xmm15,%xmm1
  movapd %xmm10,%xmm7
  movapd %xmm11,%xmm6
  movapd %xmm12,%xmm5
  unpcklpd %xmm1,%xmm1
  unpcklpd %xmm7,%xmm7
  add    %r15,%r10
  movapd %xmm13,%xmm4
  unpcklpd %xmm9,%xmm9
  add    %r15,%r9
  xor    %r8d,%r8d
  unpcklpd %xmm8,%xmm8
  unpcklpd %xmm6,%xmm6
  unpcklpd %xmm5,%xmm5
  movaps %xmm1,-0x68(%rsp)
  unpcklpd %xmm4,%xmm4
  nopl   (%rax)
L134:
  movupd (%r15,%r8,1),%xmm1
  movupd (%rcx,%r8,1),%xmm2
  mulpd  %xmm9,%xmm1
  mulpd  %xmm8,%xmm2
  addpd  %xmm1,%xmm2
  movupd (%rdx,%r8,1),%xmm1
  mulpd  %xmm7,%xmm1
  addpd  %xmm2,%xmm1
  movupd (%rdi,%r8,1),%xmm2
  mulpd  %xmm6,%xmm2
  addpd  %xmm1,%xmm2
  movupd (%r10,%r8,1),%xmm1
  mulpd  %xmm5,%xmm1
  addpd  %xmm2,%xmm1
  movupd (%rsi,%r8,1),%xmm2
  mulpd  %xmm4,%xmm2
  addpd  %xmm1,%xmm2
  movupd (%r9,%r8,1),%xmm1
  mulpd  -0x78(%rsp),%xmm1
  addpd  %xmm2,%xmm1
  movupd (%r12,%r8,1),%xmm2
  mulpd  -0x68(%rsp),%xmm2
  addpd  %xmm2,%xmm1
  movupd (%r14,%r8,1),%xmm2
  addpd  %xmm2,%xmm1
  movups %xmm1,(%r14,%r8,1)
  add    $0x10,%r8
  cmp    %rax,%r8
  jne    L134
  mov    -0x38(%rsp),%rax
  mov    0x28(%rsp),%r9d
  cmp    %r9d,-0x44(%rsp)
  je     L16
  mulsd  (%r15,%r13,1),%xmm0
  mulsd  (%rcx,%r13,1),%xmm3
  mulsd  (%rdx,%r13,1),%xmm10
  mulsd  (%rdi,%r13,1),%xmm11
  mov    -0x50(%rsp),%rdi
  mulsd  (%r11,%r13,1),%xmm12
  mulsd  (%rsi,%r13,1),%xmm13
  addsd  %xmm3,%xmm0
  mulsd  (%rax,%r13,1),%xmm14
  mulsd  (%rax,%rdi,1),%xmm15
  addsd  %xmm10,%xmm0
  addsd  %xmm0,%xmm11
  addsd  %xmm11,%xmm12
  addsd  %xmm12,%xmm13
  addsd  %xmm13,%xmm14
  addsd  %xmm15,%xmm14
  addsd  (%rbx),%xmm14
  movsd  %xmm14,(%rbx)
L16:
  mov    -0x30(%rsp),%rax
  add    $0x40,%rbp
  add    %rax,-0x40(%rsp)
  add    %rax,%r15
  add    %rax,%r12
  cmp    -0x28(%rsp),%rbp
  jne    L54
  mov    0x4c(%rsp),%eax
  lea    0x8(,%rax,8),%eax
  mov    %eax,-0x50(%rsp)
  mov    0x50(%rsp),%rax
  add    $0x1,%rax
  imul   -0x30(%rsp),%rax
  add    %rax,0x20(%rsp)
L193:
  mov    0x2c(%rsp),%ecx
  mov    -0x50(%rsp),%ebx
  lea    -0x3(%rcx),%eax
  cmp    %ebx,%eax
  jle    L216
  sub    $0x4,%ecx
  mov    0x30(%rsp),%rdi
  movslq 0x48(%rsp),%rax
  movslq %ebx,%rdx
  sub    %ebx,%ecx
  mov    -0x44(%rsp),%r13d
  mov    0x20(%rsp),%r12
  shr    $0x2,%ecx
  lea    (%rdi,%rdx,8),%r10
  mov    %rax,%rsi
  mov    %ecx,%ebx
  mov    %ecx,-0x28(%rsp)
  mov    %r13d,%ecx
  lea    0x0(,%rax,8),%r15
  lea    (%rdx,%rbx,4),%rdx
  mov    %rbx,-0x20(%rsp)
  and    $0xfffffffe,%ecx
  shl    $0x4,%rsi
  lea    0x20(%rdi,%rdx,8),%rbx
  mov    %r13d,%edx
  lea    -0x1(%r13),%edi
  mov    %ecx,-0x68(%rsp)
  shr    %edx
  mov    %rbx,-0x78(%rsp)
  lea    0x8(%r15),%rbx
  add    %r12,%rsi
  lea    0x0(,%rax,4),%ebp
  mov    %rbx,-0x58(%rsp)
  lea    -0x1(%rdx),%ebx
  movslq %ecx,%rdx
  mov    %rdi,-0x30(%rsp)
  add    %rdx,%rax
  movslq %ebp,%rbp
  add    $0x1,%rbx
  lea    0x0(,%rdx,8),%rcx
  shl    $0x3,%rbp
  shl    $0x4,%rbx
  mov    %rsi,%r8
  lea    0x0(,%rax,8),%r11
  lea    (%r14,%rcx,1),%r9
  mov    %r12,%rax
  nopl   (%rax)
L138:
  movsd  (%r10),%xmm3
  movsd  0x8(%r10),%xmm1
  lea    (%r15,%rax,1),%rdi
  movsd  0x10(%r10),%xmm4
  movsd  0x18(%r10),%xmm5
  test   %r13d,%r13d
  jle    L132
  lea    0x8(%r8),%rdx
  cmp    %rdx,%r14
  setne  %dl
  cmp    $0x1,%r13d
  setne  %r12b
  and    %r12d,%edx
  mov    -0x58(%rsp),%r12
  add    %r8,%r12
  cmp    %r12,%r14
  setne  %r12b
  and    %r12d,%edx
  lea    0x8(%rax),%r12
  cmp    %r12,%r14
  setne  %r12b
  test   %r12b,%dl
  je     L42
  mov    -0x58(%rsp),%r12
  lea    (%r12,%rax,1),%rdx
  cmp    %rdx,%r14
  je     L42
  movapd %xmm3,%xmm9
  movapd %xmm1,%xmm8
  movapd %xmm4,%xmm7
  xor    %edx,%edx
  movapd %xmm5,%xmm6
  unpcklpd %xmm9,%xmm9
  unpcklpd %xmm8,%xmm8
  unpcklpd %xmm7,%xmm7
  lea    (%r15,%r8,1),%r12
  unpcklpd %xmm6,%xmm6
  nopl   0x0(%rax,%rax,1)
L85:
  movupd (%rdi,%rdx,1),%xmm0
  movupd (%rax,%rdx,1),%xmm2
  mulpd  %xmm8,%xmm0
  mulpd  %xmm9,%xmm2
  addpd  %xmm2,%xmm0
  movupd (%r8,%rdx,1),%xmm2
  mulpd  %xmm7,%xmm2
  addpd  %xmm2,%xmm0
  movupd (%r12,%rdx,1),%xmm2
  mulpd  %xmm6,%xmm2
  addpd  %xmm2,%xmm0
  movupd (%r14,%rdx,1),%xmm2
  addpd  %xmm2,%xmm0
  movups %xmm0,(%r14,%rdx,1)
  add    $0x10,%rdx
  cmp    %rdx,%rbx
  jne    L85
  cmp    %r13d,-0x68(%rsp)
  je     L132
  mulsd  (%rdi,%rcx,1),%xmm1
  mulsd  (%rax,%rcx,1),%xmm3
  mulsd  (%rsi,%rcx,1),%xmm4
  mulsd  (%rsi,%r11,1),%xmm5
  movapd %xmm1,%xmm0
  addsd  %xmm3,%xmm0
  addsd  %xmm4,%xmm0
  addsd  %xmm5,%xmm0
  addsd  (%r9),%xmm0
  movsd  %xmm0,(%r9)
L132:
  add    %rbp,%rax
  add    $0x20,%r10
  add    %rbp,%rsi
  add    %rbp,%r8
  cmp    -0x78(%rsp),%r10
  jne    L138
  mov    -0x50(%rsp),%eax
  mov    -0x28(%rsp),%ebx
  lea    0x4(%rax,%rbx,4),%eax
  mov    %eax,-0x50(%rsp)
  mov    -0x20(%rsp),%rax
  add    $0x1,%rax
  imul   %rbp,%rax
  add    %rax,0x20(%rsp)
L216:
  mov    0x2c(%rsp),%ebp
  mov    -0x50(%rsp),%eax
  cmp    %eax,%ebp
  jle    L128
  movslq 0x48(%rsp),%r8
  mov    -0x44(%rsp),%ebx
  shl    $0x3,%r8
  test   %ebx,%ebx
  jle    L128
  movslq %eax,%rdi
  mov    %ebx,%r10d
  mov    %ebx,%eax
  mov    0x20(%rsp),%r11
  shr    %eax
  and    $0xfffffffe,%r10d
  mov    0x30(%rsp),%r12
  movslq %ebx,%r9
  lea    -0x1(%rax),%ecx
  mov    %r10d,%esi
  add    $0x1,%rcx
  shl    $0x3,%rsi
  shl    $0x4,%rcx
  lea    (%r14,%rsi,1),%rdx
L165:
  lea    0x8(%r11),%rax
  movsd  (%r12,%rdi,8),%xmm1
  cmp    %rax,%r14
  je     L49
  cmp    $0x1,%ebx
  je     L49
  movapd %xmm1,%xmm2
  xor    %eax,%eax
  unpcklpd %xmm2,%xmm2
  nopl   0x0(%rax,%rax,1)
L151:
  movupd (%r11,%rax,1),%xmm0
  movupd (%r14,%rax,1),%xmm6
  mulpd  %xmm2,%xmm0
  addpd  %xmm6,%xmm0
  movups %xmm0,(%r14,%rax,1)
  add    $0x10,%rax
  cmp    %rax,%rcx
  jne    L151
  cmp    %ebx,%r10d
  je     L262
  movsd  (%r11,%rsi,1),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  (%rdx),%xmm0
  movsd  %xmm0,(%rdx)
L262:
  add    $0x1,%rdi
  add    %r8,%r11
  cmp    %edi,%ebp
  jg     L165
L128:
  add    $0x60,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L277:
  mov    0x38(%rsp),%r10
  mov    %rbx,-0x78(%rsp)
  xor    %r8d,%r8d
  mov    0x40(%rsp),%rbx
  add    %r15,%r10
  nopl   (%rax)
L129:
  movsd  (%r15,%r8,8),%xmm1
  movsd  (%rcx,%r8,8),%xmm2
  mov    %r8,%r9
  mulsd  %xmm3,%xmm2
  mulsd  %xmm0,%xmm1
  addsd  %xmm2,%xmm1
  movsd  (%rdx,%r8,8),%xmm2
  mulsd  %xmm10,%xmm2
  addsd  %xmm2,%xmm1
  movsd  (%rdi,%r8,8),%xmm2
  mulsd  %xmm11,%xmm2
  addsd  %xmm2,%xmm1
  movsd  (%r11,%r8,8),%xmm2
  mulsd  %xmm12,%xmm2
  addsd  %xmm2,%xmm1
  movsd  (%rsi,%r8,8),%xmm2
  mulsd  %xmm13,%xmm2
  addsd  %xmm2,%xmm1
  movsd  (%rax,%r8,8),%xmm2
  mulsd  %xmm14,%xmm2
  addsd  %xmm2,%xmm1
  movsd  (%r10,%r8,8),%xmm2
  mulsd  %xmm15,%xmm2
  addsd  %xmm2,%xmm1
  addsd  (%r14,%r8,8),%xmm1
  movsd  %xmm1,(%r14,%r8,8)
  add    $0x1,%r8
  cmp    %rbx,%r9
  jne    L129
  mov    -0x78(%rsp),%rbx
  jmp    L16
  nopl   0x0(%rax,%rax,1)
L42:
  mov    %rcx,-0x38(%rsp)
  lea    (%r15,%rsi,1),%r12
  mov    -0x30(%rsp),%rcx
  xor    %edx,%edx
  mov    %r15,-0x40(%rsp)
  mov    %r12,%r15
  nopl   0x0(%rax,%rax,1)
L32:
  movsd  (%rdi,%rdx,8),%xmm0
  movsd  (%rax,%rdx,8),%xmm2
  mov    %rdx,%r12
  mulsd  %xmm3,%xmm2
  mulsd  %xmm1,%xmm0
  addsd  %xmm2,%xmm0
  movsd  (%rsi,%rdx,8),%xmm2
  mulsd  %xmm4,%xmm2
  addsd  %xmm2,%xmm0
  movsd  (%r15,%rdx,8),%xmm2
  mulsd  %xmm5,%xmm2
  addsd  %xmm2,%xmm0
  addsd  (%r14,%rdx,8),%xmm0
  movsd  %xmm0,(%r14,%rdx,8)
  add    $0x1,%rdx
  cmp    %r12,%rcx
  jne    L32
  mov    -0x40(%rsp),%r15
  mov    -0x38(%rsp),%rcx
  jmp    L132
L49:
  xor    %eax,%eax
  nopw   0x0(%rax,%rax,1)
L161:
  movsd  (%r11,%rax,8),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  (%r14,%rax,8),%xmm0
  movsd  %xmm0,(%r14,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%r9
  jne    L161
  jmp    L262
L109:
  movl   $0x0,-0x50(%rsp)
  jmp    L193
```

assistant: ```c
void dmatvec(ldm, nrow, ncol, M, vec, Mxvec) int ldm;
int nrow;
int ncol;
double *M;
double *vec;
double *Mxvec;
{
  register double vi0, vi1, vi2, vi3, vi4, vi5, vi6, vi7;
  register double *Mki0, *Mki1, *Mki2, *Mki3, *Mki4, *Mki5, *Mki6, *Mki7;
  register int firstcol = 0;
  double *M0;
  int k;
  M0 = &M[0];
  while (firstcol < ncol - 7) {
    Mki0 = M0;
    Mki1 = Mki0 + ldm;
    Mki2 = Mki1 + ldm;
    Mki3 = Mki2 + ldm;
    Mki4 = Mki3 + ldm;
    Mki5 = Mki4 + ldm;
    Mki6 = Mki5 + ldm;
    Mki7 = Mki6 + ldm;
    vi0 = vec[firstcol++];
    vi1 = vec[firstcol++];
    vi2 = vec[firstcol++];
    vi3 = vec[firstcol++];
    vi4 = vec[firstcol++];
    vi5 = vec[firstcol++];
    vi6 = vec[firstcol++];
    vi7 = vec[firstcol++];
    for (k = 0; k < nrow; k++)
      Mxvec[k] += vi0 * Mki0[k] + vi1 * Mki1[k] + vi2 * Mki2[k] + vi3 * Mki3[k] + vi4 * Mki4[k] + vi5 * Mki5[k] + vi6 * Mki6[k] + vi7 * Mki7[k];
    M0 += 8 * ldm;
  }
  while (firstcol < ncol - 3) {
    Mki0 = M0;
    Mki1 = Mki0 + ldm;
    Mki2 = Mki1 + ldm;
    Mki3 = Mki2 + ldm;
    vi0 = vec[firstcol++];
    vi1 = vec[firstcol++];
    vi2 = vec[firstcol++];
    vi3 = vec[firstcol++];
    for (k = 0; k < nrow; k++)
      Mxvec[k] += vi0 * Mki0[k] + vi1 * Mki1[k] + vi2 * Mki2[k] + vi3 * Mki3[k];
    M0 += 4 * ldm;
  }
  while (firstcol < ncol) {
    Mki0 = M0;
    vi0 = vec[firstcol++];
    for (k = 0; k < nrow; k++)
      Mxvec[k] += vi0 * Mki0[k];
    M0 += ldm;
  }
}
```