user: Please show me the C implementation that matches the assembly code listed below:

```asm
<stdlib_strided_dsumpw>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x90,%rsp
  mov    %edi,-0x74(%rbp)
  mov    %rsi,-0x80(%rbp)
  mov    %edx,-0x78(%rbp)
  cmpl   $0x0,-0x74(%rbp)
  jg     L167
  pxor   %xmm0,%xmm0
  jmp    L25
L167:
  cmpl   $0x1,-0x74(%rbp)
  je     L172
  cmpl   $0x0,-0x78(%rbp)
  jne    L146
L172:
  mov    -0x80(%rbp),%rax
  movsd  (%rax),%xmm0
  jmp    L25
L146:
  cmpl   $0x0,-0x78(%rbp)
  jns    L97
  mov    $0x1,%eax
  sub    -0x74(%rbp),%eax
  mov    -0x78(%rbp),%edx
  imul   %edx,%eax
  mov    %eax,-0x68(%rbp)
  jmp    L35
L97:
  movl   $0x0,-0x68(%rbp)
L35:
  cmpl   $0x7,-0x74(%rbp)
  jg     L66
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movl   $0x0,-0x64(%rbp)
  jmp    L132
L29:
  mov    -0x68(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x48(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    -0x78(%rbp),%eax
  add    %eax,-0x68(%rbp)
  addl   $0x1,-0x64(%rbp)
L132:
  mov    -0x64(%rbp),%eax
  cmp    -0x74(%rbp),%eax
  jl     L29
  movsd  -0x48(%rbp),%xmm0
  jmp    L25
L66:
  cmpl   $0x80,-0x74(%rbp)
  jg     L144
  mov    -0x68(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x68(%rbp),%edx
  mov    -0x78(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0x78(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x68(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x78(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x68(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0x78(%rbp),%eax
  lea    0x0(,%rax,4),%edx
  mov    -0x68(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x78(%rbp),%edx
  mov    %edx,%eax
  shl    $0x2,%eax
  add    %eax,%edx
  mov    -0x68(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x78(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  add    %eax,%eax
  mov    %eax,%edx
  mov    -0x68(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x78(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  sub    %edx,%eax
  mov    %eax,%ecx
  mov    -0x68(%rbp),%eax
  add    %ecx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x78(%rbp),%eax
  shl    $0x3,%eax
  add    %eax,-0x68(%rbp)
  mov    -0x74(%rbp),%eax
  cltd
  shr    $0x1d,%edx
  add    %edx,%eax
  and    $0x7,%eax
  sub    %edx,%eax
  mov    %eax,-0x5c(%rbp)
  movl   $0x8,-0x64(%rbp)
  jmp    L139
L94:
  mov    -0x68(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x40(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x68(%rbp),%edx
  mov    -0x78(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x38(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0x78(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x68(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x30(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x78(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x68(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x28(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0x78(%rbp),%eax
  lea    0x0(,%rax,4),%edx
  mov    -0x68(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x20(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x78(%rbp),%edx
  mov    %edx,%eax
  shl    $0x2,%eax
  add    %eax,%edx
  mov    -0x68(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x18(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x78(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  add    %eax,%eax
  mov    %eax,%edx
  mov    -0x68(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x10(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x78(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  sub    %edx,%eax
  mov    %eax,%ecx
  mov    -0x68(%rbp),%eax
  add    %ecx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x78(%rbp),%eax
  shl    $0x3,%eax
  add    %eax,-0x68(%rbp)
  addl   $0x8,-0x64(%rbp)
L139:
  mov    -0x74(%rbp),%eax
  sub    -0x5c(%rbp),%eax
  cmp    %eax,-0x64(%rbp)
  jl     L94
  movsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x38(%rbp),%xmm1
  movsd  -0x30(%rbp),%xmm0
  addsd  -0x28(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x20(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  -0x18(%rbp),%xmm2
  movsd  -0x10(%rbp),%xmm0
  addsd  -0x8(%rbp),%xmm0
  addsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  jmp    L118
L96:
  mov    -0x68(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x48(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    -0x78(%rbp),%eax
  add    %eax,-0x68(%rbp)
  addl   $0x1,-0x64(%rbp)
L118:
  mov    -0x64(%rbp),%eax
  cmp    -0x74(%rbp),%eax
  jl     L96
  movsd  -0x48(%rbp),%xmm0
  jmp    L25
L144:
  mov    -0x74(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  mov    %eax,-0x60(%rbp)
  mov    -0x60(%rbp),%eax
  cltd
  shr    $0x1d,%edx
  add    %edx,%eax
  and    $0x7,%eax
  sub    %edx,%eax
  sub    %eax,-0x60(%rbp)
  cmpl   $0x0,-0x78(%rbp)
  jns    L89
  mov    -0x60(%rbp),%eax
  sub    -0x74(%rbp),%eax
  imul   -0x78(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x58(%rbp)
  mov    -0x80(%rbp),%rax
  mov    %rax,-0x50(%rbp)
  jmp    L23
L89:
  mov    -0x80(%rbp),%rax
  mov    %rax,-0x58(%rbp)
  mov    -0x60(%rbp),%eax
  imul   -0x78(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x50(%rbp)
L23:
  mov    -0x78(%rbp),%edx
  mov    -0x58(%rbp),%rcx
  mov    -0x60(%rbp),%eax
  mov    %rcx,%rsi
  mov    %eax,%edi
  call   <stdlib_strided_dsumpw>
  movsd  %xmm0,-0x88(%rbp)
  mov    -0x74(%rbp),%eax
  sub    -0x60(%rbp),%eax
  mov    -0x78(%rbp),%edx
  mov    -0x50(%rbp),%rcx
  mov    %rcx,%rsi
  mov    %eax,%edi
  call   <stdlib_strided_dsumpw>
  addsd  -0x88(%rbp),%xmm0
L25:
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

assistant: ```c
double stdlib_strided_dsumpw(const int64_t N, const double *X, const int64_t stride) {
  double *xp1;
  double *xp2;
  double sum;
  int64_t ix;
  int64_t M;
  int64_t n;
  int64_t i;
  double s0;
  double s1;
  double s2;
  double s3;
  double s4;
  double s5;
  double s6;
  double s7;
  if (N <= 0) {
    return 0.0;
  }
  if (N == 1 || stride == 0) {
    return X[0];
  }
  if (stride < 0) {
    ix = (1 - N) * stride;
  } else {
    ix = 0;
  }
  if (N < 8) {
    sum = 0.0;
    for (i = 0; i < N; i++) {
      sum += X[ix];
      ix += stride;
    }
    return sum;
  }
  if (N <= 128) {
    s0 = X[ix];
    s1 = X[ix + stride];
    s2 = X[ix + (2 * stride)];
    s3 = X[ix + (3 * stride)];
    s4 = X[ix + (4 * stride)];
    s5 = X[ix + (5 * stride)];
    s6 = X[ix + (6 * stride)];
    s7 = X[ix + (7 * stride)];
    ix += 8 * stride;
    M = N % 8;
    for (i = 8; i < N - M; i += 8) {
      s0 += X[ix];
      s1 += X[ix + stride];
      s2 += X[ix + (2 * stride)];
      s3 += X[ix + (3 * stride)];
      s4 += X[ix + (4 * stride)];
      s5 += X[ix + (5 * stride)];
      s6 += X[ix + (6 * stride)];
      s7 += X[ix + (7 * stride)];
      ix += 8 * stride;
    }
    sum = ((s0 + s1) + (s2 + s3)) + ((s4 + s5) + (s6 + s7));
    for (; i < N; i++) {
      sum += X[ix];
      ix += stride;
    }
    return sum;
  }
  n = N / 2;
  n -= n % 8;
  if (stride < 0) {
    xp1 = (double *)X + ((n - N) * stride);
    xp2 = (double *)X;
  } else {
    xp1 = (double *)X;
    xp2 = (double *)X + (n * stride);
  }
  return stdlib_strided_dsumpw(n, xp1, stride) + stdlib_strided_dsumpw(N - n, xp2, stride);
}
```