user: Please show me the C implementation that matches the assembly code listed below:

```asm
<_inner_mv_daxpy2>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  mov    %rdi,-0x10(%rbp)
  mov    %rsi,-0x18(%rbp)
  mov    %rdx,-0x20(%rbp)
  mov    %ecx,-0x24(%rbp)
  mov    %r8,-0x30(%rbp)
  mov    %r9,-0x38(%rbp)
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm6
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm7
  mov    $0x0,%ebx
  jmp    L29
L18:
  mov    -0x10(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  movapd %xmm6,%xmm2
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x10(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  movapd %xmm7,%xmm3
  mulsd  %xmm3,%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x18(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x20(%rbp)
  mov    -0x10(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  mov    -0x10(%rbp),%rax
  add    $0x8,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm3,%xmm0
  mov    -0x18(%rbp),%rax
  add    $0x8,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x20(%rbp)
  mov    -0x10(%rbp),%rax
  add    $0x10,%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  mov    -0x10(%rbp),%rax
  add    $0x10,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  add    $0x10,%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm3,%xmm0
  mov    -0x18(%rbp),%rax
  add    $0x10,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x20(%rbp)
  mov    -0x10(%rbp),%rax
  add    $0x18,%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  mov    -0x10(%rbp),%rax
  add    $0x18,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  add    $0x18,%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm3,%xmm0
  mov    -0x18(%rbp),%rax
  add    $0x18,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x20(%rbp)
  addq   $0x20,-0x10(%rbp)
  addq   $0x20,-0x18(%rbp)
  add    $0x4,%ebx
L29:
  mov    0x10(%rbp),%eax
  sub    $0x3,%eax
  cmp    %eax,%ebx
  jl     L18
  cmp    0x10(%rbp),%ebx
  je     L24
  mov    0x10(%rbp),%eax
  sub    $0x1,%eax
  cmp    %eax,%ebx
  jge    L26
  mov    -0x10(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  movapd %xmm6,%xmm4
  mulsd  %xmm4,%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x10(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  movapd %xmm7,%xmm5
  mulsd  %xmm5,%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x18(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x20(%rbp)
  mov    -0x10(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm4,%xmm0
  mov    -0x10(%rbp),%rax
  add    $0x8,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm5,%xmm0
  mov    -0x18(%rbp),%rax
  add    $0x8,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x20(%rbp)
  addq   $0x10,-0x10(%rbp)
  addq   $0x10,-0x18(%rbp)
  add    $0x2,%ebx
L26:
  cmp    0x10(%rbp),%ebx
  jge    L36
  mov    -0x10(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm6,%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x10(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm7,%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x18(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x20(%rbp)
  add    $0x1,%ebx
  jmp    L36
L24:
  nop
L36:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

assistant: ```c
void _inner_mv_daxpy2(double *a0, double *a1, const double *x, int incX, const double *y0, const double *y1, double alpha, int nRE) {
  register int i;
  register double cf0, cf1;
  cf0 = alpha * y0[0];
  cf1 = alpha * y1[0];
  for (i = 0; i < nRE - 3; i += 4) {
    a0[0] += x[0] * cf0;
    a1[0] += x[0] * cf1;
    x += incX;
    a0[1] += x[0] * cf0;
    a1[1] += x[0] * cf1;
    x += incX;
    a0[2] += x[0] * cf0;
    a1[2] += x[0] * cf1;
    x += incX;
    a0[3] += x[0] * cf0;
    a1[3] += x[0] * cf1;
    x += incX;
    a0 += 4;
    a1 += 4;
  }
  if (i == nRE)
    return;
  if (i < nRE - 1) {
    a0[0] += x[0] * cf0;
    a1[0] += x[0] * cf1;
    x += incX;
    a0[1] += x[0] * cf0;
    a1[1] += x[0] * cf1;
    x += incX;
    a0 += 2;
    a1 += 2;
    i += 2;
  }
  if (i < nRE) {
    a0[0] += x[0] * cf0;
    a1[0] += x[0] * cf1;
    x += incX;
    i++;
  }
}
```