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

```asm
<cdflib_cdft>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xb0,%rsp
  mov    %rdi,-0x78(%rbp)
  mov    %rsi,-0x80(%rbp)
  mov    %rdx,-0x88(%rbp)
  mov    %rcx,-0x90(%rbp)
  mov    %r8,-0x98(%rbp)
  mov    %r9,-0xa0(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0xa8(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  D279(%rip),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  D67(%rip),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  mov    -0xa0(%rbp),%rax
  movl   $0x0,(%rax)
  mov    -0xa8(%rbp),%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x78(%rbp),%rax
  mov    (%rax),%eax
  test   %eax,%eax
  jle    L394
  mov    -0x78(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x3,%eax
  jle    L437
L394:
  mov    -0x78(%rbp),%rax
  mov    (%rax),%eax
  test   %eax,%eax
  jg     L399
  mov    -0xa8(%rbp),%rax
  movsd  D148(%rip),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L285
L399:
  nop
  mov    -0xa8(%rbp),%rax
  movsd  D116(%rip),%xmm0
  movsd  %xmm0,(%rax)
L285:
  mov    -0xa0(%rbp),%rax
  movl   $0xffffffff,(%rax)
  jmp    L238
L437:
  nop
  mov    -0x78(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x1,%eax
  je     L363
  mov    -0x80(%rbp),%rax
  movsd  (%rax),%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  setae  %al
  xor    $0x1,%eax
  test   %al,%al
  je     L358
  mov    -0x80(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  D148(%rip),%xmm1
  comisd %xmm1,%xmm0
  seta   %al
  xor    $0x1,%eax
  test   %al,%al
  jne    L161
L358:
  mov    -0x80(%rbp),%rax
  movsd  (%rax),%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  setae  %al
  xor    $0x1,%eax
  test   %al,%al
  jne    L274
  mov    -0xa8(%rbp),%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  jmp    L468
L274:
  nop
  mov    -0xa8(%rbp),%rax
  movsd  D148(%rip),%xmm0
  movsd  %xmm0,(%rax)
L468:
  mov    -0xa0(%rbp),%rax
  movl   $0xfffffffe,(%rax)
  jmp    L238
L363:
  nop
  jmp    L357
L161:
  nop
L357:
  mov    -0x78(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x1,%eax
  je     L259
  mov    -0x88(%rbp),%rax
  movsd  (%rax),%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  setae  %al
  xor    $0x1,%eax
  test   %al,%al
  je     L33
  mov    -0x88(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  D148(%rip),%xmm1
  comisd %xmm1,%xmm0
  seta   %al
  xor    $0x1,%eax
  test   %al,%al
  jne    L135
L33:
  mov    -0x88(%rbp),%rax
  movsd  (%rax),%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  setae  %al
  xor    $0x1,%eax
  test   %al,%al
  jne    L235
  mov    -0xa8(%rbp),%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  jmp    L83
L235:
  nop
  mov    -0xa8(%rbp),%rax
  movsd  D148(%rip),%xmm0
  movsd  %xmm0,(%rax)
L83:
  mov    -0xa0(%rbp),%rax
  movl   $0xfffffffd,(%rax)
  jmp    L238
L259:
  nop
  jmp    L360
L135:
  nop
L360:
  mov    -0x78(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x3,%eax
  je     L308
  mov    -0x98(%rbp),%rax
  movsd  (%rax),%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  setae  %al
  xor    $0x1,%eax
  test   %al,%al
  jne    L175
  mov    -0xa8(%rbp),%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xa0(%rbp),%rax
  movl   $0xfffffffb,(%rax)
  jmp    L238
L308:
  nop
  jmp    L29
L175:
  nop
L29:
  mov    -0x78(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x1,%eax
  je     L326
  mov    -0x80(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x88(%rbp),%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  movsd  D279(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  D279(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movq   D106(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0xb0(%rbp)
  mov    D165(%rip),%eax
  mov    %eax,%edi
  mov    $0x0,%eax
  call   <dpmpar@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D116(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0xb0(%rbp),%xmm2
  comisd %xmm0,%xmm2
  seta   %al
  xor    $0x1,%eax
  test   %al,%al
  jne    L200
  pxor   %xmm0,%xmm0
  comisd -0x10(%rbp),%xmm0
  seta   %al
  xor    $0x1,%eax
  test   %al,%al
  jne    L499
  mov    -0xa8(%rbp),%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  jmp    L319
L499:
  nop
  mov    -0xa8(%rbp),%rax
  movsd  D148(%rip),%xmm0
  movsd  %xmm0,(%rax)
L319:
  mov    -0xa0(%rbp),%rax
  movl   $0x3,(%rax)
  jmp    L238
L326:
  nop
  jmp    L121
L200:
  nop
L121:
  mov    -0x78(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x1,%eax
  je     L207
  mov    -0x80(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x88(%rbp),%rax
  movsd  (%rax),%xmm0
  comisd %xmm1,%xmm0
  setae  %al
  movzbl %al,%eax
  mov    %rax,-0x18(%rbp)
L207:
  mov    -0x78(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x1,%eax
  jne    L26
  mov    -0x88(%rbp),%rcx
  mov    -0x80(%rbp),%rdx
  mov    -0x98(%rbp),%rsi
  mov    -0x90(%rbp),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <cdflib_cumt@plt>
  mov    -0xa0(%rbp),%rax
  movl   $0x0,(%rax)
  jmp    L380
L26:
  mov    -0x78(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x2,%eax
  jne    L141
  mov    -0x98(%rbp),%rdx
  mov    -0x88(%rbp),%rcx
  mov    -0x80(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <cdflib_dt1@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mov    -0x90(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  D265(%rip),%xmm0
  movq   D62(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  D265(%rip),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  D160(%rip),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  D270(%rip),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  lea    -0x28(%rbp),%r9
  lea    -0x68(%rbp),%r8
  lea    -0x70(%rbp),%rcx
  lea    -0x70(%rbp),%rdx
  lea    -0x30(%rbp),%rsi
  lea    -0x38(%rbp),%rax
  sub    $0x8,%rsp
  lea    -0x20(%rbp),%rdi
  push   %rdi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <cdflib_dstinv@plt>
  add    $0x10,%rsp
  mov    -0xa0(%rbp),%rax
  movl   $0x0,(%rax)
  lea    -0x48(%rbp),%rdi
  lea    -0x40(%rbp),%rcx
  lea    -0x60(%rbp),%rdx
  mov    -0x90(%rbp),%rsi
  mov    -0xa0(%rbp),%rax
  mov    %rdi,%r8
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <cdflib_dinvr@plt>
L503:
  mov    -0xa0(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x1,%eax
  jne    L106
  lea    -0x50(%rbp),%rcx
  lea    -0x58(%rbp),%rdx
  mov    -0x98(%rbp),%rsi
  mov    -0x90(%rbp),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <cdflib_cumt@plt>
  cmpq   $0x0,-0x18(%rbp)
  je     L222
  movsd  -0x58(%rbp),%xmm0
  mov    -0x80(%rbp),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  jmp    L177
L222:
  nop
  movsd  -0x50(%rbp),%xmm0
  mov    -0x88(%rbp),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
L177:
  lea    -0x48(%rbp),%rdi
  lea    -0x40(%rbp),%rcx
  lea    -0x60(%rbp),%rdx
  mov    -0x90(%rbp),%rsi
  mov    -0xa0(%rbp),%rax
  mov    %rdi,%r8
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <cdflib_dinvr@plt>
  jmp    L503
L106:
  nop
  mov    -0xa0(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0xffffffff,%eax
  jne    L409
  mov    -0x40(%rbp),%rax
  test   %rax,%rax
  je     L389
  mov    -0xa0(%rbp),%rax
  movl   $0x1,(%rax)
  movsd  D265(%rip),%xmm0
  movq   D62(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mov    -0xa8(%rbp),%rax
  movsd  %xmm0,(%rax)
  jmp    L380
L389:
  nop
  mov    -0xa0(%rbp),%rax
  movl   $0x2,(%rax)
  movsd  D265(%rip),%xmm0
  mov    -0xa8(%rbp),%rax
  movsd  %xmm0,(%rax)
  jmp    L380
L141:
  mov    -0x78(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x3,%eax
  jne    L67
  mov    -0x98(%rbp),%rax
  movsd  D67(%rip),%xmm0
  movsd  %xmm0,(%rax)
  movsd  D50(%rip),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  D15(%rip),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  D160(%rip),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  D270(%rip),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  lea    -0x28(%rbp),%r9
  lea    -0x68(%rbp),%r8
  lea    -0x70(%rbp),%rcx
  lea    -0x70(%rbp),%rdx
  lea    -0x30(%rbp),%rsi
  lea    -0x38(%rbp),%rax
  sub    $0x8,%rsp
  lea    -0x20(%rbp),%rdi
  push   %rdi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <cdflib_dstinv@plt>
  add    $0x10,%rsp
  mov    -0xa0(%rbp),%rax
  movl   $0x0,(%rax)
  lea    -0x48(%rbp),%rdi
  lea    -0x40(%rbp),%rcx
  lea    -0x60(%rbp),%rdx
  mov    -0x98(%rbp),%rsi
  mov    -0xa0(%rbp),%rax
  mov    %rdi,%r8
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <cdflib_dinvr@plt>
L443:
  mov    -0xa0(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x1,%eax
  jne    L515
  lea    -0x50(%rbp),%rcx
  lea    -0x58(%rbp),%rdx
  mov    -0x98(%rbp),%rsi
  mov    -0x90(%rbp),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <cdflib_cumt@plt>
  cmpq   $0x0,-0x18(%rbp)
  je     L204
  movsd  -0x58(%rbp),%xmm0
  mov    -0x80(%rbp),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  jmp    L268
L204:
  nop
  movsd  -0x50(%rbp),%xmm0
  mov    -0x88(%rbp),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
L268:
  lea    -0x48(%rbp),%rdi
  lea    -0x40(%rbp),%rcx
  lea    -0x60(%rbp),%rdx
  mov    -0x98(%rbp),%rsi
  mov    -0xa0(%rbp),%rax
  mov    %rdi,%r8
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <cdflib_dinvr@plt>
  jmp    L443
L515:
  nop
  mov    -0xa0(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0xffffffff,%eax
  jne    L343
  mov    -0x40(%rbp),%rax
  test   %rax,%rax
  je     L78
  mov    -0xa0(%rbp),%rax
  movl   $0x1,(%rax)
  movsd  D50(%rip),%xmm0
  mov    -0xa8(%rbp),%rax
  movsd  %xmm0,(%rax)
  jmp    L380
L78:
  nop
  mov    -0xa0(%rbp),%rax
  movl   $0x2,(%rax)
  movsd  D15(%rip),%xmm0
  mov    -0xa8(%rbp),%rax
  movsd  %xmm0,(%rax)
  jmp    L380
L409:
  nop
  jmp    L380
L67:
  nop
  jmp    L380
L343:
  nop
L380:
  nop
L238:
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L170
  call   <__stack_chk_fail@plt>
L170:
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: tag: D279
datatype: f64
name: extract_symbolic_data
arguments: tag: D67
datatype: f64
name: extract_symbolic_data
arguments: tag: D148
datatype: f64
name: extract_symbolic_data
arguments: tag: D116
datatype: f64
name: extract_symbolic_data
arguments: tag: D106
datatype: qword
name: extract_symbolic_data
arguments: tag: D165
datatype: byte[4]
name: extract_symbolic_data
arguments: tag: D265
datatype: f64
name: extract_symbolic_data
arguments: tag: D62
datatype: f32[2]
name: extract_symbolic_data
arguments: tag: D160
datatype: byte[8]
name: extract_symbolic_data
arguments: tag: D270
datatype: byte[8]
name: extract_symbolic_data
arguments: tag: D50
datatype: byte[8]
name: extract_symbolic_data
arguments: tag: D15
datatype: f64

tool: D279: 0.5

tool: D67: 5

tool: D148: 1

tool: D116: 3

tool: D106: 9223372036854775807

tool: {"D165": [0x01,0x00,0x00,0x00]}

tool: D265: 1e+30

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

tool: {"D160": [0x1f,0xb8,0xd4,0x4a,0x7a,0xee,0x8d,0x35]}

tool: {"D270": [0x3a,0x8c,0x30,0xe2,0x8e,0x79,0x45,0x3e]}

tool: {"D50": [0x59,0xf3,0xf8,0xc2,0x1f,0x6e,0xa5,0x01]}

tool: D15: 10000000000

assistant: ```c
void cdflib_cdft(const int *which, double *p, double *q, double *t, double *df, int *status, double *bound) {
  static const double tol = 1.0e-8;
  static const double atol = 1.0e-50;
  static const double zero = 1.0e-300;
  static const double inf = 1.0e30;
  static const double maxdf = 1.0e10;
  static const int K1 = 1;
  double K4 = 0.5e0;
  double K5 = 5.0e0;
  double fx, cum, ccum, pq;
  unsigned long qhi, qleft, qporq;
  double T2, T3, T6, T7, T8, T9, T10, T11;
  *status = 0;
  *bound = 0.0;
  if (!(*which < 1 || *which > 3)) {
    goto S30;
  }
  if (!(*which < 1)) {
    goto S10;
  }
  *bound = 1.0e0;
  goto S20;
S10:
  *bound = 3.0e0;
S20:
  *status = -1;
  return;
S30:
  if (*which == 1) {
    goto S70;
  }
  if (!(*p <= 0.0e0 || *p > 1.0e0)) {
    goto S60;
  }
  if (!(*p <= 0.0e0)) {
    goto S40;
  }
  *bound = 0.0e0;
  goto S50;
S40:
  *bound = 1.0e0;
S50:
  *status = -2;
  return;
S70:
S60:
  if (*which == 1) {
    goto S110;
  }
  if (!(*q <= 0.0e0 || *q > 1.0e0)) {
    goto S100;
  }
  if (!(*q <= 0.0e0)) {
    goto S80;
  }
  *bound = 0.0e0;
  goto S90;
S80:
  *bound = 1.0e0;
S90:
  *status = -3;
  return;
S110:
S100:
  if (*which == 3) {
    goto S130;
  }
  if (!(*df <= 0.0e0)) {
    goto S120;
  }
  *bound = 0.0e0;
  *status = -5;
  return;
S130:
S120:
  if (*which == 1) {
    goto S170;
  }
  pq = *p + *q;
  if (!(fabs(pq - 0.5e0 - 0.5e0) > 3.0e0 * dpmpar(K1))) {
    goto S160;
  }
  if (!(pq < 0.0e0)) {
    goto S140;
  }
  *bound = 0.0e0;
  goto S150;
S140:
  *bound = 1.0e0;
S150:
  *status = 3;
  return;
S170:
S160:
  if (!(*which == 1)) {
    qporq = *p <= *q;
  }
  if (1 == *which) {
    cdflib_cumt(t, df, p, q);
    *status = 0;
  } else if (2 == *which) {
    *t = cdflib_dt1(p, q, df);
    T2 = -inf;
    T3 = inf;
    T6 = atol;
    T7 = tol;
    cdflib_dstinv(&T2, &T3, &K4, &K4, &K5, &T6, &T7);
    *status = 0;
    cdflib_dinvr(status, t, &fx, &qleft, &qhi);
  S180:
    if (!(*status == 1)) {
      goto S210;
    }
    cdflib_cumt(t, df, &cum, &ccum);
    if (!qporq) {
      goto S190;
    }
    fx = cum - *p;
    goto S200;
  S190:
    fx = ccum - *q;
  S200:
    cdflib_dinvr(status, t, &fx, &qleft, &qhi);
    goto S180;
  S210:
    if (!(*status == -1)) {
      goto S240;
    }
    if (!qleft) {
      goto S220;
    }
    *status = 1;
    *bound = -inf;
    goto S230;
  S220:
    *status = 2;
    *bound = inf;
  S240:
  S230:;
  } else if (3 == *which) {
    *df = 5.0e0;
    T8 = zero;
    T9 = maxdf;
    T10 = atol;
    T11 = tol;
    cdflib_dstinv(&T8, &T9, &K4, &K4, &K5, &T10, &T11);
    *status = 0;
    cdflib_dinvr(status, df, &fx, &qleft, &qhi);
  S250:
    if (!(*status == 1)) {
      goto S280;
    }
    cdflib_cumt(t, df, &cum, &ccum);
    if (!qporq) {
      goto S260;
    }
    fx = cum - *p;
    goto S270;
  S260:
    fx = ccum - *q;
  S270:
    cdflib_dinvr(status, df, &fx, &qleft, &qhi);
    goto S250;
  S280:
    if (!(*status == -1)) {
      goto S310;
    }
    if (!qleft) {
      goto S290;
    }
    *status = 1;
    *bound = zero;
    goto S300;
  S290:
    *status = 2;
    *bound = maxdf;
  S300:;
  }
S310:
  return;
}
```