user: Please show me the C implementation that matches the assembly code listed below:

```asm
<bessi0>:
  endbr64
  sub    $0x28,%rsp
  movss  %xmm0,0x8(%rsp)
  call   <fabs@plt>
  movaps %xmm0,%xmm3
  movss  D171(%rip),%xmm0
  comiss %xmm3,%xmm0
  jbe    L2
  movss  0x8(%rsp),%xmm1
  movsd  D110(%rip),%xmm0
  cvtss2sd %xmm1,%xmm1
  divsd  D103(%rip),%xmm1
  mulsd  %xmm1,%xmm1
  mulsd  %xmm1,%xmm0
  addsd  D168(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D148(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D147(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D82(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D32(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D133(%rip),%xmm0
  add    $0x28,%rsp
  cvtsd2ss %xmm0,%xmm0
  ret
  nopw   0x0(%rax,%rax,1)
L2:
  movsd  D103(%rip),%xmm2
  pxor   %xmm0,%xmm0
  movss  %xmm3,0x14(%rsp)
  cvtss2sd %xmm3,%xmm0
  divsd  %xmm0,%xmm2
  movaps %xmm3,%xmm0
  movsd  %xmm2,0x18(%rsp)
  call   <exp@plt>
  movss  0x14(%rsp),%xmm3
  movsd  %xmm0,0x8(%rsp)
  movaps %xmm3,%xmm0
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm2
  movsd  0x8(%rsp),%xmm4
  movapd %xmm0,%xmm1
  movsd  D161(%rip),%xmm0
  divsd  %xmm1,%xmm4
  mulsd  %xmm2,%xmm0
  subsd  D87(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D105(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D95(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D55(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D11(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D49(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D166(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D7(%rip),%xmm0
  add    $0x28,%rsp
  mulsd  %xmm4,%xmm0
  cvtsd2ss %xmm0,%xmm0
  ret
```

tool_calls: name: read_data
arguments: data_tag: D171
type_infered: f32
name: read_data
arguments: data_tag: D110
type_infered: byte[8]
name: read_data
arguments: data_tag: D103
type_infered: byte[8]
name: read_data
arguments: data_tag: D168
type_infered: byte[8]
name: read_data
arguments: data_tag: D148
type_infered: f64
name: read_data
arguments: data_tag: D147
type_infered: f64
name: read_data
arguments: data_tag: D82
type_infered: f64
name: read_data
arguments: data_tag: D32
type_infered: f64
name: read_data
arguments: data_tag: D133
type_infered: f64
name: read_data
arguments: data_tag: D161
type_infered: byte[8]
name: read_data
arguments: data_tag: D87
type_infered: byte[8]
name: read_data
arguments: data_tag: D105
type_infered: byte[8]
name: read_data
arguments: data_tag: D95
type_infered: byte[8]
name: read_data
arguments: data_tag: D55
type_infered: byte[8]
name: read_data
arguments: data_tag: D11
type_infered: byte[8]
name: read_data
arguments: data_tag: D49
type_infered: byte[8]
name: read_data
arguments: data_tag: D166
type_infered: byte[8]
name: read_data
arguments: data_tag: D7
type_infered: f64

tool: D171: 3.75

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

tool: {"D103": [0x00,0x00,0x00,0x00,0x00,0x00,0x0e,0x40]}

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

tool: D148: 0.2659732

tool: D147: 1.2067492

tool: D82: 3.0899424

tool: D32: 3.5156229

tool: D133: 1

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

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

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

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

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

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

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

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

tool: D7: 0.39894228

assistant: ```c
float bessi0(float x) {
  float ax, ans;
  double y;
  if ((ax = fabs(x)) < 3.75) {
    y = x / 3.75;
    y *= y;
    ans = 1.0 + y * (3.5156229 + y * (3.0899424 + y * (1.2067492 + y * (0.2659732 + y * (0.360768e-1 + y * 0.45813e-2)))));
  } else {
    y = 3.75 / ax;
    ans = (exp(ax) / sqrt(ax)) * (0.39894228 + y * (0.1328592e-1 + y * (0.225319e-2 + y * (-0.157565e-2 + y * (0.916281e-2 + y * (-0.2057706e-1 + y * (0.2635537e-1 + y * (-0.1647633e-1 + y * 0.392377e-2))))))));
  }
  return ans;
}
```