user: Translate the below assembly code into C programming language:

```asm
<cdflib_cumfnc>:
  endbr64
  push   %r15
  pxor   %xmm6,%xmm6
  push   %r14
  push   %r13
  mov    %r9,%r13
  push   %r12
  mov    %r8,%r12
  push   %rbp
  push   %rbx
  sub    $0xd8,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0xc8(%rsp)
  xor    %eax,%eax
  comisd (%rdi),%xmm6
  jb     L49
  mov    D77(%rip),%rax
  movq   $0x0,(%r8)
  mov    %rax,(%r9)
L98:
  mov    0xc8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L127
  add    $0xd8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L49:
  movsd  (%rcx),%xmm9
  mov    %rdi,%r14
  mov    %rsi,%r10
  mov    %rdx,%r15
  movsd  D150(%rip),%xmm0
  comisd %xmm9,%xmm0
  ja     L21
  mulsd  D50(%rip),%xmm9
  cvttsd2si %xmm9,%ebx
  test   %ebx,%ebx
  je     L196
  lea    0x1(%rbx),%ebp
  pxor   %xmm12,%xmm12
  pxor   %xmm7,%xmm7
  movsd  D101(%rip),%xmm4
  cvtsi2sd %ebp,%xmm12
  movsd  D77(%rip),%xmm2
  cvtsi2sd %ebx,%xmm7
  movsd  %xmm4,0x58(%rsp)
  movapd %xmm12,%xmm0
  movsd  %xmm7,0x8(%rsp)
L154:
  movsd  %xmm0,0x98(%rsp)
  movapd %xmm9,%xmm0
  mov    %r10,0x20(%rsp)
  movsd  %xmm2,0x30(%rsp)
  movsd  %xmm12,0x28(%rsp)
  movsd  %xmm9,0x10(%rsp)
  call   <log@plt>
  mulsd  0x8(%rsp),%xmm0
  movsd  0x10(%rsp),%xmm9
  xor    %eax,%eax
  lea    0x98(%rsp),%rdi
  movsd  %xmm9,0x18(%rsp)
  subsd  %xmm9,%xmm0
  movsd  %xmm0,0x10(%rsp)
  call   <cdflib_gamma_log@plt>
  pxor   %xmm1,%xmm1
  movsd  0x10(%rsp),%xmm0
  cvtsi2sd %eax,%xmm1
  subsd  %xmm1,%xmm0
  call   <exp@plt>
  mov    0x20(%rsp),%r10
  movsd  (%r14),%xmm3
  movapd %xmm0,%xmm4
  movsd  (%r15),%xmm1
  movsd  0x18(%rsp),%xmm9
  movsd  (%r10),%xmm0
  movsd  0x28(%rsp),%xmm12
  movapd %xmm1,%xmm6
  movapd %xmm1,%xmm5
  movsd  0x30(%rsp),%xmm2
  mulsd  %xmm0,%xmm3
  addsd  %xmm3,%xmm6
  divsd  %xmm6,%xmm5
  comisd D50(%rip),%xmm5
  jbe    L134
  divsd  %xmm6,%xmm3
  movapd %xmm2,%xmm5
  subsd  %xmm3,%xmm5
  movsd  %xmm5,0x78(%rsp)
L92:
  mulsd  D50(%rip),%xmm1
  movsd  %xmm12,0x48(%rsp)
  lea    0x70(%rsp),%rdx
  lea    0x78(%rsp),%rcx
  mulsd  D50(%rip),%xmm0
  movsd  %xmm9,0x40(%rsp)
  lea    0xa8(%rsp),%rsi
  lea    0xa0(%rsp),%rdi
  addsd  0x8(%rsp),%xmm0
  movsd  %xmm3,0x70(%rsp)
  sub    $0x8,%rsp
  movsd  %xmm2,0x58(%rsp)
  mov    %r10,0x20(%rsp)
  movsd  %xmm0,0xa8(%rsp)
  movsd  %xmm1,0xb0(%rsp)
  movsd  %xmm4,0x18(%rsp)
  lea    0x6c(%rsp),%rax
  push   %rax
  xor    %eax,%eax
  lea    0x78(%rsp),%r9
  lea    0xa0(%rsp),%r8
  call   <cdflib_beta_inc@plt>
  mov    0x28(%rsp),%r10
  xor    %eax,%eax
  movsd  D50(%rip),%xmm0
  mulsd  (%r15),%xmm0
  movsd  0x20(%rsp),%xmm4
  lea    0xc0(%rsp),%rdi
  lea    0x98(%rsp),%r14
  movsd  D50(%rip),%xmm3
  mulsd  (%r10),%xmm3
  addsd  0x18(%rsp),%xmm3
  movapd %xmm4,%xmm1
  movsd  0xa0(%rsp),%xmm6
  movsd  %xmm4,0x48(%rsp)
  movsd  0x60(%rsp),%xmm2
  mulsd  %xmm6,%xmm1
  movsd  %xmm6,0x40(%rsp)
  movsd  %xmm0,0x98(%rsp)
  addsd  %xmm3,%xmm0
  movsd  %xmm2,0x28(%rsp)
  movsd  %xmm3,0x90(%rsp)
  movsd  %xmm0,0xc0(%rsp)
  movapd %xmm3,%xmm0
  addsd  %xmm2,%xmm0
  movsd  %xmm1,0x38(%rsp)
  movsd  %xmm3,0x30(%rsp)
  movsd  %xmm0,0xc8(%rsp)
  call   <cdflib_gamma_log@plt>
  lea    0xc8(%rsp),%rdi
  mov    %eax,%r15d
  xor    %eax,%eax
  call   <cdflib_gamma_log@plt>
  mov    %r14,%rdi
  sub    %eax,%r15d
  xor    %eax,%eax
  call   <cdflib_gamma_log@plt>
  movsd  0x80(%rsp),%xmm14
  mov    %eax,0x20(%rsp)
  movapd %xmm14,%xmm0
  movsd  %xmm14,0x60(%rsp)
  call   <log@plt>
  movsd  %xmm0,0x18(%rsp)
  movsd  0x88(%rsp),%xmm0
  call   <log@plt>
  mov    0x20(%rsp),%eax
  movsd  0x18(%rsp),%xmm5
  movsd  0x30(%rsp),%xmm3
  movapd %xmm0,%xmm7
  pxor   %xmm0,%xmm0
  movsd  0x98(%rsp),%xmm13
  sub    %eax,%r15d
  mulsd  %xmm3,%xmm5
  cvtsi2sd %r15d,%xmm0
  movsd  %xmm13,0x18(%rsp)
  mulsd  %xmm13,%xmm7
  movsd  %xmm3,0x20(%rsp)
  addsd  %xmm5,%xmm0
  addsd  %xmm7,%xmm0
  call   <exp@plt>
  pop    %rax
  pop    %rdx
  movsd  0x28(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm13
  movapd %xmm0,%xmm7
  movsd  D76(%rip),%xmm11
  movsd  0x18(%rsp),%xmm2
  movsd  0x30(%rsp),%xmm6
  movsd  0x38(%rsp),%xmm4
  comisd %xmm1,%xmm11
  movsd  0x40(%rsp),%xmm9
  movsd  0x48(%rsp),%xmm12
  ja     L170
  movsd  D78(%rip),%xmm15
  movsd  0x50(%rsp),%xmm14
  movapd %xmm4,%xmm10
  xor    %eax,%eax
  movsd  0x90(%rsp),%xmm8
  movsd  0x10(%rsp),%xmm3
  movapd %xmm8,%xmm0
  mulsd  %xmm4,%xmm0
  jmp    L20
  nopw   0x0(%rax,%rax,1)
L182:
  test   %ebx,%ebx
  jle    L175
  pxor   %xmm0,%xmm0
  subsd  %xmm2,%xmm3
  movapd %xmm13,%xmm5
  mov    $0x1,%eax
  cvtsi2sd %ebx,%xmm0
  sub    $0x1,%ebx
  addsd  %xmm3,%xmm5
  divsd  %xmm9,%xmm0
  mulsd  %xmm14,%xmm5
  mulsd  %xmm0,%xmm10
  movapd %xmm3,%xmm0
  addsd  %xmm2,%xmm0
  divsd  %xmm5,%xmm0
  mulsd  %xmm0,%xmm7
  addsd  %xmm7,%xmm8
  movapd %xmm8,%xmm0
  mulsd  %xmm10,%xmm0
  addsd  %xmm0,%xmm1
  comisd %xmm1,%xmm11
  ja     L176
L20:
  movapd %xmm1,%xmm5
  mulsd  %xmm15,%xmm5
  comisd %xmm0,%xmm5
  jbe    L182
L175:
  test   %al,%al
  je     L170
L176:
  movsd  %xmm8,0x90(%rsp)
L170:
  movsd  0x80(%rsp),%xmm0
  pxor   %xmm7,%xmm7
  subsd  %xmm2,%xmm0
  addsd  %xmm13,%xmm0
  ucomisd %xmm7,%xmm0
  jp     L145
  jne    L145
  lea    0x80(%rsp),%rdi
  xor    %eax,%eax
  movsd  %xmm2,0x50(%rsp)
  movsd  %xmm12,0x48(%rsp)
  movsd  %xmm1,0x40(%rsp)
  movsd  %xmm9,0x38(%rsp)
  movsd  %xmm4,0x30(%rsp)
  movsd  %xmm6,0x28(%rsp)
  call   <cdflib_gamma_log@plt>
  mov    %r14,%rdi
  mov    %eax,%ebx
  xor    %eax,%eax
  call   <cdflib_gamma_log@plt>
  movsd  0x70(%rsp),%xmm7
  neg    %ebx
  movsd  0x80(%rsp),%xmm3
  mov    %eax,%r14d
  movapd %xmm7,%xmm0
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm7,0x20(%rsp)
  call   <log@plt>
  movsd  %xmm0,0x8(%rsp)
  movsd  0x78(%rsp),%xmm0
  call   <log@plt>
  movsd  0x88(%rsp),%xmm8
L206:
  movsd  0x10(%rsp),%xmm3
  mulsd  %xmm8,%xmm0
  sub    %r14d,%ebx
  pxor   %xmm5,%xmm5
  movsd  0x50(%rsp),%xmm2
  cvtsi2sd %ebx,%xmm5
  movapd %xmm3,%xmm10
  movsd  %xmm3,0x18(%rsp)
  subsd  %xmm2,%xmm10
  mulsd  0x8(%rsp),%xmm10
  movsd  %xmm2,0x10(%rsp)
  movsd  %xmm8,0x8(%rsp)
  addsd  %xmm10,%xmm5
  addsd  %xmm5,%xmm0
  call   <exp@plt>
  mov    D76(%rip),%rax
  movsd  0x8(%rsp),%xmm8
  movsd  0x10(%rsp),%xmm2
  movsd  0x18(%rsp),%xmm3
  movapd %xmm0,%xmm5
  movsd  0x20(%rsp),%xmm7
  movsd  0x28(%rsp),%xmm6
  movq   %rax,%xmm11
  movsd  0x30(%rsp),%xmm4
  movsd  0x38(%rsp),%xmm9
  movsd  0x40(%rsp),%xmm1
  movsd  0x48(%rsp),%xmm12
  movsd  D78(%rip),%xmm10
  movsd  0x58(%rsp),%xmm13
  jmp    L77
  nopl   0x0(%rax,%rax,1)
L2:
  movapd %xmm1,%xmm12
  mulsd  %xmm10,%xmm12
  comisd %xmm0,%xmm12
  ja     L100
  pxor   %xmm12,%xmm12
  cvtsi2sd %ebp,%xmm12
L77:
  movapd %xmm9,%xmm0
  addsd  %xmm2,%xmm3
  add    $0x1,%ebp
  divsd  %xmm12,%xmm0
  movapd %xmm3,%xmm12
  subsd  %xmm2,%xmm12
  mulsd  %xmm0,%xmm4
  movapd %xmm3,%xmm0
  addsd  %xmm8,%xmm0
  subsd  %xmm13,%xmm0
  mulsd  %xmm7,%xmm0
  divsd  %xmm12,%xmm0
  mulsd  %xmm0,%xmm5
  movapd %xmm4,%xmm0
  subsd  %xmm5,%xmm6
  mulsd  %xmm6,%xmm0
  addsd  %xmm0,%xmm1
  comisd %xmm1,%xmm11
  jbe    L2
L100:
  movsd  D50(%rip),%xmm0
  movsd  %xmm1,(%r12)
  subsd  %xmm1,%xmm0
  addsd  D50(%rip),%xmm0
  movsd  %xmm0,0x0(%r13)
  jmp    L98
  nopl   0x0(%rax)
L21:
  mov    0xc8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L127
  add    $0xd8,%rsp
  mov    %r12,%rcx
  mov    %r9,%r8
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    10b0 <cdflib_cumf@plt>
  nopl   (%rax)
L196:
  movsd  D101(%rip),%xmm4
  mov    $0x2,%ebp
  mov    $0x1,%ebx
  movsd  D77(%rip),%xmm2
  movapd %xmm4,%xmm12
  movsd  %xmm2,0x8(%rsp)
  movapd %xmm4,%xmm0
  movsd  %xmm4,0x58(%rsp)
  jmp    L154
  nopl   0x0(%rax)
L145:
  lea    0xc0(%rsp),%rdi
  xor    %eax,%eax
  movsd  %xmm2,0x50(%rsp)
  movsd  %xmm12,0x48(%rsp)
  movsd  %xmm1,0x40(%rsp)
  movsd  %xmm9,0x38(%rsp)
  movsd  %xmm4,0x30(%rsp)
  movsd  %xmm6,0x28(%rsp)
  movsd  %xmm0,0xc0(%rsp)
  call   <cdflib_gamma_log@plt>
  lea    0x80(%rsp),%rdi
  mov    %eax,%ebx
  xor    %eax,%eax
  call   <cdflib_gamma_log@plt>
  mov    %r14,%rdi
  sub    %eax,%ebx
  xor    %eax,%eax
  call   <cdflib_gamma_log@plt>
  movsd  0x70(%rsp),%xmm7
  movsd  0x80(%rsp),%xmm3
  mov    %eax,%r14d
  movapd %xmm7,%xmm0
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm7,0x20(%rsp)
  call   <log@plt>
  movsd  %xmm0,0x8(%rsp)
  movsd  0x78(%rsp),%xmm0
  call   <log@plt>
  movsd  0x88(%rsp),%xmm8
  jmp    L206
  nopw   0x0(%rax,%rax,1)
L134:
  movapd %xmm2,%xmm3
  movsd  %xmm5,0x78(%rsp)
  subsd  %xmm5,%xmm3
  jmp    L92
L127:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: label: D77
datatype: f64
name: unpack_data
arguments: label: D150
datatype: byte[8]
name: unpack_data
arguments: label: D50
datatype: f64[2]
name: unpack_data
arguments: label: D101
datatype: f64
name: unpack_data
arguments: label: D76
datatype: byte[8]
name: unpack_data
arguments: label: D78
datatype: byte[8]

tool: D77: 1

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

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

tool: D101: 2

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

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

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;
}
```