user: Translate the below assembly code into C programming language:

```asm
<bessel>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x90,%rsp
  mov    %rdi,-0x78(%rbp)
  mov    %esi,-0x7c(%rbp)
  movsd  %xmm0,-0x88(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  pxor   %xmm0,%xmm0
  comisd -0x88(%rbp),%xmm0
  jb     L41
  mov    -0x78(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  jmp    L243
L41:
  mov    -0x7c(%rbp),%eax
  sub    $0x1,%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  -0x88(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jb     L257
  mov    -0x7c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  jmp    L243
L257:
  mov    -0x88(%rbp),%rax
  movq   %rax,%xmm0
  call   <floor@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x50(%rbp)
  movsd  -0x88(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,-0x5c(%rbp)
  mov    -0x5c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  mov    -0x5c(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x60(%rbp)
  mov    -0x60(%rbp),%eax
  cmp    -0x7c(%rbp),%eax
  jge    L74
  movsd  -0x88(%rbp),%xmm0
  subsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0x60(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x5c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  mulsd  -0x50(%rbp),%xmm0
  movsd  -0x58(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  cmpl   $0x0,-0x5c(%rbp)
  jle    L0
  mov    -0x5c(%rbp),%eax
  add    $0x2,%eax
  cmp    %eax,-0x7c(%rbp)
  jle    L0
  mov    -0x5c(%rbp),%eax
  sub    $0x2,%eax
  mov    %eax,-0x60(%rbp)
  movl   $0x0,-0x64(%rbp)
  jmp    L90
L222:
  cmpl   $0x0,-0x60(%rbp)
  js     L184
  mov    -0x60(%rbp),%eax
  add    $0x1,%eax
  cmp    %eax,-0x7c(%rbp)
  jg     L116
L184:
  mov    -0x64(%rbp),%eax
  cltq
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x40(%rbp,%rax,8)
  jmp    L142
L116:
  mov    -0x60(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x60(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x64(%rbp),%eax
  cltq
  movsd  %xmm0,-0x40(%rbp,%rax,8)
L142:
  addl   $0x1,-0x60(%rbp)
  addl   $0x1,-0x64(%rbp)
L90:
  cmpl   $0x4,-0x64(%rbp)
  jle    L222
  movl   $0x0,-0x64(%rbp)
  jmp    L151
L158:
  mov    -0x64(%rbp),%eax
  add    $0x1,%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm0
  mov    -0x64(%rbp),%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x64(%rbp),%eax
  cltq
  movsd  %xmm0,-0x40(%rbp,%rax,8)
  addl   $0x1,-0x64(%rbp)
L151:
  cmpl   $0x3,-0x64(%rbp)
  jle    L158
  movsd  -0x50(%rbp),%xmm1
  movsd  D65(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x50(%rbp),%xmm0
  movsd  D45(%rip),%xmm2
  subsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x38(%rbp),%xmm1
  movsd  -0x30(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x48(%rbp),%xmm0
  movsd  -0x58(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  mov    -0x5c(%rbp),%eax
  add    $0x2,%eax
  cmp    %eax,-0x7c(%rbp)
  jle    L48
  movl   $0x0,-0x64(%rbp)
  jmp    L214
L172:
  mov    -0x64(%rbp),%eax
  add    $0x1,%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm0
  mov    -0x64(%rbp),%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x64(%rbp),%eax
  cltq
  movsd  %xmm0,-0x40(%rbp,%rax,8)
  addl   $0x1,-0x64(%rbp)
L214:
  cmpl   $0x2,-0x64(%rbp)
  jle    L172
  movsd  -0x48(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  D54(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x50(%rbp),%xmm0
  movsd  D38(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x48(%rbp),%xmm1
  movsd  -0x38(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x58(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  cmpl   $0x1,-0x5c(%rbp)
  jle    L247
  mov    -0x5c(%rbp),%eax
  add    $0x2,%eax
  cmp    %eax,-0x7c(%rbp)
  jle    L247
  movl   $0x0,-0x64(%rbp)
  jmp    L231
L203:
  mov    -0x64(%rbp),%eax
  add    $0x1,%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm0
  mov    -0x64(%rbp),%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x64(%rbp),%eax
  cltq
  movsd  %xmm0,-0x40(%rbp,%rax,8)
  addl   $0x1,-0x64(%rbp)
L231:
  cmpl   $0x1,-0x64(%rbp)
  jle    L203
  movsd  -0x48(%rbp),%xmm1
  movsd  D28(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x50(%rbp),%xmm2
  movsd  D45(%rip),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x50(%rbp),%xmm0
  movsd  D19(%rip),%xmm2
  subsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x40(%rbp),%xmm1
  movsd  -0x38(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x48(%rbp),%xmm0
  movsd  -0x58(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  jmp    L243
L74:
  nop
  jmp    L243
L0:
  nop
  jmp    L243
L48:
  nop
  jmp    L243
L247:
  nop
L243:
  movsd  -0x58(%rbp),%xmm0
  movq   %xmm0,%rax
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L185
  call   <__stack_chk_fail@plt>
L185:
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: unpack_data
arguments: data_label: D65
data_type: f64
name: unpack_data
arguments: data_label: D45
data_type: f64
name: unpack_data
arguments: data_label: D54
data_type: f64
name: unpack_data
arguments: data_label: D38
data_type: f64
name: unpack_data
arguments: data_label: D28
data_type: f64
name: unpack_data
arguments: data_label: D19
data_type: f64

tool: D65: 0.25

tool: D45: 1

tool: D54: 3

tool: D38: 0.5

tool: D28: 0.125

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