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

```asm
<bessi0>:
  endbr64
  movsd  D167(%rip),%xmm2
  movapd %xmm0,%xmm1
  andpd  D52(%rip),%xmm1
  comisd %xmm1,%xmm2
  jbe    L5
  movapd %xmm0,%xmm1
  movsd  D62(%rip),%xmm0
  divsd  %xmm2,%xmm1
  mulsd  %xmm1,%xmm1
  mulsd  %xmm1,%xmm0
  addsd  D1(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D151(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D37(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D22(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D178(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D74(%rip),%xmm0
  ret
  nopl   0x0(%rax,%rax,1)
L5:
  divsd  %xmm1,%xmm2
  sub    $0x18,%rsp
  movapd %xmm1,%xmm0
  movsd  %xmm1,(%rsp)
  movsd  %xmm2,0x8(%rsp)
  call   <exp@plt>
  movsd  (%rsp),%xmm1
  movsd  0x8(%rsp),%xmm2
  movapd %xmm0,%xmm3
  pxor   %xmm0,%xmm0
  ucomisd %xmm1,%xmm0
  ja     L22
  sqrtsd %xmm1,%xmm1
L20:
  movsd  D116(%rip),%xmm0
  divsd  %xmm1,%xmm3
  mulsd  %xmm2,%xmm0
  subsd  D122(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D97(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D0(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D67(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D179(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D124(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D55(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D53(%rip),%xmm0
  add    $0x18,%rsp
  mulsd  %xmm3,%xmm0
  ret
L22:
  movapd %xmm1,%xmm0
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm3,(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm2
  movsd  (%rsp),%xmm3
  movapd %xmm0,%xmm1
  jmp    L20
```

tool_calls: name: read_data
arguments: data_tag: D167
data_type: f64
name: read_data
arguments: data_tag: D52
data_type: qword[2]
name: read_data
arguments: data_tag: D62
data_type: byte[8]
name: read_data
arguments: data_tag: D1
data_type: byte[8]
name: read_data
arguments: data_tag: D151
data_type: f64
name: read_data
arguments: data_tag: D37
data_type: f64
name: read_data
arguments: data_tag: D22
data_type: f64
name: read_data
arguments: data_tag: D178
data_type: f64
name: read_data
arguments: data_tag: D74
data_type: f64
name: read_data
arguments: data_tag: D116
data_type: byte[8]
name: read_data
arguments: data_tag: D122
data_type: byte[8]
name: read_data
arguments: data_tag: D97
data_type: byte[8]
name: read_data
arguments: data_tag: D0
data_type: byte[8]
name: read_data
arguments: data_tag: D67
data_type: byte[8]
name: read_data
arguments: data_tag: D179
data_type: byte[8]
name: read_data
arguments: data_tag: D124
data_type: byte[8]
name: read_data
arguments: data_tag: D55
data_type: byte[8]
name: read_data
arguments: data_tag: D53
data_type: f64

tool: D167: 3.75

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

tool: {"D62": [0x74,0x48,0xc5,0x5a,0xd7,0xc3,0x72,0x3f]}

tool: {"D1": [0xfd,0xe5,0x49,0x88,0xa8,0x78,0xa2,0x3f]}

tool: D151: 0.2659732

tool: D37: 1.2067492

tool: D22: 3.0899424

tool: D178: 3.5156229

tool: D74: 1

tool: {"D116": [0x88,0x7f,0x3b,0xfd,0x5e,0x12,0x70,0x3f]}

tool: {"D122": [0x55,0x4c,0x08,0xca,0x2b,0xdf,0x90,0x3f]}

tool: {"D97": [0xe3,0x5a,0xe5,0xf0,0xe6,0xfc,0x9a,0x3f]}

tool: {"D0": [0x16,0xc7,0xfd,0x1e,0x27,0x12,0x95,0x3f]}

tool: {"D67": [0x54,0xd1,0x50,0x8a,0xf3,0xc3,0x82,0x3f]}

tool: {"D179": [0x4a,0x89,0x13,0x4e,0xc1,0xd0,0x59,0x3f]}

tool: {"D124": [0xdb,0xd5,0x92,0x2b,0x48,0x75,0x62,0x3f]}

tool: {"D55": [0x1c,0x9d,0x2d,0xff,0xa5,0x35,0x8b,0x3f]}

tool: D53: 0.39894228

assistant: ```c
double bessi0(double x) {
  double ax, ans;
  static double P[7] = {1.0, 3.5156229, 3.0899424, 1.2067492, 0.2659732, 0.360768e-1, 0.45813e-2};
  static double Q[9] = {0.39894228, 0.1328592e-1, 0.225319e-2, -0.157565e-2, 0.916281e-2, -0.2057706e-1, 0.2635537e-1, -0.1647633e-1, 0.392377e-2};
  double y;
  if ((ax = fabs(x)) < 3.75) {
    y = x / 3.75;
    y *= y;
    ans = P[0] + y * (P[1] + y * (P[2] + y * (P[3] + y * (P[4] + y * (P[5] + y * P[6])))));
  } else {
    y = 3.75 / ax;
    ans = (exp(ax) / sqrt(ax)) * (Q[0] + y * (Q[1] + y * (Q[2] + y * (Q[3] + y * (Q[4] + y * (Q[5] + y * (Q[6] + y * (Q[7] + y * Q[8]))))))));
  }
  return ans;
}
```