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

```asm
<bessj1>:
  endbr64
  movsd  D233(%rip),%xmm2
  movapd %xmm0,%xmm3
  movapd %xmm0,%xmm1
  andpd  D129(%rip),%xmm3
  comisd %xmm3,%xmm2
  jbe    L39
  movapd %xmm0,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  D223(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D155(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D236(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D251(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D256(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D222(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D228(%rip),%xmm1
  addsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D246(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D113(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D128(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D11(%rip),%xmm1
  divsd  %xmm1,%xmm0
  ret
  nopw   0x0(%rax,%rax,1)
L39:
  divsd  %xmm3,%xmm2
  push   %rbx
  pxor   %xmm6,%xmm6
  movapd %xmm3,%xmm5
  sub    $0x40,%rsp
  movsd  D181(%rip),%xmm0
  subsd  D173(%rip),%xmm5
  movapd %xmm2,%xmm4
  mulsd  %xmm2,%xmm4
  mulsd  %xmm4,%xmm0
  addsd  D46(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  subsd  D207(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  D121(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  D101(%rip),%xmm0
  movq   %xmm0,%rbx
  movsd  D146(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  subsd  D27(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  D104(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  subsd  D96(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  D144(%rip),%xmm0
  movsd  %xmm0,0x10(%rsp)
  movsd  D201(%rip),%xmm0
  divsd  %xmm3,%xmm0
  ucomisd %xmm0,%xmm6
  ja     L2
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x8(%rsp)
L20:
  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     L34
  add    $0x40,%rsp
  pop    %rbx
  ret
  nopl   0x0(%rax)
L34:
  xorpd  D84(%rip),%xmm0
  add    $0x40,%rsp
  pop    %rbx
  ret
L2:
  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    L20
```

tool_calls: name: retrieve_data
arguments: label: D233
datatype: f64
name: retrieve_data
arguments: label: D129
datatype: qword[2]
name: retrieve_data
arguments: label: D223
datatype: byte[8]
name: retrieve_data
arguments: label: D155
datatype: f64
name: retrieve_data
arguments: label: D236
datatype: f64
name: retrieve_data
arguments: label: D251
datatype: f64
name: retrieve_data
arguments: label: D256
datatype: f64
name: retrieve_data
arguments: label: D222
datatype: f64
name: retrieve_data
arguments: label: D228
datatype: f64
name: retrieve_data
arguments: label: D246
datatype: f64
name: retrieve_data
arguments: label: D113
datatype: f64
name: retrieve_data
arguments: label: D128
datatype: f64
name: retrieve_data
arguments: label: D11
datatype: f64
name: retrieve_data
arguments: label: D181
datatype: byte[8]
name: retrieve_data
arguments: label: D173
datatype: f64
name: retrieve_data
arguments: label: D46
datatype: byte[8]
name: retrieve_data
arguments: label: D207
datatype: byte[8]
name: retrieve_data
arguments: label: D121
datatype: byte[8]
name: retrieve_data
arguments: label: D101
datatype: f64
name: retrieve_data
arguments: label: D146
datatype: byte[8]
name: retrieve_data
arguments: label: D27
datatype: byte[8]
name: retrieve_data
arguments: label: D104
datatype: byte[8]
name: retrieve_data
arguments: label: D96
datatype: byte[8]
name: retrieve_data
arguments: label: D144
datatype: f64
name: retrieve_data
arguments: label: D201
datatype: f64
name: retrieve_data
arguments: label: D84
datatype: f32[4]

tool: D233: 8

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

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

tool: D155: 15704.4826

tool: D236: 2972611.439

tool: D251: 242396853.1

tool: D256: 7895059235

tool: D222: 72362614232

tool: D228: 376.9991397

tool: D246: 99447.43394

tool: D113: 18583304.74

tool: D128: 2300535178

tool: D11: 144725228442

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

tool: D173: 2.356194491

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

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

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

tool: D101: 1

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

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

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

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

tool: D144: 0.04687499995

tool: D201: 0.636619772

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