user: Translate the below assembly code into C programming language:

```asm
<_inner_vec2_daxpy>:
  endbr64
  push   %r14
  movslq %r9d,%r9
  mov    %rdi,%r11
  mov    %rcx,%rax
  push   %r13
  push   %r12
  push   %rbp
  mov    %rdx,%rbp
  push   %rbx
  movsd  (%r8),%xmm2
  movslq %esi,%rbx
  mov    0x30(%rsp),%r12d
  mulsd  %xmm0,%xmm2
  mulsd  (%r8,%r9,8),%xmm0
  cmp    $0x3,%r12d
  jle    L29
  lea    -0x4(%r12),%r10d
  movslq %ebx,%rdx
  mov    %rbp,%rsi
  shr    $0x2,%r10d
  mov    %rdx,%r8
  lea    0x0(,%rdx,8),%rdi
  shl    $0x4,%rdx
  lea    0x1(%r10),%r14d
  shl    $0x5,%r8
  add    %r11,%rdx
  mov    %r14,%r13
  mov    %r14,%r10
  shl    $0x5,%r13
  lea    (%rcx,%r13,1),%r9
  mov    %r11,%rcx
  nopl   0x0(%rax,%rax,1)
L18:
  movsd  (%rsi),%xmm1
  movsd  (%rax),%xmm3
  add    $0x20,%rax
  add    $0x20,%rsi
  mulsd  %xmm0,%xmm3
  mulsd  %xmm2,%xmm1
  addsd  %xmm3,%xmm1
  addsd  (%rcx),%xmm1
  movsd  %xmm1,(%rcx)
  movsd  -0x18(%rsi),%xmm1
  movsd  -0x18(%rax),%xmm3
  mulsd  %xmm2,%xmm1
  mulsd  %xmm0,%xmm3
  addsd  %xmm3,%xmm1
  addsd  (%rcx,%rdi,1),%xmm1
  movsd  %xmm1,(%rcx,%rdi,1)
  movsd  -0x10(%rsi),%xmm1
  movsd  -0x10(%rax),%xmm3
  add    %r8,%rcx
  mulsd  %xmm0,%xmm3
  mulsd  %xmm2,%xmm1
  addsd  %xmm3,%xmm1
  addsd  (%rdx),%xmm1
  movsd  %xmm1,(%rdx)
  movsd  -0x8(%rsi),%xmm1
  movsd  -0x8(%rax),%xmm3
  mulsd  %xmm2,%xmm1
  mulsd  %xmm0,%xmm3
  addsd  %xmm3,%xmm1
  addsd  (%rdx,%rdi,1),%xmm1
  movsd  %xmm1,(%rdx,%rdi,1)
  add    %r8,%rdx
  cmp    %r9,%rax
  jne    L18
  imul   %r14,%r8
  add    %r13,%rbp
  shl    $0x2,%r10d
  add    %r8,%r11
L28:
  cmp    %r10d,%r12d
  je     L56
  lea    -0x1(%r12),%eax
  cmp    %r10d,%eax
  jle    L34
  shl    $0x3,%rbx
  add    $0x10,%rbp
  add    $0x10,%r9
  add    $0x2,%r10d
  movsd  -0x10(%rbp),%xmm1
  movsd  -0x10(%r9),%xmm3
  mulsd  %xmm2,%xmm1
  mulsd  %xmm0,%xmm3
  addsd  %xmm3,%xmm1
  addsd  (%r11),%xmm1
  movsd  %xmm1,(%r11)
  movsd  -0x8(%rbp),%xmm1
  movsd  -0x8(%r9),%xmm3
  add    %rbx,%r11
  mulsd  %xmm2,%xmm1
  mulsd  %xmm0,%xmm3
  addsd  %xmm3,%xmm1
  addsd  (%r11),%xmm1
  movsd  %xmm1,(%r11)
  add    %rbx,%r11
L34:
  cmp    %r12d,%r10d
  jge    L56
  mulsd  0x0(%rbp),%xmm2
  mulsd  (%r9),%xmm0
  addsd  %xmm0,%xmm2
  addsd  (%r11),%xmm2
  movsd  %xmm2,(%r11)
L56:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  nopl   0x0(%rax)
L29:
  mov    %rcx,%r9
  xor    %r10d,%r10d
  jmp    L28
```

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