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

```asm
<Mmv>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %edi,%eax
  mov    %esi,%r15d
  movsd  %xmm0,0x8(%rsp)
  mov    %rdx,%rbp
  mov    %rdx,0x28(%rsp)
  mov    %ecx,%edi
  mov    %ecx,0x34(%rsp)
  mov    %r8,%rbx
  movsd  %xmm1,0x18(%rsp)
  mov    %r9,0x38(%rsp)
  lea    0x3(%rax),%edx
  test   %eax,%eax
  cmovns %eax,%edx
  sar    $0x2,%edx
  mov    %edx,0x10(%rsp)
  cltd
  shr    $0x1e,%edx
  lea    (%rax,%rdx,1),%ecx
  and    $0x3,%ecx
  sub    %edx,%ecx
  mov    %ecx,0x30(%rsp)
  lea    0x3(%rsi),%r10d
  test   %esi,%esi
  cmovns %esi,%r10d
  sar    $0x2,%r10d
  mov    %esi,%edx
  sar    $0x1f,%edx
  shr    $0x1e,%edx
  lea    (%rsi,%rdx,1),%ecx
  and    $0x3,%ecx
  mov    %ecx,%r14d
  sub    %edx,%r14d
  cmp    $0x3,%eax
  jle    L71
  movslq %edi,%r12
  shl    $0x3,%r12
  lea    0x0(,%r10,4),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    %rdx,0x20(%rsp)
  lea    -0x1(%r14),%edx
  lea    0x1(%rax,%rdx,1),%r11
  shl    $0x3,%r11
  mov    $0x0,%r13d
  movapd %xmm1,%xmm9
L0:
  mov    %r12,%rcx
  add    0x0(%rbp),%rcx
  mov    %r12,%rsi
  add    0x8(%rbp),%rsi
  mov    %r12,%rdi
  add    0x10(%rbp),%rdi
  mov    %r12,%r8
  add    0x18(%rbp),%r8
  cmp    $0x3,%r15d
  jle    L78
  mov    %rbx,%rdx
  mov    $0x0,%eax
  pxor   %xmm1,%xmm1
  movapd %xmm1,%xmm2
  movapd %xmm1,%xmm4
  movapd %xmm1,%xmm5
L79:
  movsd  (%rdx),%xmm0
  movsd  0x8(%rdx),%xmm8
  movsd  0x10(%rdx),%xmm7
  movsd  0x18(%rdx),%xmm6
  movapd %xmm0,%xmm3
  mulsd  (%rcx,%rax,8),%xmm3
  addsd  %xmm3,%xmm5
  movapd %xmm8,%xmm3
  mulsd  0x8(%rcx,%rax,8),%xmm3
  addsd  %xmm3,%xmm5
  movapd %xmm7,%xmm3
  mulsd  0x10(%rcx,%rax,8),%xmm3
  addsd  %xmm3,%xmm5
  movapd %xmm6,%xmm3
  mulsd  0x18(%rcx,%rax,8),%xmm3
  addsd  %xmm3,%xmm5
  movapd %xmm0,%xmm3
  mulsd  (%rsi,%rax,8),%xmm3
  addsd  %xmm3,%xmm4
  movapd %xmm8,%xmm3
  mulsd  0x8(%rsi,%rax,8),%xmm3
  addsd  %xmm3,%xmm4
  movapd %xmm7,%xmm3
  mulsd  0x10(%rsi,%rax,8),%xmm3
  addsd  %xmm3,%xmm4
  movapd %xmm6,%xmm3
  mulsd  0x18(%rsi,%rax,8),%xmm3
  addsd  %xmm3,%xmm4
  movapd %xmm0,%xmm3
  mulsd  (%rdi,%rax,8),%xmm3
  addsd  %xmm3,%xmm2
  mulsd  0x8(%rdi,%rax,8),%xmm8
  addsd  %xmm8,%xmm2
  movapd %xmm7,%xmm3
  mulsd  0x10(%rdi,%rax,8),%xmm3
  addsd  %xmm3,%xmm2
  movapd %xmm6,%xmm3
  mulsd  0x18(%rdi,%rax,8),%xmm3
  addsd  %xmm3,%xmm2
  mulsd  (%r8,%rax,8),%xmm0
  addsd  %xmm0,%xmm1
  movapd %xmm7,%xmm0
  mulsd  0x8(%r8,%rax,8),%xmm0
  addsd  %xmm0,%xmm1
  mulsd  0x10(%r8,%rax,8),%xmm7
  addsd  %xmm7,%xmm1
  mulsd  0x18(%r8,%rax,8),%xmm6
  addsd  %xmm6,%xmm1
  add    $0x1,%rax
  add    $0x20,%rdx
  cmp    %eax,%r10d
  jg     L79
L70:
  test   %r14d,%r14d
  jle    L72
  mov    0x20(%rsp),%rax
L63:
  movsd  (%rbx,%rax,1),%xmm0
  movapd %xmm0,%xmm3
  mulsd  (%rcx,%rax,1),%xmm3
  addsd  %xmm3,%xmm5
  movapd %xmm0,%xmm3
  mulsd  (%rsi,%rax,1),%xmm3
  addsd  %xmm3,%xmm4
  movapd %xmm0,%xmm3
  mulsd  (%rdi,%rax,1),%xmm3
  addsd  %xmm3,%xmm2
  mulsd  (%r8,%rax,1),%xmm0
  addsd  %xmm0,%xmm1
  add    $0x8,%rax
  cmp    %r11,%rax
  jne    L63
L72:
  movapd %xmm9,%xmm0
  mulsd  (%r9),%xmm0
  movsd  0x8(%rsp),%xmm6
  mulsd  %xmm6,%xmm5
  addsd  %xmm5,%xmm0
  movsd  %xmm0,(%r9)
  movapd %xmm9,%xmm0
  mulsd  0x8(%r9),%xmm0
  mulsd  %xmm6,%xmm4
  addsd  %xmm4,%xmm0
  movsd  %xmm0,0x8(%r9)
  movapd %xmm9,%xmm0
  mulsd  0x10(%r9),%xmm0
  mulsd  %xmm6,%xmm2
  addsd  %xmm2,%xmm0
  movsd  %xmm0,0x10(%r9)
  movapd %xmm9,%xmm0
  mulsd  0x18(%r9),%xmm0
  mulsd  %xmm6,%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x18(%r9)
  add    $0x1,%r13d
  add    $0x20,%rbp
  add    $0x20,%r9
  cmp    %r13d,0x10(%rsp)
  jg     L0
L71:
  mov    0x30(%rsp),%eax
  test   %eax,%eax
  jle    L15
  movslq 0x34(%rsp),%r14
  shl    $0x3,%r14
  mov    0x38(%rsp),%rdx
  mov    %rdx,%rbp
  mov    0x10(%rsp),%r13d
  shl    $0x2,%r13d
  movslq %r13d,%r13
  mov    0x28(%rsp),%rcx
  lea    (%rcx,%r13,8),%r12
  lea    -0x1(%rax),%eax
  lea    0x8(%rdx,%rax,8),%rax
  mov    %rax,0x20(%rsp)
L37:
  movsd  0x18(%rsp),%xmm7
  mulsd  0x0(%rbp),%xmm7
  movsd  %xmm7,0x10(%rsp)
  mov    %r14,%rsi
  add    (%r12),%rsi
  mov    %rbx,%rdx
  mov    %r15d,%edi
  mov    $0x0,%eax
  call   <Mdot@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  0x8(%rsp),%xmm0
  addsd  0x10(%rsp),%xmm0
  movsd  %xmm0,0x0(%rbp,%r13,8)
  add    $0x8,%rbp
  add    $0x8,%r12
  cmp    0x20(%rsp),%rbp
  jne    L37
L15:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L78:
  pxor   %xmm1,%xmm1
  movapd %xmm1,%xmm2
  movapd %xmm1,%xmm4
  movapd %xmm1,%xmm5
  jmp    L70
```

assistant: ```c
void Mmv(m, n, alpha, A, j0, x, beta, y) int m, n, j0;
double alpha, beta;
double **A, *x, *y;
{
  register int i, j, m4, n4;
  register double sum0, sum1, sum2, sum3, tmp0, tmp1, tmp2, tmp3;
  register double *dp0, *dp1, *dp2, *dp3;
  m4 = n4 = 0;
  m4 = m / 4;
  m = m % 4;
  n4 = n / 4;
  n = n % 4;
  for (i = 0; i < m4; i++) {
    sum0 = sum1 = sum2 = sum3 = 0.0;
    dp0 = &(A[4 * i][j0]);
    dp1 = &(A[4 * i + 1][j0]);
    dp2 = &(A[4 * i + 2][j0]);
    dp3 = &(A[4 * i + 3][j0]);
    for (j = 0; j < n4; j++) {
      tmp0 = x[4 * j];
      tmp1 = x[4 * j + 1];
      tmp2 = x[4 * j + 2];
      tmp3 = x[4 * j + 3];
      sum0 = sum0 + dp0[j] * tmp0 + dp0[j + 1] * tmp1 + dp0[j + 2] * tmp2 + dp0[j + 3] * tmp3;
      sum1 = sum1 + dp1[j] * tmp0 + dp1[j + 1] * tmp1 + dp1[j + 2] * tmp2 + dp1[j + 3] * tmp3;
      sum2 = sum2 + dp2[j] * tmp0 + dp2[j + 1] * tmp1 + dp2[j + 2] * tmp2 + dp2[j + 3] * tmp3;
      sum3 = sum3 + dp3[j] * tmp0 + dp3[j + 1] * tmp2 + dp3[j + 2] * tmp2 + dp3[j + 3] * tmp3;
    }
    for (j = 0; j < n; j++) {
      sum0 += dp0[4 * n4 + j] * x[4 * n4 + j];
      sum1 += dp1[4 * n4 + j] * x[4 * n4 + j];
      sum2 += dp2[4 * n4 + j] * x[4 * n4 + j];
      sum3 += dp3[4 * n4 + j] * x[4 * n4 + j];
    }
    y[4 * i] = beta * y[4 * i] + alpha * sum0;
    y[4 * i + 1] = beta * y[4 * i + 1] + alpha * sum1;
    y[4 * i + 2] = beta * y[4 * i + 2] + alpha * sum2;
    y[4 * i + 3] = beta * y[4 * i + 3] + alpha * sum3;
  }
  for (i = 0; i < m; i++)
    y[4 * m4 + i] = beta * y[i] + alpha * Mdot(4 * n4 + n, &(A[4 * m4 + i][j0]), x);
}
```