user: 请从给定的汇编指令推导出C源代码:

```asm
<PointChi2>:
  endbr64
  movapd %xmm0,%xmm7
  movsd  D124(%rip),%xmm0
  comisd %xmm7,%xmm0
  ja     L154
  push   %r14
  push   %rbx
  sub    $0x58,%rsp
  movq   %xmm1,%rbx
  movsd  D103(%rip),%xmm0
  movsd  %xmm7,0x38(%rsp)
  comisd D642(%rip),%xmm7
  ja     L119
  pxor   %xmm1,%xmm1
  movq   %rbx,%xmm5
  comisd %xmm5,%xmm1
  jae    L119
  movq   %rbx,%xmm4
  mulsd  D564(%rip),%xmm4
  movsd  %xmm4,0x10(%rsp)
  movapd %xmm4,%xmm0
  call   <LnGamma@plt>
  movsd  %xmm0,0x18(%rsp)
  movsd  0x10(%rsp),%xmm4
  subsd  D525(%rip),%xmm4
  movsd  %xmm4,0x20(%rsp)
  movsd  0x38(%rsp),%xmm0
  call   <log@plt>
  mulsd  D466(%rip),%xmm0
  movq   %rbx,%xmm5
  comisd %xmm0,%xmm5
  jb     L150
  comisd D140(%rip),%xmm5
  ja     L109
  movsd  D525(%rip),%xmm0
  subsd  0x38(%rsp),%xmm0
  call   <log@plt>
  movsd  0x20(%rsp),%xmm5
  mulsd  D497(%rip),%xmm5
  movsd  %xmm5,0x40(%rsp)
  movsd  D216(%rip),%xmm5
  movsd  %xmm5,0x8(%rsp)
  addsd  0x18(%rsp),%xmm0
  movsd  %xmm0,0x48(%rsp)
L53:
  movsd  D699(%rip),%xmm0
  movsd  0x8(%rsp),%xmm5
  addsd  %xmm5,%xmm0
  mulsd  %xmm5,%xmm0
  addsd  D525(%rip),%xmm0
  movsd  %xmm0,0x28(%rsp)
  movapd %xmm5,%xmm0
  addsd  D469(%rip),%xmm0
  mulsd  %xmm5,%xmm0
  addsd  D511(%rip),%xmm0
  movapd %xmm0,%xmm7
  mulsd  %xmm5,%xmm7
  movsd  %xmm7,0x30(%rsp)
  movsd  D564(%rip),%xmm2
  mulsd  %xmm5,%xmm2
  movapd %xmm2,%xmm0
  addsd  0x48(%rsp),%xmm0
  addsd  0x40(%rsp),%xmm0
  call   <exp@plt>
  movsd  0x30(%rsp),%xmm7
  mulsd  %xmm7,%xmm0
  movsd  0x28(%rsp),%xmm4
  divsd  %xmm4,%xmm0
  movsd  D525(%rip),%xmm2
  subsd  %xmm0,%xmm2
  movsd  0x8(%rsp),%xmm5
  movapd %xmm5,%xmm0
  addsd  %xmm5,%xmm0
  addsd  D699(%rip),%xmm0
  divsd  %xmm4,%xmm0
  subsd  D564(%rip),%xmm0
  movapd %xmm5,%xmm1
  mulsd  D252(%rip),%xmm1
  addsd  D624(%rip),%xmm1
  mulsd  %xmm5,%xmm1
  addsd  D511(%rip),%xmm1
  divsd  %xmm7,%xmm1
  subsd  %xmm1,%xmm0
  divsd  %xmm0,%xmm2
  movapd %xmm5,%xmm0
  subsd  %xmm2,%xmm5
  movsd  %xmm5,0x8(%rsp)
  divsd  %xmm5,%xmm0
  subsd  D525(%rip),%xmm0
  andpd  D440(%rip),%xmm0
  subsd  D185(%rip),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm0,%xmm1
  jb     L53
  jmp    L21
L150:
  movsd  0x10(%rsp),%xmm0
  mulsd  D497(%rip),%xmm0
  addsd  0x18(%rsp),%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm1
  movsd  0x10(%rsp),%xmm5
  movapd %xmm5,%xmm0
  mulsd  0x38(%rsp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D525(%rip),%xmm1
  divsd  %xmm5,%xmm1
  call   <pow@plt>
  movsd  %xmm0,0x8(%rsp)
  movapd %xmm0,%xmm1
  subsd  D391(%rip),%xmm1
  pxor   %xmm2,%xmm2
  comisd %xmm1,%xmm2
  jbe    L21
  jmp    L119
L109:
  movsd  0x38(%rsp),%xmm0
  call   <PointNormal@plt>
  movq   %xmm0,%r14
  movsd  D39(%rip),%xmm0
  movq   %rbx,%xmm5
  divsd  %xmm5,%xmm0
  movapd %xmm0,%xmm4
  movsd  %xmm0,0x8(%rsp)
  pxor   %xmm0,%xmm0
  ucomisd %xmm4,%xmm0
  ja     L107
  sqrtsd %xmm4,%xmm4
  movapd %xmm4,%xmm0
L77:
  movq   %r14,%xmm5
  mulsd  %xmm5,%xmm0
  addsd  D525(%rip),%xmm0
  subsd  0x8(%rsp),%xmm0
  movsd  D252(%rip),%xmm1
  call   <pow@plt>
  movapd %xmm0,%xmm4
  movq   %rbx,%xmm5
  mulsd  %xmm5,%xmm4
  movsd  %xmm4,0x8(%rsp)
  movq   %rbx,%xmm0
  mulsd  D217(%rip),%xmm0
  addsd  D468(%rip),%xmm0
  comisd %xmm0,%xmm4
  jbe    L21
  movsd  D525(%rip),%xmm0
  subsd  0x38(%rsp),%xmm0
  call   <log@plt>
  movq   %xmm0,%rbx
  movsd  0x8(%rsp),%xmm0
  mulsd  D564(%rip),%xmm0
  call   <log@plt>
  movapd %xmm0,%xmm1
  mulsd  0x20(%rsp),%xmm1
  movq   %rbx,%xmm0
  subsd  %xmm1,%xmm0
  addsd  0x18(%rsp),%xmm0
  mulsd  D604(%rip),%xmm0
  movsd  %xmm0,0x8(%rsp)
  jmp    L21
L107:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L77
L83:
  movsd  %xmm0,0x8(%rsp)
L21:
  movsd  D564(%rip),%xmm4
  mulsd  0x8(%rsp),%xmm4
  movsd  %xmm4,0x28(%rsp)
  movsd  0x18(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm1
  movapd %xmm4,%xmm0
  call   <IncompleteGamma@plt>
  pxor   %xmm4,%xmm4
  comisd %xmm0,%xmm4
  ja     L71
  movsd  0x38(%rsp),%xmm4
  subsd  %xmm0,%xmm4
  movsd  %xmm4,0x30(%rsp)
  movsd  0x8(%rsp),%xmm0
  call   <log@plt>
  movapd %xmm0,%xmm1
  movsd  0x10(%rsp),%xmm7
  mulsd  D497(%rip),%xmm7
  movapd %xmm7,%xmm0
  addsd  0x18(%rsp),%xmm0
  addsd  0x28(%rsp),%xmm0
  mulsd  0x20(%rsp),%xmm1
  subsd  %xmm1,%xmm0
  call   <exp@plt>
  movsd  0x30(%rsp),%xmm4
  mulsd  %xmm0,%xmm4
  movapd %xmm4,%xmm5
  divsd  0x8(%rsp),%xmm5
  movapd %xmm4,%xmm0
  mulsd  D564(%rip),%xmm0
  movsd  0x20(%rsp),%xmm2
  movapd %xmm2,%xmm9
  mulsd  %xmm5,%xmm9
  movapd %xmm0,%xmm3
  subsd  %xmm9,%xmm3
  movapd %xmm3,%xmm1
  mulsd  D539(%rip),%xmm1
  addsd  D375(%rip),%xmm1
  mulsd  %xmm3,%xmm1
  movsd  D480(%rip),%xmm15
  addsd  %xmm15,%xmm1
  mulsd  %xmm3,%xmm1
  addsd  D580(%rip),%xmm1
  mulsd  %xmm3,%xmm1
  addsd  D517(%rip),%xmm1
  mulsd  %xmm3,%xmm1
  movsd  D434(%rip),%xmm11
  addsd  %xmm11,%xmm1
  movsd  D228(%rip),%xmm10
  divsd  %xmm10,%xmm1
  mulsd  %xmm1,%xmm0
  movsd  D525(%rip),%xmm8
  addsd  %xmm8,%xmm0
  movapd %xmm3,%xmm6
  mulsd  D223(%rip),%xmm6
  addsd  D358(%rip),%xmm6
  mulsd  %xmm3,%xmm6
  addsd  D473(%rip),%xmm6
  mulsd  %xmm3,%xmm6
  addsd  D583(%rip),%xmm6
  mulsd  %xmm3,%xmm6
  addsd  %xmm10,%xmm6
  movsd  D144(%rip),%xmm14
  divsd  %xmm14,%xmm6
  movapd %xmm6,%xmm10
  movapd %xmm3,%xmm6
  mulsd  D704(%rip),%xmm6
  addsd  D151(%rip),%xmm6
  mulsd  %xmm3,%xmm6
  addsd  D110(%rip),%xmm6
  mulsd  %xmm3,%xmm6
  addsd  %xmm6,%xmm11
  divsd  %xmm14,%xmm11
  movapd %xmm3,%xmm6
  mulsd  D606(%rip),%xmm6
  addsd  D673(%rip),%xmm6
  mulsd  %xmm3,%xmm6
  addsd  D15(%rip),%xmm6
  mulsd  %xmm2,%xmm6
  movapd %xmm3,%xmm12
  mulsd  D112(%rip),%xmm12
  addsd  D386(%rip),%xmm12
  mulsd  %xmm3,%xmm12
  addsd  D371(%rip),%xmm12
  addsd  %xmm12,%xmm6
  movsd  D18(%rip),%xmm12
  divsd  %xmm12,%xmm6
  movapd %xmm3,%xmm13
  mulsd  D439(%rip),%xmm13
  addsd  D416(%rip),%xmm13
  mulsd  %xmm2,%xmm13
  movapd %xmm13,%xmm7
  mulsd  D83(%rip),%xmm3
  addsd  %xmm15,%xmm3
  addsd  %xmm3,%xmm7
  divsd  %xmm14,%xmm7
  movapd %xmm2,%xmm3
  mulsd  D640(%rip),%xmm2
  addsd  D595(%rip),%xmm2
  mulsd  %xmm3,%xmm2
  addsd  D607(%rip),%xmm2
  divsd  %xmm12,%xmm2
  mulsd  %xmm5,%xmm2
  subsd  %xmm2,%xmm7
  mulsd  %xmm5,%xmm7
  subsd  %xmm7,%xmm6
  movapd %xmm6,%xmm2
  mulsd  %xmm5,%xmm2
  subsd  %xmm2,%xmm11
  mulsd  %xmm5,%xmm11
  subsd  %xmm11,%xmm10
  mulsd  %xmm5,%xmm10
  subsd  %xmm10,%xmm1
  mulsd  %xmm9,%xmm1
  subsd  %xmm1,%xmm0
  mulsd  %xmm4,%xmm0
  movsd  0x8(%rsp),%xmm4
  addsd  %xmm4,%xmm0
  divsd  %xmm0,%xmm4
  movapd %xmm4,%xmm1
  subsd  %xmm8,%xmm1
  andpd  D440(%rip),%xmm1
  comisd D391(%rip),%xmm1
  ja     L83
  jmp    L119
L71:
  lea    D496(%rip),%rdi
  call   <printf@plt>
  movsd  D103(%rip),%xmm0
L119:
  add    $0x58,%rsp
  pop    %rbx
  pop    %r14
  ret
L154:
  movsd  D103(%rip),%xmm0
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D124
data_type: byte[8]
name: extract_symbolic_data
arguments: data_tag: D103
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D642
data_type: byte[16]
name: extract_symbolic_data
arguments: data_tag: D564
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D525
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D466
data_type: byte[8]
name: extract_symbolic_data
arguments: data_tag: D140
data_type: byte[16]
name: extract_symbolic_data
arguments: data_tag: D497
data_type: byte[8]
name: extract_symbolic_data
arguments: data_tag: D216
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D699
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D469
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D511
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D252
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D624
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D440
data_type: qword[2]
name: extract_symbolic_data
arguments: data_tag: D185
data_type: byte[8]
name: extract_symbolic_data
arguments: data_tag: D391
data_type: byte[16]
name: extract_symbolic_data
arguments: data_tag: D39
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D217
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D468
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D604
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D539
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D375
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D480
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D580
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D517
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D434
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D228
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D223
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D358
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D473
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D583
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D144
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D704
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D151
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D110
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D606
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D673
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D15
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D112
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D386
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D371
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D18
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D439
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D416
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D83
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D640
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D595
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D607
data_type: f32[2]
name: extract_symbolic_data
arguments: data_tag: D496
data_type: string

tool: {"D124": [0x8d,0xed,0xb5,0xa0,0xf7,0xc6,0xc0,0x3e]}

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

tool: {"D642": [0xd3,0x17,0x42,0xce,0xfb,0xff,0xef,0x3f,0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0x3f]}

tool: D564: 0.5

tool: {"D525": [0.0,1.875]}

tool: {"D466": [0xd7,0xa3,0x70,0x3d,0x0a,0xd7,0xf3,0xbf]}

tool: {"D140": [0x7b,0x14,0xae,0x47,0xe1,0x7a,0xd4,0x3f,0xcc,0xfc,0xf1,0xfe,0x42,0x2e,0xe6,0x3f]}

tool: {"D497": [0xcc,0xfc,0xf1,0xfe,0x42,0x2e,0xe6,0x3f]}

tool: D216: 0.4

tool: D699: 4.67

tool: D469: 6.66

tool: D511: 6.73

tool: D252: 3

tool: D624: 13.32

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

tool: {"D185": [0x7b,0x14,0xae,0x47,0xe1,0x7a,0x84,0x3f]}

tool: {"D391": [0x8d,0xed,0xb5,0xa0,0xf7,0xc6,0xa0,0x3e,0xb3,0xd1,0x39,0x3f,0xc5,0x71,0xcc,0x3f]}

tool: D39: 0.222222

tool: D217: 2.2

tool: {"D468": [0.0,2.375]}

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

tool: {"D539": [0.0,3.21875]}

tool: {"D375": [0.0,3.2734375]}

tool: {"D480": [0.0,3.328125]}

tool: {"D580": [0.0,3.41015625]}

tool: {"D517": [0.0,3.5234375]}

tool: {"D434": [0.0,3.66015625]}

tool: {"D228": [0.0,3.91015625]}

tool: {"D223": [0.0,4.6240234375]}

tool: {"D358": [0.0,4.55712890625]}

tool: {"D473": [0.0,4.443359375]}

tool: {"D583": [0.0,4.2177734375]}

tool: {"D144": [0.0,5.115234375]}

tool: {"D704": [0.0,4.41015625]}

tool: {"D151": [0.0,4.1904296875]}

tool: {"D110": [0.0,3.951171875]}

tool: {"D606": [0.0,4.849609375]}

tool: {"D673": [0.0,4.3681640625]}

tool: {"D15": [0.0,3.787109375]}

tool: {"D112": [0.0,4.5771484375]}

tool: {"D386": [0.0,4.15625]}

tool: {"D371": [0.0,3.7421875]}

tool: {"D18": [0.0,5.615234375]}

tool: {"D439": [0.0,4.091796875]}

tool: {"D416": [0.0,3.591796875]}

tool: {"D83": [0.0,3.7578125]}

tool: {"D640": [0.0,3.49609375]}

tool: {"D595": [0.0,3.837890625]}

tool: {"D607": [0.0,3.46875]}

tool: D496: IncompleteGamma 


assistant: ```c
double PointChi2(double prob, double v) {
  double e = .5e-6, aa = .6931471805, p = prob, g;
  double xx, c, ch, a = 0, q = 0, p1 = 0, p2 = 0, t = 0, x = 0, b = 0, s1, s2, s3, s4, s5, s6;
  if (p < .000002 || p > .999998 || v <= 0)
    return (-1);
  g = LnGamma(v / 2);
  xx = v / 2;
  c = xx - 1;
  if (v >= -1.24 * log(p))
    goto l1;
  ch = pow((p * xx * exp(g + xx * aa)), 1 / xx);
  if (ch - e < 0)
    return (ch);
  goto l4;
l1:
  if (v > .32)
    goto l3;
  ch = 0.4;
  a = log(1 - p);
l2:
  q = ch;
  p1 = 1 + ch * (4.67 + ch);
  p2 = ch * (6.73 + ch * (6.66 + ch));
  t = -0.5 + (4.67 + 2 * ch) / p1 - (6.73 + ch * (13.32 + 3 * ch)) / p2;
  ch -= (1 - exp(a + g + .5 * ch + c * aa) * p2 / p1) / t;
  if (fabs(q / ch - 1) - .01 <= 0)
    goto l4;
  else
    goto l2;
l3:
  x = PointNormal(p);
  p1 = 0.222222 / v;
  ch = v * pow((x * sqrt(p1) + 1 - p1), 3.0);
  if (ch > 2.2 * v + 6)
    ch = -2 * (log(1 - p) - c * log(.5 * ch) + g);
l4:
  q = ch;
  p1 = .5 * ch;
  if ((t = IncompleteGamma(p1, xx, g)) < 0.0) {
    printf("IncompleteGamma \n");
    return (-1);
  }
  p2 = p - t;
  t = p2 * exp(xx * aa + g + p1 - c * log(ch));
  b = t / ch;
  a = 0.5 * t - b * c;
  s1 = (210 + a * (140 + a * (105 + a * (84 + a * (70 + 60 * a))))) / 420;
  s2 = (420 + a * (735 + a * (966 + a * (1141 + 1278 * a)))) / 2520;
  s3 = (210 + a * (462 + a * (707 + 932 * a))) / 2520;
  s4 = (252 + a * (672 + 1182 * a) + c * (294 + a * (889 + 1740 * a))) / 5040;
  s5 = (84 + 264 * a + c * (175 + 606 * a)) / 2520;
  s6 = (120 + c * (346 + 127 * c)) / 5040;
  ch += t * (1 + 0.5 * t * s1 - b * c * (s1 - b * (s2 - b * (s3 - b * (s4 - b * (s5 - b * s6))))));
  if (fabs(q / ch - 1) > e)
    goto l4;
  return (ch);
}
```