user: What is the c source code of the assembly code below?

```asm
<PointChi2>:
  endbr64
  movapd %xmm0,%xmm6
  movsd  D707(%rip),%xmm0
  comisd %xmm6,%xmm0
  ja     L55
  sub    $0x58,%rsp
  comisd D444(%rip),%xmm6
  movsd  D83(%rip),%xmm0
  movsd  %xmm6,0x38(%rsp)
  ja     L36
  pxor   %xmm5,%xmm5
  comisd %xmm1,%xmm5
  jae    L36
  movsd  D378(%rip),%xmm5
  movsd  %xmm1,0x20(%rsp)
  mulsd  %xmm1,%xmm5
  movapd %xmm5,%xmm0
  movsd  %xmm5,0x8(%rsp)
  call   <LnGamma@plt>
  movsd  0x8(%rsp),%xmm5
  subsd  D709(%rip),%xmm5
  movsd  %xmm0,0x10(%rsp)
  movsd  0x38(%rsp),%xmm0
  movsd  %xmm5,0x18(%rsp)
  call   <log@plt>
  mulsd  D359(%rip),%xmm0
  movsd  0x20(%rsp),%xmm2
  comisd %xmm0,%xmm2
  jb     L20
  comisd D503(%rip),%xmm2
  jbe    L21
  movsd  0x38(%rsp),%xmm0
  movsd  %xmm2,0x20(%rsp)
  call   <PointNormal@plt>
  movsd  0x20(%rsp),%xmm2
  pxor   %xmm5,%xmm5
  movsd  D402(%rip),%xmm3
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm3
  ucomisd %xmm3,%xmm5
  ja     L64
  movapd %xmm3,%xmm0
  sqrtsd %xmm0,%xmm0
L16:
  mulsd  %xmm1,%xmm0
  movsd  D99(%rip),%xmm1
  movsd  %xmm2,0x20(%rsp)
  addsd  D709(%rip),%xmm0
  subsd  %xmm3,%xmm0
  call   <pow@plt>
  movsd  0x20(%rsp),%xmm2
  movapd %xmm0,%xmm3
  mulsd  %xmm2,%xmm3
  mulsd  D588(%rip),%xmm2
  addsd  D143(%rip),%xmm2
  comisd %xmm2,%xmm3
  jbe    L28
  movsd  D709(%rip),%xmm0
  subsd  0x38(%rsp),%xmm0
  movsd  %xmm3,0x28(%rsp)
  call   <log@plt>
  movsd  0x28(%rsp),%xmm3
  movsd  %xmm0,0x20(%rsp)
  movsd  D378(%rip),%xmm0
  mulsd  %xmm3,%xmm0
  call   <log@plt>
  mulsd  0x18(%rsp),%xmm0
  movsd  0x20(%rsp),%xmm1
  movsd  D342(%rip),%xmm3
  subsd  %xmm0,%xmm1
  addsd  0x10(%rsp),%xmm1
  mulsd  %xmm1,%xmm3
  jmp    L28
  nopl   0x0(%rax)
L131:
  movsd  0x38(%rsp),%xmm6
  movsd  %xmm3,0x30(%rsp)
  subsd  %xmm0,%xmm6
  movapd %xmm3,%xmm0
  movsd  %xmm6,0x20(%rsp)
  call   <log@plt>
  movapd %xmm0,%xmm1
  movsd  D565(%rip),%xmm0
  mulsd  0x18(%rsp),%xmm1
  mulsd  0x8(%rsp),%xmm0
  addsd  0x10(%rsp),%xmm0
  addsd  0x28(%rsp),%xmm0
  subsd  %xmm1,%xmm0
  call   <exp@plt>
  movsd  0x20(%rsp),%xmm5
  movsd  0x30(%rsp),%xmm3
  movsd  0x18(%rsp),%xmm2
  movsd  D302(%rip),%xmm1
  movsd  D467(%rip),%xmm8
  mulsd  %xmm0,%xmm5
  movsd  D378(%rip),%xmm0
  movsd  D81(%rip),%xmm15
  movapd %xmm2,%xmm9
  movsd  D515(%rip),%xmm10
  movsd  D290(%rip),%xmm11
  movsd  D367(%rip),%xmm12
  movsd  D182(%rip),%xmm13
  movsd  D440(%rip),%xmm14
  mulsd  %xmm5,%xmm0
  movapd %xmm5,%xmm6
  divsd  %xmm3,%xmm6
  movapd %xmm0,%xmm4
  mulsd  %xmm6,%xmm9
  subsd  %xmm9,%xmm4
  mulsd  %xmm4,%xmm1
  addsd  D480(%rip),%xmm1
  mulsd  %xmm4,%xmm8
  addsd  D304(%rip),%xmm8
  mulsd  %xmm4,%xmm12
  addsd  D472(%rip),%xmm12
  mulsd  %xmm4,%xmm13
  addsd  D666(%rip),%xmm13
  mulsd  %xmm4,%xmm1
  mulsd  %xmm4,%xmm8
  addsd  D541(%rip),%xmm8
  mulsd  %xmm4,%xmm12
  addsd  D656(%rip),%xmm12
  mulsd  %xmm2,%xmm13
  addsd  %xmm15,%xmm1
  mulsd  %xmm4,%xmm8
  addsd  D248(%rip),%xmm8
  mulsd  %xmm4,%xmm1
  movapd %xmm13,%xmm7
  addsd  D13(%rip),%xmm1
  mulsd  %xmm4,%xmm8
  mulsd  %xmm4,%xmm1
  addsd  D90(%rip),%xmm1
  addsd  %xmm10,%xmm8
  mulsd  %xmm4,%xmm1
  addsd  %xmm11,%xmm1
  divsd  %xmm10,%xmm1
  movapd %xmm8,%xmm10
  movsd  D546(%rip),%xmm8
  mulsd  %xmm4,%xmm8
  addsd  D598(%rip),%xmm8
  mulsd  %xmm4,%xmm8
  addsd  D294(%rip),%xmm8
  mulsd  %xmm4,%xmm8
  addsd  %xmm8,%xmm11
  movsd  D728(%rip),%xmm8
  mulsd  %xmm4,%xmm8
  addsd  D453(%rip),%xmm8
  divsd  %xmm14,%xmm11
  mulsd  %xmm4,%xmm8
  mulsd  D350(%rip),%xmm4
  addsd  D86(%rip),%xmm8
  mulsd  %xmm2,%xmm8
  addsd  %xmm15,%xmm4
  mulsd  %xmm1,%xmm0
  addsd  D709(%rip),%xmm0
  addsd  %xmm4,%xmm7
  movapd %xmm2,%xmm4
  addsd  %xmm12,%xmm8
  movsd  D400(%rip),%xmm12
  mulsd  D37(%rip),%xmm2
  addsd  D705(%rip),%xmm2
  mulsd  %xmm4,%xmm2
  addsd  D547(%rip),%xmm2
  divsd  %xmm12,%xmm2
  divsd  %xmm14,%xmm7
  mulsd  %xmm6,%xmm2
  divsd  %xmm12,%xmm8
  subsd  %xmm2,%xmm7
  mulsd  %xmm6,%xmm7
  divsd  %xmm14,%xmm10
  movapd %xmm8,%xmm2
  subsd  %xmm7,%xmm2
  mulsd  %xmm6,%xmm2
  subsd  %xmm2,%xmm11
  mulsd  %xmm6,%xmm11
  subsd  %xmm11,%xmm10
  mulsd  %xmm6,%xmm10
  subsd  %xmm10,%xmm1
  mulsd  %xmm9,%xmm1
  subsd  %xmm1,%xmm0
  mulsd  %xmm5,%xmm0
  addsd  %xmm3,%xmm0
  divsd  %xmm0,%xmm3
  subsd  D709(%rip),%xmm3
  andpd  D652(%rip),%xmm3
  comisd D419(%rip),%xmm3
  jbe    L36
  movapd %xmm0,%xmm3
L28:
  movsd  D378(%rip),%xmm5
  movsd  0x10(%rsp),%xmm2
  movsd  %xmm3,0x20(%rsp)
  movsd  0x8(%rsp),%xmm1
  mulsd  %xmm3,%xmm5
  movapd %xmm5,%xmm0
  movsd  %xmm5,0x28(%rsp)
  call   <IncompleteGamma@plt>
  pxor   %xmm5,%xmm5
  movsd  0x20(%rsp),%xmm3
  comisd %xmm0,%xmm5
  jbe    L131
  lea    D602(%rip),%rdi
  call   <printf@plt>
  movsd  D83(%rip),%xmm0
L36:
  add    $0x58,%rsp
  ret
  nopl   0x0(%rax)
L55:
  movsd  D83(%rip),%xmm0
  ret
  nopl   0x0(%rax)
L20:
  movsd  0x8(%rsp),%xmm0
  mulsd  D565(%rip),%xmm0
  addsd  0x10(%rsp),%xmm0
  call   <exp@plt>
  movsd  0x8(%rsp),%xmm5
  movapd %xmm0,%xmm1
  movsd  0x38(%rsp),%xmm0
  mulsd  %xmm5,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D709(%rip),%xmm1
  divsd  %xmm5,%xmm1
  call   <pow@plt>
  pxor   %xmm5,%xmm5
  movapd %xmm0,%xmm1
  subsd  D419(%rip),%xmm1
  movapd %xmm0,%xmm3
  comisd %xmm1,%xmm5
  jbe    L28
  jmp    L36
  nopw   0x0(%rax,%rax,1)
L21:
  movsd  D709(%rip),%xmm0
  subsd  0x38(%rsp),%xmm0
  call   <log@plt>
  movsd  0x18(%rsp),%xmm6
  mulsd  D565(%rip),%xmm6
  addsd  0x10(%rsp),%xmm0
  movsd  D153(%rip),%xmm3
  movsd  %xmm0,0x48(%rsp)
  movsd  %xmm6,0x40(%rsp)
  nopl   0x0(%rax,%rax,1)
L53:
  movsd  D604(%rip),%xmm2
  movsd  D578(%rip),%xmm1
  movsd  %xmm3,0x20(%rsp)
  movsd  D378(%rip),%xmm0
  addsd  %xmm3,%xmm2
  addsd  %xmm3,%xmm1
  mulsd  %xmm3,%xmm0
  addsd  0x48(%rsp),%xmm0
  addsd  0x40(%rsp),%xmm0
  mulsd  %xmm3,%xmm2
  addsd  D498(%rip),%xmm2
  mulsd  %xmm3,%xmm1
  addsd  D709(%rip),%xmm1
  mulsd  %xmm3,%xmm2
  movsd  %xmm1,0x30(%rsp)
  movsd  %xmm2,0x28(%rsp)
  call   <exp@plt>
  movsd  0x28(%rsp),%xmm2
  movsd  0x30(%rsp),%xmm1
  pxor   %xmm6,%xmm6
  movsd  0x20(%rsp),%xmm3
  movsd  D709(%rip),%xmm4
  mulsd  %xmm2,%xmm0
  divsd  %xmm1,%xmm0
  subsd  %xmm0,%xmm4
  movapd %xmm3,%xmm0
  addsd  %xmm3,%xmm0
  addsd  D578(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  D99(%rip),%xmm1
  mulsd  %xmm3,%xmm1
  addsd  D520(%rip),%xmm1
  mulsd  %xmm3,%xmm1
  addsd  D498(%rip),%xmm1
  divsd  %xmm2,%xmm1
  subsd  D378(%rip),%xmm0
  subsd  %xmm1,%xmm0
  divsd  %xmm0,%xmm4
  movapd %xmm3,%xmm0
  subsd  %xmm4,%xmm3
  divsd  %xmm3,%xmm0
  subsd  D709(%rip),%xmm0
  andpd  D652(%rip),%xmm0
  subsd  D27(%rip),%xmm0
  comisd %xmm0,%xmm6
  jb     L53
  jmp    L28
L64:
  movsd  %xmm0,0x28(%rsp)
  movapd %xmm3,%xmm0
  movsd  %xmm2,0x30(%rsp)
  movsd  %xmm3,0x20(%rsp)
  call   <sqrt@plt>
  movsd  0x30(%rsp),%xmm2
  movsd  0x28(%rsp),%xmm1
  movsd  0x20(%rsp),%xmm3
  jmp    L16
```

tool_calls: name: retrieve_data
arguments: label: D707
datatype: byte[8]
name: retrieve_data
arguments: label: D444
datatype: byte[16]
name: retrieve_data
arguments: label: D83
datatype: f32[2]
name: retrieve_data
arguments: label: D378
datatype: f64
name: retrieve_data
arguments: label: D709
datatype: f32[2]
name: retrieve_data
arguments: label: D359
datatype: byte[8]
name: retrieve_data
arguments: label: D503
datatype: byte[16]
name: retrieve_data
arguments: label: D402
datatype: f64
name: retrieve_data
arguments: label: D99
datatype: f64
name: retrieve_data
arguments: label: D588
datatype: f64
name: retrieve_data
arguments: label: D143
datatype: f32[2]
name: retrieve_data
arguments: label: D342
datatype: f32[2]
name: retrieve_data
arguments: label: D565
datatype: byte[8]
name: retrieve_data
arguments: label: D302
datatype: f32[2]
name: retrieve_data
arguments: label: D467
datatype: f32[2]
name: retrieve_data
arguments: label: D81
datatype: f32[2]
name: retrieve_data
arguments: label: D515
datatype: f32[2]
name: retrieve_data
arguments: label: D290
datatype: f32[2]
name: retrieve_data
arguments: label: D367
datatype: f32[2]
name: retrieve_data
arguments: label: D182
datatype: f32[2]
name: retrieve_data
arguments: label: D440
datatype: f32[2]
name: retrieve_data
arguments: label: D480
datatype: f32[2]
name: retrieve_data
arguments: label: D304
datatype: f32[2]
name: retrieve_data
arguments: label: D472
datatype: f32[2]
name: retrieve_data
arguments: label: D666
datatype: f32[2]
name: retrieve_data
arguments: label: D541
datatype: f32[2]
name: retrieve_data
arguments: label: D656
datatype: f32[2]
name: retrieve_data
arguments: label: D248
datatype: f32[2]
name: retrieve_data
arguments: label: D13
datatype: f32[2]
name: retrieve_data
arguments: label: D90
datatype: f32[2]
name: retrieve_data
arguments: label: D546
datatype: f32[2]
name: retrieve_data
arguments: label: D598
datatype: f32[2]
name: retrieve_data
arguments: label: D294
datatype: f32[2]
name: retrieve_data
arguments: label: D728
datatype: f32[2]
name: retrieve_data
arguments: label: D453
datatype: f32[2]
name: retrieve_data
arguments: label: D350
datatype: f32[2]
name: retrieve_data
arguments: label: D86
datatype: f32[2]
name: retrieve_data
arguments: label: D400
datatype: f32[2]
name: retrieve_data
arguments: label: D37
datatype: f32[2]
name: retrieve_data
arguments: label: D705
datatype: f32[2]
name: retrieve_data
arguments: label: D547
datatype: f32[2]
name: retrieve_data
arguments: label: D652
datatype: qword[2]
name: retrieve_data
arguments: label: D419
datatype: byte[16]
name: retrieve_data
arguments: label: D602
datatype: string
name: retrieve_data
arguments: label: D153
datatype: f64
name: retrieve_data
arguments: label: D604
datatype: f64
name: retrieve_data
arguments: label: D578
datatype: f64
name: retrieve_data
arguments: label: D498
datatype: f64
name: retrieve_data
arguments: label: D520
datatype: f64
name: retrieve_data
arguments: label: D27
datatype: byte[8]

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

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

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

tool: D378: 0.5

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

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

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

tool: D402: 0.222222

tool: D99: 3

tool: D588: 2.2

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

tool: D602: IncompleteGamma 


tool: D153: 0.4

tool: D604: 6.66

tool: D578: 4.67

tool: D498: 6.73

tool: D520: 13.32

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

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);
}
```