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

```asm
<besselpoly>:
  endbr64
  pxor   %xmm3,%xmm3
  movapd %xmm0,%xmm4
  ucomisd %xmm3,%xmm0
  jp     L81
  jne    L81
  ucomisd %xmm3,%xmm2
  jp     L118
  jne    L118
  movsd  D20(%rip),%xmm0
  addsd  %xmm0,%xmm1
  divsd  %xmm1,%xmm0
  ret
  nopl   (%rax)
L81:
  push   %rbx
  sub    $0x50,%rsp
  comisd %xmm2,%xmm3
  ja     L18
  movsd  D20(%rip),%xmm6
  xor    %ebx,%ebx
  movq   D37(%rip),%xmm8
  movq   D62(%rip),%xmm9
L69:
  movapd %xmm4,%xmm0
  movaps %xmm9,0x40(%rsp)
  movaps %xmm8,0x30(%rsp)
  movsd  %xmm1,0x20(%rsp)
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm6,0x28(%rsp)
  movsd  %xmm2,(%rsp)
  call   <log@plt>
  movsd  (%rsp),%xmm2
  mulsd  %xmm2,%xmm0
  call   <exp@plt>
  movsd  (%rsp),%xmm2
  movsd  0x28(%rsp),%xmm6
  movsd  %xmm0,0x10(%rsp)
  movapd %xmm2,%xmm0
  movsd  %xmm2,0x8(%rsp)
  addsd  %xmm6,%xmm0
  movsd  %xmm6,(%rsp)
  call   <Gamma@plt>
  movsd  0x8(%rsp),%xmm2
  movsd  (%rsp),%xmm6
  xor    %eax,%eax
  movsd  0x20(%rsp),%xmm1
  movsd  0x18(%rsp),%xmm4
  pxor   %xmm3,%xmm3
  movapd 0x40(%rsp),%xmm9
  movsd  0x10(%rsp),%xmm5
  movapd %xmm3,%xmm11
  addsd  %xmm2,%xmm1
  movapd %xmm4,%xmm7
  movsd  D55(%rip),%xmm10
  movapd 0x30(%rsp),%xmm8
  xorpd  %xmm9,%xmm7
  mulsd  %xmm4,%xmm7
  addsd  %xmm6,%xmm1
  mulsd  %xmm1,%xmm0
  divsd  %xmm0,%xmm5
  movapd %xmm3,%xmm0
  jmp    L112
  nopl   0x0(%rax,%rax,1)
L9:
  cmp    $0x3e8,%eax
  je     L141
L112:
  movapd %xmm11,%xmm3
  lea    (%rax,%rax,1),%edx
  pxor   %xmm4,%xmm4
  add    $0x1,%eax
  addsd  %xmm2,%xmm3
  pxor   %xmm11,%xmm11
  movapd %xmm7,%xmm13
  cvtsi2sd %edx,%xmm4
  movapd %xmm5,%xmm12
  addsd  %xmm5,%xmm0
  cvtsi2sd %eax,%xmm11
  addsd  %xmm6,%xmm3
  addsd  %xmm1,%xmm4
  mulsd  %xmm11,%xmm3
  mulsd  %xmm4,%xmm13
  addsd  %xmm10,%xmm4
  mulsd  %xmm4,%xmm3
  movapd %xmm13,%xmm4
  divsd  %xmm3,%xmm4
  mulsd  %xmm4,%xmm5
  movapd %xmm5,%xmm3
  subsd  %xmm12,%xmm3
  divsd  %xmm5,%xmm3
  andpd  %xmm8,%xmm3
  comisd D106(%rip),%xmm3
  ja     L9
L141:
  test   %ebx,%ebx
  je     L140
  xorpd  %xmm9,%xmm0
L140:
  add    $0x50,%rsp
  pop    %rbx
  ret
  nopl   0x0(%rax)
L18:
  movsd  D56(%rip),%xmm6
  movapd %xmm2,%xmm0
  movapd %xmm2,%xmm7
  movq   D37(%rip),%xmm8
  andpd  %xmm8,%xmm0
  movapd %xmm8,%xmm5
  ucomisd %xmm0,%xmm6
  movsd  D20(%rip),%xmm6
  ja     L34
L145:
  ucomisd %xmm2,%xmm7
  jp     L116
  jne    L116
  movq   D62(%rip),%xmm9
  xorpd  %xmm9,%xmm2
  cvttsd2si %xmm2,%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  lea    (%rax,%rdx,1),%ebx
  and    $0x1,%ebx
  sub    %edx,%ebx
  jmp    L69
  nopl   (%rax)
L118:
  movapd %xmm3,%xmm0
  ret
  nopl   (%rax)
L34:
  cvttsd2si %xmm2,%rax
  pxor   %xmm0,%xmm0
  andnpd %xmm2,%xmm5
  cvtsi2sd %rax,%xmm0
  movapd %xmm0,%xmm7
  cmpnlesd %xmm2,%xmm7
  andpd  %xmm6,%xmm7
  subsd  %xmm7,%xmm0
  movapd %xmm0,%xmm7
  orpd   %xmm5,%xmm7
  jmp    L145
  nopl   (%rax)
L116:
  movq   D62(%rip),%xmm9
  xor    %ebx,%ebx
  jmp    L69
```

tool_calls: name: parse_data
arguments: data_label: D20
datatype: f64
name: parse_data
arguments: data_label: D37
datatype: qword
name: parse_data
arguments: data_label: D62
datatype: f32[2]
name: parse_data
arguments: data_label: D55
datatype: f32[2]
name: parse_data
arguments: data_label: D106
datatype: byte[16]
name: parse_data
arguments: data_label: D56
datatype: f32[2]

tool: D20: 1

tool: D37: 9223372036854775807

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

tool: {"D55": [0.0,2.0]}

tool: {"D106": [0x97,0xd4,0x46,0x46,0xf5,0x0e,0x67,0x3c,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0x7f]}

tool: {"D56": [0.0,176.0]}

assistant: ```c
double besselpoly(double a, double lambda, double nu) {
  int m, factor = 0;
  double Sm, relerr, Sol;
  double sum = 0.0;
  if (a == 0.0) {
    if (nu == 0.0)
      return 1.0 / (lambda + 1);
    else
      return 0.0;
  }
  if ((nu < 0) && (floor(nu) == nu)) {
    nu = -nu;
    factor = ((int)nu) % 2;
  }
  Sm = exp(nu * log(a)) / (Gamma(nu + 1) * (lambda + nu + 1));
  m = 0;
  do {
    sum += Sm;
    Sol = Sm;
    Sm *= -a * a * (lambda + nu + 1 + 2 * m) / ((nu + m + 1) * (m + 1) * (lambda + nu + 1 + 2 * m + 2));
    m++;
    relerr = fabs((Sm - Sol) / Sm);
  } while (relerr > 1.0e-17 && m < 1000);
  if (!factor)
    return sum;
  else
    return -sum;
}
```