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

```asm
<cdflib_beta_asym>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x368,%rsp
  mov    %rdi,0x78(%rsp)
  mov    %rsi,0x70(%rsp)
  mov    %rdx,%rbp
  mov    %rcx,0x48(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x358(%rsp)
  xor    %eax,%eax
  movsd  (%rdi),%xmm5
  movsd  %xmm5,0x8(%rsp)
  mov    (%rsi),%rbx
  movq   %rbx,%xmm7
  comisd %xmm7,%xmm5
  jb     L168
  divsd  %xmm5,%xmm7
  movsd  %xmm7,0x30(%rsp)
  movsd  D208(%rip),%xmm1
  addsd  %xmm1,%xmm7
  movapd %xmm7,%xmm0
  divsd  %xmm7,%xmm1
  movsd  %xmm1,0x28(%rsp)
  movq   %rbx,%xmm1
  subsd  %xmm5,%xmm1
  divsd  %xmm5,%xmm1
  movsd  %xmm1,0x20(%rsp)
  movq   %rbx,%xmm5
  mulsd  %xmm5,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L93
  sqrtsd %xmm0,%xmm0
L52:
  movsd  D208(%rip),%xmm1
  divsd  %xmm0,%xmm1
  movsd  %xmm1,0x38(%rsp)
L183:
  movsd  0x0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  0x8(%rsp),%xmm1
  xorpd  D233(%rip),%xmm1
  movsd  %xmm1,0x90(%rsp)
  movq   %rbx,%xmm5
  divsd  %xmm5,%xmm0
  movsd  %xmm0,0x98(%rsp)
  lea    0x90(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_rlog1@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  0x8(%rsp),%xmm0
  movsd  %xmm0,0x8(%rsp)
  mov    0x70(%rsp),%rax
  mov    (%rax),%rbx
  lea    0x98(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_rlog1@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movq   %rbx,%xmm7
  mulsd  %xmm7,%xmm0
  addsd  0x8(%rsp),%xmm0
  movsd  %xmm0,0x8(%rsp)
  xorpd  D233(%rip),%xmm0
  call   <exp@plt>
  movsd  %xmm0,0x68(%rsp)
  ucomisd D67(%rip),%xmm0
  jp     L249
  pxor   %xmm0,%xmm0
  jne    L249
L232:
  mov    0x358(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L39
  add    $0x368,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L93:
  call   <sqrt@plt>
  jmp    L52
L168:
  movsd  0x8(%rsp),%xmm7
  movapd %xmm7,%xmm5
  movq   %rbx,%xmm3
  divsd  %xmm3,%xmm5
  movsd  %xmm5,0x30(%rsp)
  movsd  D208(%rip),%xmm1
  movapd %xmm5,%xmm0
  addsd  %xmm1,%xmm0
  divsd  %xmm0,%xmm1
  movsd  %xmm1,0x28(%rsp)
  movq   %rbx,%xmm1
  subsd  %xmm7,%xmm1
  divsd  %xmm3,%xmm1
  movsd  %xmm1,0x20(%rsp)
  mulsd  %xmm7,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L129
  sqrtsd %xmm0,%xmm0
L184:
  movsd  D208(%rip),%xmm1
  divsd  %xmm0,%xmm1
  movsd  %xmm1,0x38(%rsp)
  jmp    L183
L129:
  call   <sqrt@plt>
  jmp    L184
L249:
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm5
  ucomisd %xmm5,%xmm0
  ja     L72
  sqrtsd %xmm5,%xmm5
  movapd %xmm5,%xmm0
L128:
  movsd  %xmm0,0x88(%rsp)
  divsd  D108(%rip),%xmm0
  movsd  D90(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,0x18(%rsp)
  movsd  0x8(%rsp),%xmm5
  addsd  %xmm5,%xmm5
  movsd  %xmm5,0x40(%rsp)
  movsd  0x20(%rsp),%xmm0
  mulsd  D193(%rip),%xmm0
  movsd  %xmm0,0xa0(%rsp)
  movapd %xmm0,%xmm7
  mulsd  %xmm1,%xmm7
  movq   %xmm7,%rbx
  xorpd  D233(%rip),%xmm7
  movsd  %xmm7,0x200(%rsp)
  mov    %rbx,0x2b0(%rsp)
  lea    0x88(%rsp),%rsi
  lea    D256(%rip),%rdi
  mov    $0x0,%eax
  call   <cdflib_error_fc@plt>
  pxor   %xmm7,%xmm7
  cvtsi2sd %eax,%xmm7
  mulsd  D36(%rip),%xmm7
  movq   %rbx,%xmm3
  movsd  0x38(%rsp),%xmm8
  mulsd  %xmm8,%xmm3
  movapd %xmm3,%xmm0
  mulsd  D108(%rip),%xmm0
  movapd %xmm0,%xmm3
  addsd  %xmm7,%xmm3
  movsd  %xmm3,0x10(%rsp)
  movsd  0x30(%rsp),%xmm6
  movapd %xmm6,%xmm3
  mulsd  %xmm6,%xmm3
  movsd  %xmm3,0x60(%rsp)
  movsd  0x28(%rsp),%xmm3
  addsd  %xmm3,%xmm3
  movsd  %xmm3,0x50(%rsp)
  movsd  0x20(%rsp),%xmm5
  addsd  %xmm5,%xmm5
  movsd  %xmm5,0x58(%rsp)
  lea    0xa8(%rsp),%rax
  mov    %rax,0x20(%rsp)
  lea    0x2b0(%rsp),%r15
  movsd  0x40(%rsp),%xmm11
  mov    $0x2,%r13d
  movsd  D208(%rip),%xmm10
  mov    D108(%rip),%rax
  movq   %rax,%xmm5
  movapd %xmm10,%xmm9
  movq   D233(%rip),%xmm15
  lea    0xa0(%rsp),%rdi
  lea    0x208(%rsp),%rax
  mov    %rax,0x28(%rsp)
  movapd %xmm15,%xmm13
  movq   D188(%rip),%xmm12
  jmp    L225
L72:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L128
L227:
  mov    %r9,%rdx
  mov    $0x1,%eax
  pxor   %xmm1,%xmm1
L206:
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  %xmm2,%xmm0
  mov    %ecx,%r10d
  sub    %eax,%r10d
  pxor   %xmm4,%xmm4
  cvtsi2sd %r10d,%xmm4
  subsd  %xmm4,%xmm0
  mulsd  -0x8(%rdi,%rax,8),%xmm0
  mulsd  -0x10(%rdx),%xmm0
  addsd  %xmm0,%xmm1
  add    $0x1,%rax
  sub    $0x8,%rdx
  cmp    %rsi,%rax
  jne    L206
L119:
  movapd %xmm2,%xmm0
  mulsd  -0x8(%rdi,%rsi,8),%xmm0
  pxor   %xmm4,%xmm4
  cvtsi2sd %r8d,%xmm4
  divsd  %xmm4,%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%r9)
  add    $0x1,%rsi
  add    $0x8,%r9
  cmp    %r11,%rsi
  je     L148
L85:
  mov    %esi,%ecx
  mov    %esi,%r8d
  pxor   %xmm1,%xmm1
  cmp    $0x1,%esi
  jg     L227
  jmp    L119
L117:
  movapd %xmm6,%xmm1
  subsd  D208(%rip),%xmm1
  mulsd  %xmm7,%xmm1
  movsd  D108(%rip),%xmm0
  movsd  0x18(%rsp),%xmm3
  movapd %xmm3,%xmm7
  mulsd  %xmm0,%xmm7
  addsd  %xmm1,%xmm7
  mulsd  %xmm6,%xmm5
  mulsd  %xmm11,%xmm0
  addsd  %xmm0,%xmm5
  movsd  0x40(%rsp),%xmm6
  mulsd  %xmm6,%xmm3
  movsd  %xmm3,0x18(%rsp)
  mulsd  %xmm6,%xmm11
  movsd  0x38(%rsp),%xmm6
  mulsd  %xmm6,%xmm8
  movapd %xmm8,%xmm0
  mulsd  0x8(%r15),%xmm0
  mulsd  %xmm7,%xmm0
  mulsd  %xmm6,%xmm8
  movapd %xmm8,%xmm1
  mulsd  0x10(%r15),%xmm1
  mulsd  %xmm5,%xmm1
  movapd %xmm0,%xmm2
  addsd  %xmm1,%xmm2
  movapd %xmm2,%xmm3
  addsd  0x10(%rsp),%xmm3
  movsd  %xmm3,0x10(%rsp)
  andpd  %xmm12,%xmm0
  andpd  %xmm12,%xmm1
  addsd  %xmm1,%xmm0
  mov    0x48(%rsp),%rax
  mulsd  (%rax),%xmm3
  comisd %xmm0,%xmm3
  jae    L208
  add    $0x2,%r13d
  addq   $0x10,0x20(%rsp)
  add    $0x10,%r15
  cmp    $0x16,%r13d
  je     L208
L225:
  mulsd  0x60(%rsp),%xmm9
  pxor   %xmm6,%xmm6
  cvtsi2sd %r13d,%xmm6
  movsd  0x30(%rsp),%xmm1
  mulsd  %xmm9,%xmm1
  addsd  D208(%rip),%xmm1
  mulsd  0x50(%rsp),%xmm1
  movapd %xmm6,%xmm0
  addsd  D29(%rip),%xmm0
  divsd  %xmm0,%xmm1
  mov    0x20(%rsp),%rax
  movsd  %xmm1,(%rax)
  addsd  %xmm9,%xmm10
  movsd  0x58(%rsp),%xmm1
  mulsd  %xmm10,%xmm1
  movapd %xmm6,%xmm0
  addsd  D27(%rip),%xmm0
  divsd  %xmm0,%xmm1
  movsd  %xmm1,0x8(%rax)
  movsd  0xa0(%rsp),%xmm14
  mov    %r13d,0x8(%rsp)
  lea    -0x2(%r13),%ebp
  mov    %r15,%r12
  mov    %r13d,%ebx
  jmp    L12
L208:
  mov    0x70(%rsp),%rsi
  mov    0x78(%rsp),%rdi
  mov    $0x0,%eax
  call   <cdflib_bcorr@plt>
  neg    %eax
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm1
  movsd  0x68(%rsp),%xmm0
  mulsd  D87(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  0x10(%rsp),%xmm0
  jmp    L232
L148:
  lea    -0x1(%rbx),%esi
  movslq %esi,%rax
  movsd  0x150(%rsp,%rax,8),%xmm2
  divsd  %xmm3,%xmm2
  movsd  %xmm2,0x200(%rsp,%rax,8)
  lea    0x200(%rsp),%rax
  mov    0x28(%rsp),%rcx
  lea    (%rcx,%r14,8),%rcx
  mov    %r12,%rdx
  pxor   %xmm1,%xmm1
L170:
  movsd  (%rdx),%xmm0
  mulsd  (%rax),%xmm0
  addsd  %xmm0,%xmm1
  sub    $0x8,%rdx
  add    $0x8,%rax
  cmp    %rcx,%rax
  jne    L170
L132:
  movslq %esi,%rsi
  addsd  %xmm1,%xmm2
  xorpd  %xmm13,%xmm2
  movsd  %xmm2,0x2b0(%rsp,%rsi,8)
  add    $0x1,%ebx
  add    $0x8,%r12
  add    $0x1,%ebp
  cmp    0x8(%rsp),%ebp
  je     L117
L12:
  pxor   %xmm3,%xmm3
  cvtsi2sd %ebx,%xmm3
  addsd  D208(%rip),%xmm3
  movsd  D90(%rip),%xmm2
  mulsd  %xmm3,%xmm2
  xorpd  %xmm15,%xmm2
  movapd %xmm14,%xmm0
  mulsd  %xmm2,%xmm0
  movsd  %xmm0,0x150(%rsp)
  cmp    $0x1,%ebx
  jle    L149
  lea    0x160(%rsp),%r9
  mov    %ebp,%r14d
  lea    0x3(%r14),%r11
  mov    $0x2,%esi
  jmp    L85
L149:
  lea    -0x1(%rbx),%esi
  movslq %esi,%rax
  movsd  0x150(%rsp,%rax,8),%xmm2
  divsd  %xmm3,%xmm2
  movsd  %xmm2,0x200(%rsp,%rax,8)
  pxor   %xmm1,%xmm1
  jmp    L132
L39:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: tag: D208
type_infered: f64
name: retrieve_data
arguments: tag: D233
type_infered: f32[4]
name: retrieve_data
arguments: tag: D67
type_infered: f32[2]
name: retrieve_data
arguments: tag: D108
type_infered: byte[8]
name: retrieve_data
arguments: tag: D90
type_infered: f64
name: retrieve_data
arguments: tag: D193
type_infered: byte[8]
name: retrieve_data
arguments: tag: D256
type_infered: string
name: retrieve_data
arguments: tag: D36
type_infered: byte[8]
name: retrieve_data
arguments: tag: D188
type_infered: qword
name: retrieve_data
arguments: tag: D29
type_infered: f64
name: retrieve_data
arguments: tag: D27
type_infered: f64
name: retrieve_data
arguments: tag: D87
type_infered: f64

tool: D208: 1

tool: {"D233": [0.0,-0.0,0.0,0.0]}

tool: {"D67": [0.0,0.0]}

tool: {"D108": [0xd1,0x3b,0x7f,0x66,0x9e,0xa0,0xd6,0x3f]}

tool: D90: 0.5

tool: {"D193": [0x55,0x55,0x55,0x55,0x55,0x55,0xe5,0x3f]}

tool: D256: 

tool: {"D36": [0x7e,0xef,0xb4,0x91,0xf8,0x5b,0xdc,0x3f]}

tool: D188: 9223372036854775807

tool: D29: 2

tool: D27: 3

tool: D87: 1.12837916709551

assistant: ```c
double cdflib_beta_asym(const double *a, const double *b, const double *lambda, const double *eps) {
  static const double e0 = 1.12837916709551e0;
  static const double e1 = .353553390593274e0;
  static const int num = 20;
  static const int K3 = 1;
  double value;
  double bsum, dsum, f, h, h2, hn, j0, j1, r, r0, r1, s, sum, t, t0, t1, u, w, w0, z, z0, z2, zn, znm1;
  int i, im1, imj, j, m, mm1, mmj, n, np1;
  double a0[21], b0[21], c[21], d[21], T1, T2;
  value = 0.0e0;
  if (*a >= *b) {
    goto S10;
  }
  h = *a / *b;
  r0 = 1.0e0 / (1.0e0 + h);
  r1 = (*b - *a) / *b;
  w0 = 1.0e0 / sqrt(*a * (1.0e0 + h));
  goto S20;
S10:
  h = *b / *a;
  r0 = 1.0e0 / (1.0e0 + h);
  r1 = (*b - *a) / *a;
  w0 = 1.0e0 / sqrt(*b * (1.0e0 + h));
S20:
  T1 = -(*lambda / *a);
  T2 = *lambda / *b;
  f = *a * cdflib_rlog1(&T1) + *b * cdflib_rlog1(&T2);
  t = exp(-f);
  if (t == 0.0e0) {
    return value;
  }
  z0 = sqrt(f);
  z = 0.5e0 * (z0 / e1);
  z2 = f + f;
  a0[0] = 2.0e0 / 3.0e0 * r1;
  c[0] = -(0.5e0 * a0[0]);
  d[0] = -c[0];
  j0 = 0.5e0 / e0 * cdflib_error_fc(&K3, &z0);
  j1 = e1;
  sum = j0 + d[0] * w0 * j1;
  s = 1.0e0;
  h2 = h * h;
  hn = 1.0e0;
  w = w0;
  znm1 = z;
  zn = z2;
  for (n = 2; n <= num; n += 2) {
    hn = h2 * hn;
    a0[n - 1] = 2.0e0 * r0 * (1.0e0 + h * hn) / ((double)n + 2.0e0);
    np1 = n + 1;
    s += hn;
    a0[np1 - 1] = 2.0e0 * r1 * s / ((double)n + 3.0e0);
    for (i = n; i <= np1; i++) {
      r = -(0.5e0 * ((double)i + 1.0e0));
      b0[0] = r * a0[0];
      for (m = 2; m <= i; m++) {
        bsum = 0.0e0;
        mm1 = m - 1;
        for (j = 1; j <= mm1; j++) {
          mmj = m - j;
          bsum += (((double)j * r - (double)mmj) * a0[j - 1] * b0[mmj - 1]);
        }
        b0[m - 1] = r * a0[m - 1] + bsum / (double)m;
      }
      c[i - 1] = b0[i - 1] / ((double)i + 1.0e0);
      dsum = 0.0e0;
      im1 = i - 1;
      for (j = 1; j <= im1; j++) {
        imj = i - j;
        dsum += (d[imj - 1] * c[j - 1]);
      }
      d[i - 1] = -(dsum + c[i - 1]);
    }
    j0 = e1 * znm1 + ((double)n - 1.0e0) * j0;
    j1 = e1 * zn + (double)n * j1;
    znm1 = z2 * znm1;
    zn = z2 * zn;
    w = w0 * w;
    t0 = d[n - 1] * w * j0;
    w = w0 * w;
    t1 = d[np1 - 1] * w * j1;
    sum += (t0 + t1);
    if (fabs(t0) + fabs(t1) <= *eps * sum) {
      goto S80;
    }
  }
S80:
  u = exp(-cdflib_bcorr(a, b));
  value = e0 * t * u * sum;
  return value;
}
```