user: 请将以下的汇编代码转换成C语言代码：

```asm
<_dmvec_daxpy>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rdi,%r14
  mov    %rsi,0x40(%rsp)
  mov    %rdx,0x38(%rsp)
  movsd  %xmm0,0x8(%rsp)
  mov    %ecx,%r15d
  mov    %r8d,0x4c(%rsp)
  mov    %r9d,0x14(%rsp)
  mov    0x98(%rsp),%eax
  sub    $0x3,%eax
  mov    %eax,0x48(%rsp)
  test   %eax,%eax
  jle    L57
  mov    %rdx,%rbp
  movslq %r8d,%rbx
  lea    0x0(,%rbx,8),%rax
  mov    %rax,0x18(%rsp)
  lea    (%r9,%r9,1),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    %rdx,0x30(%rsp)
  lea    0x0(,%r8,4),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rdi
  mov    %rdi,0x20(%rsp)
  shl    $0x4,%rbx
  add    %rsi,%rbx
  shl    $0x4,%rax
  mov    %rax,0x28(%rsp)
  mov    %rsi,%r12
  mov    $0x0,%r13d
  mov    %r14,(%rsp)
  mov    0x90(%rsp),%r14d
L71:
  mov    0x18(%rsp),%rax
  lea    (%r12,%rax,1),%rcx
  sub    $0x8,%rsp
  push   %r14
  movsd  0x18(%rsp),%xmm0
  mov    0x24(%rsp),%r9d
  mov    %rbp,%r8
  mov    %r12,%rdx
  mov    %r15d,%esi
  mov    0x10(%rsp),%rdi
  call   <_inner_vec2_daxpy@plt>
  mov    0x28(%rsp),%rax
  lea    (%rax,%rbx,1),%rcx
  mov    %r14d,(%rsp)
  movsd  0x18(%rsp),%xmm0
  mov    0x24(%rsp),%r9d
  mov    0x40(%rsp),%rax
  lea    0x0(%rbp,%rax,1),%r8
  mov    %rbx,%rdx
  mov    %r15d,%esi
  mov    0x10(%rsp),%rdi
  call   <_inner_vec2_daxpy@plt>
  add    0x38(%rsp),%rbp
  mov    0x30(%rsp),%rax
  add    %rax,%r12
  add    $0x4,%r13d
  add    %rax,%rbx
  add    $0x10,%rsp
  cmp    0x48(%rsp),%r13d
  jl     L71
  mov    (%rsp),%r14
  mov    0x98(%rsp),%eax
  lea    -0x4(%rax),%ebx
  shr    $0x2,%ebx
  add    $0x1,%ebx
  mov    %ebx,%edx
  mov    0x28(%rsp),%rax
  imul   %rdx,%rax
  add    %rax,0x38(%rsp)
  mov    0x20(%rsp),%rax
  imul   %rdx,%rax
  add    %rax,0x40(%rsp)
  shl    $0x2,%ebx
L42:
  cmp    %ebx,0x98(%rsp)
  je     L66
  mov    0x98(%rsp),%eax
  sub    $0x1,%eax
  cmp    %ebx,%eax
  jg     L28
L46:
  cmp    0x98(%rsp),%ebx
  jl     L10
L66:
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L57:
  mov    $0x0,%ebx
  jmp    L42
L28:
  movslq 0x4c(%rsp),%rax
  lea    0x0(,%rax,8),%r12
  mov    0x40(%rsp),%rdx
  lea    (%rdx,%r12,1),%rbp
  sub    $0x8,%rsp
  mov    0x98(%rsp),%eax
  push   %rax
  movsd  0x18(%rsp),%xmm0
  mov    0x24(%rsp),%r13d
  mov    %r13d,%r9d
  mov    0x48(%rsp),%r8
  mov    %rbp,%rcx
  mov    %r15d,%esi
  mov    %r14,%rdi
  call   <_inner_vec2_daxpy@plt>
  lea    0x0(%r13,%r13,1),%eax
  cltq
  mov    0x48(%rsp),%rdi
  lea    (%rdi,%rax,8),%rax
  mov    %rax,0x48(%rsp)
  lea    0x0(%rbp,%r12,1),%rax
  mov    %rax,0x50(%rsp)
  add    $0x2,%ebx
  add    $0x10,%rsp
  jmp    L46
L10:
  mov    0x90(%rsp),%r9d
  movsd  0x8(%rsp),%xmm0
  mov    0x14(%rsp),%r8d
  mov    0x38(%rsp),%rcx
  mov    0x40(%rsp),%rdx
  mov    %r15d,%esi
  mov    %r14,%rdi
  call   <_inner_vec_daxpy@plt>
  jmp    L66
```

assistant: ```c
void _dmvec_daxpy(double *Yc, const double *Aroot, const double *Xc, double alpha, int incY, int ldA, int incX, int nRE, int nC) {
  register int j, k;
  register const double *a0, *a1, *a2, *a3;
  const double *Ac;
  Ac = Aroot;
  for (j = 0; j < nC - 3; j += 4) {
    a0 = Ac;
    a1 = a0 + ldA;
    a2 = a1 + ldA;
    a3 = a2 + ldA;
    _inner_vec2_daxpy(Yc, incY, a0, a1, Xc, incX, alpha, nRE);
    Xc += 2 * incX;
    _inner_vec2_daxpy(Yc, incY, a2, a3, Xc, incX, alpha, nRE);
    Xc += 2 * incX;
    Ac += 4 * ldA;
  }
  if (j == nC)
    return;
  if (j < nC - 1) {
    a0 = Ac;
    a1 = a0 + ldA;
    _inner_vec2_daxpy(Yc, incY, a0, a1, Xc, incX, alpha, nRE);
    Xc += 2 * incX;
    Ac += 2 * ldA;
    j += 2;
  }
  if (j < nC) {
    a0 = Ac;
    _inner_vec_daxpy(Yc, incY, a0, Xc, incX, alpha, nRE);
    Xc += incX;
    Ac += ldA;
    j++;
  }
}
```