user: 这段汇编代码用C语言表示是什么样的?

```asm
<invgammp>:
  endbr64
  push   %r14
  push   %rbx
  sub    $0x48,%rsp
  movsd  %xmm0,0x18(%rsp)
  movsd  %xmm1,0x8(%rsp)
  movapd %xmm1,%xmm0
  call   <gammln@plt>
  movsd  %xmm0,0x28(%rsp)
  pxor   %xmm0,%xmm0
  comisd 0x8(%rsp),%xmm0
  jae    L21
L69:
  movsd  0x18(%rsp),%xmm7
  comisd D169(%rip),%xmm7
  jae    L43
  pxor   %xmm0,%xmm0
  comisd 0x18(%rsp),%xmm0
  jae    L175
  movsd  D169(%rip),%xmm0
  movsd  0x8(%rsp),%xmm2
  movapd %xmm2,%xmm5
  subsd  %xmm0,%xmm5
  movsd  %xmm5,0x10(%rsp)
  comisd %xmm0,%xmm2
  ja     L198
  movsd  0x8(%rsp),%xmm5
  movapd %xmm5,%xmm0
  mulsd  D298(%rip),%xmm0
  addsd  D9(%rip),%xmm0
  mulsd  %xmm5,%xmm0
  movsd  D169(%rip),%xmm1
  subsd  %xmm0,%xmm1
  movsd  0x18(%rsp),%xmm7
  comisd %xmm7,%xmm1
  jbe    L96
  movsd  D169(%rip),%xmm0
  divsd  %xmm5,%xmm0
  cvttsd2si %xmm0,%edi
  divsd  %xmm1,%xmm7
  movapd %xmm7,%xmm0
  call   <pow@plt>
  movsd  %xmm0,(%rsp)
  pxor   %xmm5,%xmm5
  movsd  %xmm5,0x38(%rsp)
  movsd  %xmm5,0x30(%rsp)
  jmp    L70
L21:
  lea    D110(%rip),%rdi
  call   <error@plt>
  jmp    L69
L43:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D112(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  addsd  0x8(%rsp),%xmm1
  call   <dmax@plt>
  movsd  %xmm0,(%rsp)
L24:
  movsd  (%rsp),%xmm0
  add    $0x48,%rsp
  pop    %rbx
  pop    %r14
  ret
L198:
  movapd %xmm5,%xmm0
  call   <log@plt>
  movapd %xmm0,%xmm6
  movsd  %xmm0,0x38(%rsp)
  subsd  D169(%rip),%xmm6
  movapd %xmm6,%xmm0
  mulsd  0x10(%rsp),%xmm0
  subsd  0x28(%rsp),%xmm0
  call   <exp@plt>
  movsd  %xmm0,0x30(%rsp)
  movsd  D262(%rip),%xmm0
  movsd  0x18(%rsp),%xmm2
  comisd %xmm2,%xmm0
  ja     L35
  movsd  D169(%rip),%xmm0
  subsd  %xmm2,%xmm0
  call   <log@plt>
  mulsd  D380(%rip),%xmm0
  call   <sqrt@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movapd %xmm1,%xmm6
  mulsd  D286(%rip),%xmm6
  movapd %xmm6,%xmm2
  addsd  D297(%rip),%xmm2
  movapd %xmm1,%xmm0
  mulsd  D97(%rip),%xmm0
  addsd  D33(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D169(%rip),%xmm0
  divsd  %xmm0,%xmm2
  movapd %xmm2,%xmm6
  subsd  %xmm1,%xmm6
  movq   %xmm6,%rbx
  movsd  D262(%rip),%xmm0
  comisd 0x18(%rsp),%xmm0
  ja     L73
L131:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  movsd  D169(%rip),%xmm0
  movsd  0x8(%rsp),%xmm1
  mulsd  D194(%rip),%xmm1
  movapd %xmm0,%xmm2
  divsd  %xmm1,%xmm2
  subsd  %xmm2,%xmm0
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
  mulsd  D160(%rip),%xmm2
  movq   %rbx,%xmm5
  divsd  %xmm2,%xmm5
  subsd  %xmm5,%xmm0
  mov    $0x3,%edi
  call   <pow@plt>
  mulsd  0x8(%rsp),%xmm0
  movapd %xmm0,%xmm1
  movsd  D40(%rip),%xmm0
  call   <dmax@plt>
  movsd  %xmm0,(%rsp)
L70:
  mov    $0xc,%ebx
  jmp    L142
L96:
  movsd  0x18(%rsp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  D169(%rip),%xmm2
  subsd  %xmm1,%xmm2
  divsd  %xmm2,%xmm0
  movsd  D169(%rip),%xmm1
  subsd  %xmm0,%xmm1
  movapd %xmm1,%xmm0
  call   <log@plt>
  movsd  D169(%rip),%xmm7
  subsd  %xmm0,%xmm7
  movsd  %xmm7,(%rsp)
  pxor   %xmm7,%xmm7
  movsd  %xmm7,0x38(%rsp)
  movsd  %xmm7,0x30(%rsp)
  jmp    L70
L67:
  movsd  (%rsp),%xmm0
  call   <log@plt>
  mulsd  0x10(%rsp),%xmm0
  subsd  (%rsp),%xmm0
  subsd  0x28(%rsp),%xmm0
  call   <exp@plt>
  jmp    L118
L192:
  andpd  D271(%rip),%xmm0
  movsd  (%rsp),%xmm1
  mulsd  D11(%rip),%xmm1
  comisd %xmm0,%xmm1
  ja     L24
  sub    $0x1,%ebx
  je     L24
L142:
  pxor   %xmm7,%xmm7
  movsd  (%rsp),%xmm1
  comisd %xmm1,%xmm7
  jae    L44
  movsd  0x8(%rsp),%xmm0
  call   <gammp@plt>
  movapd %xmm0,%xmm5
  subsd  0x18(%rsp),%xmm5
  movq   %xmm5,%r14
  movsd  0x8(%rsp),%xmm6
  comisd D169(%rip),%xmm6
  jbe    L67
  movsd  (%rsp),%xmm0
  call   <log@plt>
  subsd  0x38(%rsp),%xmm0
  movsd  0x10(%rsp),%xmm6
  mulsd  %xmm6,%xmm0
  movsd  (%rsp),%xmm1
  subsd  %xmm6,%xmm1
  subsd  %xmm1,%xmm0
  call   <exp@plt>
  mulsd  0x30(%rsp),%xmm0
L118:
  movq   %r14,%xmm4
  divsd  %xmm0,%xmm4
  movsd  0x10(%rsp),%xmm0
  divsd  (%rsp),%xmm0
  subsd  D169(%rip),%xmm0
  movsd  %xmm4,0x20(%rsp)
  mulsd  %xmm4,%xmm0
  mov    $0x1,%edi
  call   <dmin@plt>
  mulsd  D262(%rip),%xmm0
  movsd  D169(%rip),%xmm1
  subsd  %xmm0,%xmm1
  movsd  0x20(%rsp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  (%rsp),%xmm3
  subsd  %xmm0,%xmm3
  movsd  %xmm3,(%rsp)
  pxor   %xmm1,%xmm1
  comisd %xmm3,%xmm1
  jb     L192
  movapd %xmm3,%xmm1
  addsd  %xmm0,%xmm1
  mulsd  D262(%rip),%xmm1
  movsd  %xmm1,(%rsp)
  jmp    L192
L175:
  pxor   %xmm2,%xmm2
  movsd  %xmm2,(%rsp)
  jmp    L24
L44:
  pxor   %xmm5,%xmm5
  movsd  %xmm5,(%rsp)
  jmp    L24
L35:
  movsd  0x18(%rsp),%xmm0
  call   <log@plt>
  mulsd  D380(%rip),%xmm0
  call   <sqrt@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movapd %xmm1,%xmm2
  mulsd  D286(%rip),%xmm2
  addsd  D297(%rip),%xmm2
  movapd %xmm1,%xmm0
  mulsd  D97(%rip),%xmm0
  addsd  D33(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D169(%rip),%xmm0
  divsd  %xmm0,%xmm2
  subsd  %xmm1,%xmm2
  movq   %xmm2,%rbx
L73:
  movq   %rbx,%xmm2
  xorpd  D100(%rip),%xmm2
  movq   %xmm2,%rbx
  jmp    L131
```

tool_calls: name: parse_data
arguments: label: D169
type_guessed: f64[2]
name: parse_data
arguments: label: D298
type_guessed: f64
name: parse_data
arguments: label: D9
type_guessed: f64
name: parse_data
arguments: label: D110
type_guessed: string
name: parse_data
arguments: label: D112
type_guessed: f64
name: parse_data
arguments: label: D262
type_guessed: f64
name: parse_data
arguments: label: D380
type_guessed: f32[2]
name: parse_data
arguments: label: D286
type_guessed: f64
name: parse_data
arguments: label: D297
type_guessed: f64
name: parse_data
arguments: label: D97
type_guessed: f64
name: parse_data
arguments: label: D33
type_guessed: f64
name: parse_data
arguments: label: D194
type_guessed: f64
name: parse_data
arguments: label: D160
type_guessed: f64
name: parse_data
arguments: label: D40
type_guessed: byte[8]
name: parse_data
arguments: label: D271
type_guessed: qword[2]
name: parse_data
arguments: label: D11
type_guessed: byte[8]
name: parse_data
arguments: label: D100
type_guessed: f32[4]

tool: {"D169": [1.0,100.0]}

tool: D298: 0.12

tool: D9: 0.253

tool: D110: a must be pos in invgammap

tool: D112: 100

tool: D262: 0.5

tool: {"D380": [0.0,-2.0]}

tool: D286: 0.27061

tool: D297: 2.30753

tool: D97: 0.04481

tool: D33: 0.99229

tool: D194: 9

tool: D160: 3

tool: {"D40": [0xfc,0xa9,0xf1,0xd2,0x4d,0x62,0x50,0x3f]}

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

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

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

assistant: ```c
double invgammp(const double p, const double a) {
  int j;
  double x, err, t, u, pp, a1 = a - 1;
  double afac = 0;
  double lna1 = 0;
  double eps = 1.e-8;
  double gln = gammln(a);
  if (a <= 0.)
    error("a must be pos in invgammap");
  if (p >= 1.)
    return dmax(100., a + 100. * sqrt(a));
  if (p <= 0.)
    return 0.0;
  if (a > 1.) {
    lna1 = log(a1);
    afac = exp(a1 * (lna1 - 1.) - gln);
    pp = (p < 0.5) ? p : 1. - p;
    t = sqrt(-2. * log(pp));
    x = (2.30753 + t * 0.27061) / (1. + t * (0.99229 + t * 0.04481)) - t;
    if (p < 0.5)
      x = -x;
    x = dmax(1.e-3, a * pow(1. - 1. / (9. * a) - x / (3. * sqrt(a)), 3));
  } else {
    t = 1.0 - a * (0.253 + a * 0.12);
    if (p < t)
      x = pow(p / t, 1. / a);
    else
      x = 1. - log(1. - (p - t) / (1. - t));
  }
  for (j = 0; j < 12; j++) {
    if (x <= 0.0)
      return 0.0;
    err = gammp(a, x) - p;
    if (a > 1.)
      t = afac * exp(-(x - a1) + a1 * (log(x) - lna1));
    else
      t = exp(-x + a1 * log(x) - gln);
    u = err / t;
    x -= (t = u / (1. - 0.5 * dmin(1., u * ((a - 1.) / x - 1))));
    if (x <= 0.)
      x = 0.5 * (x + t);
    if (fabs(t) < eps * x)
      break;
  }
  return x;
}
```