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

```asm
<stdlib_strided_dnanvarianceyc>:
  endbr64
  movapd %xmm0,%xmm9
  test   %edi,%edi
  jle    L80
  cmp    $0x1,%edi
  je     L216
  test   %edx,%edx
  je     L216
  js     L177
  pxor   %xmm3,%xmm3
  cvtsi2sd %edi,%xmm3
  cmp    $0x1,%edx
  jne    L24
  pxor   %xmm8,%xmm8
  movsd  D27(%rip),%xmm1
  mov    %rsi,%rcx
  xor    %eax,%eax
  movapd %xmm8,%xmm0
  jmp    L30
  nopw   0x0(%rax,%rax,1)
L76:
  addsd  %xmm1,%xmm0
  add    $0x1,%eax
  add    $0x8,%rcx
  comisd %xmm0,%xmm3
  jbe    L92
L30:
  movsd  (%rcx),%xmm2
  ucomisd %xmm2,%xmm2
  jp     L76
L92:
  ucomisd %xmm0,%xmm3
  jp     L19
  je     L233
L19:
  movsd  D27(%rip),%xmm5
  add    %edx,%eax
  addsd  %xmm5,%xmm0
  comisd %xmm0,%xmm3
  jbe    L255
  cmp    $0x1,%edx
  jne    L121
  cltq
  movapd %xmm8,%xmm7
  movapd %xmm5,%xmm4
  lea    (%rsi,%rax,8),%rax
  nop
L102:
  movsd  (%rax),%xmm1
  ucomisd %xmm1,%xmm1
  jp     L179
  addsd  %xmm5,%xmm4
  movapd %xmm5,%xmm6
  addsd  %xmm1,%xmm2
  mulsd  %xmm4,%xmm1
  movapd %xmm4,%xmm10
  subsd  %xmm5,%xmm10
  mulsd  %xmm4,%xmm10
  subsd  %xmm2,%xmm1
  divsd  %xmm10,%xmm6
  mulsd  %xmm1,%xmm6
  mulsd  %xmm6,%xmm1
  addsd  %xmm1,%xmm7
L179:
  addsd  %xmm5,%xmm0
  add    $0x8,%rax
  comisd %xmm0,%xmm3
  ja     L102
L61:
  subsd  %xmm9,%xmm4
  comisd %xmm4,%xmm8
  jae    L233
  divsd  %xmm4,%xmm7
  movapd %xmm7,%xmm0
L7:
  ret
  nopl   0x0(%rax)
L216:
  movsd  (%rsi),%xmm0
  ucomisd %xmm0,%xmm0
  jp     L80
  pxor   %xmm1,%xmm1
  pxor   %xmm8,%xmm8
  cvtsi2sd %edi,%xmm1
  movapd %xmm8,%xmm0
  subsd  %xmm9,%xmm1
  comisd %xmm8,%xmm1
  ja     L7
L233:
  movapd %xmm8,%xmm0
  divsd  %xmm8,%xmm0
  ret
  nopl   0x0(%rax,%rax,1)
L80:
  pxor   %xmm8,%xmm8
  movapd %xmm8,%xmm0
  divsd  %xmm8,%xmm0
  ret
L177:
  mov    $0x1,%eax
  pxor   %xmm3,%xmm3
  sub    %edi,%eax
  cvtsi2sd %edi,%xmm3
  imul   %edx,%eax
L144:
  movslq %edx,%rdi
  movslq %eax,%rcx
  pxor   %xmm8,%xmm8
  movsd  D27(%rip),%xmm1
  shl    $0x3,%rdi
  lea    (%rsi,%rcx,8),%rcx
  movapd %xmm8,%xmm0
  jmp    L213
  nopl   0x0(%rax)
L200:
  addsd  %xmm1,%xmm0
  add    %edx,%eax
  add    %rdi,%rcx
  comisd %xmm0,%xmm3
  jbe    L92
L213:
  movsd  (%rcx),%xmm2
  ucomisd %xmm2,%xmm2
  jp     L200
  jmp    L92
  nopw   0x0(%rax,%rax,1)
L121:
  movslq %edx,%rdx
  movapd %xmm8,%xmm7
  movapd %xmm5,%xmm4
  cltq
  shl    $0x3,%rdx
  lea    (%rsi,%rax,8),%rax
  cs nopw 0x0(%rax,%rax,1)
L140:
  movsd  (%rax),%xmm1
  ucomisd %xmm1,%xmm1
  jp     L134
  addsd  %xmm5,%xmm4
  movapd %xmm5,%xmm6
  addsd  %xmm1,%xmm2
  mulsd  %xmm4,%xmm1
  movapd %xmm4,%xmm10
  subsd  %xmm5,%xmm10
  mulsd  %xmm4,%xmm10
  subsd  %xmm2,%xmm1
  divsd  %xmm10,%xmm6
  mulsd  %xmm1,%xmm6
  mulsd  %xmm6,%xmm1
  addsd  %xmm1,%xmm7
L134:
  addsd  %xmm5,%xmm0
  add    %rdx,%rax
  comisd %xmm0,%xmm3
  ja     L140
  jmp    L61
L255:
  movapd %xmm8,%xmm7
  movapd %xmm5,%xmm4
  jmp    L61
L24:
  xor    %eax,%eax
  jmp    L144
```

tool_calls: name: disassemble_data
arguments: tag: D27
type_infered: f64

tool: D27: 1

assistant: ```c
double stdlib_strided_dnanvarianceyc(const int64_t N, const double correction, const double *X, const int64_t stride) {
  double sum;
  int64_t ix;
  double nc;
  double n;
  double S;
  double v;
  double d;
  double i;
  if (N <= 0) {
    return 0.0 / 0.0;
  }
  if (N == 1 || stride == 0) {
    v = X[0];
    if (v == v && (double)N - correction > 0.0) {
      return 0.0;
    }
    return 0.0 / 0.0;
  }
  if (stride < 0) {
    ix = (1 - N) * stride;
  } else {
    ix = 0;
  }
  for (i = 0; i < N; i++) {
    v = X[ix];
    if (v == v) {
      break;
    }
    ix += stride;
  }
  if (i == N) {
    return 0.0 / 0.0;
  }
  ix += stride;
  sum = v;
  S = 0.0;
  n = 1.0;
  i += 1;
  for (; i < N; i++) {
    v = X[ix];
    if (v == v) {
      n += 1.0;
      sum += v;
      d = (n * v) - sum;
      S += (1.0 / (n * (n - 1.0))) * d * d;
    }
    ix += stride;
  }
  nc = n - correction;
  if (nc <= 0.0) {
    return 0.0 / 0.0;
  }
  return S / nc;
}
```