user: What is the c source code of the assembly code below?

```asm
<_inner_vec4_daxpy>:
  endbr64
  push   %r15
  mov    %rdx,%rax
  push   %r14
  push   %r13
  mov    %rcx,%r13
  push   %r12
  mov    %rdi,%r12
  push   %rbp
  push   %rbx
  mov    0x38(%rsp),%rdx
  mov    %esi,%ebx
  movslq 0x40(%rsp),%rdi
  mov    0x48(%rsp),%r14d
  movsd  (%rdx),%xmm4
  mov    %esi,-0x4(%rsp)
  lea    (%rdx,%rdi,8),%rsi
  lea    (%rdi,%rdi,2),%ecx
  movsd  (%rsi),%xmm3
  movsd  (%rsi,%rdi,8),%xmm2
  mulsd  %xmm0,%xmm4
  movslq %ecx,%rcx
  mulsd  %xmm0,%xmm3
  mulsd  %xmm0,%xmm2
  mulsd  (%rdx,%rcx,8),%xmm0
  cmp    $0x3,%r14d
  jle    L1
  lea    -0x4(%r14),%ecx
  movslq %ebx,%rdx
  mov    %r9,%r10
  mov    %r8,%rdi
  shr    $0x2,%ecx
  mov    %rdx,%rbx
  lea    0x0(,%rdx,8),%r11
  shl    $0x4,%rdx
  lea    0x1(%rcx),%esi
  shl    $0x5,%rbx
  add    %r12,%rdx
  mov    %r12,%rcx
  mov    %rsi,%r15
  mov    %esi,-0x8(%rsp)
  shl    $0x5,%r15
  mov    %rsi,-0x10(%rsp)
  mov    %r13,%rsi
  lea    (%rax,%r15,1),%rbp
  nopl   0x0(%rax)
L58:
  movsd  (%rax),%xmm5
  movsd  (%rsi),%xmm1
  add    $0x20,%rax
  add    $0x20,%rsi
  add    $0x20,%rdi
  add    $0x20,%r10
  mulsd  %xmm4,%xmm5
  mulsd  %xmm3,%xmm1
  addsd  %xmm5,%xmm1
  movsd  -0x20(%rdi),%xmm5
  mulsd  %xmm2,%xmm5
  addsd  %xmm1,%xmm5
  movsd  -0x20(%r10),%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm5,%xmm1
  addsd  (%rcx),%xmm1
  movsd  %xmm1,(%rcx)
  movsd  -0x18(%rax),%xmm5
  movsd  -0x18(%rsi),%xmm1
  mulsd  %xmm4,%xmm5
  mulsd  %xmm3,%xmm1
  addsd  %xmm5,%xmm1
  movsd  -0x18(%rdi),%xmm5
  mulsd  %xmm2,%xmm5
  addsd  %xmm1,%xmm5
  movsd  -0x18(%r10),%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm5,%xmm1
  addsd  (%rcx,%r11,1),%xmm1
  movsd  %xmm1,(%rcx,%r11,1)
  movsd  -0x10(%rax),%xmm5
  movsd  -0x10(%rsi),%xmm1
  add    %rbx,%rcx
  mulsd  %xmm4,%xmm5
  mulsd  %xmm3,%xmm1
  addsd  %xmm5,%xmm1
  movsd  -0x10(%rdi),%xmm5
  mulsd  %xmm2,%xmm5
  addsd  %xmm1,%xmm5
  movsd  -0x10(%r10),%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm5,%xmm1
  addsd  (%rdx),%xmm1
  movsd  %xmm1,(%rdx)
  movsd  -0x8(%rax),%xmm5
  movsd  -0x8(%rsi),%xmm1
  mulsd  %xmm4,%xmm5
  mulsd  %xmm3,%xmm1
  addsd  %xmm5,%xmm1
  movsd  -0x8(%rdi),%xmm5
  mulsd  %xmm2,%xmm5
  addsd  %xmm1,%xmm5
  movsd  -0x8(%r10),%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm5,%xmm1
  addsd  (%rdx,%r11,1),%xmm1
  movsd  %xmm1,(%rdx,%r11,1)
  add    %rbx,%rdx
  cmp    %rax,%rbp
  jne    L58
  mov    -0x10(%rsp),%rax
  add    %r15,%r13
  add    %r15,%r8
  add    %r15,%r9
  imul   %rbx,%rax
  add    %rax,%r12
  mov    -0x8(%rsp),%eax
  shl    $0x2,%eax
L14:
  cmp    %eax,%r14d
  je     L2
  lea    -0x1(%r14),%edx
  cmp    %eax,%edx
  jle    L13
  movsd  0x0(%rbp),%xmm5
  add    $0x10,%r13
  add    $0x10,%rbp
  add    $0x2,%eax
  movsd  -0x10(%r13),%xmm1
  movslq -0x4(%rsp),%rdx
  add    $0x10,%r8
  add    $0x10,%r9
  mulsd  %xmm4,%xmm5
  mulsd  %xmm3,%xmm1
  shl    $0x3,%rdx
  addsd  %xmm5,%xmm1
  movsd  -0x10(%r8),%xmm5
  mulsd  %xmm2,%xmm5
  addsd  %xmm1,%xmm5
  movsd  -0x10(%r9),%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm5,%xmm1
  addsd  (%r12),%xmm1
  movsd  %xmm1,(%r12)
  movsd  -0x8(%rbp),%xmm5
  movsd  -0x8(%r13),%xmm1
  add    %rdx,%r12
  mulsd  %xmm4,%xmm5
  mulsd  %xmm3,%xmm1
  addsd  %xmm5,%xmm1
  movsd  -0x8(%r8),%xmm5
  mulsd  %xmm2,%xmm5
  addsd  %xmm1,%xmm5
  movsd  -0x8(%r9),%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm5,%xmm1
  addsd  (%r12),%xmm1
  movsd  %xmm1,(%r12)
  add    %rdx,%r12
L13:
  cmp    %r14d,%eax
  jge    L2
  mulsd  0x0(%rbp),%xmm4
  mulsd  0x0(%r13),%xmm3
  mulsd  (%r8),%xmm2
  mulsd  (%r9),%xmm0
  addsd  %xmm4,%xmm3
  addsd  %xmm3,%xmm2
  addsd  %xmm2,%xmm0
  addsd  (%r12),%xmm0
  movsd  %xmm0,(%r12)
L2:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nop
L1:
  mov    %rax,%rbp
  xor    %eax,%eax
  jmp    L14
```

assistant: ```c
void _inner_vec4_daxpy(double *y0, int incY, const double *a0, const double *a1, const double *a2, const double *a3, const double *x0, int incX, double alpha, int nRE) {
  register int i, ix;
  register double cf0, cf1, cf2, cf3;
  register double ytmp, t0, t1;
  ix = 0;
  cf0 = alpha * x0[ix];
  ix += incX;
  cf1 = alpha * x0[ix];
  ix += incX;
  cf2 = alpha * x0[ix];
  ix += incX;
  cf3 = alpha * x0[ix];
  for (i = 0; i < nRE - 3; i += 4) {
    ytmp = a0[0] * cf0;
    t0 = a1[0] * cf1;
    ytmp += t0;
    t1 = a2[0] * cf2;
    ytmp += t1;
    t0 = a3[0] * cf3;
    ytmp += t0;
    y0[0] += ytmp;
    y0 += incY;
    ytmp = a0[1] * cf0;
    t0 = a1[1] * cf1;
    ytmp += t0;
    t1 = a2[1] * cf2;
    ytmp += t1;
    t0 = a3[1] * cf3;
    ytmp += t0;
    y0[0] += ytmp;
    y0 += incY;
    ytmp = a0[2] * cf0;
    t0 = a1[2] * cf1;
    ytmp += t0;
    t1 = a2[2] * cf2;
    ytmp += t1;
    t0 = a3[2] * cf3;
    ytmp += t0;
    y0[0] += ytmp;
    y0 += incY;
    ytmp = a0[3] * cf0;
    t0 = a1[3] * cf1;
    ytmp += t0;
    t1 = a2[3] * cf2;
    ytmp += t1;
    t0 = a3[3] * cf3;
    ytmp += t0;
    y0[0] += ytmp;
    y0 += incY;
    a0 += 4;
    a1 += 4;
    a2 += 4;
    a3 += 4;
  }
  if (i == nRE)
    return;
  if (i < nRE - 1) {
    ytmp = a0[0] * cf0;
    t0 = a1[0] * cf1;
    ytmp += t0;
    t1 = a2[0] * cf2;
    ytmp += t1;
    t0 = a3[0] * cf3;
    ytmp += t0;
    y0[0] += ytmp;
    y0 += incY;
    ytmp = a0[1] * cf0;
    t0 = a1[1] * cf1;
    ytmp += t0;
    t1 = a2[1] * cf2;
    ytmp += t1;
    t0 = a3[1] * cf3;
    ytmp += t0;
    y0[0] += ytmp;
    y0 += incY;
    a0 += 2;
    a1 += 2;
    a2 += 2;
    a3 += 2;
    i += 2;
  }
  if (i < nRE) {
    ytmp = a0[0] * cf0;
    t0 = a1[0] * cf1;
    ytmp += t0;
    t1 = a2[0] * cf2;
    ytmp += t1;
    t0 = a3[0] * cf3;
    ytmp += t0;
    y0[0] += ytmp;
  }
}
```