user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<cdflib_beta_up>:
  endbr64
  push   %r15
  mov    %r8,%r15
  push   %r14
  mov    %rcx,%r14
  push   %r13
  mov    %rsi,%r13
  push   %r12
  mov    %rdx,%r12
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  mov    %r9,%rbx
  sub    $0x38,%rsp
  movsd  (%rdi),%xmm0
  movsd  (%rsi),%xmm4
  mov    %fs:0x28,%rax
  mov    %rax,0x28(%rsp)
  xor    %eax,%eax
  movsd  D60(%rip),%xmm5
  cmpl   $0x1,(%r8)
  movl   $0x0,0x24(%rsp)
  mov    D60(%rip),%rax
  addsd  %xmm0,%xmm4
  addsd  %xmm0,%xmm5
  movq   %rax,%xmm3
  je     L70
  comisd %xmm0,%xmm3
  ja     L70
  movsd  D66(%rip),%xmm0
  mulsd  %xmm5,%xmm0
  comisd %xmm4,%xmm0
  jbe    L82
L70:
  mov    %r12,%rcx
  mov    %r13,%rdx
  lea    0x24(%rsp),%rdi
  mov    %r14,%r8
  mov    %rbp,%rsi
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm4,0x8(%rsp)
  call   <cdflib_beta_rcomp1@plt>
  movsd  0x0(%rbp),%xmm1
  mov    $0x0,%edx
  mov    (%r15),%ecx
  movapd %xmm0,%xmm6
  divsd  %xmm1,%xmm6
  ucomisd D7(%rip),%xmm6
  setnp  %al
  cmovne %edx,%eax
  test   %al,%al
  jne    L202
  cmp    $0x1,%ecx
  je     L202
  movsd  0x0(%r13),%xmm0
  movsd  0x8(%rsp),%xmm4
  lea    -0x1(%rcx),%esi
  movsd  D60(%rip),%xmm7
  movsd  0x10(%rsp),%xmm5
  movsd  0x18(%rsp),%xmm3
  comisd %xmm0,%xmm7
  jae    L115
  movsd  (%r14),%xmm2
  comisd D45(%rip),%xmm2
  jbe    L205
  subsd  %xmm7,%xmm0
  mulsd  (%r12),%xmm0
  divsd  %xmm2,%xmm0
  subsd  %xmm1,%xmm0
  comisd %xmm0,%xmm7
  ja     L115
  pxor   %xmm1,%xmm1
  cvtsi2sd %esi,%xmm1
  comisd %xmm0,%xmm1
  jbe    L205
  cvttsd2si %xmm0,%edx
  test   %edx,%edx
  jg     L74
  movapd %xmm3,%xmm1
  nopl   0x0(%rax)
L75:
  cmp    %edx,%esi
  je     L42
  lea    0x1(%rdx),%eax
L36:
  cmp    %eax,%esi
  jl     L42
  pxor   %xmm0,%xmm0
  movapd %xmm4,%xmm2
  movsd  (%r12),%xmm7
  movsd  (%rbx),%xmm8
  cvtsi2sd %edx,%xmm0
  addsd  %xmm0,%xmm2
  addsd  %xmm5,%xmm0
  divsd  %xmm0,%xmm2
  mulsd  %xmm7,%xmm2
  mulsd  %xmm3,%xmm2
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  mulsd  %xmm8,%xmm0
  comisd %xmm2,%xmm0
  jb     L145
  jmp    L42
  nop
L13:
  pxor   %xmm3,%xmm3
  movapd %xmm4,%xmm0
  cvtsi2sd %edx,%xmm3
  addsd  %xmm3,%xmm0
  addsd  %xmm5,%xmm3
  divsd  %xmm3,%xmm0
  mulsd  %xmm7,%xmm0
  mulsd  %xmm0,%xmm2
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  mulsd  %xmm8,%xmm0
  comisd %xmm2,%xmm0
  jae    L42
L145:
  mov    %eax,%edx
  add    $0x1,%eax
  cmp    %eax,%ecx
  jne    L13
L42:
  mulsd  %xmm1,%xmm6
L202:
  mov    0x28(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L89
  add    $0x38,%rsp
  movapd %xmm6,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L82:
  lea    D59(%rip),%rdi
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm4,0x8(%rsp)
  call   <cdflib_exparg@plt>
  mov    %rax,%rdi
  call   <fabs@plt>
  lea    D24(%rip),%rdi
  mov    %eax,0x24(%rsp)
  call   <cdflib_exparg@plt>
  mov    0x24(%rsp),%edx
  movsd  0x8(%rsp),%xmm4
  movsd  0x10(%rsp),%xmm5
  cmp    %eax,%edx
  jle    L2
  mov    %eax,0x24(%rsp)
  mov    %eax,%edx
L2:
  pxor   %xmm0,%xmm0
  movsd  %xmm5,0x10(%rsp)
  cvtsi2sd %edx,%xmm0
  xorpd  D69(%rip),%xmm0
  movsd  %xmm4,0x8(%rsp)
  call   <exp@plt>
  movsd  0x10(%rsp),%xmm5
  movsd  0x8(%rsp),%xmm4
  movapd %xmm0,%xmm3
  jmp    L70
L205:
  mov    %esi,%edx
  movapd %xmm3,%xmm1
  test   %esi,%esi
  jle    L42
L74:
  movsd  (%r12),%xmm7
  movapd %xmm3,%xmm1
  xor    %eax,%eax
  nopw   0x0(%rax,%rax,1)
L26:
  pxor   %xmm2,%xmm2
  movapd %xmm4,%xmm0
  cvtsi2sd %eax,%xmm2
  add    $0x1,%eax
  addsd  %xmm2,%xmm0
  addsd  %xmm5,%xmm2
  divsd  %xmm2,%xmm0
  mulsd  %xmm7,%xmm0
  mulsd  %xmm0,%xmm3
  addsd  %xmm3,%xmm1
  cmp    %eax,%edx
  jne    L26
  jmp    L75
L115:
  movapd %xmm3,%xmm1
  mov    $0x1,%eax
  xor    %edx,%edx
  jmp    L36
L89:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: label: D60
datatype: f64
name: parse_data
arguments: label: D66
datatype: f64
name: parse_data
arguments: label: D7
datatype: f64
name: parse_data
arguments: label: D45
datatype: byte[16]
name: parse_data
arguments: label: D59
datatype: string
name: parse_data
arguments: label: D24
datatype: string
name: parse_data
arguments: label: D69
datatype: byte[16]

tool: D60: 1

tool: D66: 1.1

tool: D7: 0

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

tool: D59: 

tool: 

tool: {"D69": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

assistant: ```c
double cdflib_beta_up(const double *a, const double *b, const double *x, const double *y, const int *n, const double *eps) {
  static const int K1 = 1;
  static const int K2 = 0;
  double bup, ap1, apb, d, l, r, t, w;
  int i, k, kp1, mu, nm1;
  apb = *a + *b;
  ap1 = *a + 1.0e0;
  mu = 0;
  d = 1.0e0;
  if (*n == 1 || *a < 1.0e0) {
    goto S10;
  }
  if (apb < 1.1e0 * ap1) {
    goto S10;
  }
  mu = (int)fabs(cdflib_exparg(&K1));
  k = (int)cdflib_exparg(&K2);
  if (k < mu) {
    mu = k;
  }
  t = mu;
  d = exp(-t);
S10:
  bup = cdflib_beta_rcomp1(&mu, a, b, x, y) / *a;
  if (*n == 1 || bup == 0.0e0) {
    return bup;
  }
  nm1 = *n - 1;
  w = d;
  k = 0;
  if (*b <= 1.0e0) {
    goto S50;
  }
  if (*y > 1.e-4) {
    goto S20;
  }
  k = nm1;
  goto S30;
S20:
  r = (*b - 1.0e0) * *x / *y - *a;
  if (r < 1.0e0) {
    goto S50;
  }
  t = (double)nm1;
  k = nm1;
  if (r < t) {
    k = (int)r;
  }
S30:
  for (i = 1; i <= k; i++) {
    l = i - 1;
    d = (apb + l) / (ap1 + l) * *x * d;
    w = w + d;
  }
  if (k == nm1) {
    goto S70;
  }
S50:
  kp1 = k + 1;
  for (i = kp1; i <= nm1; i++) {
    l = i - 1;
    d = (apb + l) / (ap1 + l) * *x * d;
    w = w + d;
    if (d <= *eps * w) {
      goto S70;
    }
  }
S70:
  bup *= w;
  return bup;
}
```