user: 请从给定的汇编指令推导出C源代码：

```asm
<stdlib_strided_dcuminabs>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %r8d,0x2c(%rsp)
  test   %edi,%edi
  jle    L16
  mov    %edx,%r14d
  mov    %rcx,%r13
  mov    %edi,%ebx
  mov    %rsi,%r15
  mov    %rsi,%rcx
  xor    %edx,%edx
  test   %r14d,%r14d
  js     L65
  mov    0x2c(%rsp),%esi
  mov    %r13,%rax
  mov    %esi,%r12d
  test   %esi,%esi
  js     L12
L37:
  movsd  (%rcx),%xmm0
  andpd  D1(%rip),%xmm0
  mov    %edx,0x10(%rsp)
  movsd  %xmm0,(%rax)
  movsd  %xmm0,0x8(%rsp)
  call   <stdlib_base_is_nan@plt>
  movsd  0x8(%rsp),%xmm0
  test   %rax,%rax
  je     L105
  mov    $0x1,%r14d
L26:
  movsd  %xmm0,0x8(%rsp)
  call   <stdlib_base_is_nan@plt>
  test   %rax,%rax
  je     L16
  cmp    %ebx,%r14d
  jge    L16
  movslq 0x2c(%rsp),%rbp
  movslq %r12d,%r12
  movsd  0x8(%rsp),%xmm0
  lea    0x0(%r13,%r12,8),%rax
  shl    $0x3,%rbp
  cs nopw 0x0(%rax,%rax,1)
L76:
  add    $0x1,%r14d
  movsd  %xmm0,(%rax)
  add    %rbp,%rax
  cmp    %r14d,%ebx
  jne    L76
L16:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nop
L105:
  cmp    $0x1,%ebx
  je     L124
  movslq 0x2c(%rsp),%rax
  movslq 0x10(%rsp),%rdx
  movslq %r12d,%rbp
  movapd %xmm0,%xmm1
  mov    %rax,0x20(%rsp)
  movslq %r14d,%rax
  mov    $0x1,%r14d
  lea    0x0(,%rax,8),%rdi
  add    %rax,%rdx
  mov    %rdi,0x18(%rsp)
  lea    (%r15,%rdx,8),%r15
  nopl   0x0(%rax)
L118:
  movsd  (%r15),%xmm0
  andpd  D1(%rip),%xmm0
  movsd  %xmm1,0x10(%rsp)
  mov    %ebp,%r12d
  movsd  %xmm0,0x8(%rsp)
  call   <stdlib_base_is_nan@plt>
  movsd  0x8(%rsp),%xmm0
  test   %rax,%rax
  jne    L26
  movsd  0x10(%rsp),%xmm1
  add    $0x1,%r14d
  add    0x18(%rsp),%r15
  minsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13,%rbp,8)
  movapd %xmm0,%xmm1
  add    0x20(%rsp),%rbp
  cmp    %r14d,%ebx
  jne    L118
L40:
  add    $0x38,%rsp
  movapd %xmm1,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1050 <stdlib_base_is_nan@plt>
  nop
L12:
  mov    $0x1,%r12d
  sub    %ebx,%r12d
  imul   %esi,%r12d
  movslq %r12d,%rax
  add    %esi,%r12d
  lea    0x0(%r13,%rax,8),%rax
  jmp    L37
  nopl   (%rax)
L65:
  mov    $0x1,%edx
  sub    %edi,%edx
  imul   %r14d,%edx
  movslq %edx,%rax
  lea    (%rsi,%rax,8),%rcx
  mov    0x2c(%rsp),%esi
  mov    %r13,%rax
  mov    %esi,%r12d
  test   %esi,%esi
  jns    L37
  jmp    L12
L124:
  movapd %xmm0,%xmm1
  jmp    L40
```

tool_calls: name: disassemble_data
arguments: data_label: D1
data_type: qword[2]

tool: {"D1": [0x7fffffffffffffff,0x0000000000000000]}

assistant: ```c
void stdlib_strided_dcuminabs(const int64_t N, const double *X, const int64_t strideX, double *Y, const int64_t strideY) {
  double min;
  int64_t ix;
  int64_t iy;
  int64_t i;
  double v;
  if (N <= 0) {
    return;
  }
  if (strideX < 0) {
    ix = (1 - N) * strideX;
  } else {
    ix = 0;
  }
  if (strideY < 0) {
    iy = (1 - N) * strideY;
  } else {
    iy = 0;
  }
  min = fabs(X[ix]);
  Y[iy] = min;
  iy += strideY;
  i = 1;
  if (!stdlib_base_is_nan(min)) {
    for (; i < N; i++) {
      ix += strideX;
      v = fabs(X[ix]);
      if (stdlib_base_is_nan(v)) {
        min = v;
        break;
      }
      if (v < min) {
        min = v;
      }
      Y[iy] = min;
      iy += strideY;
    }
  }
  if (stdlib_base_is_nan(min)) {
    for (; i < N; i++) {
      Y[iy] = min;
      iy += strideY;
    }
  }
  return;
}
```