user: 请从给定的汇编指令推导出C源代码:

```asm
<cdflib_beta_asym>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %rdx,%rbp
  push   %rbx
  mov    %rcx,%rbx
  sub    $0x388,%rsp
  movsd  (%rdi),%xmm1
  movsd  (%rsi),%xmm2
  mov    %rdi,0x88(%rsp)
  mov    %rsi,0x80(%rsp)
  movapd %xmm2,%xmm3
  subsd  %xmm1,%xmm3
  mov    %fs:0x28,%rax
  mov    %rax,0x378(%rsp)
  xor    %eax,%eax
  comisd %xmm2,%xmm1
  jb     L11
  divsd  %xmm1,%xmm3
  movapd %xmm2,%xmm6
  pxor   %xmm7,%xmm7
  movsd  D131(%rip),%xmm12
  divsd  %xmm1,%xmm6
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm6,0x38(%rsp)
  addsd  %xmm12,%xmm6
  movapd %xmm6,%xmm0
  movapd %xmm12,%xmm6
  divsd  %xmm0,%xmm6
  mulsd  %xmm2,%xmm0
  ucomisd %xmm0,%xmm7
  movsd  %xmm6,0x28(%rsp)
  ja     L161
L61:
  sqrtsd %xmm0,%xmm0
L139:
  movapd %xmm12,%xmm6
  lea    0xb0(%rsp),%rdi
  xor    %eax,%eax
  movsd  %xmm7,0x20(%rsp)
  divsd  %xmm0,%xmm6
  movsd  0x0(%rbp),%xmm0
  movsd  %xmm1,(%rsp)
  movsd  %xmm12,0x18(%rsp)
  movapd %xmm0,%xmm3
  divsd  %xmm2,%xmm0
  movsd  %xmm6,0x40(%rsp)
  divsd  %xmm1,%xmm3
  movsd  %xmm0,0xb8(%rsp)
  xorpd  D218(%rip),%xmm3
  movsd  %xmm3,0xb0(%rsp)
  call   <cdflib_rlog1@plt>
  pxor   %xmm0,%xmm0
  movsd  (%rsp),%xmm1
  lea    0xb8(%rsp),%rdi
  cvtsi2sd %eax,%xmm0
  mov    0x80(%rsp),%rax
  movsd  (%rax),%xmm6
  xor    %eax,%eax
  mulsd  %xmm1,%xmm0
  movsd  %xmm6,0x8(%rsp)
  movsd  %xmm0,(%rsp)
  call   <cdflib_rlog1@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mulsd  0x8(%rsp),%xmm1
  addsd  (%rsp),%xmm1
  movapd %xmm1,%xmm0
  xorpd  D218(%rip),%xmm0
  movsd  %xmm1,(%rsp)
  call   <exp@plt>
  movsd  0x20(%rsp),%xmm7
  movsd  (%rsp),%xmm1
  movsd  0x18(%rsp),%xmm12
  movsd  %xmm0,0x78(%rsp)
  ucomisd %xmm7,%xmm0
  jp     L70
  movapd %xmm7,%xmm0
  jne    L70
L180:
  mov    0x378(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L163
  add    $0x388,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L11:
  divsd  %xmm2,%xmm3
  movapd %xmm1,%xmm6
  pxor   %xmm7,%xmm7
  movsd  D131(%rip),%xmm12
  divsd  %xmm2,%xmm6
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm6,0x38(%rsp)
  addsd  %xmm12,%xmm6
  movapd %xmm6,%xmm0
  movapd %xmm12,%xmm6
  divsd  %xmm0,%xmm6
  mulsd  %xmm1,%xmm0
  ucomisd %xmm0,%xmm7
  movsd  %xmm6,0x28(%rsp)
  jbe    L61
L161:
  movsd  %xmm7,0x20(%rsp)
  movsd  %xmm12,0x18(%rsp)
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm1,(%rsp)
  call   <sqrt@plt>
  movsd  0x20(%rsp),%xmm7
  movsd  0x8(%rsp),%xmm2
  movsd  0x18(%rsp),%xmm12
  movsd  (%rsp),%xmm1
  jmp    L139
  nopl   0x0(%rax)
L70:
  ucomisd %xmm1,%xmm7
  ja     L85
  movapd %xmm1,%xmm0
  sqrtsd %xmm0,%xmm0
L52:
  movsd  0x10(%rsp),%xmm8
  movsd  D151(%rip),%xmm6
  addsd  %xmm1,%xmm1
  xor    %eax,%eax
  mulsd  D192(%rip),%xmm8
  lea    0xa8(%rsp),%rsi
  movsd  %xmm0,0xa8(%rsp)
  divsd  D228(%rip),%xmm0
  mulsd  D151(%rip),%xmm0
  lea    D206(%rip),%rdi
  movsd  %xmm7,0x98(%rsp)
  lea    0x2e0(%rsp),%r13
  movsd  %xmm12,0x18(%rsp)
  lea    0x170(%rsp),%r15
  lea    0x220(%rsp),%rbp
  movsd  %xmm1,0x48(%rsp)
  movsd  %xmm8,0xc0(%rsp)
  mulsd  %xmm8,%xmm6
  movsd  %xmm0,0x30(%rsp)
  movsd  %xmm8,0x90(%rsp)
  movapd %xmm6,%xmm5
  xorpd  D218(%rip),%xmm5
  movsd  %xmm6,0x2d0(%rsp)
  movsd  %xmm6,0x8(%rsp)
  movsd  %xmm5,(%rsp)
  movsd  %xmm5,0x220(%rsp)
  call   <cdflib_error_fc@plt>
  pxor   %xmm0,%xmm0
  movsd  D180(%rip),%xmm5
  movsd  0x40(%rsp),%xmm11
  cvtsi2sd %eax,%xmm0
  movsd  0x8(%rsp),%xmm6
  movsd  0x38(%rsp),%xmm4
  lea    0xc0(%rsp),%r8
  movsd  0x10(%rsp),%xmm7
  mov    D228(%rip),%rax
  movsd  0x18(%rsp),%xmm12
  movsd  0x48(%rsp),%xmm15
  movsd  0x90(%rsp),%xmm8
  addsd  %xmm7,%xmm7
  mulsd  %xmm0,%xmm5
  movapd %xmm11,%xmm0
  movapd %xmm12,%xmm14
  mulsd  %xmm6,%xmm0
  movapd %xmm12,%xmm13
  mulsd  D228(%rip),%xmm0
  movsd  %xmm7,0x58(%rsp)
  movsd  (%rbx),%xmm7
  mov    $0x2,%ebx
  movsd  %xmm7,0x68(%rsp)
  movsd  0x98(%rsp),%xmm7
  movsd  %xmm5,0x20(%rsp)
  addsd  %xmm5,%xmm0
  movapd %xmm4,%xmm5
  mulsd  %xmm4,%xmm5
  movsd  %xmm0,0x8(%rsp)
  movsd  %xmm5,0x60(%rsp)
  movsd  0x28(%rsp),%xmm5
  mov    %rax,0x28(%rsp)
  lea    0x180(%rsp),%rax
  mov    %rax,0x18(%rsp)
  lea    0x228(%rsp),%rax
  addsd  %xmm5,%xmm5
  mov    %rax,0x10(%rsp)
  lea    0x2d0(%rsp),%rax
  mov    %rax,0x70(%rsp)
  movsd  %xmm5,0x50(%rsp)
L45:
  pxor   %xmm10,%xmm10
  lea    0x1(%rbx),%r10
  lea    0x3(%rbx),%r14
  mov    %r13,%r11
  mulsd  0x60(%rsp),%xmm13
  movsd  0x38(%rsp),%xmm1
  cvtsi2sd %ebx,%xmm10
  movsd  D201(%rip),%xmm0
  addsd  %xmm10,%xmm0
  mulsd  %xmm13,%xmm1
  addsd  %xmm13,%xmm14
  addsd  %xmm12,%xmm1
  mulsd  0x50(%rsp),%xmm1
  divsd  %xmm0,%xmm1
  movsd  D178(%rip),%xmm0
  addsd  %xmm10,%xmm0
  movsd  %xmm1,-0x8(%r8,%rbx,8)
  movsd  0x58(%rsp),%xmm1
  mulsd  %xmm14,%xmm1
  divsd  %xmm0,%xmm1
  movsd  %xmm1,(%r8,%rbx,8)
L33:
  lea    -0x1(%r10),%eax
  pxor   %xmm9,%xmm9
  movapd %xmm8,%xmm0
  movsd  D151(%rip),%xmm5
  cvtsi2sd %eax,%xmm9
  mov    0x18(%rsp),%r9
  mov    %r10d,%r12d
  mov    $0x1,%edi
  mov    $0x2,%esi
  addsd  %xmm12,%xmm9
  mulsd  %xmm9,%xmm5
  xorpd  D218(%rip),%xmm5
  mulsd  %xmm5,%xmm0
  movsd  %xmm0,0x170(%rsp)
  nopl   0x0(%rax,%rax,1)
L175:
  mov    %r9,%rdx
  movapd %xmm7,%xmm3
  movapd %xmm0,%xmm4
  mov    $0x1,%eax
  movapd %xmm8,%xmm1
  jmp    L86
  nopw   0x0(%rax,%rax,1)
L99:
  movsd  -0x10(%rdx),%xmm4
  movsd  (%r8,%rax,8),%xmm1
  add    $0x1,%rax
L86:
  pxor   %xmm0,%xmm0
  mov    %esi,%ecx
  pxor   %xmm2,%xmm2
  sub    $0x8,%rdx
  cvtsi2sd %eax,%xmm0
  sub    %eax,%ecx
  cvtsi2sd %ecx,%xmm2
  mulsd  %xmm5,%xmm0
  subsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  %xmm4,%xmm0
  addsd  %xmm0,%xmm3
  cmp    %rax,%rdi
  jne    L99
  pxor   %xmm1,%xmm1
  movsd  (%r8,%rdi,8),%xmm0
  add    $0x8,%r9
  cvtsi2sd %esi,%xmm1
  add    $0x1,%esi
  mulsd  %xmm5,%xmm0
  divsd  %xmm1,%xmm3
  addsd  %xmm3,%xmm0
  movsd  %xmm0,-0x10(%r9)
  cmp    %r10d,%esi
  je     L200
  movsd  -0x10(%r9),%xmm0
  add    $0x1,%rdi
  jmp    L175
  nopl   0x0(%rax,%rax,1)
L200:
  movsd  -0x10(%r15,%r10,8),%xmm2
  movsd  (%rsp),%xmm0
  lea    -0x3(%r12),%edx
  mov    %rbp,%rax
  mov    0x10(%rsp),%rdi
  movapd %xmm7,%xmm1
  divsd  %xmm9,%xmm2
  movapd %xmm0,%xmm3
  movapd %xmm6,%xmm0
  lea    (%rdi,%rdx,8),%rcx
  mov    %r11,%rdx
  movsd  %xmm2,-0x10(%rbp,%r10,8)
  jmp    L174
  nopw   0x0(%rax,%rax,1)
L42:
  movsd  -0x10(%rdx),%xmm0
  movsd  (%rax),%xmm3
L174:
  mulsd  %xmm3,%xmm0
  add    $0x8,%rax
  sub    $0x8,%rdx
  addsd  %xmm0,%xmm1
  cmp    %rax,%rcx
  jne    L42
  addsd  %xmm1,%xmm2
  xorpd  D218(%rip),%xmm2
  add    $0x1,%r10
  add    $0x8,%r11
  movsd  %xmm2,-0x10(%r11)
  cmp    %r10,%r14
  je     L159
  movsd  -0x10(%r11),%xmm6
  jmp    L33
L159:
  movsd  0x30(%rsp),%xmm6
  movapd %xmm10,%xmm0
  mov    0x70(%rsp),%rax
  movsd  D228(%rip),%xmm1
  movsd  0x48(%rsp),%xmm2
  subsd  %xmm12,%xmm0
  mulsd  0x20(%rsp),%xmm0
  mulsd  %xmm6,%xmm1
  mulsd  %xmm2,%xmm6
  mulsd  0x28(%rsp),%xmm10
  movapd %xmm1,%xmm5
  addsd  %xmm0,%xmm5
  movsd  D228(%rip),%xmm0
  movsd  %xmm6,0x30(%rsp)
  movsd  0x40(%rsp),%xmm6
  mulsd  %xmm15,%xmm0
  mulsd  %xmm6,%xmm11
  movsd  %xmm5,0x20(%rsp)
  mulsd  %xmm2,%xmm15
  movapd %xmm0,%xmm4
  movsd  -0x8(%rax,%rbx,8),%xmm0
  addsd  %xmm10,%xmm4
  mulsd  %xmm11,%xmm0
  mulsd  %xmm6,%xmm11
  movsd  0x0(%r13),%xmm6
  movsd  %xmm4,0x28(%rsp)
  movapd %xmm6,%xmm1
  mulsd  %xmm5,%xmm0
  movsd  0x8(%rsp),%xmm5
  mulsd  %xmm11,%xmm1
  movapd %xmm0,%xmm2
  andpd  D92(%rip),%xmm0
  mulsd  %xmm4,%xmm1
  addsd  %xmm1,%xmm2
  andpd  D92(%rip),%xmm1
  addsd  %xmm1,%xmm0
  addsd  %xmm2,%xmm5
  movsd  %xmm5,0x8(%rsp)
  mulsd  0x68(%rsp),%xmm5
  comisd %xmm0,%xmm5
  jae    L80
  add    $0x2,%rbx
  cmp    $0x16,%rbx
  je     L80
  add    $0x10,%r13
  jmp    L45
L80:
  mov    0x80(%rsp),%rsi
  mov    0x88(%rsp),%rdi
  xor    %eax,%eax
  call   <cdflib_bcorr@plt>
  pxor   %xmm0,%xmm0
  neg    %eax
  cvtsi2sd %eax,%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm1
  movsd  0x78(%rsp),%xmm0
  mulsd  D91(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  0x8(%rsp),%xmm0
  jmp    L180
L85:
  movapd %xmm1,%xmm0
  movsd  %xmm7,0x18(%rsp)
  movsd  %xmm12,0x8(%rsp)
  movsd  %xmm1,(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm7
  movsd  (%rsp),%xmm1
  movsd  0x8(%rsp),%xmm12
  jmp    L52
L163:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: label: D131
data_type: f64
name: unpack_data
arguments: label: D218
data_type: f32[4]
name: unpack_data
arguments: label: D151
data_type: f64
name: unpack_data
arguments: label: D192
data_type: byte[8]
name: unpack_data
arguments: label: D228
data_type: byte[8]
name: unpack_data
arguments: label: D206
data_type: string
name: unpack_data
arguments: label: D180
data_type: byte[8]
name: unpack_data
arguments: label: D201
data_type: f64
name: unpack_data
arguments: label: D178
data_type: f64
name: unpack_data
arguments: label: D92
data_type: qword[2]
name: unpack_data
arguments: label: D91
data_type: f64

tool: D131: 1

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

tool: D151: 0.5

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

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

tool: D206: 

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

tool: D201: 2

tool: D178: 3

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

tool: D91: 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;
}
```