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

```asm
<PointChi2>:
  endbr64
  movapd %xmm0,%xmm6
  movsd  D75(%rip),%xmm0
  comisd %xmm6,%xmm0
  ja     L29
  sub    $0x58,%rsp
  comisd D450(%rip),%xmm6
  movsd  D94(%rip),%xmm0
  movsd  %xmm6,0x38(%rsp)
  ja     L58
  pxor   %xmm5,%xmm5
  comisd %xmm1,%xmm5
  jae    L58
  movsd  D723(%rip),%xmm5
  movsd  %xmm1,0x20(%rsp)
  mulsd  %xmm1,%xmm5
  movapd %xmm5,%xmm0
  movsd  %xmm5,0x8(%rsp)
  call   <LnGamma@plt>
  movsd  0x8(%rsp),%xmm5
  subsd  D699(%rip),%xmm5
  movsd  %xmm0,0x10(%rsp)
  movsd  0x38(%rsp),%xmm0
  movsd  %xmm5,0x18(%rsp)
  call   <log@plt>
  mulsd  D167(%rip),%xmm0
  movsd  0x20(%rsp),%xmm2
  comisd %xmm0,%xmm2
  jb     L15
  comisd D411(%rip),%xmm2
  jbe    L41
  movsd  0x38(%rsp),%xmm0
  movsd  %xmm2,0x20(%rsp)
  call   <PointNormal@plt>
  movsd  0x20(%rsp),%xmm2
  pxor   %xmm5,%xmm5
  movsd  D28(%rip),%xmm3
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm3
  ucomisd %xmm3,%xmm5
  ja     L135
  movapd %xmm3,%xmm0
  sqrtsd %xmm0,%xmm0
L110:
  mulsd  %xmm1,%xmm0
  movsd  D499(%rip),%xmm1
  movsd  %xmm2,0x20(%rsp)
  addsd  D699(%rip),%xmm0
  subsd  %xmm3,%xmm0
  call   <pow@plt>
  movsd  0x20(%rsp),%xmm2
  movapd %xmm0,%xmm3
  mulsd  %xmm2,%xmm3
  mulsd  D342(%rip),%xmm2
  addsd  D183(%rip),%xmm2
  comisd %xmm2,%xmm3
  jbe    L122
  movsd  D699(%rip),%xmm0
  subsd  0x38(%rsp),%xmm0
  movsd  %xmm3,0x28(%rsp)
  call   <log@plt>
  movsd  0x28(%rsp),%xmm3
  movsd  %xmm0,0x20(%rsp)
  movsd  D723(%rip),%xmm0
  mulsd  %xmm3,%xmm0
  call   <log@plt>
  mulsd  0x18(%rsp),%xmm0
  movsd  0x20(%rsp),%xmm1
  movsd  D78(%rip),%xmm3
  subsd  %xmm0,%xmm1
  addsd  0x10(%rsp),%xmm1
  mulsd  %xmm1,%xmm3
  jmp    L122
  nopl   0x0(%rax)
L114:
  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  D545(%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  D82(%rip),%xmm1
  movsd  D286(%rip),%xmm8
  mulsd  %xmm0,%xmm5
  movsd  D723(%rip),%xmm0
  movsd  D693(%rip),%xmm15
  movapd %xmm2,%xmm9
  movsd  D177(%rip),%xmm10
  movsd  D329(%rip),%xmm11
  movsd  D616(%rip),%xmm12
  movsd  D93(%rip),%xmm13
  movsd  D81(%rip),%xmm14
  mulsd  %xmm5,%xmm0
  movapd %xmm5,%xmm6
  divsd  %xmm3,%xmm6
  movapd %xmm0,%xmm4
  mulsd  %xmm6,%xmm9
  subsd  %xmm9,%xmm4
  mulsd  %xmm4,%xmm1
  addsd  D429(%rip),%xmm1
  mulsd  %xmm4,%xmm8
  addsd  D112(%rip),%xmm8
  mulsd  %xmm4,%xmm12
  addsd  D636(%rip),%xmm12
  mulsd  %xmm4,%xmm13
  addsd  D379(%rip),%xmm13
  mulsd  %xmm4,%xmm1
  mulsd  %xmm4,%xmm8
  addsd  D422(%rip),%xmm8
  mulsd  %xmm4,%xmm12
  addsd  D486(%rip),%xmm12
  mulsd  %xmm2,%xmm13
  addsd  %xmm15,%xmm1
  mulsd  %xmm4,%xmm8
  addsd  D216(%rip),%xmm8
  mulsd  %xmm4,%xmm1
  movapd %xmm13,%xmm7
  addsd  D340(%rip),%xmm1
  mulsd  %xmm4,%xmm8
  mulsd  %xmm4,%xmm1
  addsd  D333(%rip),%xmm1
  addsd  %xmm10,%xmm8
  mulsd  %xmm4,%xmm1
  addsd  %xmm11,%xmm1
  divsd  %xmm10,%xmm1
  movapd %xmm8,%xmm10
  movsd  D272(%rip),%xmm8
  mulsd  %xmm4,%xmm8
  addsd  D41(%rip),%xmm8
  mulsd  %xmm4,%xmm8
  addsd  D228(%rip),%xmm8
  mulsd  %xmm4,%xmm8
  addsd  %xmm8,%xmm11
  movsd  D364(%rip),%xmm8
  mulsd  %xmm4,%xmm8
  addsd  D24(%rip),%xmm8
  divsd  %xmm14,%xmm11
  mulsd  %xmm4,%xmm8
  mulsd  D555(%rip),%xmm4
  addsd  D420(%rip),%xmm8
  mulsd  %xmm2,%xmm8
  addsd  %xmm15,%xmm4
  mulsd  %xmm1,%xmm0
  addsd  D699(%rip),%xmm0
  addsd  %xmm4,%xmm7
  movapd %xmm2,%xmm4
  addsd  %xmm12,%xmm8
  movsd  D708(%rip),%xmm12
  mulsd  D40(%rip),%xmm2
  addsd  D322(%rip),%xmm2
  mulsd  %xmm4,%xmm2
  addsd  D423(%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  D699(%rip),%xmm3
  andpd  D716(%rip),%xmm3
  comisd D331(%rip),%xmm3
  jbe    L58
  movapd %xmm0,%xmm3
L122:
  movsd  D723(%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    L114
  lea    D483(%rip),%rdi
  call   <printf@plt>
  movsd  D94(%rip),%xmm0
L58:
  add    $0x58,%rsp
  ret
  nopl   0x0(%rax)
L29:
  movsd  D94(%rip),%xmm0
  ret
  nopl   0x0(%rax)
L15:
  movsd  0x8(%rsp),%xmm0
  mulsd  D545(%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  D699(%rip),%xmm1
  divsd  %xmm5,%xmm1
  call   <pow@plt>
  pxor   %xmm5,%xmm5
  movapd %xmm0,%xmm1
  subsd  D331(%rip),%xmm1
  movapd %xmm0,%xmm3
  comisd %xmm1,%xmm5
  jbe    L122
  jmp    L58
  nopw   0x0(%rax,%rax,1)
L41:
  movsd  D699(%rip),%xmm0
  subsd  0x38(%rsp),%xmm0
  call   <log@plt>
  movsd  0x18(%rsp),%xmm6
  mulsd  D545(%rip),%xmm6
  addsd  0x10(%rsp),%xmm0
  movsd  D107(%rip),%xmm3
  movsd  %xmm0,0x48(%rsp)
  movsd  %xmm6,0x40(%rsp)
  nopl   0x0(%rax,%rax,1)
L1:
  movsd  D278(%rip),%xmm2
  movsd  D356(%rip),%xmm1
  movsd  %xmm3,0x20(%rsp)
  movsd  D723(%rip),%xmm0
  addsd  %xmm3,%xmm2
  addsd  %xmm3,%xmm1
  mulsd  %xmm3,%xmm0
  addsd  0x48(%rsp),%xmm0
  addsd  0x40(%rsp),%xmm0
  mulsd  %xmm3,%xmm2
  addsd  D568(%rip),%xmm2
  mulsd  %xmm3,%xmm1
  addsd  D699(%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  D699(%rip),%xmm4
  mulsd  %xmm2,%xmm0
  divsd  %xmm1,%xmm0
  subsd  %xmm0,%xmm4
  movapd %xmm3,%xmm0
  addsd  %xmm3,%xmm0
  addsd  D356(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  D499(%rip),%xmm1
  mulsd  %xmm3,%xmm1
  addsd  D262(%rip),%xmm1
  mulsd  %xmm3,%xmm1
  addsd  D568(%rip),%xmm1
  divsd  %xmm2,%xmm1
  subsd  D723(%rip),%xmm0
  subsd  %xmm1,%xmm0
  divsd  %xmm0,%xmm4
  movapd %xmm3,%xmm0
  subsd  %xmm4,%xmm3
  divsd  %xmm3,%xmm0
  subsd  D699(%rip),%xmm0
  andpd  D716(%rip),%xmm0
  subsd  D224(%rip),%xmm0
  comisd %xmm0,%xmm6
  jb     L1
  jmp    L122
L135:
  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    L110
```

tool_calls: name: read_data
arguments: label: D75
data_type: byte[8]
name: read_data
arguments: label: D450
data_type: byte[16]
name: read_data
arguments: label: D94
data_type: f32[2]
name: read_data
arguments: label: D723
data_type: f64
name: read_data
arguments: label: D699
data_type: f32[2]
name: read_data
arguments: label: D167
data_type: byte[8]
name: read_data
arguments: label: D411
data_type: byte[16]
name: read_data
arguments: label: D28
data_type: f64
name: read_data
arguments: label: D499
data_type: f64
name: read_data
arguments: label: D342
data_type: f64
name: read_data
arguments: label: D183
data_type: f32[2]
name: read_data
arguments: label: D78
data_type: f32[2]
name: read_data
arguments: label: D545
data_type: byte[8]
name: read_data
arguments: label: D82
data_type: f32[2]
name: read_data
arguments: label: D286
data_type: f32[2]
name: read_data
arguments: label: D693
data_type: f32[2]
name: read_data
arguments: label: D177
data_type: f32[2]
name: read_data
arguments: label: D329
data_type: f32[2]
name: read_data
arguments: label: D616
data_type: f32[2]
name: read_data
arguments: label: D93
data_type: f32[2]
name: read_data
arguments: label: D81
data_type: f32[2]
name: read_data
arguments: label: D429
data_type: f32[2]
name: read_data
arguments: label: D112
data_type: f32[2]
name: read_data
arguments: label: D636
data_type: f32[2]
name: read_data
arguments: label: D379
data_type: f32[2]
name: read_data
arguments: label: D422
data_type: f32[2]
name: read_data
arguments: label: D486
data_type: f32[2]
name: read_data
arguments: label: D216
data_type: f32[2]
name: read_data
arguments: label: D340
data_type: f32[2]
name: read_data
arguments: label: D333
data_type: f32[2]
name: read_data
arguments: label: D272
data_type: f32[2]
name: read_data
arguments: label: D41
data_type: f32[2]
name: read_data
arguments: label: D228
data_type: f32[2]
name: read_data
arguments: label: D364
data_type: f32[2]
name: read_data
arguments: label: D24
data_type: f32[2]
name: read_data
arguments: label: D555
data_type: f32[2]
name: read_data
arguments: label: D420
data_type: f32[2]
name: read_data
arguments: label: D708
data_type: f32[2]
name: read_data
arguments: label: D40
data_type: f32[2]
name: read_data
arguments: label: D322
data_type: f32[2]
name: read_data
arguments: label: D423
data_type: f32[2]
name: read_data
arguments: label: D716
data_type: qword[2]
name: read_data
arguments: label: D331
data_type: byte[16]
name: read_data
arguments: label: D483
data_type: string
name: read_data
arguments: label: D107
data_type: f64
name: read_data
arguments: label: D278
data_type: f64
name: read_data
arguments: label: D356
data_type: f64
name: read_data
arguments: label: D568
data_type: f64
name: read_data
arguments: label: D262
data_type: f64
name: read_data
arguments: label: D224
data_type: byte[8]

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

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

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

tool: D723: 0.5

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

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

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

tool: D28: 0.222222

tool: D499: 3

tool: D342: 2.2

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

tool: D483: IncompleteGamma 


tool: D107: 0.4

tool: D278: 6.66

tool: D356: 4.67

tool: D568: 6.73

tool: D262: 13.32

tool: {"D224": [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);
}
```