user: What would the C source look like for the assembly code provided here?

```asm
<Mmv>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x48,%rsp
  mov    %edi,-0x34(%rbp)
  mov    %esi,-0x38(%rbp)
  movsd  %xmm0,-0x40(%rbp)
  mov    %rdx,-0x48(%rbp)
  mov    %ecx,-0x4c(%rbp)
  mov    %r8,-0x58(%rbp)
  movsd  %xmm1,-0x60(%rbp)
  mov    %r9,-0x68(%rbp)
  mov    -0x34(%rbp),%eax
  lea    0x3(%rax),%edx
  test   %eax,%eax
  cmovs  %edx,%eax
  sar    $0x2,%eax
  mov    %eax,%r14d
  mov    -0x34(%rbp),%eax
  cltd
  shr    $0x1e,%edx
  add    %edx,%eax
  and    $0x3,%eax
  sub    %edx,%eax
  mov    %eax,-0x34(%rbp)
  mov    -0x38(%rbp),%eax
  lea    0x3(%rax),%edx
  test   %eax,%eax
  cmovs  %edx,%eax
  sar    $0x2,%eax
  mov    %eax,%r13d
  mov    -0x38(%rbp),%eax
  cltd
  shr    $0x1e,%edx
  add    %edx,%eax
  and    $0x3,%eax
  sub    %edx,%eax
  mov    %eax,-0x38(%rbp)
  mov    $0x0,%r12d
  jmp    L48
L9:
  pxor   %xmm3,%xmm3
  movapd %xmm3,%xmm10
  movapd %xmm3,%xmm9
  movapd %xmm3,%xmm8
  lea    0x0(,%r12,4),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    (%rdx,%rax,1),%r15
  lea    0x0(,%r12,4),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    %rax,%rsi
  lea    0x0(,%r12,4),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    %rax,%rcx
  lea    0x0(,%r12,4),%eax
  cltq
  add    $0x3,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    %rax,%rdi
  mov    $0x0,%ebx
  jmp    L70
L16:
  lea    0x0(,%rbx,4),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm4
  lea    0x0(,%rbx,4),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm6
  lea    0x0(,%rbx,4),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  lea    0x0(,%rbx,4),%eax
  cltq
  add    $0x3,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm5
  movslq %ebx,%rax
  shl    $0x3,%rax
  add    %r15,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm4,%xmm0
  movapd %xmm3,%xmm1
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x1,%rax
  shl    $0x3,%rax
  add    %r15,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm6,%xmm0
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x2,%rax
  shl    $0x3,%rax
  add    %r15,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x3,%rax
  shl    $0x3,%rax
  add    %r15,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm5,%xmm0
  addsd  %xmm0,%xmm1
  movapd %xmm1,%xmm3
  movslq %ebx,%rax
  shl    $0x3,%rax
  add    %rsi,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm4,%xmm0
  movapd %xmm8,%xmm1
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x1,%rax
  shl    $0x3,%rax
  add    %rsi,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm6,%xmm0
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x2,%rax
  shl    $0x3,%rax
  add    %rsi,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x3,%rax
  shl    $0x3,%rax
  add    %rsi,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm5,%xmm0
  addsd  %xmm0,%xmm1
  movapd %xmm1,%xmm8
  movslq %ebx,%rax
  shl    $0x3,%rax
  add    %rcx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm4,%xmm0
  movapd %xmm9,%xmm1
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x1,%rax
  shl    $0x3,%rax
  add    %rcx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm6,%xmm0
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x2,%rax
  shl    $0x3,%rax
  add    %rcx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x3,%rax
  shl    $0x3,%rax
  add    %rcx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm5,%xmm0
  addsd  %xmm0,%xmm1
  movapd %xmm1,%xmm9
  movslq %ebx,%rax
  shl    $0x3,%rax
  add    %rdi,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm4,%xmm0
  movapd %xmm10,%xmm1
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x1,%rax
  shl    $0x3,%rax
  add    %rdi,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x2,%rax
  shl    $0x3,%rax
  add    %rdi,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movslq %ebx,%rax
  add    $0x3,%rax
  shl    $0x3,%rax
  add    %rdi,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm5,%xmm0
  addsd  %xmm0,%xmm1
  movapd %xmm1,%xmm10
  add    $0x1,%ebx
L70:
  cmp    %r13d,%ebx
  jl     L16
  mov    $0x0,%ebx
  jmp    L23
L71:
  lea    0x0(,%r13,4),%eax
  add    %ebx,%eax
  cltq
  shl    $0x3,%rax
  add    %r15,%rax
  movsd  (%rax),%xmm1
  lea    0x0(,%r13,4),%eax
  add    %ebx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm3,%xmm7
  addsd  %xmm7,%xmm0
  movapd %xmm0,%xmm3
  lea    0x0(,%r13,4),%eax
  add    %ebx,%eax
  cltq
  shl    $0x3,%rax
  add    %rsi,%rax
  movsd  (%rax),%xmm1
  lea    0x0(,%r13,4),%eax
  add    %ebx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm8,%xmm7
  addsd  %xmm7,%xmm0
  movapd %xmm0,%xmm8
  lea    0x0(,%r13,4),%eax
  add    %ebx,%eax
  cltq
  shl    $0x3,%rax
  add    %rcx,%rax
  movsd  (%rax),%xmm1
  lea    0x0(,%r13,4),%eax
  add    %ebx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm9,%xmm7
  addsd  %xmm7,%xmm0
  movapd %xmm0,%xmm9
  lea    0x0(,%r13,4),%eax
  add    %ebx,%eax
  cltq
  shl    $0x3,%rax
  add    %rdi,%rax
  movsd  (%rax),%xmm1
  lea    0x0(,%r13,4),%eax
  add    %ebx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm10,%xmm7
  addsd  %xmm7,%xmm0
  movapd %xmm0,%xmm10
  add    $0x1,%ebx
L23:
  cmp    -0x38(%rbp),%ebx
  jl     L71
  lea    0x0(,%r12,4),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  movapd %xmm3,%xmm0
  mulsd  -0x40(%rbp),%xmm0
  lea    0x0(,%r12,4),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  lea    0x0(,%r12,4),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  movapd %xmm8,%xmm0
  mulsd  -0x40(%rbp),%xmm0
  lea    0x0(,%r12,4),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  lea    0x0(,%r12,4),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  movapd %xmm9,%xmm0
  mulsd  -0x40(%rbp),%xmm0
  lea    0x0(,%r12,4),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  lea    0x0(,%r12,4),%eax
  cltq
  add    $0x3,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  movapd %xmm10,%xmm0
  mulsd  -0x40(%rbp),%xmm0
  lea    0x0(,%r12,4),%eax
  cltq
  add    $0x3,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  add    $0x1,%r12d
L48:
  cmp    %r14d,%r12d
  jl     L9
  mov    $0x0,%r12d
  jmp    L22
L35:
  movslq %r12d,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  lea    0x0(,%r14,4),%eax
  add    %r12d,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    (%rdx,%rax,1),%rsi
  lea    0x0(,%r13,4),%edx
  mov    -0x38(%rbp),%eax
  lea    (%rdx,%rax,1),%ecx
  mov    -0x58(%rbp),%rax
  mov    %rax,%rdx
  mov    %ecx,%edi
  mov    $0x0,%eax
  call   <Mdot@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  -0x40(%rbp),%xmm0
  lea    0x0(,%r14,4),%eax
  add    %r12d,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x70(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  add    $0x1,%r12d
L22:
  cmp    -0x34(%rbp),%r12d
  jl     L35
  nop
  nop
  add    $0x48,%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
```

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);
}
```