user: What is the c source code of the assembly code below:

```asm
<bessel>:
  endbr64
  sub    $0x48,%rsp
  mov    %fs:0x28,%rdx
  mov    %rdx,0x38(%rsp)
  xor    %edx,%edx
  pxor   %xmm1,%xmm1
  comisd %xmm0,%xmm1
  jae    L105
  mov    %esi,%eax
  lea    -0x1(%rsi),%edx
  pxor   %xmm1,%xmm1
  cvtsi2sd %edx,%xmm1
  comisd %xmm1,%xmm0
  jae    L74
  cvttsd2si %xmm0,%esi
  movslq %esi,%rcx
  lea    0x0(,%rcx,8),%r8
  movsd  (%rdi,%rcx,8),%xmm2
  lea    0x1(%rsi),%edx
  cmp    %edx,%eax
  jle    L30
  movapd %xmm0,%xmm4
  movsd  D3(%rip),%xmm3
  movapd %xmm0,%xmm1
  andpd  %xmm3,%xmm1
  movsd  D9(%rip),%xmm5
  ucomisd %xmm1,%xmm5
  jbe    L109
  cvttsd2si %xmm0,%rdx
  pxor   %xmm1,%xmm1
  cvtsi2sd %rdx,%xmm1
  movapd %xmm1,%xmm4
  cmpnlesd %xmm0,%xmm4
  movsd  D38(%rip),%xmm5
  andpd  %xmm5,%xmm4
  subsd  %xmm4,%xmm1
  andnpd %xmm0,%xmm3
  movapd %xmm1,%xmm4
  orpd   %xmm3,%xmm4
L109:
  subsd  %xmm4,%xmm0
  movsd  0x8(%rdi,%r8,1),%xmm1
  subsd  %xmm2,%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm1,%xmm2
  test   %esi,%esi
  jle    L30
  lea    0x2(%rsi),%edx
  cmp    %eax,%edx
  jge    L30
  lea    -0x1(%rsi),%edx
  movslq %edx,%rdx
  lea    0x4(%rcx),%r8
  neg    %rcx
  lea    (%rsp,%rcx,8),%rcx
  jmp    L116
L105:
  movsd  (%rdi),%xmm2
L30:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L43
  movapd %xmm2,%xmm0
  add    $0x48,%rsp
  ret
L74:
  cltq
  movsd  -0x8(%rdi,%rax,8),%xmm2
  jmp    L30
L28:
  movsd  %xmm1,0x8(%rcx,%rdx,8)
  add    $0x1,%rdx
  cmp    %r8,%rdx
  je     L96
L116:
  pxor   %xmm1,%xmm1
  mov    %edx,%r9d
  sub    $0x1,%r9d
  js     L28
  cmp    %edx,%eax
  jle    L28
  movsd  (%rdi,%rdx,8),%xmm1
  subsd  -0x8(%rdi,%rdx,8),%xmm1
  jmp    L28
L96:
  mov    %rsp,%rax
  lea    0x20(%rsp),%rdx
L29:
  movsd  0x8(%rax),%xmm1
  subsd  (%rax),%xmm1
  movsd  %xmm1,(%rax)
  add    $0x8,%rax
  cmp    %rdx,%rax
  jne    L29
  movapd %xmm0,%xmm3
  subsd  D38(%rip),%xmm3
  movapd %xmm0,%xmm1
  mulsd  D52(%rip),%xmm1
  mulsd  %xmm3,%xmm1
  movsd  0x8(%rsp),%xmm8
  movsd  0x10(%rsp),%xmm7
  movapd %xmm8,%xmm3
  addsd  %xmm7,%xmm3
  mulsd  %xmm1,%xmm3
  movapd %xmm3,%xmm4
  addsd  %xmm2,%xmm4
  movapd %xmm8,%xmm6
  subsd  (%rsp),%xmm6
  movapd %xmm7,%xmm5
  subsd  %xmm8,%xmm5
  movsd  %xmm5,0x8(%rsp)
  movsd  0x18(%rsp),%xmm2
  movapd %xmm2,%xmm3
  subsd  %xmm7,%xmm3
  addsd  %xmm1,%xmm1
  divsd  D15(%rip),%xmm1
  movapd %xmm0,%xmm2
  subsd  D74(%rip),%xmm2
  mulsd  %xmm1,%xmm2
  mulsd  %xmm5,%xmm2
  addsd  %xmm4,%xmm2
  cmp    $0x1,%esi
  jle    L30
  mulsd  D11(%rip),%xmm1
  movapd %xmm0,%xmm4
  addsd  D38(%rip),%xmm4
  mulsd  %xmm4,%xmm1
  subsd  %xmm5,%xmm3
  subsd  %xmm6,%xmm5
  addsd  %xmm5,%xmm3
  subsd  D53(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  %xmm0,%xmm3
  addsd  %xmm3,%xmm2
  jmp    L30
L43:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: disassemble_data
arguments: tag: D3
datatype: byte[8]
name: disassemble_data
arguments: tag: D9
datatype: byte[8]
name: disassemble_data
arguments: tag: D38
datatype: f64
name: disassemble_data
arguments: tag: D52
datatype: f64
name: disassemble_data
arguments: tag: D15
datatype: f64
name: disassemble_data
arguments: tag: D74
datatype: f64
name: disassemble_data
arguments: tag: D11
datatype: f64
name: disassemble_data
arguments: tag: D53
datatype: f64

tool: {"D3": [0xff,0xff,0xff,0xff,0xff,0xff,0xff,0x7f]}

tool: {"D9": [0x00,0x00,0x00,0x00,0x00,0x00,0x30,0x43]}

tool: D38: 1

tool: D52: 0.25

tool: D15: 3

tool: D74: 0.5

tool: D11: 0.125

tool: D53: 2

assistant: ```c
double bessel(double *v, int n, double t) {
  int i, iy, k;
  double ans, p, B, d[6];
  if (t <= 0) {
    ans = v[0];
    goto done;
  }
  if (t >= n - 1) {
    ans = v[n - 1];
    goto done;
  }
  p = floor(t);
  iy = (int)t;
  ans = v[iy];
  k = iy + 1;
  if (k >= n)
    goto done;
  p = t - p;
  ans += p * (v[k] - v[iy]);
  if ((iy - 1 < 0) || (iy + 2 >= n))
    goto done;
  k = iy - 2;
  for (i = 0; i < 5; i++) {
    if ((k < 0) || (k + 1 >= n))
      d[i] = 0;
    else
      d[i] = v[k + 1] - v[k];
    k += 1;
  }
  for (i = 0; i < 4; i++)
    d[i] = d[i + 1] - d[i];
  B = 0.25 * p * (p - 1.0);
  ans += B * (d[1] + d[2]);
  if (iy + 2 >= n)
    goto done;
  for (i = 0; i < 3; i++)
    d[i] = d[i + 1] - d[i];
  B = 2.0 * B / 3.0;
  ans += (p - 0.5) * B * d[1];
  if ((iy - 2 < 0) || (iy + 3 > n))
    goto done;
  for (i = 0; i < 2; i++)
    d[i] = d[i + 1] - d[i];
  B = 0.125 * B * (p + 1.0) * (p - 2.0);
  ans += B * (d[0] + d[1]);
done:
  return ans;
}
```