user: Translate the below assembly code into C programming language:

```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     L15
  comisd %xmm0,%xmm3
  ja     L15
  movsd  D47(%rip),%xmm0
  mulsd  %xmm5,%xmm0
  comisd %xmm4,%xmm0
  jbe    L139
L15:
  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 D48(%rip),%xmm6
  setnp  %al
  cmovne %edx,%eax
  test   %al,%al
  jne    L150
  cmp    $0x1,%ecx
  je     L150
  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    L203
  movsd  (%r14),%xmm2
  comisd D29(%rip),%xmm2
  jbe    L188
  subsd  %xmm7,%xmm0
  mulsd  (%r12),%xmm0
  divsd  %xmm2,%xmm0
  subsd  %xmm1,%xmm0
  comisd %xmm0,%xmm7
  ja     L203
  pxor   %xmm1,%xmm1
  cvtsi2sd %esi,%xmm1
  comisd %xmm0,%xmm1
  jbe    L188
  cvttsd2si %xmm0,%edx
  test   %edx,%edx
  jg     L229
  movapd %xmm3,%xmm1
  nopl   0x0(%rax)
L8:
  cmp    %edx,%esi
  je     L33
  lea    0x1(%rdx),%eax
L204:
  cmp    %eax,%esi
  jl     L33
  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     L0
  jmp    L33
  nop
L156:
  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    L33
L0:
  mov    %eax,%edx
  add    $0x1,%eax
  cmp    %eax,%ecx
  jne    L156
L33:
  mulsd  %xmm1,%xmm6
L150:
  mov    0x28(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L26
  add    $0x38,%rsp
  movapd %xmm6,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L139:
  lea    D86(%rip),%rdi
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm4,0x8(%rsp)
  call   <cdflib_exparg@plt>
  mov    %rax,%rdi
  call   <fabs@plt>
  lea    D62(%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    L135
  mov    %eax,0x24(%rsp)
  mov    %eax,%edx
L135:
  pxor   %xmm0,%xmm0
  movsd  %xmm5,0x10(%rsp)
  cvtsi2sd %edx,%xmm0
  xorpd  D65(%rip),%xmm0
  movsd  %xmm4,0x8(%rsp)
  call   <exp@plt>
  movsd  0x10(%rsp),%xmm5
  movsd  0x8(%rsp),%xmm4
  movapd %xmm0,%xmm3
  jmp    L15
L188:
  mov    %esi,%edx
  movapd %xmm3,%xmm1
  test   %esi,%esi
  jle    L33
L229:
  movsd  (%r12),%xmm7
  movapd %xmm3,%xmm1
  xor    %eax,%eax
  nopw   0x0(%rax,%rax,1)
L71:
  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    L71
  jmp    L8
L203:
  movapd %xmm3,%xmm1
  mov    $0x1,%eax
  xor    %edx,%edx
  jmp    L204
L26:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_label: D60
type: f64
name: read_data
arguments: data_label: D47
type: f64
name: read_data
arguments: data_label: D48
type: f64
name: read_data
arguments: data_label: D29
type: byte[16]
name: read_data
arguments: data_label: D86
type: string
name: read_data
arguments: data_label: D62
type: string
name: read_data
arguments: data_label: D65
type: byte[16]

tool: D60: 1

tool: D47: 1.1

tool: D48: 0

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

tool: D86: 

tool: 

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