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

```asm
<cdflib_beta_frac>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xd0,%rsp
  mov    %rdi,-0xa8(%rbp)
  mov    %rsi,-0xb0(%rbp)
  mov    %rdx,-0xb8(%rbp)
  mov    %rcx,-0xc0(%rbp)
  mov    %r8,-0xc8(%rbp)
  mov    %r9,-0xd0(%rbp)
  mov    -0xc0(%rbp),%rcx
  mov    -0xb8(%rbp),%rdx
  mov    -0xb0(%rbp),%rsi
  mov    -0xa8(%rbp),%rax
  mov    %rax,%rdi
  call   <cdflib_beta_cdflib_rcomp@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x58(%rbp)
  pxor   %xmm0,%xmm0
  ucomisd -0x58(%rbp),%xmm0
  jp     L14
  pxor   %xmm0,%xmm0
  ucomisd -0x58(%rbp),%xmm0
  jne    L14
  movsd  -0x58(%rbp),%xmm0
  jmp    L7
L14:
  mov    -0xc8(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D104(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0xb0(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0xa8(%rbp),%rax
  movsd  (%rax),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    -0xa8(%rbp),%rax
  movsd  (%rax),%xmm2
  movsd  D104(%rip),%xmm0
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  D104(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0xc0(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D104(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movsd  D104(%rip),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  mov    -0xa8(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D104(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  D104(%rip),%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0x50(%rbp),%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0x40(%rbp),%xmm0
  divsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x68(%rbp)
L3:
  movsd  -0x78(%rbp),%xmm1
  movsd  D104(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  mov    -0xa8(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  -0x78(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0xb0(%rbp),%rax
  movsd  (%rax),%xmm0
  subsd  -0x78(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x78(%rbp),%xmm1
  mov    -0xb8(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0xa8(%rbp),%rax
  movsd  (%rax),%xmm0
  divsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x70(%rbp),%xmm0
  addsd  -0x48(%rbp),%xmm0
  mulsd  -0x70(%rbp),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x20(%rbp),%xmm1
  mov    -0xb8(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  -0x28(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x30(%rbp),%xmm1
  movsd  D104(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x40(%rbp),%xmm1
  addsd  -0x30(%rbp),%xmm1
  addsd  -0x30(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x28(%rbp),%xmm0
  divsd  -0x60(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x78(%rbp),%xmm1
  movsd  -0x78(%rbp),%xmm0
  mulsd  -0x38(%rbp),%xmm0
  addsd  -0x50(%rbp),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x30(%rbp),%xmm1
  movsd  D104(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  -0x60(%rbp),%xmm1
  movsd  D98(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0x18(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x98(%rbp),%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  -0x90(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0x18(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x88(%rbp),%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  -0x80(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0x68(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x90(%rbp),%xmm0
  divsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0x68(%rbp),%xmm0
  subsd  -0x8(%rbp),%xmm0
  movq   D55(%rip),%xmm1
  andpd  %xmm0,%xmm1
  mov    -0xd0(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  -0x68(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jae    L47
  movsd  -0x98(%rbp),%xmm0
  divsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0x88(%rbp),%xmm0
  divsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0x68(%rbp),%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  D104(%rip),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  jmp    L3
L47:
  nop
  movsd  -0x58(%rbp),%xmm0
  mulsd  -0x68(%rbp),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0x58(%rbp),%xmm0
L7:
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: unpack_data
arguments: label: D104
type_guessed: f64
name: unpack_data
arguments: label: D98
type_guessed: f64
name: unpack_data
arguments: label: D55
type_guessed: qword

tool: D104: 1

tool: D98: 2

tool: D55: 9223372036854775807

assistant: ```c
double cdflib_beta_frac(const double *a, const double *b, const double *x, const double *y, const double *lambda, const double *eps) {
  double bfrac, alpha, an, anp1, beta, bn, bnp1, c, c0, c1, e, n, p, r, r0, s, t, w, yp1;
  bfrac = cdflib_beta_cdflib_rcomp(a, b, x, y);
  if (bfrac == 0.0e0) {
    return bfrac;
  }
  c = 1.0e0 + *lambda;
  c0 = *b / *a;
  c1 = 1.0e0 + 1.0e0 / *a;
  yp1 = *y + 1.0e0;
  n = 0.0e0;
  p = 1.0e0;
  s = *a + 1.0e0;
  an = 0.0e0;
  bn = anp1 = 1.0e0;
  bnp1 = c / c1;
  r = c1 / c;
S10:
  n = n + 1.0e0;
  t = n / *a;
  w = n * (*b - n) * *x;
  e = *a / s;
  alpha = p * (p + c0) * e * e * (w * *x);
  e = (1.0e0 + t) / (c1 + t + t);
  beta = n + w / s + e * (c + n * yp1);
  p = 1.0e0 + t;
  s += 2.0e0;
  t = alpha * an + beta * anp1;
  an = anp1;
  anp1 = t;
  t = alpha * bn + beta * bnp1;
  bn = bnp1;
  bnp1 = t;
  r0 = r;
  r = anp1 / bnp1;
  if (fabs(r - r0) <= (*eps) * r) {
    goto S20;
  }
  an /= bnp1;
  bn /= bnp1;
  anp1 = r;
  bnp1 = 1.0e0;
  goto S10;
S20:
  bfrac = bfrac * r;
  return bfrac;
}
```