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

```asm
<sph_bessel_asym>:
  endbr64
  movsd  D47(%rip),%xmm1
  mov    %rsi,%rcx
  mov    %rdx,%rsi
  movsd  %xmm1,(%rcx)
  test   %edi,%edi
  js     L35
  mov    %edi,%edx
  movsd  D31(%rip),%xmm3
  lea    0x8(%rcx),%rax
  movapd %xmm1,%xmm2
  lea    0x10(%rcx,%rdx,8),%rdx
  nopl   0x0(%rax)
L64:
  mulsd  %xmm0,%xmm1
  addsd  %xmm3,%xmm2
  add    $0x8,%rax
  divsd  %xmm2,%xmm1
  movsd  %xmm1,-0x8(%rax)
  cmp    %rdx,%rax
  jne    L64
L35:
  movq   $0x0,(%rsi)
  test   %edi,%edi
  jle    L83
  mov    %rsi,%rax
  sub    %rcx,%rax
  cmp    $0x20,%rax
  jbe    L44
  lea    -0x1(%rdi),%eax
  cmp    $0x2,%eax
  jbe    L44
  mov    %edi,%r8d
  lea    0x10(%rcx),%rax
  lea    0x8(%rsi),%rdx
  movdqa D5(%rip),%xmm4
  shr    $0x2,%r8d
  movdqa D33(%rip),%xmm6
  movdqa D46(%rip),%xmm5
  sub    $0x1,%r8d
  shl    $0x5,%r8
  lea    0x30(%rcx,%r8,1),%r8
  nopl   0x0(%rax)
L97:
  movdqa %xmm4,%xmm0
  movupd (%rax),%xmm7
  add    $0x20,%rax
  paddd  %xmm6,%xmm4
  movdqa %xmm0,%xmm3
  movupd -0x10(%rax),%xmm8
  movdqa %xmm0,%xmm2
  add    $0x20,%rdx
  paddd  %xmm5,%xmm3
  pslld  $0x1,%xmm2
  pshufd $0xee,%xmm0,%xmm1
  cvtdq2pd %xmm0,%xmm0
  cvtdq2pd %xmm1,%xmm1
  mulpd  %xmm7,%xmm1
  pshufd $0xee,%xmm3,%xmm7
  paddd  %xmm5,%xmm2
  cvtdq2pd %xmm7,%xmm7
  mulpd  %xmm8,%xmm7
  cvtdq2pd %xmm3,%xmm3
  subpd  %xmm7,%xmm1
  pshufd $0xee,%xmm2,%xmm7
  cvtdq2pd %xmm2,%xmm2
  cvtdq2pd %xmm7,%xmm7
  divpd  %xmm7,%xmm1
  movupd -0x30(%rax),%xmm7
  mulpd  %xmm7,%xmm0
  movupd -0x20(%rax),%xmm7
  mulpd  %xmm7,%xmm3
  subpd  %xmm3,%xmm0
  divpd  %xmm2,%xmm0
  movups %xmm1,-0x10(%rdx)
  movups %xmm0,-0x20(%rdx)
  cmp    %r8,%rax
  jne    L97
  mov    %edi,%eax
  and    $0xfffffffc,%eax
  lea    0x1(%rax),%r10d
  cmp    %eax,%edi
  je     L83
  lea    0x2(%rax),%r11d
  pxor   %xmm0,%xmm0
  pxor   %xmm1,%xmm1
  movslq %r10d,%r9
  cvtsi2sd %r11d,%xmm0
  movsd  0x8(%rcx,%r9,8),%xmm2
  lea    (%r10,%r10,1),%edx
  cvtsi2sd %r10d,%xmm1
  mulsd  -0x8(%rcx,%r9,8),%xmm1
  lea    0x1(%rdx),%r10d
  lea    0x0(,%r9,8),%r8
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  pxor   %xmm2,%xmm2
  cvtsi2sd %r10d,%xmm2
  divsd  %xmm2,%xmm1
  movsd  %xmm1,(%rsi,%r9,8)
  cmp    %edi,%r11d
  jg     L83
  lea    0x3(%rax),%r10d
  mulsd  (%rcx,%r9,8),%xmm0
  pxor   %xmm1,%xmm1
  movsd  0x10(%rcx,%r8,1),%xmm2
  cvtsi2sd %r10d,%xmm1
  lea    0x3(%rdx),%r9d
  mulsd  %xmm1,%xmm2
  subsd  %xmm2,%xmm0
  pxor   %xmm2,%xmm2
  cvtsi2sd %r9d,%xmm2
  divsd  %xmm2,%xmm0
  movsd  %xmm0,0x8(%rsi,%r8,1)
  cmp    %r10d,%edi
  jl     L83
  add    $0x4,%eax
  pxor   %xmm0,%xmm0
  add    $0x5,%edx
  mulsd  0x8(%rcx,%r8,1),%xmm1
  cvtsi2sd %eax,%xmm0
  mulsd  0x18(%rcx,%r8,1),%xmm0
  subsd  %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2sd %edx,%xmm0
  divsd  %xmm0,%xmm1
  movsd  %xmm1,0x10(%rsi,%r8,1)
  ret
  nopl   (%rax)
L44:
  lea    0x2(%rdi),%r8d
  mov    $0x3,%edx
  mov    $0x2,%eax
  xchg   %ax,%ax
L84:
  lea    -0x1(%rax),%edi
  pxor   %xmm0,%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mulsd  (%rcx,%rax,8),%xmm1
  cvtsi2sd %edi,%xmm0
  mulsd  -0x10(%rcx,%rax,8),%xmm0
  subsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sd %edx,%xmm1
  add    $0x2,%edx
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rsi,%rax,8)
  add    $0x1,%rax
  cmp    %r8,%rax
  jne    L84
L83:
  ret
```

tool_calls: name: parse_data
arguments: data_tag: D47
type_infered: f64
name: parse_data
arguments: data_tag: D31
type_infered: f64
name: parse_data
arguments: data_tag: D5
type_infered: byte[16]
name: parse_data
arguments: data_tag: D33
type_infered: byte[16]
name: parse_data
arguments: data_tag: D46
type_infered: byte[16]

tool: D47: 1

tool: D31: 2

tool: {"D5": [0x01,0x00,0x00,0x00,0x02,0x00,0x00,0x00,0x03,0x00,0x00,0x00,0x04,0x00,0x00,0x00]}

tool: {"D33": [0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00]}

tool: {"D46": [0x01,0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x01,0x00,0x00,0x00]}

assistant: ```c
void sph_bessel_asym(double x, int l, double *sb, double *dsb) {
  int i;
  double f = 1.0;
  sb[0] = 1.0;
  for (i = 1; i <= l + 1; i++) {
    f += 2.0;
    sb[i] = sb[i - 1] * x / f;
  }
  dsb[0] = 0.0;
  for (i = 1; i <= l; i++) {
    dsb[i] = ((double)i * sb[i - 1] - (double)(i + 1) * sb[i + 1]) / (double)(2 * i + 1);
  }
}
```