user: What would the C source look like for the assembly code provided here?

```asm
<__ip__>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r13
  push   %r12
  push   %rbx
  mov    %rdi,%rcx
  mov    %rsi,%rax
  mov    %edx,-0x1c(%rbp)
  pxor   %xmm3,%xmm3
  movapd %xmm3,%xmm7
  mov    -0x1c(%rbp),%edx
  sar    $0x3,%edx
  mov    %edx,%r12d
  mov    $0x0,%ebx
  mov    -0x1c(%rbp),%edx
  and    $0x7,%edx
  mov    %edx,%r13d
  jmp    L10
L21:
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm2
  movapd %xmm7,%xmm4
  addsd  %xmm4,%xmm2
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm3
  addsd  %xmm2,%xmm3
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm4
  addsd  %xmm3,%xmm4
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm5
  addsd  %xmm4,%xmm5
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm6
  addsd  %xmm5,%xmm6
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm7
  addsd  %xmm6,%xmm7
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm7,%xmm2
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  %xmm2,%xmm0
  movapd %xmm0,%xmm7
  add    $0x1,%ebx
L10:
  mov    %r12d,%edx
  lea    -0x1(%rdx),%r12d
  test   %edx,%edx
  jne    L21
  cmp    $0x7,%r13d
  ja     L24
  mov    %r13d,%edx
  lea    0x0(,%rdx,4),%rsi
  lea    D6(%rip),%rdx
  mov    (%rsi,%rdx,1),%edx
  movslq %edx,%rdx
  lea    D6(%rip),%rsi
  add    %rsi,%rdx
  notrack jmp *%rdx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm7,%xmm5
  addsd  %xmm5,%xmm0
  movapd %xmm0,%xmm7
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm7,%xmm6
  addsd  %xmm6,%xmm0
  movapd %xmm0,%xmm7
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  %xmm7,%xmm0
  movapd %xmm0,%xmm7
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm7,%xmm3
  addsd  %xmm3,%xmm0
  movapd %xmm0,%xmm7
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm7,%xmm4
  addsd  %xmm4,%xmm0
  movapd %xmm0,%xmm7
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm7,%xmm5
  addsd  %xmm5,%xmm0
  movapd %xmm0,%xmm7
  add    $0x1,%ebx
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rdx),%xmm1
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm7,%xmm6
  addsd  %xmm6,%xmm0
  movapd %xmm0,%xmm7
  add    $0x1,%ebx
L24:
  movq   %xmm7,%rax
  movq   %rax,%xmm0
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %rbp
  ret
```

assistant: ```c
double __ip__(dp1, dp2, len)
register double *dp1, *dp2;
int len;
{
  register int j, k;
  register int i;
  register double sum;
  sum = 0.0;
  i = len >> 3;
  j = 0;
  k = len & 7;
  while (i--) {
    sum += dp1[j] * dp2[j];
    j++;
    sum += dp1[j] * dp2[j];
    j++;
    sum += dp1[j] * dp2[j];
    j++;
    sum += dp1[j] * dp2[j];
    j++;
    sum += dp1[j] * dp2[j];
    j++;
    sum += dp1[j] * dp2[j];
    j++;
    sum += dp1[j] * dp2[j];
    j++;
    sum += dp1[j] * dp2[j];
    j++;
  }
  switch (k) {
  case 7:
    sum += dp1[j] * dp2[j];
    j++;
  case 6:
    sum += dp1[j] * dp2[j];
    j++;
  case 5:
    sum += dp1[j] * dp2[j];
    j++;
  case 4:
    sum += dp1[j] * dp2[j];
    j++;
  case 3:
    sum += dp1[j] * dp2[j];
    j++;
  case 2:
    sum += dp1[j] * dp2[j];
    j++;
  case 1:
    sum += dp1[j] * dp2[j];
    j++;
  }
  return sum;
}
```