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

```asm
<bessel>:
  endbr64
  sub    $0x48,%rsp
  pxor   %xmm7,%xmm7
  mov    %rdi,%rdx
  mov    %fs:0x28,%rcx
  mov    %rcx,0x38(%rsp)
  xor    %ecx,%ecx
  comisd %xmm0,%xmm7
  jae    L117
  movslq %esi,%rax
  pxor   %xmm2,%xmm2
  lea    -0x1(%rax),%ecx
  cvtsi2sd %ecx,%xmm2
  comisd %xmm2,%xmm0
  jae    L37
  cvttsd2si %xmm0,%ecx
  movslq %ecx,%rsi
  movsd  (%rdx,%rsi,8),%xmm1
  lea    0x0(,%rsi,8),%rdi
  lea    0x1(%rcx),%esi
  cmp    %esi,%eax
  jle    L56
  movsd  D23(%rip),%xmm3
  movsd  D60(%rip),%xmm5
  movapd %xmm0,%xmm2
  movapd %xmm0,%xmm4
  andpd  %xmm3,%xmm2
  ucomisd %xmm2,%xmm5
  ja     L132
L168:
  movsd  0x8(%rdx,%rdi,1),%xmm3
  subsd  %xmm4,%xmm0
  movapd %xmm3,%xmm2
  subsd  %xmm1,%xmm2
  movapd %xmm2,%xmm6
  mulsd  %xmm0,%xmm6
  addsd  %xmm1,%xmm6
  test   %ecx,%ecx
  jle    L106
  lea    0x2(%rcx),%r8d
  cmp    %eax,%r8d
  jge    L106
  movapd %xmm7,%xmm4
  cmp    $0x1,%ecx
  je     L20
  cmp    %ecx,%eax
  jge    L13
L20:
  movsd  %xmm4,(%rsp)
  movapd %xmm7,%xmm4
  cmp    %ecx,%eax
  jle    L76
  movapd %xmm1,%xmm4
  subsd  -0x8(%rdx,%rdi,1),%xmm4
L76:
  movsd  %xmm4,0x8(%rsp)
  movsd  %xmm2,0x10(%rsp)
  cmp    $0xffffffff,%esi
  je     L176
  movsd  0x10(%rdx,%rdi,1),%xmm1
  lea    0x3(%rcx),%esi
  subsd  %xmm3,%xmm1
  movsd  %xmm1,0x18(%rsp)
  cmp    %esi,%eax
  jg     L150
L123:
  movapd %xmm0,%xmm4
  movsd  D21(%rip),%xmm1
  movsd  D44(%rip),%xmm5
  movsd  %xmm7,0x20(%rsp)
  movupd 0x8(%rsp),%xmm3
  subpd  (%rsp),%xmm3
  mulsd  %xmm0,%xmm1
  subsd  %xmm5,%xmm4
  movupd 0x18(%rsp),%xmm2
  subpd  0x10(%rsp),%xmm2
  movaps %xmm3,(%rsp)
  movaps %xmm2,0x10(%rsp)
  movapd %xmm2,%xmm7
  mulsd  %xmm4,%xmm1
  movapd %xmm3,%xmm4
  unpckhpd %xmm4,%xmm4
  addsd  %xmm2,%xmm4
  unpckhpd %xmm2,%xmm2
  mulsd  %xmm1,%xmm4
  addsd  %xmm1,%xmm1
  addsd  %xmm6,%xmm4
  movupd 0x8(%rsp),%xmm6
  subpd  %xmm3,%xmm6
  movapd %xmm2,%xmm3
  movapd %xmm1,%xmm2
  divsd  D24(%rip),%xmm2
  movapd %xmm0,%xmm1
  subsd  D22(%rip),%xmm1
  subsd  %xmm7,%xmm3
  movapd %xmm6,%xmm7
  mulsd  %xmm2,%xmm1
  unpckhpd %xmm7,%xmm7
  mulsd  %xmm7,%xmm1
  addsd  %xmm4,%xmm1
  cmp    $0x1,%ecx
  je     L56
  mulsd  D54(%rip),%xmm2
  addsd  %xmm0,%xmm5
  subsd  D88(%rip),%xmm0
  movapd %xmm7,%xmm4
  subsd  %xmm7,%xmm3
  subsd  %xmm6,%xmm4
  mulsd  %xmm5,%xmm2
  addsd  %xmm4,%xmm3
  mulsd  %xmm2,%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm0,%xmm1
  jmp    L56
L37:
  movsd  -0x8(%rdi,%rax,8),%xmm1
L56:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L61
  movapd %xmm1,%xmm0
  add    $0x48,%rsp
  ret
  nopl   0x0(%rax,%rax,1)
L117:
  movsd  (%rdi),%xmm1
  jmp    L56
  xchg   %ax,%ax
L132:
  cvttsd2si %xmm0,%r8
  pxor   %xmm2,%xmm2
  movsd  D44(%rip),%xmm5
  andnpd %xmm0,%xmm3
  cvtsi2sd %r8,%xmm2
  movapd %xmm2,%xmm4
  cmpnlesd %xmm0,%xmm4
  andpd  %xmm5,%xmm4
  subsd  %xmm4,%xmm2
  movapd %xmm2,%xmm4
  orpd   %xmm3,%xmm4
  jmp    L168
  nopl   0x0(%rax,%rax,1)
L106:
  movapd %xmm6,%xmm1
  jmp    L56
  cs nopw 0x0(%rax,%rax,1)
L176:
  movq   $0x0,0x18(%rsp)
L150:
  movsd  0x18(%rdx,%rdi,1),%xmm1
  movapd %xmm1,%xmm7
  subsd  0x10(%rdx,%rdi,1),%xmm7
  jmp    L123
  xchg   %ax,%ax
L13:
  movsd  -0x8(%rdx,%rdi,1),%xmm4
  subsd  -0x10(%rdx,%rdi,1),%xmm4
  jmp    L20
L61:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: disassemble_data
arguments: tag: D23
data_type: byte[8]
name: disassemble_data
arguments: tag: D60
data_type: byte[8]
name: disassemble_data
arguments: tag: D21
data_type: f64
name: disassemble_data
arguments: tag: D44
data_type: f64
name: disassemble_data
arguments: tag: D24
data_type: f64
name: disassemble_data
arguments: tag: D22
data_type: f64
name: disassemble_data
arguments: tag: D54
data_type: f64
name: disassemble_data
arguments: tag: D88
data_type: f64

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

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

tool: D21: 0.25

tool: D44: 1

tool: D24: 3

tool: D22: 0.5

tool: D54: 0.125

tool: D88: 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;
}
```