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

```asm
<stdlib_strided_dnansumkbn2>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x78,%rsp
  mov    %rdi,-0x68(%rbp)
  mov    %rsi,-0x70(%rbp)
  mov    %rdx,-0x78(%rbp)
  cmpq   $0x0,-0x68(%rbp)
  jg     L63
  pxor   %xmm0,%xmm0
  jmp    L42
L63:
  cmpq   $0x1,-0x68(%rbp)
  je     L3
  cmpq   $0x0,-0x78(%rbp)
  jne    L74
L3:
  mov    -0x70(%rbp),%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <stdlib_base_is_nan@plt>
  test   %rax,%rax
  je     L27
  pxor   %xmm0,%xmm0
  jmp    L42
L27:
  mov    -0x70(%rbp),%rax
  movsd  (%rax),%xmm0
  jmp    L42
L74:
  cmpq   $0x0,-0x78(%rbp)
  jns    L118
  mov    $0x1,%eax
  sub    -0x68(%rbp),%rax
  mov    -0x78(%rbp),%rdx
  imul   %rdx,%rax
  mov    %rax,-0x48(%rbp)
  jmp    L22
L118:
  movq   $0x0,-0x48(%rbp)
L22:
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movq   $0x0,-0x40(%rbp)
  jmp    L80
L44:
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x20(%rbp),%rax
  movq   %rax,%xmm0
  call   <stdlib_base_is_nan@plt>
  test   %rax,%rax
  jne    L64
  movsd  -0x58(%rbp),%xmm0
  addsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x58(%rbp),%rax
  movq   %rax,%xmm0
  call   <fabs@plt>
  mov    %rax,%rbx
  mov    -0x20(%rbp),%rax
  movq   %rax,%xmm0
  call   <fabs@plt>
  cmp    %rax,%rbx
  jl     L43
  movsd  -0x58(%rbp),%xmm0
  subsd  -0x18(%rbp),%xmm0
  movsd  -0x20(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  jmp    L132
L43:
  movsd  -0x20(%rbp),%xmm0
  subsd  -0x18(%rbp),%xmm0
  movsd  -0x58(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
L132:
  movsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0x38(%rbp),%xmm0
  addsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x38(%rbp),%rax
  movq   %rax,%xmm0
  call   <fabs@plt>
  mov    %rax,%rbx
  mov    -0x28(%rbp),%rax
  movq   %rax,%xmm0
  call   <fabs@plt>
  cmp    %rax,%rbx
  jl     L140
  movsd  -0x38(%rbp),%xmm0
  subsd  -0x18(%rbp),%xmm0
  movsd  -0x28(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  jmp    L158
L140:
  movsd  -0x28(%rbp),%xmm0
  subsd  -0x18(%rbp),%xmm0
  movsd  -0x38(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
L158:
  movsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x50(%rbp),%xmm0
  addsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
L64:
  mov    -0x78(%rbp),%rax
  add    %rax,-0x48(%rbp)
  addq   $0x1,-0x40(%rbp)
L80:
  mov    -0x40(%rbp),%rax
  cmp    -0x68(%rbp),%rax
  jl     L44
  movsd  -0x58(%rbp),%xmm0
  addsd  -0x38(%rbp),%xmm0
  addsd  -0x50(%rbp),%xmm0
L42:
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

assistant: ```c
double stdlib_strided_dnansumkbn2(const int64_t N, const double *X, const int64_t stride) {
  double sum;
  double ccs;
  int64_t ix;
  int64_t i;
  double cs;
  double cc;
  double v;
  double t;
  double c;
  if (N <= 0) {
    return 0.0;
  }
  if (N == 1 || stride == 0) {
    if (stdlib_base_is_nan(X[0])) {
      return 0.0;
    }
    return X[0];
  }
  if (stride < 0) {
    ix = (1 - N) * stride;
  } else {
    ix = 0;
  }
  sum = 0.0;
  ccs = 0.0;
  cs = 0.0;
  for (i = 0; i < N; i++) {
    v = X[ix];
    if (!stdlib_base_is_nan(v)) {
      t = sum + v;
      if (fabs(sum) >= fabs(v)) {
        c = (sum - t) + v;
      } else {
        c = (v - t) + sum;
      }
      sum = t;
      t = cs + c;
      if (fabs(cs) >= fabs(c)) {
        cc = (cs - t) + c;
      } else {
        cc = (c - t) + cs;
      }
      cs = t;
      ccs += cc;
    }
    ix += stride;
  }
  return sum + cs + ccs;
}
```