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

```asm
<stdlib_strided_dmskrange>:
  endbr64
  test   %rdi,%rdi
  jle    L116
  push   %r15
  mov    %r8,%r15
  push   %r14
  mov    %rdx,%r14
  push   %r13
  mov    %rdi,%r13
  push   %r12
  xor    %r12d,%r12d
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  test   %rdx,%rdx
  js     L82
  xor    %ebx,%ebx
  test   %r15,%r15
  js     L170
L159:
  xor    %ebp,%ebp
  nopl   (%rax)
L20:
  cmpq   $0x0,(%rcx,%rbx,8)
  je     L166
  add    $0x1,%rbp
  add    %r14,%r12
  add    %r15,%rbx
  cmp    %rbp,%r13
  jne    L20
L77:
  add    $0x28,%rsp
  pxor   %xmm0,%xmm0
  divsd  %xmm0,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L170:
  mov    $0x1,%ebx
  sub    %r13,%rbx
  imul   %r15,%rbx
  jmp    L159
  cs nopw 0x0(%rax,%rax,1)
L82:
  mov    $0x1,%r12d
  xor    %ebx,%ebx
  sub    %rdi,%r12
  imul   %rdx,%r12
  test   %r15,%r15
  jns    L159
  jmp    L170
  cs nopw 0x0(%rax,%rax,1)
L166:
  mov    %rcx,0x8(%rsp)
  cmp    %rbp,%r13
  je     L77
  movsd  (%rsi,%r12,8),%xmm3
  mov    %rsi,(%rsp)
  movapd %xmm3,%xmm0
  movsd  %xmm3,0x18(%rsp)
  call   <stdlib_base_is_nan@plt>
  test   %rax,%rax
  jne    L77
  add    $0x1,%rbp
  cmp    %rbp,%r13
  jle    L150
  mov    0x8(%rsp),%rcx
  mov    (%rsp),%rsi
  lea    0x0(,%r15,8),%rdx
  add    %rbx,%r15
  movsd  0x18(%rsp),%xmm3
  lea    (%rcx,%r15,8),%rbx
  lea    0x0(,%r14,8),%r15
  add    %r12,%r14
  lea    (%rsi,%r14,8),%r12
  movsd  %xmm3,0x10(%rsp)
  jmp    L24
  nopw   0x0(%rax,%rax,1)
L87:
  maxsd  0x18(%rsp),%xmm0
  movsd  %xmm0,0x18(%rsp)
L156:
  add    $0x1,%rbp
  add    %rdx,%rbx
  add    %r15,%r12
  cmp    %rbp,%r13
  je     L12
L24:
  cmpq   $0x0,(%rbx)
  jne    L156
  movsd  (%r12),%xmm0
  mov    %rdx,0x8(%rsp)
  movsd  %xmm0,(%rsp)
  call   <stdlib_base_is_nan@plt>
  test   %rax,%rax
  jne    L77
  movsd  0x10(%rsp),%xmm1
  movsd  (%rsp),%xmm0
  mov    0x8(%rsp),%rdx
  comisd %xmm0,%xmm1
  jbe    L87
  add    $0x1,%rbp
  movsd  %xmm0,0x10(%rsp)
  add    %rdx,%rbx
  add    %r15,%r12
  cmp    %rbp,%r13
  jne    L24
L12:
  movsd  0x18(%rsp),%xmm0
  subsd  0x10(%rsp),%xmm0
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L116:
  pxor   %xmm0,%xmm0
  divsd  %xmm0,%xmm0
  ret
L150:
  movsd  0x18(%rsp),%xmm5
  movsd  %xmm5,0x10(%rsp)
  jmp    L12
```

assistant: ```c
double stdlib_strided_dmskrange(const int64_t N, const double *X, const int64_t strideX, const uint8_t *Mask, const int64_t strideMask) {
  double max;
  double min;
  int64_t ix;
  int64_t im;
  int64_t i;
  double v;
  if (N <= 0) {
    return 0.0 / 0.0;
  }
  if (strideX < 0) {
    ix = (1 - N) * strideX;
  } else {
    ix = 0;
  }
  if (strideMask < 0) {
    im = (1 - N) * strideMask;
  } else {
    im = 0;
  }
  for (i = 0; i < N; i++) {
    if (Mask[im] == 0) {
      break;
    }
    ix += strideX;
    im += strideMask;
  }
  if (i == N) {
    return 0.0 / 0.0;
  }
  min = X[ix];
  if (stdlib_base_is_nan(min)) {
    return 0.0 / 0.0;
  }
  max = min;
  i += 1;
  for (; i < N; i++) {
    ix += strideX;
    im += strideMask;
    if (Mask[im]) {
      continue;
    }
    v = X[ix];
    if (stdlib_base_is_nan(v)) {
      return 0.0 / 0.0;
    }
    if (v < min) {
      min = v;
    } else if (v > max) {
      max = v;
    }
  }
  return max - min;
}
```