user: Demonstrate the C code that generates the assembly output below:

```asm
<bessel>:
  endbr64
  sub    $0x48,%rsp
  pxor   %xmm1,%xmm1
  mov    %fs:0x28,%rax
  mov    %rax,0x38(%rsp)
  xor    %eax,%eax
  comisd %xmm0,%xmm1
  jae    L20
  lea    -0x1(%rsi),%eax
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
  comisd %xmm2,%xmm0
  jae    L54
  cvttsd2si %xmm0,%r8d
  movslq %r8d,%rdx
  lea    0x1(%r8),%eax
  movsd  (%rdi,%rdx,8),%xmm3
  lea    0x0(,%rdx,8),%r9
  cmp    %eax,%esi
  jle    L23
  movsd  D88(%rip),%xmm4
  movsd  D60(%rip),%xmm6
  movapd %xmm0,%xmm2
  movapd %xmm0,%xmm5
  andpd  %xmm4,%xmm2
  ucomisd %xmm2,%xmm6
  ja     L65
L134:
  movsd  0x8(%rdi,%r9,1),%xmm2
  subsd  %xmm5,%xmm0
  subsd  %xmm3,%xmm2
  mulsd  %xmm0,%xmm2
  addsd  %xmm2,%xmm3
  test   %r8d,%r8d
  jle    L23
  lea    0x2(%r8),%eax
  cmp    %esi,%eax
  jge    L23
  lea    0x3(%rdx),%rcx
  mov    %rsp,%rdx
  lea    -0x2(%r8),%eax
  mov    %rdx,%r10
  cltq
  sub    %r9,%r10
L3:
  movapd %xmm1,%xmm2
  cmp    $0xffffffff,%eax
  je     L103
  lea    0x1(%rax),%r9d
  cmp    %r9d,%esi
  jle    L103
  movsd  0x8(%rdi,%rax,8),%xmm2
  subsd  (%rdi,%rax,8),%xmm2
L103:
  movsd  %xmm2,0x10(%r10,%rax,8)
  add    $0x1,%rax
  cmp    %rcx,%rax
  jne    L3
  mov    %rdx,%rax
  lea    0x20(%rsp),%rdx
L38:
  movsd  0x8(%rax),%xmm1
  subsd  (%rax),%xmm1
  add    $0x8,%rax
  movsd  %xmm1,-0x8(%rax)
  cmp    %rdx,%rax
  jne    L38
  movsd  D46(%rip),%xmm1
  movapd %xmm0,%xmm2
  movsd  D80(%rip),%xmm5
  movsd  0x8(%rsp),%xmm4
  movsd  0x10(%rsp),%xmm8
  mulsd  %xmm0,%xmm1
  subsd  %xmm5,%xmm2
  movapd %xmm8,%xmm6
  movapd %xmm4,%xmm7
  subsd  (%rsp),%xmm7
  subsd  %xmm4,%xmm6
  mulsd  %xmm2,%xmm1
  movapd %xmm4,%xmm2
  addsd  %xmm8,%xmm2
  mulsd  %xmm1,%xmm2
  addsd  %xmm1,%xmm1
  divsd  D51(%rip),%xmm1
  addsd  %xmm3,%xmm2
  movsd  0x18(%rsp),%xmm3
  movapd %xmm3,%xmm4
  movapd %xmm0,%xmm3
  subsd  D33(%rip),%xmm3
  subsd  %xmm8,%xmm4
  mulsd  %xmm1,%xmm3
  mulsd  %xmm6,%xmm3
  addsd  %xmm2,%xmm3
  cmp    $0x1,%r8d
  je     L23
  mulsd  D42(%rip),%xmm1
  addsd  %xmm0,%xmm5
  subsd  D78(%rip),%xmm0
  movapd %xmm6,%xmm2
  subsd  %xmm6,%xmm4
  subsd  %xmm7,%xmm2
  mulsd  %xmm5,%xmm1
  addsd  %xmm4,%xmm2
  mulsd  %xmm1,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm3
  jmp    L23
  nopl   0x0(%rax)
L54:
  movslq %esi,%rsi
  movsd  -0x8(%rdi,%rsi,8),%xmm3
L23:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L26
  movapd %xmm3,%xmm0
  add    $0x48,%rsp
  ret
  nopw   0x0(%rax,%rax,1)
L20:
  movsd  (%rdi),%xmm3
  jmp    L23
  xchg   %ax,%ax
L65:
  cvttsd2si %xmm0,%rax
  pxor   %xmm2,%xmm2
  movsd  D80(%rip),%xmm6
  andnpd %xmm0,%xmm4
  cvtsi2sd %rax,%xmm2
  movapd %xmm2,%xmm5
  cmpnlesd %xmm0,%xmm5
  andpd  %xmm6,%xmm5
  subsd  %xmm5,%xmm2
  movapd %xmm2,%xmm5
  orpd   %xmm4,%xmm5
  jmp    L134
L26:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: data_tag: D88
type_infered: byte[8]
name: parse_data
arguments: data_tag: D60
type_infered: byte[8]
name: parse_data
arguments: data_tag: D46
type_infered: f64
name: parse_data
arguments: data_tag: D80
type_infered: f64
name: parse_data
arguments: data_tag: D51
type_infered: f64
name: parse_data
arguments: data_tag: D33
type_infered: f64
name: parse_data
arguments: data_tag: D42
type_infered: f64
name: parse_data
arguments: data_tag: D78
type_infered: f64

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

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

tool: D46: 0.25

tool: D80: 1

tool: D51: 3

tool: D33: 0.5

tool: D42: 0.125

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