user: Converte this assembly snippet into its equivalent C code:

```asm
<_dmvec_daxpy>:
  endbr64
  push   %r15
  mov    %rdi,%r15
  push   %r14
  mov    %ecx,%r14d
  push   %r13
  mov    %r9d,%r13d
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    0x90(%rsp),%eax
  mov    %rsi,0x40(%rsp)
  mov    %rdx,0x38(%rsp)
  mov    %eax,0xc(%rsp)
  mov    0x98(%rsp),%eax
  mov    %r8d,0x4c(%rsp)
  mov    %eax,0x48(%rsp)
  sub    $0x3,%eax
  mov    %eax,0x2c(%rsp)
  movsd  %xmm0,(%rsp)
  test   %eax,%eax
  jle    L15
  movslq %r8d,%rbx
  mov    %rsi,%r12
  mov    %rdx,%rbp
  xor    %r10d,%r10d
  lea    0x0(,%rbx,8),%rax
  lea    0x0(,%r8,4),%r11d
  shl    $0x4,%rbx
  mov    %rax,0x10(%rsp)
  lea    (%r9,%r9,1),%eax
  movslq %r11d,%r11
  add    %rsi,%rbx
  cltq
  shl    $0x3,%r11
  lea    0x0(,%rax,8),%rdi
  shl    $0x4,%rax
  mov    %rdi,0x30(%rsp)
  mov    %rax,0x20(%rsp)
  nop
L68:
  mov    0x10(%rsp),%rax
  mov    %r10d,0x28(%rsp)
  mov    %rbp,%r8
  mov    %r12,%rdx
  mov    %r11,0x18(%rsp)
  sub    $0x8,%rsp
  mov    %r13d,%r9d
  mov    %r15,%rdi
  mov    0x14(%rsp),%esi
  lea    (%r12,%rax,1),%rcx
  push   %rsi
  movsd  0x10(%rsp),%xmm0
  mov    %r14d,%esi
  call   <_inner_vec2_daxpy@plt>
  mov    0x20(%rsp),%rax
  mov    %rbx,%rdx
  mov    %r13d,%r9d
  mov    0x1c(%rsp),%esi
  movsd  0x10(%rsp),%xmm0
  mov    %r15,%rdi
  lea    (%rax,%rbx,1),%rcx
  mov    0x40(%rsp),%rax
  mov    %esi,(%rsp)
  mov    %r14d,%esi
  lea    0x0(%rbp,%rax,1),%r8
  call   <_inner_vec2_daxpy@plt>
  mov    0x28(%rsp),%r11
  mov    0x38(%rsp),%r10d
  add    0x30(%rsp),%rbp
  pop    %rcx
  add    %r11,%r12
  add    $0x4,%r10d
  add    %r11,%rbx
  pop    %rsi
  cmp    0x2c(%rsp),%r10d
  jl     L68
  mov    0x48(%rsp),%eax
  lea    -0x4(%rax),%ebx
  mov    0x20(%rsp),%rax
  shr    $0x2,%ebx
  lea    0x1(%rbx),%edx
  imul   %rdx,%rax
  mov    %rdx,%rbx
  add    %rax,0x38(%rsp)
  imul   %r11,%rdx
  add    %rdx,0x40(%rsp)
  shl    $0x2,%ebx
L11:
  mov    0x48(%rsp),%eax
  cmp    %ebx,%eax
  je     L30
  sub    $0x1,%eax
  cmp    %ebx,%eax
  jg     L31
  cmp    0x48(%rsp),%ebx
  jl     L18
L30:
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L31:
  movslq 0x4c(%rsp),%rbp
  mov    0x40(%rsp),%rdx
  mov    %r15,%rdi
  mov    %r13d,%r9d
  sub    $0x8,%rsp
  mov    %r14d,%esi
  add    $0x2,%ebx
  mov    0x14(%rsp),%eax
  shl    $0x3,%rbp
  lea    (%rdx,%rbp,1),%r12
  push   %rax
  movsd  0x10(%rsp),%xmm0
  mov    %r12,%rcx
  mov    0x48(%rsp),%r8
  call   <_inner_vec2_daxpy@plt>
  mov    0x48(%rsp),%rdi
  lea    0x0(%r13,%r13,1),%eax
  cltq
  lea    (%rdi,%rax,8),%rax
  mov    %rax,0x48(%rsp)
  lea    (%r12,%rbp,1),%rax
  mov    %rax,0x50(%rsp)
  pop    %rax
  pop    %rdx
  cmp    0x48(%rsp),%ebx
  jge    L30
L18:
  mov    0xc(%rsp),%r9d
  movsd  (%rsp),%xmm0
  mov    %r13d,%r8d
  mov    %r14d,%esi
  mov    0x38(%rsp),%rcx
  mov    0x40(%rsp),%rdx
  add    $0x58,%rsp
  mov    %r15,%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1070 <_inner_vec_daxpy@plt>
L15:
  xor    %ebx,%ebx
  jmp    L11
```

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