user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<cdflib_cumfnc>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xb8,%rsp
  mov    %rdx,(%rsp)
  mov    %r8,%r13
  mov    %r9,%r14
  mov    %fs:0x28,%rax
  mov    %rax,0xa8(%rsp)
  xor    %eax,%eax
  pxor   %xmm0,%xmm0
  comisd (%rdi),%xmm0
  jb     L85
  movq   $0x0,(%r8)
  mov    D8(%rip),%rax
  mov    %rax,(%r9)
L72:
  mov    0xa8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L191
  add    $0xb8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L85:
  mov    %rdi,%r12
  mov    %rsi,%r15
  movsd  (%rcx),%xmm0
  movsd  D73(%rip),%xmm1
  comisd %xmm0,%xmm1
  ja     L99
  mulsd  D163(%rip),%xmm0
  movapd %xmm0,%xmm6
  cvttsd2si %xmm0,%ebx
  test   %ebx,%ebx
  mov    $0x1,%eax
  cmove  %eax,%ebx
  lea    0x1(%rbx),%ebp
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebp,%xmm0
  movsd  %xmm0,0x78(%rsp)
  pxor   %xmm5,%xmm5
  cvtsi2sd %ebx,%xmm5
  movsd  %xmm5,0x18(%rsp)
  movsd  %xmm6,0x10(%rsp)
  movapd %xmm6,%xmm0
  call   <log@plt>
  mulsd  0x18(%rsp),%xmm0
  movapd %xmm0,%xmm5
  subsd  0x10(%rsp),%xmm5
  movsd  %xmm5,0x8(%rsp)
  lea    0x78(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_gamma_log@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  0x8(%rsp),%xmm0
  subsd  %xmm1,%xmm0
  call   <exp@plt>
  movsd  %xmm0,0x8(%rsp)
  movsd  (%r15),%xmm0
  movapd %xmm0,%xmm2
  mulsd  (%r12),%xmm2
  mov    (%rsp),%rax
  movsd  (%rax),%xmm1
  movapd %xmm1,%xmm4
  addsd  %xmm2,%xmm4
  movapd %xmm1,%xmm3
  divsd  %xmm4,%xmm3
  movsd  %xmm3,0x58(%rsp)
  comisd D163(%rip),%xmm3
  ja     L45
  movsd  D8(%rip),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0x50(%rsp)
  jmp    L144
L99:
  mov    %r9,%r8
  mov    %r13,%rcx
  mov    (%rsp),%rdx
  mov    $0x0,%eax
  call   <cdflib_cumf@plt>
  jmp    L72
L45:
  divsd  %xmm4,%xmm2
  movsd  %xmm2,0x50(%rsp)
  movsd  D8(%rip),%xmm3
  subsd  %xmm2,%xmm3
  movsd  %xmm3,0x58(%rsp)
L144:
  mulsd  D163(%rip),%xmm0
  addsd  0x18(%rsp),%xmm0
  movsd  %xmm0,0x80(%rsp)
  mulsd  D163(%rip),%xmm1
  movsd  %xmm1,0x88(%rsp)
  lea    0x58(%rsp),%rcx
  lea    0x50(%rsp),%rdx
  lea    0x88(%rsp),%rsi
  lea    0x80(%rsp),%rdi
  sub    $0x8,%rsp
  lea    0x4c(%rsp),%rax
  push   %rax
  lea    0x58(%rsp),%r9
  lea    0x80(%rsp),%r8
  mov    $0x0,%eax
  call   <cdflib_beta_inc@plt>
  movsd  D163(%rip),%xmm0
  mulsd  (%r15),%xmm0
  addsd  0x28(%rsp),%xmm0
  movapd %xmm0,%xmm7
  movsd  %xmm0,0x70(%rsp)
  mov    0x10(%rsp),%rax
  movsd  D163(%rip),%xmm0
  mulsd  (%rax),%xmm0
  movsd  %xmm0,0x78(%rsp)
  movsd  0x80(%rsp),%xmm6
  movsd  %xmm6,0x28(%rsp)
  movsd  0x18(%rsp),%xmm5
  mulsd  %xmm6,%xmm5
  movsd  %xmm5,0x10(%rsp)
  addsd  %xmm7,%xmm0
  movsd  %xmm0,0xa0(%rsp)
  movsd  %xmm7,0x30(%rsp)
  movapd %xmm7,%xmm0
  addsd  D8(%rip),%xmm0
  movsd  %xmm0,0xa8(%rsp)
  lea    0xa0(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_gamma_log@plt>
  mov    %eax,%r12d
  lea    0xa8(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_gamma_log@plt>
  sub    %eax,%r12d
  lea    0x78(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_gamma_log@plt>
  mov    %eax,0x4c(%rsp)
  movsd  0x60(%rsp),%xmm2
  movsd  %xmm2,0x40(%rsp)
  movapd %xmm2,%xmm0
  call   <log@plt>
  movq   %xmm0,%r15
  movsd  0x68(%rsp),%xmm0
  call   <log@plt>
  movsd  0x78(%rsp),%xmm4
  movsd  %xmm4,0x38(%rsp)
  sub    0x4c(%rsp),%r12d
  pxor   %xmm1,%xmm1
  cvtsi2sd %r12d,%xmm1
  movq   %r15,%xmm2
  mulsd  0x30(%rsp),%xmm2
  addsd  %xmm2,%xmm1
  mulsd  %xmm4,%xmm0
  addsd  %xmm1,%xmm0
  call   <exp@plt>
  add    $0x10,%rsp
  movsd  D100(%rip),%xmm1
  comisd (%rsp),%xmm1
  ja     L181
  movsd  0x70(%rsp),%xmm2
  movsd  0x8(%rsp),%xmm3
  mov    $0x0,%eax
  movsd  D139(%rip),%xmm6
  movsd  D8(%rip),%xmm5
  mov    $0x1,%edx
  movapd %xmm1,%xmm7
L194:
  movsd  (%rsp),%xmm4
  mulsd  %xmm6,%xmm4
  movapd %xmm2,%xmm1
  mulsd  %xmm3,%xmm1
  comisd %xmm1,%xmm4
  ja     L26
  test   %ebx,%ebx
  jle    L6
  pxor   %xmm1,%xmm1
  cvtsi2sd %ebx,%xmm1
  divsd  0x10(%rsp),%xmm1
  mulsd  %xmm1,%xmm3
  sub    $0x1,%ebx
  movsd  0x20(%rsp),%xmm10
  subsd  %xmm5,%xmm10
  movsd  %xmm10,0x20(%rsp)
  movapd %xmm10,%xmm1
  addsd  %xmm5,%xmm1
  movsd  0x28(%rsp),%xmm4
  addsd  %xmm10,%xmm4
  mulsd  0x30(%rsp),%xmm4
  divsd  %xmm4,%xmm1
  mulsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm2
  movapd %xmm2,%xmm1
  mulsd  %xmm3,%xmm1
  addsd  (%rsp),%xmm1
  movsd  %xmm1,(%rsp)
  mov    %edx,%eax
  comisd %xmm1,%xmm7
  jbe    L194
  movsd  %xmm2,0x70(%rsp)
L181:
  movsd  0x60(%rsp),%xmm0
  subsd  D8(%rip),%xmm0
  addsd  0x28(%rsp),%xmm0
  ucomisd D159(%rip),%xmm0
  jp     L143
  jne    L143
  lea    0x60(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_gamma_log@plt>
  mov    %eax,%ebx
  lea    0x68(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_gamma_log@plt>
  mov    %eax,%r15d
  movsd  0x60(%rsp),%xmm0
  subsd  D8(%rip),%xmm0
  movsd  %xmm0,0x20(%rsp)
  movsd  0x50(%rsp),%xmm0
  call   <log@plt>
  movq   %xmm0,%r12
  movsd  0x58(%rsp),%xmm0
  call   <log@plt>
  neg    %ebx
  sub    %r15d,%ebx
  pxor   %xmm1,%xmm1
  cvtsi2sd %ebx,%xmm1
  movq   %r12,%xmm2
  mulsd  0x20(%rsp),%xmm2
  addsd  %xmm2,%xmm1
  mulsd  0x68(%rsp),%xmm0
  addsd  %xmm1,%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm2
L3:
  movsd  0x60(%rsp),%xmm1
  movsd  0x68(%rsp),%xmm8
  movsd  0x50(%rsp),%xmm7
  movsd  D8(%rip),%xmm3
  movsd  D104(%rip),%xmm6
  movsd  D100(%rip),%xmm5
  movsd  D139(%rip),%xmm9
L135:
  pxor   %xmm4,%xmm4
  cvtsi2sd %ebp,%xmm4
  movsd  0x10(%rsp),%xmm0
  divsd  %xmm4,%xmm0
  mulsd  0x8(%rsp),%xmm0
  movapd %xmm0,%xmm11
  movsd  %xmm0,0x8(%rsp)
  add    $0x1,%ebp
  addsd  %xmm3,%xmm1
  movapd %xmm1,%xmm0
  addsd  %xmm8,%xmm0
  subsd  %xmm6,%xmm0
  mulsd  %xmm7,%xmm0
  movapd %xmm1,%xmm4
  subsd  %xmm3,%xmm4
  divsd  %xmm4,%xmm0
  mulsd  %xmm0,%xmm2
  movsd  0x18(%rsp),%xmm4
  subsd  %xmm2,%xmm4
  movsd  %xmm4,0x18(%rsp)
  movapd %xmm11,%xmm0
  mulsd  %xmm4,%xmm0
  movsd  (%rsp),%xmm4
  addsd  %xmm0,%xmm4
  movsd  %xmm4,(%rsp)
  comisd %xmm4,%xmm5
  ja     L93
  movsd  (%rsp),%xmm4
  mulsd  %xmm9,%xmm4
  comisd %xmm0,%xmm4
  jbe    L135
L93:
  movsd  (%rsp),%xmm7
  movsd  %xmm7,0x0(%r13)
  movsd  D163(%rip),%xmm1
  movapd %xmm1,%xmm0
  subsd  %xmm7,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%r14)
  jmp    L72
L26:
  test   %al,%al
  je     L181
  movsd  %xmm2,0x70(%rsp)
  jmp    L181
L6:
  test   %al,%al
  je     L181
  movsd  %xmm2,0x70(%rsp)
  jmp    L181
L143:
  movsd  %xmm0,0xa0(%rsp)
  lea    0xa0(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_gamma_log@plt>
  mov    %eax,%ebx
  lea    0x60(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_gamma_log@plt>
  sub    %eax,%ebx
  lea    0x68(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_gamma_log@plt>
  mov    %eax,%r15d
  movsd  0x60(%rsp),%xmm0
  subsd  D8(%rip),%xmm0
  movsd  %xmm0,0x20(%rsp)
  movsd  0x50(%rsp),%xmm0
  call   <log@plt>
  movq   %xmm0,%r12
  movsd  0x58(%rsp),%xmm0
  call   <log@plt>
  sub    %r15d,%ebx
  pxor   %xmm1,%xmm1
  cvtsi2sd %ebx,%xmm1
  movq   %r12,%xmm2
  mulsd  0x20(%rsp),%xmm2
  addsd  %xmm2,%xmm1
  mulsd  0x68(%rsp),%xmm0
  addsd  %xmm1,%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm2
  jmp    L3
L191:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: disassemble_data
arguments: data_label: D8
datatype: f64
name: disassemble_data
arguments: data_label: D73
datatype: byte[8]
name: disassemble_data
arguments: data_label: D163
datatype: f64[2]
name: disassemble_data
arguments: data_label: D100
datatype: byte[8]
name: disassemble_data
arguments: data_label: D139
datatype: byte[8]
name: disassemble_data
arguments: data_label: D159
datatype: f64
name: disassemble_data
arguments: data_label: D104
datatype: f64

tool: D8: 1

tool: {"D73": [0xbb,0xbd,0xd7,0xd9,0xdf,0x7c,0xdb,0x3d]}

tool: {"D163": [0.5,1e-20]}

tool: {"D100": [0x23,0x42,0x92,0x0c,0xa1,0x9c,0xc7,0x3b]}

tool: {"D139": [0x2d,0x43,0x1c,0xeb,0xe2,0x36,0x1a,0x3f]}

tool: D159: 0

tool: D104: 2

assistant: ```c
void cdflib_cumfnc(const double *f, const double *dfn, const double *dfd, const double *pnonc, double *cum, double *ccum) {
  static const double half = 0.5e0;
  static const double done = 1.0e0;
  static const double eps = 1.0e-4;
  double dsum, dummy, prod, xx, yy, adn, aup, b, betdn, betup, centwt, dnterm, sum, upterm, xmult, xnonc;
  int i, icent, ierr;
  double T1, T2, T3, T4, T5, T6;
  if (!(*f <= 0.0e0)) {
    goto S10;
  }
  *cum = 0.0e0;
  *ccum = 1.0e0;
  return;
S10:
  if (!(*pnonc < 1.0e-10)) {
    goto S20;
  }
  cdflib_cumf(f, dfn, dfd, cum, ccum);
  return;
S20:
  xnonc = *pnonc / 2.0e0;
  icent = (int)xnonc;
  if (icent == 0) {
    icent = 1;
  }
  T1 = (double)(icent + 1);
  centwt = exp(-xnonc + (double)icent * log(xnonc) - cdflib_gamma_log(&T1));
  prod = *dfn * *f;
  dsum = *dfd + prod;
  yy = *dfd / dsum;
  if (yy > half) {
    xx = prod / dsum;
    yy = done - xx;
  } else {
    xx = done - yy;
  }
  T2 = *dfn * half + (double)icent;
  T3 = *dfd * half;
  cdflib_beta_inc(&T2, &T3, &xx, &yy, &betdn, &dummy, &ierr);
  adn = *dfn / 2.0e0 + (double)icent;
  aup = adn;
  b = *dfd / 2.0e0;
  betup = betdn;
  sum = centwt * betdn;
  xmult = centwt;
  i = icent;
  T4 = adn + b;
  T5 = adn + 1.0e0;
  dnterm = exp(cdflib_gamma_log(&T4) - cdflib_gamma_log(&T5) - cdflib_gamma_log(&b) + adn * log(xx) + b * log(yy));
S30:
  if ((int)(sum < 1.0e-20 || (xmult * betdn) < eps * sum) || i <= 0) {
    goto S40;
  }
  xmult *= ((double)i / xnonc);
  i -= 1;
  adn -= 1.0;
  dnterm = (adn + 1.0) / ((adn + b) * xx) * dnterm;
  betdn += dnterm;
  sum += (xmult * betdn);
  goto S30;
S40:
  i = icent + 1;
  xmult = centwt;
  if (aup - 1.0 + b == 0)
    upterm = exp(-cdflib_gamma_log(&aup) - cdflib_gamma_log(&b) + (aup - 1.0) * log(xx) + b * log(yy));
  else {
    T6 = aup - 1.0 + b;
    upterm = exp(cdflib_gamma_log(&T6) - cdflib_gamma_log(&aup) - cdflib_gamma_log(&b) + (aup - 1.0) * log(xx) + b * log(yy));
  }
  goto S60;
S50:
  if ((int)(sum < 1.0e-20 || (xmult * betup) < eps * sum)) {
    goto S70;
  }
S60:
  xmult *= (xnonc / (double)i);
  i += 1;
  aup += 1.0;
  upterm = (aup + b - 2.0e0) * xx / (aup - 1.0) * upterm;
  betup -= upterm;
  sum += (xmult * betup);
  goto S50;
S70:
  *cum = sum;
  *ccum = 0.5e0 + (0.5e0 - *cum);
  return;
}
```