user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<stdlib_strided_sdsapxsumpw>:
  endbr64
  test   %edi,%edi
  jle    L61
  push   %r15
  movaps %xmm0,%xmm1
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %edi,%ebx
  sub    $0x18,%rsp
  cmp    $0x1,%edi
  je     L31
  movslq %edx,%r12
  test   %r12d,%r12d
  je     L31
  mov    $0x0,%eax
  jns    L76
  mov    $0x1,%eax
  sub    %edi,%eax
  imul   %r12d,%eax
L76:
  cmp    $0x7,%ebx
  jle    L40
  cmp    $0x80,%ebx
  jle    L59
  mov    %ebx,%ebp
  sar    %ebp
  and    $0xfffffff8,%ebp
  test   %r12d,%r12d
  js     L88
  mov    %r12d,%eax
  imul   %ebp,%eax
  cltq
  lea    (%rsi,%rax,4),%r13
L75:
  mov    %r12d,%edx
  mov    %ebp,%edi
  movaps %xmm1,%xmm0
  movss  %xmm1,0xc(%rsp)
  call   <stdlib_strided_sdsapxsumpw>
  movss  0xc(%rsp),%xmm1
  mov    %ebx,%edi
  mov    %r12d,%edx
  movss  %xmm0,0x8(%rsp)
  sub    %ebp,%edi
  mov    %r13,%rsi
  movaps %xmm1,%xmm0
  call   <stdlib_strided_sdsapxsumpw>
  addss  0x8(%rsp),%xmm0
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L31:
  addss  (%rsi),%xmm1
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  movaps %xmm1,%xmm0
  pop    %r14
  pop    %r15
  ret
  cs nopw 0x0(%rax,%rax,1)
L61:
  pxor   %xmm0,%xmm0
  ret
  nopl   (%rax)
L88:
  mov    %ebp,%eax
  mov    %rsi,%r13
  sub    %ebx,%eax
  imul   %r12d,%eax
  cltq
  lea    (%rsi,%rax,4),%rsi
  jmp    L75
  xchg   %ax,%ax
L40:
  cltq
  shl    $0x2,%r12
  pxor   %xmm2,%xmm2
  lea    (%rsi,%rax,4),%rdx
  xor    %eax,%eax
L5:
  movss  (%rdx),%xmm0
  add    $0x1,%eax
  add    %r12,%rdx
  addss  %xmm1,%xmm0
  cvtss2sd %xmm0,%xmm0
  addsd  %xmm0,%xmm2
  cmp    %eax,%ebx
  jne    L5
  add    $0x18,%rsp
  pxor   %xmm0,%xmm0
  pop    %rbx
  cvtsd2ss %xmm2,%xmm0
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L59:
  movslq %eax,%rdx
  add    %r12d,%eax
  lea    (%r12,%r12,1),%r9d
  movss  (%rsi,%rdx,4),%xmm9
  movslq %eax,%rdx
  add    %r12d,%eax
  lea    0x0(,%r12,8),%r15d
  movss  (%rsi,%rdx,4),%xmm0
  movslq %eax,%rdx
  add    %r12d,%eax
  lea    (%r9,%r12,1),%r8d
  movss  (%rsi,%rdx,4),%xmm8
  movslq %eax,%rdx
  add    %r12d,%eax
  addss  %xmm1,%xmm9
  movss  (%rsi,%rdx,4),%xmm5
  movslq %eax,%rdx
  add    %r12d,%eax
  addss  %xmm1,%xmm0
  movss  (%rsi,%rdx,4),%xmm7
  movslq %eax,%rdx
  add    %r12d,%eax
  addss  %xmm1,%xmm8
  movss  (%rsi,%rdx,4),%xmm3
  movslq %eax,%rdx
  add    %r12d,%eax
  addss  %xmm1,%xmm5
  movss  (%rsi,%rdx,4),%xmm4
  movslq %eax,%rdx
  addss  %xmm1,%xmm7
  mov    %r15d,%r10d
  movss  (%rsi,%rdx,4),%xmm6
  addss  %xmm1,%xmm3
  mov    %ebx,%edx
  add    %r12d,%eax
  addss  %xmm1,%xmm4
  and    $0xfffffff8,%edx
  mov    %eax,0x8(%rsp)
  lea    (%r8,%r8,1),%r11d
  addss  %xmm1,%xmm6
  sub    %r12d,%r10d
  cvtss2sd %xmm9,%xmm9
  lea    0x0(,%r12,4),%edi
  lea    (%rdi,%r12,1),%ebp
  cvtss2sd %xmm0,%xmm0
  cvtss2sd %xmm5,%xmm5
  cvtss2sd %xmm7,%xmm7
  cvtss2sd %xmm8,%xmm8
  cvtss2sd %xmm3,%xmm3
  cvtss2sd %xmm4,%xmm4
  cvtss2sd %xmm6,%xmm6
  cmp    $0x8,%edx
  je     L95
  movslq %r15d,%r13
  cltq
  mov    $0x8,%ecx
  movslq %r12d,%r14
  shl    $0x2,%r13
  lea    (%rsi,%rax,4),%rax
  movslq %r9d,%r9
  movslq %r8d,%r8
  movslq %edi,%rdi
  movslq %ebp,%rbp
  movslq %r11d,%r11
  movslq %r10d,%r10
  nopl   (%rax)
L111:
  movss  (%rax),%xmm2
  add    $0x8,%ecx
  addss  %xmm1,%xmm2
  cvtss2sd %xmm2,%xmm2
  addsd  %xmm2,%xmm9
  movss  (%rax,%r14,4),%xmm2
  addss  %xmm1,%xmm2
  cvtss2sd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  movss  (%rax,%r9,4),%xmm2
  addss  %xmm1,%xmm2
  cvtss2sd %xmm2,%xmm2
  addsd  %xmm2,%xmm8
  movss  (%rax,%r8,4),%xmm2
  addss  %xmm1,%xmm2
  cvtss2sd %xmm2,%xmm2
  addsd  %xmm2,%xmm5
  movss  (%rax,%rdi,4),%xmm2
  addss  %xmm1,%xmm2
  cvtss2sd %xmm2,%xmm2
  addsd  %xmm2,%xmm7
  movss  (%rax,%rbp,4),%xmm2
  addss  %xmm1,%xmm2
  cvtss2sd %xmm2,%xmm2
  addsd  %xmm2,%xmm3
  movss  (%rax,%r11,4),%xmm2
  addss  %xmm1,%xmm2
  cvtss2sd %xmm2,%xmm2
  addsd  %xmm2,%xmm4
  movss  (%rax,%r10,4),%xmm2
  add    %r13,%rax
  addss  %xmm1,%xmm2
  cvtss2sd %xmm2,%xmm2
  addsd  %xmm2,%xmm6
  cmp    %edx,%ecx
  jl     L111
  mov    0x8(%rsp),%eax
  sub    $0x9,%edx
  shr    $0x3,%edx
  add    %r15d,%eax
  imul   %edx,%r15d
  lea    0x10(,%rdx,8),%edx
  add    %r15d,%eax
  mov    %eax,0x8(%rsp)
L95:
  addsd  %xmm9,%xmm0
  addsd  %xmm8,%xmm5
  addsd  %xmm7,%xmm3
  addsd  %xmm6,%xmm4
  addsd  %xmm5,%xmm0
  addsd  %xmm4,%xmm3
  addsd  %xmm3,%xmm0
  cmp    %ebx,%edx
  jge    L17
  movslq 0x8(%rsp),%rax
  shl    $0x2,%r12
  lea    (%rsi,%rax,4),%rax
  nopw   0x0(%rax,%rax,1)
L92:
  movss  (%rax),%xmm2
  add    $0x1,%edx
  add    %r12,%rax
  addss  %xmm1,%xmm2
  cvtss2sd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  cmp    %edx,%ebx
  jne    L92
L17:
  add    $0x18,%rsp
  cvtsd2ss %xmm0,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

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