user: 请将以下的汇编代码转换成C语言代码：

```asm
<stdlib_strided_dsumpw>:
  endbr64
  test   %edi,%edi
  jle    L63
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x18,%rsp
  mov    %edi,%ebp
  mov    %edx,%ebx
  cmp    $0x1,%edi
  je     L20
  test   %edx,%edx
  je     L20
  mov    $0x0,%eax
  js     L65
L15:
  cmp    $0x7,%ebp
  jle    L26
  cmp    $0x80,%ebp
  jle    L70
  mov    %ebp,%eax
  shr    $0x1f,%eax
  add    %ebp,%eax
  mov    %eax,%r13d
  sar    %r13d
  sar    $0x1f,%eax
  shr    $0x1d,%eax
  lea    0x0(%r13,%rax,1),%edx
  and    $0x7,%edx
  sub    %eax,%edx
  sub    %edx,%r13d
  test   %ebx,%ebx
  js     L48
  mov    %ebx,%eax
  imul   %r13d,%eax
  cltq
  lea    (%rsi,%rax,8),%r14
L100:
  mov    %ebx,%edx
  mov    %r13d,%edi
  call   <stdlib_strided_dsumpw>
  movq   %xmm0,%r12
  mov    %ebp,%edi
  sub    %r13d,%edi
  mov    %ebx,%edx
  mov    %r14,%rsi
  call   <stdlib_strided_dsumpw>
  movq   %r12,%xmm4
  addsd  %xmm0,%xmm4
  movapd %xmm4,%xmm0
  jmp    L1
L20:
  movsd  (%rsi),%xmm0
L1:
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L65:
  mov    $0x1,%eax
  sub    %edi,%eax
  imul   %edx,%eax
  jmp    L15
L26:
  movslq %ebx,%rbx
  shl    $0x3,%rbx
  cltq
  lea    (%rsi,%rax,8),%rdx
  mov    $0x0,%eax
  pxor   %xmm0,%xmm0
L173:
  addsd  (%rdx),%xmm0
  add    $0x1,%eax
  add    %rbx,%rdx
  cmp    %eax,%ebp
  jne    L173
  jmp    L1
L70:
  movslq %eax,%rdx
  movsd  (%rsi,%rdx,8),%xmm0
  lea    (%rax,%rbx,1),%edx
  movslq %edx,%rax
  movsd  (%rsi,%rax,8),%xmm7
  lea    (%rbx,%rbx,1),%r10d
  add    %ebx,%edx
  movslq %edx,%rax
  movsd  (%rsi,%rax,8),%xmm3
  lea    (%r10,%rbx,1),%r9d
  add    %ebx,%edx
  movslq %edx,%rax
  movsd  (%rsi,%rax,8),%xmm6
  lea    0x0(,%rbx,4),%r8d
  add    %ebx,%edx
  movslq %edx,%rax
  movsd  (%rsi,%rax,8),%xmm1
  lea    (%r8,%rbx,1),%r13d
  add    %ebx,%edx
  movslq %edx,%rax
  movsd  (%rsi,%rax,8),%xmm5
  lea    (%r9,%r9,1),%r12d
  add    %ebx,%edx
  movslq %edx,%rax
  movsd  (%rsi,%rax,8),%xmm2
  lea    0x0(,%rbx,8),%r15d
  mov    %r15d,%r11d
  sub    %ebx,%r11d
  add    %ebx,%edx
  movslq %edx,%rax
  movsd  (%rsi,%rax,8),%xmm4
  add    %ebx,%edx
  mov    %ebp,%ecx
  sar    $0x1f,%ecx
  shr    $0x1d,%ecx
  lea    0x0(%rbp,%rcx,1),%eax
  and    $0x7,%eax
  sub    %ecx,%eax
  mov    %ebp,%edi
  sub    %eax,%edi
  cmp    $0x8,%edi
  jle    L158
  movslq %r15d,%rax
  shl    $0x3,%rax
  mov    %rax,0x8(%rsp)
  movslq %edx,%rax
  lea    (%rsi,%rax,8),%rax
  mov    $0x8,%ecx
  movslq %ebx,%r14
  movslq %r10d,%r10
  movslq %r9d,%r9
  movslq %r8d,%r8
  movslq %r13d,%r13
  movslq %r12d,%r12
  movslq %r11d,%r11
L117:
  addsd  (%rax),%xmm0
  addsd  (%rax,%r14,8),%xmm7
  addsd  (%rax,%r10,8),%xmm3
  addsd  (%rax,%r9,8),%xmm6
  addsd  (%rax,%r8,8),%xmm1
  addsd  (%rax,%r13,8),%xmm5
  addsd  (%rax,%r12,8),%xmm2
  addsd  (%rax,%r11,8),%xmm4
  add    $0x8,%ecx
  add    0x8(%rsp),%rax
  cmp    %edi,%ecx
  jl     L117
  add    %r15d,%edx
  sub    $0x9,%edi
  shr    $0x3,%edi
  imul   %edi,%r15d
  add    %r15d,%edx
  lea    0x10(,%rdi,8),%eax
L33:
  addsd  %xmm7,%xmm0
  addsd  %xmm6,%xmm3
  addsd  %xmm3,%xmm0
  addsd  %xmm5,%xmm1
  addsd  %xmm4,%xmm2
  addsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm0
  cmp    %eax,%ebp
  jle    L1
  movslq %ebx,%rbx
  shl    $0x3,%rbx
  movslq %edx,%rdx
  lea    (%rsi,%rdx,8),%rdx
L91:
  addsd  (%rdx),%xmm0
  add    $0x1,%eax
  add    %rbx,%rdx
  cmp    %eax,%ebp
  jne    L91
  jmp    L1
L158:
  mov    $0x8,%eax
  jmp    L33
L48:
  mov    %r13d,%eax
  sub    %ebp,%eax
  imul   %ebx,%eax
  cltq
  mov    %rsi,%r14
  lea    (%rsi,%rax,8),%rsi
  jmp    L100
L63:
  pxor   %xmm0,%xmm0
  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);
}
```