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

```asm
<BesselI0>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x20,%rsp
  movsd  %xmm0,-0x98(%rbp)
  pxor   %xmm0,%xmm0
  ucomisd -0x98(%rbp),%xmm0
  jp     L27
  pxor   %xmm0,%xmm0
  ucomisd -0x98(%rbp),%xmm0
  jne    L27
  movsd  D79(%rip),%xmm0
  jmp    L4
L27:
  movsd  -0x98(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0x90(%rbp),%xmm1
  movsd  D61(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D83(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x88(%rbp),%xmm1
  movsd  D174(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x80(%rbp),%xmm1
  movsd  D127(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x78(%rbp),%xmm1
  movsd  D100(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x70(%rbp),%xmm1
  movsd  D157(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x68(%rbp),%xmm1
  movsd  D173(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  movsd  D38(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x58(%rbp),%xmm1
  movsd  D44(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x50(%rbp),%xmm1
  movsd  D85(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x48(%rbp),%xmm1
  movsd  D193(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x40(%rbp),%xmm1
  movsd  D135(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x38(%rbp),%xmm1
  movsd  D183(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x30(%rbp),%xmm1
  movsd  D43(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x28(%rbp),%xmm1
  movsd  D28(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movsd  D2(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x18(%rbp),%xmm1
  movsd  D102(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x90(%rbp),%xmm0
  mulsd  -0x10(%rbp),%xmm0
  movsd  D28(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x20(%rbp),%xmm0
  movq   D70(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  divsd  -0x8(%rbp),%xmm0
L4:
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: read_data
arguments: tag: D79
datatype: f64
name: read_data
arguments: tag: D61
datatype: byte[8]
name: read_data
arguments: tag: D83
datatype: byte[8]
name: read_data
arguments: tag: D174
datatype: byte[8]
name: read_data
arguments: tag: D127
datatype: byte[8]
name: read_data
arguments: tag: D100
datatype: byte[8]
name: read_data
arguments: tag: D157
datatype: byte[8]
name: read_data
arguments: tag: D173
datatype: byte[8]
name: read_data
arguments: tag: D38
datatype: byte[8]
name: read_data
arguments: tag: D44
datatype: f64
name: read_data
arguments: tag: D85
datatype: f64
name: read_data
arguments: tag: D193
datatype: f64
name: read_data
arguments: tag: D135
datatype: f64
name: read_data
arguments: tag: D183
datatype: f64
name: read_data
arguments: tag: D43
datatype: f64
name: read_data
arguments: tag: D28
datatype: f64
name: read_data
arguments: tag: D2
datatype: f64
name: read_data
arguments: tag: D102
datatype: f64
name: read_data
arguments: tag: D70
datatype: f32[2]

tool: D79: 1

tool: {"D61": [0xb7,0x44,0xa1,0x5f,0x28,0x75,0x39,0x3b]}

tool: {"D83": [0x86,0x04,0xa7,0xc1,0x34,0x79,0xe6,0x3b]}

tool: {"D174": [0xf8,0xdd,0x0c,0x15,0x4a,0xa3,0x8b,0x3c]}

tool: {"D127": [0xd3,0xdd,0x80,0x7c,0xd3,0x7e,0x28,0x3d]}

tool: {"D100": [0x77,0x85,0x48,0xa6,0x3a,0x8b,0xc0,0x3d]}

tool: {"D157": [0xb0,0x75,0x91,0xf2,0x38,0x34,0x51,0x3e]}

tool: {"D173": [0xea,0xee,0x2d,0xfb,0x7c,0x77,0xdb,0x3e]}

tool: {"D38": [0xc0,0xcf,0xdf,0x99,0xa4,0x98,0x60,0x3f]}

tool: D44: 0.463076284721

tool: D85: 75.4337328948189

tool: D193: 8307.92541809429

tool: D135: 571661.130563785

tool: D183: 21641557.2361227

tool: D43: 356644482.244025

tool: D28: 1440482982.27235

tool: D2: 3076.46912682801

tool: D102: 3476263.32405882

tool: {"D70": [0.0,-0.0]}

assistant: ```c
double BesselI0(double arg) {
  double numerator, denominator, z, z1, z2, z3, z4, z5, z6, z7, z8, z9, z10, z11, z12, z13, z_1, z_2;
  if (arg == 0.0) {
    return 1.0;
  }
  z = arg * arg;
  z1 = z * 0.210580722890567e-22 + 0.380715242345326e-19;
  z2 = z * z1 + 0.479440257548300e-16;
  z3 = z * z2 + 0.435125971262668e-13;
  z4 = z * z3 + 0.300931127112960e-10;
  z5 = z * z4 + 0.160224679395361e-7;
  z6 = z * z5 + 0.654858370096785e-5;
  z7 = z * z6 + 0.202591084143397e-2;
  z8 = z * z7 + 0.463076284721000e0;
  z9 = z * z8 + 0.754337328948189e2;
  z10 = z * z9 + 0.830792541809429e4;
  z11 = z * z10 + 0.571661130563785e6;
  z12 = z * z11 + 0.216415572361227e8;
  z13 = z * z12 + 0.356644482244025e9;
  numerator = z * z13 + 0.144048298227235e10;
  z_1 = z - 0.307646912682801e4;
  z_2 = z * z_1 + 0.347626332405882e7;
  denominator = z * z_2 - 0.144048298227235e10;
  return -numerator / denominator;
}
```