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

```asm
<stdlib_strided_dssumors>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %edi,-0x24(%rbp)
  mov    %rsi,-0x30(%rbp)
  mov    %edx,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  cmpl   $0x0,-0x24(%rbp)
  jg     L56
  movsd  -0x8(%rbp),%xmm0
  jmp    L144
L56:
  cmpl   $0x1,-0x24(%rbp)
  je     L0
  cmpl   $0x0,-0x28(%rbp)
  jne    L42
L0:
  mov    -0x30(%rbp),%rax
  movss  (%rax),%xmm0
  cvtss2sd %xmm0,%xmm0
  jmp    L144
L42:
  cmpl   $0x1,-0x28(%rbp)
  jne    L118
  mov    -0x24(%rbp),%edx
  movslq %edx,%rax
  imul   $0x2aaaaaab,%rax,%rax
  shr    $0x20,%rax
  mov    %edx,%esi
  sar    $0x1f,%esi
  mov    %eax,%ecx
  sub    %esi,%ecx
  mov    %ecx,%eax
  add    %eax,%eax
  add    %ecx,%eax
  add    %eax,%eax
  sub    %eax,%edx
  mov    %edx,-0xc(%rbp)
  cmpl   $0x0,-0xc(%rbp)
  jle    L158
  movl   $0x0,-0x10(%rbp)
  jmp    L17
L98:
  mov    -0x10(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  cvtss2sd %xmm0,%xmm0
  movsd  -0x8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  addl   $0x1,-0x10(%rbp)
L17:
  mov    -0x10(%rbp),%eax
  cmp    -0xc(%rbp),%eax
  jl     L98
L158:
  cmpl   $0x5,-0x24(%rbp)
  jg     L126
  movsd  -0x8(%rbp),%xmm0
  jmp    L144
L126:
  mov    -0xc(%rbp),%eax
  mov    %eax,-0x10(%rbp)
  jmp    L97
L108:
  mov    -0x10(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd %xmm0,%xmm1
  mov    -0x10(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  cvtss2sd %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  mov    -0x10(%rbp),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  cvtss2sd %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  mov    -0x10(%rbp),%eax
  cltq
  add    $0x3,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  cvtss2sd %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  mov    -0x10(%rbp),%eax
  cltq
  add    $0x4,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  cvtss2sd %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  mov    -0x10(%rbp),%eax
  cltq
  add    $0x5,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  cvtss2sd %xmm0,%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  addl   $0x6,-0x10(%rbp)
L97:
  mov    -0x10(%rbp),%eax
  cmp    -0x24(%rbp),%eax
  jl     L108
  movsd  -0x8(%rbp),%xmm0
  jmp    L144
L118:
  cmpl   $0x0,-0x28(%rbp)
  jns    L170
  mov    $0x1,%eax
  sub    -0x24(%rbp),%eax
  mov    -0x28(%rbp),%edx
  imul   %edx,%eax
  mov    %eax,-0x14(%rbp)
  jmp    L121
L170:
  movl   $0x0,-0x14(%rbp)
L121:
  movl   $0x0,-0x10(%rbp)
  jmp    L145
L132:
  mov    -0x14(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  cvtss2sd %xmm0,%xmm0
  movsd  -0x8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x28(%rbp),%eax
  add    %eax,-0x14(%rbp)
  addl   $0x1,-0x10(%rbp)
L145:
  mov    -0x10(%rbp),%eax
  cmp    -0x24(%rbp),%eax
  jl     L132
  movsd  -0x8(%rbp),%xmm0
L144:
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  pop    %rbp
  ret
```

assistant: ```c
double stdlib_strided_dssumors(const int64_t N, const float *X, const int64_t stride) {
  double sum;
  int64_t ix;
  int64_t m;
  int64_t i;
  sum = 0.0;
  if (N <= 0) {
    return sum;
  }
  if (N == 1 || stride == 0) {
    return X[0];
  }
  if (stride == 1) {
    m = N % 6;
    if (m > 0) {
      for (i = 0; i < m; i++) {
        sum += (double)X[i];
      }
    }
    if (N < 6) {
      return sum;
    }
    for (i = m; i < N; i += 6) {
      sum += (double)X[i] + (double)X[i + 1] + (double)X[i + 2] + (double)X[i + 3] + (double)X[i + 4] + (double)X[i + 5];
    }
    return sum;
  }
  if (stride < 0) {
    ix = (1 - N) * stride;
  } else {
    ix = 0;
  }
  for (i = 0; i < N; i++) {
    sum += (double)X[ix];
    ix += stride;
  }
  return sum;
}
```