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

```asm
<polint>:
  endbr64
  push   %r15
  movapd %xmm0,%xmm1
  xor    %eax,%eax
  push   %r14
  push   %r13
  push   %r12
  mov    %rdi,%r12
  push   %rbp
  mov    %edx,%ebp
  push   %rbx
  mov    %rsi,%rbx
  mov    %edx,%esi
  sub    $0x58,%rsp
  subsd  0x8(%rdi),%xmm1
  movq   D25(%rip),%xmm2
  mov    $0x1,%edi
  mov    %edx,0x4c(%rsp)
  mov    %rcx,0x38(%rsp)
  andpd  %xmm2,%xmm1
  mov    %r8,0x40(%rsp)
  movsd  %xmm0,0x10(%rsp)
  movsd  %xmm1,0x8(%rsp)
  call   <vector@plt>
  mov    %ebp,%esi
  mov    $0x1,%edi
  mov    %rax,%r14
  xor    %eax,%eax
  call   <vector@plt>
  mov    %rax,%r15
  test   %ebp,%ebp
  jle    L119
  movsd  0x10(%rsp),%xmm4
  lea    0x1(%rbp),%ecx
  mov    $0x1,%eax
  mov    $0x1,%edx
  movq   D25(%rip),%xmm2
  movsd  0x8(%rsp),%xmm1
  nopl   0x0(%rax)
L80:
  movapd %xmm4,%xmm0
  subsd  (%r12,%rax,8),%xmm0
  andpd  %xmm2,%xmm0
  comisd %xmm0,%xmm1
  minsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  movsd  (%rbx,%rax,8),%xmm0
  cmova  %eax,%edx
  movsd  %xmm0,(%r14,%rax,8)
  movsd  %xmm0,(%r15,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rcx
  jne    L80
  lea    -0x1(%rdx),%eax
  movslq %edx,%rdx
  mov    0x38(%rsp),%rcx
  movsd  (%rbx,%rdx,8),%xmm0
  mov    %eax,0x48(%rsp)
  movsd  %xmm0,(%rcx)
  mov    0x4c(%rsp),%ecx
  cmp    $0x1,%ecx
  jle    L37
  lea    -0x1(%rcx),%ebp
  lea    0x8(%r12),%r13
  pxor   %xmm5,%xmm5
  cltq
  mov    %rax,0x30(%rsp)
  nopw   0x0(%rax,%rax,1)
L68:
  mov    $0x1,%ebx
  nopl   (%rax)
L53:
  movsd  (%r12,%rbx,8),%xmm0
  movsd  0x0(%r13,%rbx,8),%xmm1
  movsd  0x8(%r14,%rbx,8),%xmm2
  subsd  (%r15,%rbx,8),%xmm2
  subsd  %xmm4,%xmm0
  subsd  %xmm4,%xmm1
  movapd %xmm0,%xmm3
  subsd  %xmm1,%xmm3
  ucomisd %xmm5,%xmm3
  jp     L20
  je     L115
L20:
  divsd  %xmm3,%xmm2
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  movsd  %xmm1,(%r15,%rbx,8)
  movsd  %xmm0,(%r14,%rbx,8)
  add    $0x1,%rbx
  cmp    %ebx,%ebp
  jge    L53
  mov    0x48(%rsp),%eax
  add    %eax,%eax
  cmp    %eax,%ebp
  jle    L124
L33:
  mov    0x30(%rsp),%rax
  add    $0x8,%r13
  movsd  0x8(%r14,%rax,8),%xmm0
  mov    0x40(%rsp),%rax
  movsd  %xmm0,(%rax)
  mov    0x38(%rsp),%rax
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  sub    $0x1,%ebp
  jne    L68
L37:
  mov    0x4c(%rsp),%ebx
  mov    %r15,%rdi
  mov    $0x1,%esi
  xor    %eax,%eax
  mov    %ebx,%edx
  call   <free_vector@plt>
  add    $0x58,%rsp
  mov    %ebx,%edx
  mov    %r14,%rdi
  pop    %rbx
  mov    $0x1,%esi
  pop    %rbp
  xor    %eax,%eax
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1070 <free_vector@plt>
  nop
L115:
  xor    %eax,%eax
  lea    D10(%rip),%rdi
  movsd  %xmm4,0x28(%rsp)
  movsd  %xmm3,0x20(%rsp)
  movsd  %xmm2,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  movsd  %xmm0,0x8(%rsp)
  call   <nrerror@plt>
  movsd  0x18(%rsp),%xmm2
  pxor   %xmm5,%xmm5
  movsd  0x20(%rsp),%xmm3
  movsd  0x10(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm0
  divsd  %xmm3,%xmm2
  movsd  0x28(%rsp),%xmm4
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  movsd  %xmm1,(%r15,%rbx,8)
  movsd  %xmm0,(%r14,%rbx,8)
  add    $0x1,%rbx
  cmp    %ebx,%ebp
  jge    L53
  mov    0x48(%rsp),%eax
  add    %eax,%eax
  cmp    %eax,%ebp
  jg     L33
L124:
  mov    0x30(%rsp),%rcx
  mov    0x38(%rsp),%rsi
  add    $0x8,%r13
  mov    0x40(%rsp),%rdx
  subl   $0x1,0x48(%rsp)
  movsd  (%r15,%rcx,8),%xmm0
  movslq 0x48(%rsp),%rax
  movsd  %xmm0,(%rdx)
  addsd  (%rsi),%xmm0
  movsd  %xmm0,(%rsi)
  sub    $0x1,%ebp
  je     L37
  mov    %rax,0x30(%rsp)
  jmp    L68
L119:
  movsd  0x8(%rbx),%xmm0
  mov    0x38(%rsp),%rax
  movsd  %xmm0,(%rax)
  jmp    L37
```

tool_calls: name: disassemble_data
arguments: data_tag: D25
type: qword
name: disassemble_data
arguments: data_tag: D10
type: string

tool: D25: 9223372036854775807

tool: D10: Error in routine POLINT

assistant: ```c
void polint(xa, ya, n, x, y, dy) double xa[], ya[], x, *y, *dy;
int n;
{
  int i, m, ns = 1;
  double den, dif, dift, ho, hp, w;
  double *c, *d, *vector();
  void nrerror(), free_vector();
  dif = fabs(x - xa[1]);
  c = vector(1, n);
  d = vector(1, n);
  for (i = 1; i <= n; i++) {
    if ((dift = fabs(x - xa[i])) < dif) {
      ns = i;
      dif = dift;
    }
    c[i] = ya[i];
    d[i] = ya[i];
  }
  *y = ya[ns--];
  for (m = 1; m < n; m++) {
    for (i = 1; i <= n - m; i++) {
      ho = xa[i] - x;
      hp = xa[i + m] - x;
      w = c[i + 1] - d[i];
      if ((den = ho - hp) == 0.0)
        nrerror("Error in routine POLINT");
      den = w / den;
      d[i] = hp * den;
      c[i] = ho * den;
    }
    *y += (*dy = (2 * ns < (n - m) ? c[ns + 1] : d[ns--]));
  }
  free_vector(d, 1, n);
  free_vector(c, 1, n);
}
```