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

```asm
<bessj1>:
  endbr64
  push   %r14
  push   %rbx
  sub    $0x38,%rsp
  movapd %xmm0,%xmm3
  movsd  %xmm0,(%rsp)
  movapd %xmm0,%xmm1
  andpd  D166(%rip),%xmm1
  movsd  D185(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L4
  movapd %xmm3,%xmm2
  mulsd  %xmm3,%xmm2
  movapd %xmm2,%xmm0
  mulsd  D167(%rip),%xmm0
  addsd  D19(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D141(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D216(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D257(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D152(%rip),%xmm0
  mulsd  %xmm3,%xmm0
  movapd %xmm2,%xmm1
  addsd  D30(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D41(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D259(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D234(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D20(%rip),%xmm1
  divsd  %xmm1,%xmm0
L43:
  add    $0x38,%rsp
  pop    %rbx
  pop    %r14
  ret
L4:
  movsd  D185(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x10(%rsp)
  movapd %xmm0,%xmm2
  mulsd  %xmm0,%xmm2
  movapd %xmm1,%xmm4
  subsd  D146(%rip),%xmm4
  movq   %xmm4,%rbx
  movapd %xmm2,%xmm0
  mulsd  D262(%rip),%xmm0
  addsd  D237(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D2(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D252(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D209(%rip),%xmm0
  movq   %xmm0,%r14
  movapd %xmm2,%xmm0
  mulsd  D8(%rip),%xmm0
  subsd  D112(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D150(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D172(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D3(%rip),%xmm0
  movsd  %xmm0,0x18(%rsp)
  movsd  D12(%rip),%xmm0
  divsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L0
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x8(%rsp)
L17:
  lea    0x28(%rsp),%rdi
  lea    0x20(%rsp),%rsi
  movq   %rbx,%xmm0
  call   <sincos@plt>
  movq   %r14,%xmm0
  mulsd  0x20(%rsp),%xmm0
  movsd  0x10(%rsp),%xmm1
  mulsd  0x28(%rsp),%xmm1
  mulsd  0x18(%rsp),%xmm1
  subsd  %xmm1,%xmm0
  mulsd  0x8(%rsp),%xmm0
  pxor   %xmm1,%xmm1
  comisd (%rsp),%xmm1
  jbe    L43
  xorpd  D178(%rip),%xmm0
  jmp    L43
L0:
  call   <sqrt@plt>
  movsd  %xmm0,0x8(%rsp)
  jmp    L17
```

tool_calls: name: read_data
arguments: label: D166
type: qword[2]
name: read_data
arguments: label: D185
type: f64
name: read_data
arguments: label: D167
type: byte[8]
name: read_data
arguments: label: D19
type: f64
name: read_data
arguments: label: D141
type: f64
name: read_data
arguments: label: D216
type: f64
name: read_data
arguments: label: D257
type: f64
name: read_data
arguments: label: D152
type: f64
name: read_data
arguments: label: D30
type: f64
name: read_data
arguments: label: D41
type: f64
name: read_data
arguments: label: D259
type: f64
name: read_data
arguments: label: D234
type: f64
name: read_data
arguments: label: D20
type: f64
name: read_data
arguments: label: D146
type: f64
name: read_data
arguments: label: D262
type: byte[8]
name: read_data
arguments: label: D237
type: byte[8]
name: read_data
arguments: label: D2
type: byte[8]
name: read_data
arguments: label: D252
type: byte[8]
name: read_data
arguments: label: D209
type: f64
name: read_data
arguments: label: D8
type: byte[8]
name: read_data
arguments: label: D112
type: byte[8]
name: read_data
arguments: label: D150
type: byte[8]
name: read_data
arguments: label: D172
type: byte[8]
name: read_data
arguments: label: D3
type: f64
name: read_data
arguments: label: D12
type: f64
name: read_data
arguments: label: D178
type: f32[4]

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

tool: D185: 8

tool: {"D167": [0xa0,0x16,0x07,0xc0,0x0d,0x29,0x3e,0xc0]}

tool: D19: 15704.4826

tool: D141: 2972611.439

tool: D216: 242396853.1

tool: D257: 7895059235

tool: D152: 72362614232

tool: D30: 376.9991397

tool: D41: 99447.43394

tool: D259: 18583304.74

tool: D234: 2300535178

tool: D20: 144725228442

tool: D146: 2.356194491

tool: {"D262": [0x14,0x8d,0x4d,0x63,0xf5,0x20,0x90,0xbe]}

tool: {"D237": [0x15,0x0e,0xdf,0x00,0x7c,0x9d,0xc4,0x3e]}

tool: {"D2": [0xd3,0x5c,0xd1,0xa2,0xa0,0x6f,0x02,0x3f]}

tool: {"D252": [0xfd,0x82,0x82,0xf7,0xfa,0xff,0x5d,0x3f]}

tool: D209: 1

tool: {"D8": [0x22,0x81,0x43,0xd5,0xa7,0x65,0x7c,0x3e]}

tool: {"D112": [0xf1,0x78,0x38,0xf1,0xcf,0x9a,0xad,0x3e]}

tool: {"D150": [0x43,0xa6,0xd8,0x14,0x21,0xb8,0xe1,0x3e]}

tool: {"D172": [0x80,0x43,0x57,0x5c,0xea,0x3f,0x2a,0x3f]}

tool: D3: 0.04687499995

tool: D12: 0.636619772

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

assistant: ```c
double bessj1(double x) {
  double ax, z;
  double xx, y, ans, ans1, ans2;
  if ((ax = fabs(x)) < 8.0) {
    y = x * x;
    ans1 = x * (72362614232.0 + y * (-7895059235.0 + y * (242396853.1 + y * (-2972611.439 + y * (15704.48260 + y * (-30.16036606))))));
    ans2 = 144725228442.0 + y * (2300535178.0 + y * (18583304.74 + y * (99447.43394 + y * (376.9991397 + y * 1.0))));
    ans = ans1 / ans2;
  } else {
    z = 8.0 / ax;
    y = z * z;
    xx = ax - 2.356194491;
    ans1 = 1.0 + y * (0.183105e-2 + y * (-0.3516396496e-4 + y * (0.2457520174e-5 + y * (-0.240337019e-6))));
    ans2 = 0.04687499995 + y * (-0.2002690873e-3 + y * (0.8449199096e-5 + y * (-0.88228987e-6 + y * 0.105787412e-6)));
    ans = sqrt(0.636619772 / ax) * (cos(xx) * ans1 - z * sin(xx) * ans2);
    if (x < 0.0)
      ans = -ans;
  }
  return ans;
}
```