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

```asm
<cdflib_cdfnor>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rdx,%r12
  mov    %rcx,0x8(%rsp)
  mov    %r8,0x10(%rsp)
  mov    0x70(%rsp),%r14
  mov    0x78(%rsp),%r15
  mov    %fs:0x28,%rax
  mov    %rax,0x28(%rsp)
  xor    %eax,%eax
  movq   $0x0,(%r15)
  movl   $0x0,(%r14)
  mov    (%rdi),%eax
  lea    -0x1(%rax),%edx
  cmp    $0x3,%edx
  jbe    L83
  movsd  D37(%rip),%xmm0
  test   %eax,%eax
  jle    L208
L309:
  movsd  %xmm0,(%r15)
  movl   $0xffffffff,(%r14)
L168:
  mov    0x28(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L144
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L83:
  mov    %rdi,%rbx
  mov    %rsi,%rbp
  mov    %r9,%r13
  cmp    $0x1,%eax
  je     L132
  movsd  (%rsi),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm0,%xmm1
  jae    L286
  comisd D87(%rip),%xmm0
  jbe    L61
  comisd %xmm0,%xmm1
  jae    L235
  movsd  D87(%rip),%xmm0
L97:
  movsd  %xmm0,(%r15)
  movl   $0xfffffffe,(%r14)
  jmp    L168
L208:
  movsd  D87(%rip),%xmm0
  jmp    L309
L235:
  movapd %xmm1,%xmm0
  jmp    L97
L286:
  pxor   %xmm0,%xmm0
  jmp    L97
L61:
  movsd  (%r12),%xmm1
  pxor   %xmm2,%xmm2
  comisd %xmm1,%xmm2
  jae    L127
  comisd D87(%rip),%xmm1
  jbe    L248
  movapd %xmm2,%xmm0
  comisd %xmm1,%xmm2
  jae    L47
  movsd  D87(%rip),%xmm0
L47:
  movsd  %xmm0,(%r15)
  movl   $0xfffffffd,(%r14)
  jmp    L168
L127:
  pxor   %xmm0,%xmm0
  jmp    L47
L248:
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x18(%rsp)
  mov    $0x1,%edi
  call   <dpmpar@plt>
  movsd  D17(%rip),%xmm2
  movsd  0x18(%rsp),%xmm3
  movapd %xmm3,%xmm1
  subsd  %xmm2,%xmm1
  subsd  %xmm2,%xmm1
  andpd  D69(%rip),%xmm1
  mulsd  D84(%rip),%xmm0
  comisd %xmm0,%xmm1
  jbe    L89
  pxor   %xmm0,%xmm0
  comisd %xmm3,%xmm0
  ja     L19
  movsd  D87(%rip),%xmm0
L19:
  movsd  %xmm0,(%r15)
  movl   $0x3,(%r14)
  jmp    L168
L89:
  mov    (%rbx),%eax
  cmp    $0x4,%eax
  jne    L163
L247:
  mov    %r12,%rsi
  mov    %rbp,%rdi
  call   <cdflib_dinvnr@plt>
  movapd %xmm0,%xmm1
  mov    0x8(%rsp),%rax
  movsd  (%rax),%xmm0
  mov    0x10(%rsp),%rax
  subsd  (%rax),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13)
  jmp    L168
L163:
  movsd  0x0(%r13),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm0,%xmm1
  jae    L267
  cmp    $0x1,%eax
  je     L275
  cmp    $0x2,%eax
  je     L187
  cmp    $0x3,%eax
  je     L183
  cmp    $0x4,%eax
  jne    L168
  jmp    L247
L275:
  mov    0x8(%rsp),%rax
  movsd  (%rax),%xmm1
  mov    0x10(%rsp),%rax
  subsd  (%rax),%xmm1
  divsd  %xmm0,%xmm1
  movsd  %xmm1,0x20(%rsp)
  lea    0x20(%rsp),%rdi
  mov    %r12,%rdx
  mov    %rbp,%rsi
  call   <cdflib_cumnor@plt>
  jmp    L168
L187:
  mov    %r12,%rsi
  mov    %rbp,%rdi
  call   <cdflib_dinvnr@plt>
  mulsd  0x0(%r13),%xmm0
  mov    0x10(%rsp),%rax
  addsd  (%rax),%xmm0
  mov    0x8(%rsp),%rax
  movsd  %xmm0,(%rax)
  jmp    L168
L183:
  mov    %r12,%rsi
  mov    %rbp,%rdi
  call   <cdflib_dinvnr@plt>
  movapd %xmm0,%xmm1
  mulsd  0x0(%r13),%xmm1
  mov    0x8(%rsp),%rax
  movsd  (%rax),%xmm0
  subsd  %xmm1,%xmm0
  mov    0x10(%rsp),%rax
  movsd  %xmm0,(%rax)
  jmp    L168
L132:
  movsd  (%r9),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm0,%xmm1
  jb     L275
L267:
  movq   $0x0,(%r15)
  movl   $0xfffffffa,(%r14)
  jmp    L168
L144:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: label: D37
data_type: f64
name: extract_symbolic_data
arguments: label: D87
data_type: f64[2]
name: extract_symbolic_data
arguments: label: D17
data_type: f64
name: extract_symbolic_data
arguments: label: D69
data_type: qword[2]
name: extract_symbolic_data
arguments: label: D84
data_type: f64

tool: D37: 4

tool: {"D87": [1.0,4.0]}

tool: D17: 0.5

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

tool: D84: 3

assistant: ```c
void cdflib_cdfnor(const int *which, double *p, double *q, double *x, double *mean, double *sd, int *status, double *bound) {
  static const int K1 = 1;
  double z, pq;
  *status = 0;
  *bound = 0.0;
  *status = 0;
  if (!(*which < 1 || *which > 4)) {
    goto S30;
  }
  if (!(*which < 1)) {
    goto S10;
  }
  *bound = 1.0e0;
  goto S20;
S10:
  *bound = 4.0e0;
S20:
  *status = -1;
  return;
S30:
  if (*which == 1) {
    goto S70;
  }
  if (!(*p <= 0.0e0 || *p > 1.0e0)) {
    goto S60;
  }
  if (!(*p <= 0.0e0)) {
    goto S40;
  }
  *bound = 0.0e0;
  goto S50;
S40:
  *bound = 1.0e0;
S50:
  *status = -2;
  return;
S70:
S60:
  if (*which == 1) {
    goto S110;
  }
  if (!(*q <= 0.0e0 || *q > 1.0e0)) {
    goto S100;
  }
  if (!(*q <= 0.0e0)) {
    goto S80;
  }
  *bound = 0.0e0;
  goto S90;
S80:
  *bound = 1.0e0;
S90:
  *status = -3;
  return;
S110:
S100:
  if (*which == 1) {
    goto S150;
  }
  pq = *p + *q;
  if (!(fabs(pq - 0.5e0 - 0.5e0) > 3.0e0 * dpmpar(K1))) {
    goto S140;
  }
  if (!(pq < 0.0e0)) {
    goto S120;
  }
  *bound = 0.0e0;
  goto S130;
S120:
  *bound = 1.0e0;
S130:
  *status = 3;
  return;
S150:
S140:
  if (*which == 4) {
    goto S170;
  }
  if (!(*sd <= 0.0e0)) {
    goto S160;
  }
  *bound = 0.0e0;
  *status = -6;
  return;
S170:
S160:
  if (1 == *which) {
    z = (*x - *mean) / *sd;
    cdflib_cumnor(&z, p, q);
  } else if (2 == *which) {
    z = cdflib_dinvnr(p, q);
    *x = *sd * z + *mean;
  } else if (3 == *which) {
    z = cdflib_dinvnr(p, q);
    *mean = *x - *sd * z;
  } else if (4 == *which) {
    z = cdflib_dinvnr(p, q);
    *sd = (*x - *mean) / z;
  }
  return;
}
```