user: 请从给定的汇编指令推导出C源代码:

```asm
<bessj1>:
  endbr64
  movsd  D113(%rip),%xmm2
  movapd %xmm0,%xmm3
  movapd %xmm0,%xmm1
  andpd  D21(%rip),%xmm3
  comisd %xmm3,%xmm2
  jbe    L33
  movapd %xmm0,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  D199(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D65(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D24(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D114(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D10(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D49(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D167(%rip),%xmm1
  addsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D218(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D152(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D42(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D257(%rip),%xmm1
  divsd  %xmm1,%xmm0
  ret
  nopw   0x0(%rax,%rax,1)
L33:
  divsd  %xmm3,%xmm2
  push   %rbx
  pxor   %xmm6,%xmm6
  movapd %xmm3,%xmm5
  sub    $0x40,%rsp
  movsd  D32(%rip),%xmm0
  subsd  D158(%rip),%xmm5
  movapd %xmm2,%xmm4
  mulsd  %xmm2,%xmm4
  mulsd  %xmm4,%xmm0
  addsd  D189(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  subsd  D149(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  D108(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  D18(%rip),%xmm0
  movq   %xmm0,%rbx
  movsd  D181(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  subsd  D226(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  D248(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  subsd  D36(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  D159(%rip),%xmm0
  movsd  %xmm0,0x10(%rsp)
  movsd  D121(%rip),%xmm0
  divsd  %xmm3,%xmm0
  ucomisd %xmm0,%xmm6
  ja     L34
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x8(%rsp)
L1:
  movapd %xmm5,%xmm0
  lea    0x38(%rsp),%rdi
  lea    0x30(%rsp),%rsi
  movsd  %xmm1,0x20(%rsp)
  movsd  %xmm2,0x18(%rsp)
  call   <sincos@plt>
  movsd  0x18(%rsp),%xmm2
  mulsd  0x38(%rsp),%xmm2
  movq   %rbx,%xmm0
  mulsd  0x10(%rsp),%xmm2
  movsd  0x20(%rsp),%xmm1
  pxor   %xmm7,%xmm7
  mulsd  0x30(%rsp),%xmm0
  comisd %xmm1,%xmm7
  subsd  %xmm2,%xmm0
  mulsd  0x8(%rsp),%xmm0
  ja     L39
  add    $0x40,%rsp
  pop    %rbx
  ret
  nopl   0x0(%rax)
L39:
  xorpd  D170(%rip),%xmm0
  add    $0x40,%rsp
  pop    %rbx
  ret
L34:
  movsd  %xmm1,0x28(%rsp)
  movsd  %xmm5,0x20(%rsp)
  movsd  %xmm2,0x18(%rsp)
  call   <sqrt@plt>
  movsd  0x28(%rsp),%xmm1
  movsd  0x20(%rsp),%xmm5
  movsd  0x18(%rsp),%xmm2
  movsd  %xmm0,0x8(%rsp)
  jmp    L1
```

tool_calls: name: parse_data
arguments: data_label: D113
type_guessed: f64
name: parse_data
arguments: data_label: D21
type_guessed: qword[2]
name: parse_data
arguments: data_label: D199
type_guessed: byte[8]
name: parse_data
arguments: data_label: D65
type_guessed: f64
name: parse_data
arguments: data_label: D24
type_guessed: f64
name: parse_data
arguments: data_label: D114
type_guessed: f64
name: parse_data
arguments: data_label: D10
type_guessed: f64
name: parse_data
arguments: data_label: D49
type_guessed: f64
name: parse_data
arguments: data_label: D167
type_guessed: f64
name: parse_data
arguments: data_label: D218
type_guessed: f64
name: parse_data
arguments: data_label: D152
type_guessed: f64
name: parse_data
arguments: data_label: D42
type_guessed: f64
name: parse_data
arguments: data_label: D257
type_guessed: f64
name: parse_data
arguments: data_label: D32
type_guessed: byte[8]
name: parse_data
arguments: data_label: D158
type_guessed: f64
name: parse_data
arguments: data_label: D189
type_guessed: byte[8]
name: parse_data
arguments: data_label: D149
type_guessed: byte[8]
name: parse_data
arguments: data_label: D108
type_guessed: byte[8]
name: parse_data
arguments: data_label: D18
type_guessed: f64
name: parse_data
arguments: data_label: D181
type_guessed: byte[8]
name: parse_data
arguments: data_label: D226
type_guessed: byte[8]
name: parse_data
arguments: data_label: D248
type_guessed: byte[8]
name: parse_data
arguments: data_label: D36
type_guessed: byte[8]
name: parse_data
arguments: data_label: D159
type_guessed: f64
name: parse_data
arguments: data_label: D121
type_guessed: f64
name: parse_data
arguments: data_label: D170
type_guessed: f32[4]

tool: D113: 8

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

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

tool: D65: 15704.4826

tool: D24: 2972611.439

tool: D114: 242396853.1

tool: D10: 7895059235

tool: D49: 72362614232

tool: D167: 376.9991397

tool: D218: 99447.43394

tool: D152: 18583304.74

tool: D42: 2300535178

tool: D257: 144725228442

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

tool: D158: 2.356194491

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

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

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

tool: D18: 1

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

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

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

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

tool: D159: 0.04687499995

tool: D121: 0.636619772

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