user: 下面的汇编指令对应的C源代码是什么：

```asm
<cdflib_cdfgam>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %rdx,%r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xc8,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0xb8(%rsp)
  xor    %eax,%eax
  mov    D139(%rip),%rax
  mov    0x100(%rsp),%r10
  mov    0x108(%rsp),%r11
  mov    %rax,0x40(%rsp)
  mov    D207(%rip),%rax
  movl   $0x0,(%r10)
  mov    %rax,0x48(%rsp)
  mov    (%rdi),%eax
  movq   $0x0,(%r11)
  lea    -0x1(%rax),%edx
  cmp    $0x3,%edx
  jbe    L252
  movsd  D116(%rip),%xmm0
  test   %eax,%eax
  jle    L354
L372:
  movsd  %xmm0,(%r11)
  movl   $0xffffffff,(%r10)
L268:
  mov    0xb8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L418
  add    $0xc8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L354:
  movsd  D12(%rip),%xmm0
  jmp    L372
  nopw   0x0(%rax,%rax,1)
L252:
  mov    %rdi,%rbx
  mov    %rsi,%r12
  mov    %rcx,%r14
  mov    %r8,%rbp
  mov    %r9,%r15
  cmp    $0x1,%eax
  je     L250
  movsd  (%rsi),%xmm1
  pxor   %xmm3,%xmm3
  comisd %xmm1,%xmm3
  ja     L143
  movsd  D12(%rip),%xmm0
  comisd %xmm0,%xmm1
  movapd %xmm0,%xmm2
  jbe    L368
L20:
  movsd  %xmm2,(%r11)
  movl   $0xfffffffe,(%r10)
  jmp    L268
  nopl   0x0(%rax)
L143:
  pxor   %xmm2,%xmm2
  jmp    L20
  cs nopw 0x0(%rax,%rax,1)
L368:
  movsd  0x0(%r13),%xmm2
  comisd %xmm2,%xmm3
  jae    L317
  comisd %xmm0,%xmm2
  jbe    L118
L309:
  movsd  %xmm0,(%r11)
  movl   $0xfffffffd,(%r10)
  jmp    L268
  nopl   0x0(%rax,%rax,1)
L118:
  cmp    $0x2,%eax
  jne    L99
  comisd (%r8),%xmm3
  jb     L402
L72:
  movl   $0xfffffffb,(%r10)
  jmp    L268
  nopl   0x0(%rax,%rax,1)
L317:
  pxor   %xmm0,%xmm0
  jmp    L309
  cs nopw 0x0(%rax,%rax,1)
L250:
  pxor   %xmm5,%xmm5
  comisd (%rcx),%xmm5
  jbe    L273
L306:
  movl   $0xfffffffc,(%r10)
  jmp    L268
L273:
  pxor   %xmm6,%xmm6
  comisd 0x0(%rbp),%xmm6
  jae    L72
  cmp    $0x4,%eax
  jne    L402
L501:
  movsd  (%r12),%xmm1
L466:
  addsd  0x0(%r13),%xmm1
  xor    %eax,%eax
  mov    $0x1,%edi
  mov    %r11,0x18(%rsp)
  mov    %r10,0x10(%rsp)
  movsd  %xmm1,0x8(%rsp)
  call   <dpmpar@plt>
  pxor   %xmm2,%xmm2
  movsd  0x8(%rsp),%xmm1
  mov    0x10(%rsp),%r10
  cvtsi2sd %eax,%xmm2
  mulsd  D164(%rip),%xmm2
  mov    0x18(%rsp),%r11
  movapd %xmm1,%xmm0
  subsd  D139(%rip),%xmm0
  subsd  D139(%rip),%xmm0
  andpd  D155(%rip),%xmm0
  comisd %xmm2,%xmm0
  jbe    L34
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  ja     L126
  movsd  D12(%rip),%xmm0
L126:
  movsd  %xmm0,(%r11)
  movl   $0x3,(%r10)
  jmp    L268
L402:
  movsd  (%r15),%xmm0
  pxor   %xmm7,%xmm7
  comisd %xmm0,%xmm7
  jb     L425
L365:
  movl   $0xfffffffa,(%r10)
  jmp    L268
L34:
  mov    (%rbx),%eax
  cmp    $0x1,%eax
  je     L277
  movsd  (%r12),%xmm1
  movsd  0x0(%r13),%xmm2
  cmp    $0x2,%eax
  je     L218
  cmp    $0x3,%eax
  je     L202
  cmp    $0x4,%eax
  jne    L268
  mov    D26(%rip),%rax
  mov    %r13,%r8
  mov    %r12,%rcx
  mov    %rbp,%rdi
  lea    0x50(%rsp),%rsi
  lea    0x3c(%rsp),%r9
  mov    %r10,0x8(%rsp)
  mov    %rax,0xb0(%rsp)
  lea    0xb0(%rsp),%rdx
  xor    %eax,%eax
  call   <cdflib_gamma_inc_inv@plt>
  cmpl   $0x0,0x3c(%rsp)
  mov    0x8(%rsp),%r10
  js     L47
  movsd  0x50(%rsp),%xmm0
  divsd  (%r14),%xmm0
  movsd  %xmm0,(%r15)
  movl   $0x0,(%r10)
  jmp    L268
L425:
  cmp    $0x1,%eax
  jne    L501
L108:
  mulsd  (%r14),%xmm0
  lea    0x60(%rsp),%rdi
  mov    %r13,%rcx
  mov    %r12,%rdx
  movl   $0x0,(%r10)
  mov    %rbp,%rsi
  xor    %eax,%eax
  movsd  %xmm0,0x60(%rsp)
  call   <cdflib_cumgam@plt>
  jmp    L268
L218:
  mov    D26(%rip),%rax
  mov    %r13,%r8
  mov    %r12,%rcx
  mov    %rbp,%rdi
  lea    0x50(%rsp),%rsi
  lea    0x3c(%rsp),%r9
  mov    %r10,0x8(%rsp)
  mov    %rax,0x88(%rsp)
  lea    0x88(%rsp),%rdx
  xor    %eax,%eax
  call   <cdflib_gamma_inc_inv@plt>
  cmpl   $0x0,0x3c(%rsp)
  mov    0x8(%rsp),%r10
  js     L47
  movsd  0x50(%rsp),%xmm0
  divsd  (%r15),%xmm0
  movsd  %xmm0,(%r14)
  movl   $0x0,(%r10)
  jmp    L268
L419:
  mov    %rbp,%r10
L47:
  movl   $0xa,(%r10)
  jmp    L268
L202:
  mov    D207(%rip),%rax
  mov    %r11,0x20(%rsp)
  lea    0x40(%rsp),%rdx
  lea    0x98(%rsp),%rsi
  mov    %r10,0x8(%rsp)
  lea    0x90(%rsp),%rdi
  mov    %rdx,%rcx
  mov    %rax,0x0(%rbp)
  mov    D113(%rip),%rax
  movsd  (%r14),%xmm0
  mulsd  (%r15),%xmm0
  movsd  %xmm2,0x18(%rsp)
  mov    %rax,0x90(%rsp)
  mov    D27(%rip),%rax
  movsd  %xmm1,0x10(%rsp)
  mov    %rax,0x98(%rsp)
  mov    D153(%rip),%rax
  mov    %rax,0xa0(%rsp)
  mov    D182(%rip),%rax
  movsd  %xmm0,0x60(%rsp)
  mov    %rax,0xa8(%rsp)
  push   %rax
  lea    0xb0(%rsp),%rax
  push   %rax
  xor    %eax,%eax
  lea    0xb0(%rsp),%r9
  lea    0x58(%rsp),%r8
  call   <cdflib_dstinv@plt>
  mov    0x18(%rsp),%r10
  xor    %eax,%eax
  mov    %rbp,%rsi
  lea    0x88(%rsp),%r14
  lea    0x90(%rsp),%r15
  movl   $0x0,(%r10)
  mov    %r15,%rcx
  mov    %r10,%rdi
  mov    %r14,%r8
  lea    0x68(%rsp),%rbx
  mov    %rbx,%rdx
  call   <cdflib_dinvr@plt>
  mov    0x18(%rsp),%r10
  pop    %rdx
  pop    %rcx
  movsd  0x10(%rsp),%xmm1
  mov    (%r10),%eax
  movsd  0x18(%rsp),%xmm2
  mov    0x20(%rsp),%r11
  cmp    $0x1,%eax
  jne    L91
  lea    0x70(%rsp),%rax
  comisd %xmm1,%xmm2
  mov    %rax,0x18(%rsp)
  lea    0x68(%rsp),%rax
  mov    %rax,0x10(%rsp)
  lea    0x60(%rsp),%rax
  mov    %rax,0x8(%rsp)
  jb     L258
  movsd  D48(%rip),%xmm6
  mov    %r11,0x28(%rsp)
  mov    %rbp,%r13
  mov    %r10,%rbp
  movsd  %xmm6,0x20(%rsp)
L461:
  mov    0x18(%rsp),%rcx
  mov    0x10(%rsp),%rdx
  xor    %eax,%eax
  mov    %r13,%rsi
  mov    0x8(%rsp),%rdi
  call   <cdflib_cumgam@plt>
  movsd  0x68(%rsp),%xmm0
  comisd 0x20(%rsp),%xmm0
  movapd %xmm0,%xmm1
  subsd  (%r12),%xmm1
  movsd  %xmm1,0x58(%rsp)
  ja     L419
  mov    %r14,%r8
  mov    %r15,%rcx
  mov    %rbx,%rdx
  mov    %r13,%rsi
  mov    %rbp,%rdi
  xor    %eax,%eax
  call   <cdflib_dinvr@plt>
  mov    0x0(%rbp),%eax
  cmp    $0x1,%eax
  je     L461
L17:
  mov    0x28(%rsp),%r11
  mov    %rbp,%r10
L91:
  add    $0x1,%eax
  jne    L268
  cmpq   $0x0,0x80(%rsp)
  jne    L474
  mov    D27(%rip),%rax
  movl   $0x2,(%r10)
  mov    %rax,(%r11)
  jmp    L268
L418:
  call   <__stack_chk_fail@plt>
L258:
  movsd  D48(%rip),%xmm7
  mov    %r11,0x28(%rsp)
  mov    %rbp,%r12
  mov    %r10,%rbp
  movsd  %xmm7,0x20(%rsp)
L95:
  mov    0x18(%rsp),%rcx
  mov    0x10(%rsp),%rdx
  xor    %eax,%eax
  mov    %r12,%rsi
  mov    0x8(%rsp),%rdi
  call   <cdflib_cumgam@plt>
  movsd  0x70(%rsp),%xmm0
  comisd 0x20(%rsp),%xmm0
  movapd %xmm0,%xmm1
  subsd  0x0(%r13),%xmm1
  movsd  %xmm1,0x58(%rsp)
  ja     L419
  mov    %r14,%r8
  mov    %r15,%rcx
  mov    %rbx,%rdx
  mov    %r12,%rsi
  mov    %rbp,%rdi
  xor    %eax,%eax
  call   <cdflib_dinvr@plt>
  mov    0x0(%rbp),%eax
  cmp    $0x1,%eax
  je     L95
  jmp    L17
L474:
  mov    D113(%rip),%rax
  movl   $0x1,(%r10)
  mov    %rax,(%r11)
  jmp    L268
L277:
  movsd  (%r15),%xmm0
  jmp    L108
L99:
  pxor   %xmm7,%xmm7
  comisd (%rcx),%xmm7
  ja     L306
  cmp    $0x3,%eax
  jne    L273
  comisd (%r9),%xmm7
  jae    L365
  jmp    L466
```

tool_calls: name: unpack_data
arguments: data_tag: D139
data_type: f64
name: unpack_data
arguments: data_tag: D207
data_type: f64
name: unpack_data
arguments: data_tag: D116
data_type: f64
name: unpack_data
arguments: data_tag: D12
data_type: f64
name: unpack_data
arguments: data_tag: D164
data_type: f64
name: unpack_data
arguments: data_tag: D155
data_type: qword[2]
name: unpack_data
arguments: data_tag: D26
data_type: f32[2]
name: unpack_data
arguments: data_tag: D113
data_type: byte[8]
name: unpack_data
arguments: data_tag: D27
data_type: f64
name: unpack_data
arguments: data_tag: D153
data_type: byte[8]
name: unpack_data
arguments: data_tag: D182
data_type: byte[8]
name: unpack_data
arguments: data_tag: D48
data_type: f64

tool: D139: 0.5

tool: D207: 5

tool: D116: 4

tool: D12: 1

tool: D164: 3

tool: {"D155": [0x7fffffffffffffff,0x0000000000000000]}

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

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

tool: D27: 1e+300

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

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

tool: D48: 1.5

assistant: ```c
void cdflib_cdfgam(const int *which, double *p, double *q, double *x, double *shape, double *scale, 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.0e300;
  static const int K1 = 1;
  double K5 = 0.5e0;
  double K6 = 5.0e0;
  double xx, fx, xscale, cum, ccum, pq, porq = 0;
  int ierr;
  unsigned long qhi, qleft, qporq;
  double T2, T3, T4, T7, T8, T9;
  *status = 0;
  *bound = 0.0;
  if (!(*which < 1 || *which > 4)) {
    goto S30;
  }
  if (!(*which < 1)) {
    goto S10;
  }
  *bound = 1.0e0;
  goto S20;
S10:
  *bound = 4.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 == 2) {
    goto S130;
  }
  if (!(*x < 0.0e0)) {
    goto S120;
  }
  *bound = 0.0e0;
  *status = -4;
  return;
S130:
S120:
  if (*which == 3) {
    goto S150;
  }
  if (!(*shape <= 0.0e0)) {
    goto S140;
  }
  *bound = 0.0e0;
  *status = -5;
  return;
S150:
S140:
  if (*which == 4) {
    goto S170;
  }
  if (!(*scale <= 0.0e0)) {
    goto S160;
  }
  *bound = 0.0e0;
  *status = -6;
  return;
S170:
S160:
  if (*which == 1) {
    goto S210;
  }
  pq = *p + *q;
  if (!(fabs(pq - 0.5e0 - 0.5e0) > 3.0e0 * dpmpar(K1))) {
    goto S200;
  }
  if (!(pq < 0.0e0)) {
    goto S180;
  }
  *bound = 0.0e0;
  goto S190;
S180:
  *bound = 1.0e0;
S190:
  *status = 3;
  return;
S210:
S200:
  if (*which == 1) {
    goto S240;
  }
  qporq = *p <= *q;
  if (!qporq) {
    goto S220;
  }
  porq = *p;
  goto S230;
S220:
  porq = *q;
S240:
S230:
  if (1 == *which) {
    *status = 0;
    xscale = *x * *scale;
    cdflib_cumgam(&xscale, shape, p, q);
    if (porq > 1.5e0) {
      *status = 10;
    }
  } else if (2 == *which) {
    T2 = -1.0e0;
    cdflib_gamma_inc_inv(shape, &xx, &T2, p, q, &ierr);
    if (ierr < 0.0e0) {
      *status = 10;
      return;
    } else {
      *x = xx / *scale;
      *status = 0;
    }
  } else if (3 == *which) {
    *shape = 5.0e0;
    xscale = *x * *scale;
    T3 = zero;
    T4 = inf;
    T7 = atol;
    T8 = tol;
    cdflib_dstinv(&T3, &T4, &K5, &K5, &K6, &T7, &T8);
    *status = 0;
    cdflib_dinvr(status, shape, &fx, &qleft, &qhi);
  S250:
    if (!(*status == 1)) {
      goto S290;
    }
    cdflib_cumgam(&xscale, shape, &cum, &ccum);
    if (!qporq) {
      goto S260;
    }
    fx = cum - *p;
    goto S270;
  S260:
    fx = ccum - *q;
  S270:
    if (!((qporq && cum > 1.5e0) || (!qporq && ccum > 1.5e0))) {
      goto S280;
    }
    *status = 10;
    return;
  S280:
    cdflib_dinvr(status, shape, &fx, &qleft, &qhi);
    goto S250;
  S290:
    if (!(*status == -1)) {
      goto S320;
    }
    if (!qleft) {
      goto S300;
    }
    *status = 1;
    *bound = zero;
    goto S310;
  S300:
    *status = 2;
    *bound = inf;
  S320:
  S310:;
  } else if (4 == *which) {
    T9 = -1.0e0;
    cdflib_gamma_inc_inv(shape, &xx, &T9, p, q, &ierr);
    if (ierr < 0.0e0) {
      *status = 10;
      return;
    } else {
      *scale = xx / *x;
      *status = 0;
    }
  }
  return;
}
```