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

```asm
<cdflib_beta_up>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rdi,%r12
  mov    %rsi,%r13
  mov    %rdx,%rbp
  mov    %rcx,%r14
  mov    %r8,%r15
  mov    %r9,%rbx
  mov    %fs:0x28,%rax
  mov    %rax,0x28(%rsp)
  xor    %eax,%eax
  movsd  (%rdi),%xmm0
  movapd %xmm0,%xmm7
  addsd  (%rsi),%xmm7
  movsd  %xmm7,0x8(%rsp)
  movapd %xmm0,%xmm4
  addsd  D88(%rip),%xmm4
  movsd  %xmm4,0x10(%rsp)
  movl   $0x0,0x24(%rsp)
  cmpl   $0x1,(%r8)
  je     L212
  movsd  D88(%rip),%xmm1
  comisd %xmm0,%xmm1
  ja     L219
  movsd  %xmm1,0x18(%rsp)
  mulsd  D36(%rip),%xmm4
  comisd %xmm7,%xmm4
  jbe    L150
L39:
  lea    0x24(%rsp),%rdi
  mov    %r14,%r8
  mov    %rbp,%rcx
  mov    %r13,%rdx
  mov    %r12,%rsi
  call   <cdflib_beta_rcomp1@plt>
  movsd  (%r12),%xmm1
  movapd %xmm0,%xmm4
  divsd  %xmm1,%xmm4
  mov    (%r15),%ecx
  ucomisd D117(%rip),%xmm4
  setnp  %al
  mov    $0x0,%edx
  cmovne %edx,%eax
  test   %al,%al
  jne    L125
  cmp    $0x1,%ecx
  je     L125
  lea    -0x1(%rcx),%esi
  movsd  0x0(%r13),%xmm0
  movsd  D88(%rip),%xmm2
  comisd %xmm0,%xmm2
  jae    L119
  movsd  (%r14),%xmm2
  comisd D35(%rip),%xmm2
  jbe    L132
  movsd  D88(%rip),%xmm3
  subsd  %xmm3,%xmm0
  mulsd  0x0(%rbp),%xmm0
  divsd  %xmm2,%xmm0
  subsd  %xmm1,%xmm0
  comisd %xmm0,%xmm3
  ja     L94
  pxor   %xmm1,%xmm1
  cvtsi2sd %esi,%xmm1
  comisd %xmm0,%xmm1
  jbe    L132
  cvttsd2si %xmm0,%edi
  test   %edi,%edi
  jg     L90
  movsd  0x18(%rsp),%xmm1
  jmp    L71
L150:
  lea    D2(%rip),%rdi
  call   <cdflib_exparg@plt>
  mov    %rax,%rdi
  call   <fabs@plt>
  mov    %eax,0x24(%rsp)
  lea    D25(%rip),%rdi
  call   <cdflib_exparg@plt>
  cmp    %eax,0x24(%rsp)
  jle    L133
  mov    %eax,0x24(%rsp)
L133:
  pxor   %xmm0,%xmm0
  cvtsi2sdl 0x24(%rsp),%xmm0
  xorpd  D71(%rip),%xmm0
  call   <exp@plt>
  movsd  %xmm0,0x18(%rsp)
  jmp    L39
L212:
  movsd  D88(%rip),%xmm6
  movsd  %xmm6,0x18(%rsp)
  jmp    L39
L219:
  movsd  D88(%rip),%xmm6
  movsd  %xmm6,0x18(%rsp)
  jmp    L39
L119:
  movsd  0x18(%rsp),%xmm1
  mov    $0x0,%edi
  jmp    L175
L94:
  movsd  0x18(%rsp),%xmm1
  mov    $0x0,%edi
  jmp    L175
L132:
  mov    %esi,%edi
  movsd  0x18(%rsp),%xmm1
  test   %esi,%esi
  jle    L165
L90:
  movsd  0x0(%rbp),%xmm3
  mov    %edi,%edx
  movsd  0x18(%rsp),%xmm1
  mov    $0x0,%eax
L22:
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
  movsd  0x8(%rsp),%xmm0
  addsd  %xmm2,%xmm0
  addsd  0x10(%rsp),%xmm2
  divsd  %xmm2,%xmm0
  mulsd  %xmm3,%xmm0
  mulsd  0x18(%rsp),%xmm0
  movsd  %xmm0,0x18(%rsp)
  addsd  %xmm0,%xmm1
  add    $0x1,%eax
  cmp    %eax,%edx
  jne    L22
L71:
  cmp    %edi,%esi
  je     L165
L175:
  lea    0x1(%rdi),%eax
  cmp    %eax,%esi
  jl     L165
  pxor   %xmm0,%xmm0
  cvtsi2sd %edi,%xmm0
  movsd  0x0(%rbp),%xmm6
  movsd  0x8(%rsp),%xmm2
  addsd  %xmm0,%xmm2
  addsd  0x10(%rsp),%xmm0
  divsd  %xmm0,%xmm2
  mulsd  %xmm6,%xmm2
  mulsd  0x18(%rsp),%xmm2
  addsd  %xmm2,%xmm1
  movsd  (%rbx),%xmm5
  movapd %xmm1,%xmm0
  mulsd  %xmm5,%xmm0
  comisd %xmm2,%xmm0
  jae    L165
L184:
  mov    %eax,%edx
  add    $0x1,%eax
  cmp    %eax,%ecx
  je     L165
  pxor   %xmm3,%xmm3
  cvtsi2sd %edx,%xmm3
  movsd  0x8(%rsp),%xmm0
  addsd  %xmm3,%xmm0
  addsd  0x10(%rsp),%xmm3
  divsd  %xmm3,%xmm0
  mulsd  %xmm6,%xmm0
  mulsd  %xmm0,%xmm2
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  mulsd  %xmm5,%xmm0
  comisd %xmm2,%xmm0
  jb     L184
L165:
  mulsd  %xmm1,%xmm4
L125:
  mov    0x28(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L234
  movapd %xmm4,%xmm0
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L234:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D88
type_infered: f64
name: extract_symbolic_data
arguments: data_label: D36
type_infered: f64
name: extract_symbolic_data
arguments: data_label: D117
type_infered: f64
name: extract_symbolic_data
arguments: data_label: D35
type_infered: byte[16]
name: extract_symbolic_data
arguments: data_label: D2
type_infered: string
name: extract_symbolic_data
arguments: data_label: D25
type_infered: string
name: extract_symbolic_data
arguments: data_label: D71
type_infered: byte[16]

tool: D88: 1

tool: D36: 1.1

tool: D117: 0

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

tool: D2: 

tool: 

tool: {"D71": [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;
}
```