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

```asm
<betaln>:
  endbr64
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  movapd %xmm0,%xmm5
  movapd %xmm0,%xmm4
  minsd  %xmm1,%xmm4
  movsd  %xmm4,0x8(%rsp)
  maxsd  %xmm1,%xmm5
  movsd  %xmm5,0x10(%rsp)
  comisd D10(%rip),%xmm4
  jae    L148
  movsd  0x8(%rsp),%xmm7
  comisd D213(%rip),%xmm7
  jae    L89
  movsd  0x10(%rsp),%xmm2
  comisd D10(%rip),%xmm2
  jae    L169
  movsd  0x8(%rsp),%xmm0
  call   <gamln@plt>
  mov    %eax,%ebp
  movsd  0x10(%rsp),%xmm0
  call   <gamln@plt>
  mov    %eax,%ebx
  movsd  0x8(%rsp),%xmm2
  addsd  0x10(%rsp),%xmm2
  movapd %xmm2,%xmm0
  call   <gamln@plt>
  sub    %eax,%ebx
  add    %ebp,%ebx
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebx,%xmm0
L143:
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  ret
L148:
  movapd %xmm5,%xmm1
  movapd %xmm4,%xmm0
  call   <bcorr@plt>
  movq   %xmm0,%rbp
  movsd  0x8(%rsp),%xmm2
  divsd  0x10(%rsp),%xmm2
  movq   %xmm2,%rbx
  movapd %xmm2,%xmm1
  addsd  D213(%rip),%xmm1
  movapd %xmm2,%xmm3
  divsd  %xmm1,%xmm3
  movapd %xmm3,%xmm0
  call   <log@plt>
  movapd %xmm0,%xmm1
  movsd  0x8(%rsp),%xmm0
  subsd  D103(%rip),%xmm0
  xorpd  D56(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,0x8(%rsp)
  movq   %rbx,%xmm0
  call   <alnrel@plt>
  movapd %xmm0,%xmm2
  mulsd  0x10(%rsp),%xmm2
  movsd  %xmm2,0x18(%rsp)
  comisd 0x8(%rsp),%xmm2
  jae    L249
  movsd  0x10(%rsp),%xmm0
  call   <log@plt>
  mulsd  D60(%rip),%xmm0
  addsd  D65(%rip),%xmm0
  movq   %rbp,%xmm7
  addsd  %xmm7,%xmm0
  subsd  0x18(%rsp),%xmm0
  subsd  0x8(%rsp),%xmm0
  jmp    L143
L89:
  comisd D64(%rip),%xmm7
  ja     L124
  movsd  0x10(%rsp),%xmm5
  comisd D64(%rip),%xmm5
  jbe    L128
  mov    D26(%rip),%rbx
  movsd  D10(%rip),%xmm0
  comisd %xmm5,%xmm0
  jbe    L152
L159:
  movsd  0x10(%rsp),%xmm0
  subsd  D213(%rip),%xmm0
  cvttsd2si %xmm0,%edx
  test   %edx,%edx
  jle    L154
  add    $0x1,%edx
  movsd  D213(%rip),%xmm0
  mov    $0x1,%eax
  movapd %xmm0,%xmm3
L234:
  movsd  0x10(%rsp),%xmm6
  subsd  %xmm3,%xmm6
  movsd  %xmm6,0x10(%rsp)
  movsd  0x8(%rsp),%xmm2
  addsd  %xmm6,%xmm2
  divsd  %xmm2,%xmm6
  mulsd  %xmm6,%xmm0
  add    $0x1,%eax
  cmp    %eax,%edx
  jne    L234
L182:
  call   <log@plt>
  movq   %rbx,%xmm3
  addsd  %xmm0,%xmm3
  movsd  %xmm3,0x18(%rsp)
  movsd  0x8(%rsp),%xmm0
  call   <gamln@plt>
  mov    %eax,%ebx
  movsd  0x10(%rsp),%xmm0
  call   <gamln@plt>
  mov    %eax,%ebp
  movsd  0x10(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm0
  call   <gsumln@plt>
  movapd %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebp,%xmm0
  subsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sd %ebx,%xmm1
  addsd  %xmm1,%xmm0
  addsd  0x18(%rsp),%xmm0
  jmp    L143
L169:
  movsd  0x8(%rsp),%xmm0
  call   <gamln@plt>
  mov    %eax,%ebx
  movsd  0x10(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm0
  call   <algdiv@plt>
  add    %eax,%ebx
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebx,%xmm0
  jmp    L143
L124:
  movsd  0x10(%rsp),%xmm3
  comisd D77(%rip),%xmm3
  ja     L77
  movsd  0x8(%rsp),%xmm0
  subsd  D213(%rip),%xmm0
  cvttsd2si %xmm0,%edx
  test   %edx,%edx
  jle    L153
  add    $0x1,%edx
  movsd  D213(%rip),%xmm1
  mov    $0x1,%eax
  movapd %xmm1,%xmm2
L13:
  movsd  0x8(%rsp),%xmm7
  subsd  %xmm2,%xmm7
  movsd  %xmm7,0x8(%rsp)
  movapd %xmm7,%xmm0
  divsd  0x10(%rsp),%xmm0
  movapd %xmm0,%xmm3
  addsd  %xmm2,%xmm3
  divsd  %xmm3,%xmm0
  mulsd  %xmm0,%xmm1
  add    $0x1,%eax
  cmp    %edx,%eax
  jne    L13
L63:
  movapd %xmm1,%xmm0
  call   <log@plt>
  movq   %xmm0,%rbx
  movsd  D10(%rip),%xmm0
  comisd 0x10(%rsp),%xmm0
  ja     L159
  movsd  0x8(%rsp),%xmm0
  call   <gamln@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movq   %rbx,%xmm7
  addsd  %xmm0,%xmm7
  movsd  %xmm7,0x18(%rsp)
  movsd  0x10(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm0
  call   <algdiv@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  addsd  0x18(%rsp),%xmm0
  jmp    L143
L128:
  movsd  0x8(%rsp),%xmm0
  call   <gamln@plt>
  mov    %eax,%ebx
  movsd  0x10(%rsp),%xmm0
  call   <gamln@plt>
  mov    %eax,%ebp
  movsd  0x10(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm0
  call   <gsumln@plt>
  movapd %xmm0,%xmm1
  add    %ebp,%ebx
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebx,%xmm0
  subsd  %xmm1,%xmm0
  jmp    L143
L152:
  movsd  0x8(%rsp),%xmm0
  call   <gamln@plt>
  mov    %eax,%ebx
  movsd  0x10(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm0
  call   <algdiv@plt>
  add    %eax,%ebx
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebx,%xmm0
  jmp    L143
L77:
  movsd  0x8(%rsp),%xmm0
  subsd  D213(%rip),%xmm0
  cvttsd2si %xmm0,%ebp
  movsd  D213(%rip),%xmm1
  test   %ebp,%ebp
  jle    L14
  lea    0x1(%rbp),%edx
  movsd  D213(%rip),%xmm1
  mov    $0x1,%eax
  movapd %xmm1,%xmm2
L40:
  movsd  0x8(%rsp),%xmm3
  subsd  %xmm2,%xmm3
  movsd  %xmm3,0x8(%rsp)
  movapd %xmm3,%xmm0
  divsd  0x10(%rsp),%xmm0
  addsd  %xmm2,%xmm0
  divsd  %xmm0,%xmm3
  mulsd  %xmm3,%xmm1
  add    $0x1,%eax
  cmp    %eax,%edx
  jne    L40
L14:
  movapd %xmm1,%xmm0
  call   <log@plt>
  movq   %xmm0,%rbx
  movsd  0x10(%rsp),%xmm0
  call   <log@plt>
  movapd %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebp,%xmm0
  mulsd  %xmm1,%xmm0
  movq   %rbx,%xmm4
  subsd  %xmm0,%xmm4
  movsd  %xmm4,0x18(%rsp)
  movsd  0x8(%rsp),%xmm0
  call   <gamln@plt>
  mov    %eax,%ebx
  movsd  0x10(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm0
  call   <algdiv@plt>
  add    %eax,%ebx
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebx,%xmm0
  addsd  0x18(%rsp),%xmm0
  jmp    L143
L153:
  movsd  D213(%rip),%xmm1
  jmp    L63
L154:
  movsd  D213(%rip),%xmm0
  jmp    L182
L249:
  movsd  0x10(%rsp),%xmm0
  call   <log@plt>
  mulsd  D60(%rip),%xmm0
  addsd  D65(%rip),%xmm0
  movq   %rbp,%xmm4
  addsd  %xmm4,%xmm0
  subsd  0x8(%rsp),%xmm0
  subsd  0x18(%rsp),%xmm0
  jmp    L143
```

tool_calls: name: disassemble_data
arguments: data_label: D10
type_infered: f64[2]
name: disassemble_data
arguments: data_label: D213
type_infered: f64[2]
name: disassemble_data
arguments: data_label: D103
type_infered: f64
name: disassemble_data
arguments: data_label: D56
type_infered: f32[4]
name: disassemble_data
arguments: data_label: D60
type_infered: f32[2]
name: disassemble_data
arguments: data_label: D65
type_infered: byte[8]
name: disassemble_data
arguments: data_label: D64
type_infered: f64[2]
name: disassemble_data
arguments: data_label: D26
type_infered: f32[2]
name: disassemble_data
arguments: data_label: D77
type_infered: f32[4]

tool: {"D10": [8.0,0.5]}

tool: {"D213": [1.0,8.0]}

tool: D103: 0.5

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

tool: {"D60": [0.0,-1.75]}

tool: {"D65": [0xb7,0xbe,0x64,0xc8,0xf1,0x67,0xed,0x3f]}

tool: {"D64": [2.0,1000.0]}

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

tool: {"D77": [0.0,4.4765625,0.0,-1.75]}

assistant: ```c
double betaln(double a0, double b0) {
  static double e = .918938533204673;
  int i__1;
  double ret_val;
  double a, b, c__, h__;
  int i__, n;
  double u, v, w, z__;
  a = ((a0 < b0) ? a0 : b0);
  b = ((a0 > b0) ? a0 : b0);
  if (a >= 8.0) {
    goto L60;
  }
  if (a >= 1.0) {
    goto L20;
  }
  if (b >= 8.0) {
    goto L10;
  }
  ret_val = gamln(a) + (gamln(b) - gamln(a + b));
  return ret_val;
L10:
  ret_val = gamln(a) + algdiv(a, b);
  return ret_val;
L20:
  if (a > 2.0) {
    goto L30;
  }
  if (b > 2.0) {
    goto L21;
  }
  ret_val = gamln(a) + gamln(b) - gsumln(a, b);
  return ret_val;
L21:
  w = 0.0;
  if (b < 8.0) {
    goto L40;
  }
  ret_val = gamln(a) + algdiv(a, b);
  return ret_val;
L30:
  if (b > 1e3f) {
    goto L50;
  }
  n = a - 1.0;
  w = 1.0;
  i__1 = n;
  for (i__ = 1; i__ <= i__1; ++i__) {
    a += -1.0;
    h__ = a / b;
    w *= h__ / (h__ + 1.0);
  }
  w = log(w);
  if (b < 8.0) {
    goto L40;
  }
  ret_val = w + gamln(a) + algdiv(a, b);
  return ret_val;
L40:
  n = b - 1.0;
  z__ = 1.0;
  i__1 = n;
  for (i__ = 1; i__ <= i__1; ++i__) {
    b += -1.0;
    z__ *= b / (a + b);
  }
  ret_val = w + log(z__) + (gamln(a) + (gamln(b) - gsumln(a, b)));
  return ret_val;
L50:
  n = a - 1.0;
  w = 1.0;
  i__1 = n;
  for (i__ = 1; i__ <= i__1; ++i__) {
    a += -1.0;
    w *= a / (a / b + 1.0);
  }
  ret_val = log(w) - n * log(b) + (gamln(a) + algdiv(a, b));
  return ret_val;
L60:
  w = bcorr(a, b);
  h__ = a / b;
  c__ = h__ / (h__ + 1.0);
  u = -(a - 0.5) * log(c__);
  v = b * alnrel(h__);
  if (u <= v) {
    goto L61;
  }
  ret_val = log(b) * -0.5 + e + w - v - u;
  return ret_val;
L61:
  ret_val = log(b) * -0.5 + e + w - u - v;
  return ret_val;
}
```