user: Please show me the C implementation that matches the assembly code listed below:

```asm
<polint>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rdi,%r14
  mov    %rsi,%rbx
  mov    %edx,%r15d
  mov    %edx,0x2c(%rsp)
  mov    %rcx,0x18(%rsp)
  mov    %r8,0x20(%rsp)
  pxor   %xmm7,%xmm7
  cvtsd2ss %xmm0,%xmm7
  movss  %xmm7,0x10(%rsp)
  movaps %xmm7,%xmm0
  subss  0x4(%rdi),%xmm0
  andps  D9(%rip),%xmm0
  movss  %xmm0,(%rsp)
  mov    %edx,%esi
  mov    $0x1,%edi
  mov    $0x0,%eax
  call   <vector@plt>
  mov    %rax,%rbp
  mov    %r15d,%esi
  mov    $0x1,%edi
  mov    $0x0,%eax
  call   <vector@plt>
  mov    %rax,%r12
  test   %r15d,%r15d
  jle    L78
  lea    0x1(%r15),%ecx
  mov    $0x1,%eax
  mov    $0x1,%edx
  movss  D9(%rip),%xmm1
L137:
  movss  0x10(%rsp),%xmm0
  subss  (%r14,%rax,4),%xmm0
  andps  %xmm1,%xmm0
  movss  (%rsp),%xmm7
  comiss %xmm0,%xmm7
  minss  %xmm7,%xmm0
  movss  %xmm0,(%rsp)
  cmova  %eax,%edx
  movss  (%rbx,%rax,4),%xmm0
  movss  %xmm0,0x0(%rbp,%rax,4)
  movss  %xmm0,(%r12,%rax,4)
  add    $0x1,%rax
  cmp    %rcx,%rax
  jne    L137
L89:
  lea    -0x1(%rdx),%eax
  mov    %eax,0x14(%rsp)
  movslq %edx,%rdx
  movss  (%rbx,%rdx,4),%xmm0
  mov    0x18(%rsp),%rax
  movss  %xmm0,(%rax)
  mov    0x2c(%rsp),%eax
  cmp    $0x1,%eax
  jle    L91
  sub    $0x1,%eax
  mov    %eax,0x28(%rsp)
  lea    0x4(%r14),%r13
  jmp    L38
L78:
  mov    $0x1,%edx
  jmp    L89
L102:
  lea    D19(%rip),%rdi
  mov    $0x0,%eax
  call   <nrerror@plt>
L97:
  movss  0x8(%rsp),%xmm0
  divss  0xc(%rsp),%xmm0
  movss  0x4(%rsp),%xmm1
  mulss  %xmm0,%xmm1
  movss  %xmm1,(%r12,%rbx,4)
  mulss  (%rsp),%xmm0
  movss  %xmm0,0x0(%rbp,%rbx,4)
  add    $0x1,%rbx
  cmp    %ebx,%r15d
  jl     L71
L18:
  movss  (%r14,%rbx,4),%xmm0
  movaps %xmm0,%xmm2
  movss  0x10(%rsp),%xmm4
  subss  %xmm4,%xmm2
  movss  %xmm2,(%rsp)
  movss  0x0(%r13,%rbx,4),%xmm0
  movaps %xmm0,%xmm3
  subss  %xmm4,%xmm3
  movss  %xmm3,0x4(%rsp)
  movss  0x4(%rbp,%rbx,4),%xmm0
  subss  (%r12,%rbx,4),%xmm0
  movss  %xmm0,0x8(%rsp)
  subss  %xmm3,%xmm2
  movss  %xmm2,0xc(%rsp)
  ucomiss D26(%rip),%xmm2
  jp     L97
  jne    L97
  jmp    L102
L71:
  mov    0x14(%rsp),%esi
  lea    (%rsi,%rsi,1),%eax
  cmp    %r15d,%eax
  jge    L48
  movslq %esi,%rax
  movss  0x4(%rbp,%rax,4),%xmm0
L35:
  mov    0x20(%rsp),%rax
  movss  %xmm0,(%rax)
  mov    0x18(%rsp),%rax
  addss  (%rax),%xmm0
  movss  %xmm0,(%rax)
  add    $0x4,%r13
  subl   $0x1,0x28(%rsp)
  je     L91
L38:
  mov    0x28(%rsp),%eax
  mov    %eax,%r15d
  mov    $0x1,%ebx
  test   %eax,%eax
  jg     L18
  jmp    L71
L48:
  mov    0x14(%rsp),%ecx
  movslq %ecx,%rax
  movss  (%r12,%rax,4),%xmm0
  lea    -0x1(%rcx),%eax
  mov    %eax,0x14(%rsp)
  jmp    L35
L91:
  mov    0x2c(%rsp),%ebx
  mov    %ebx,%edx
  mov    $0x1,%esi
  mov    %r12,%rdi
  mov    $0x0,%eax
  call   <free_vector@plt>
  mov    %ebx,%edx
  mov    $0x1,%esi
  mov    %rbp,%rdi
  mov    $0x0,%eax
  call   <free_vector@plt>
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: parse_data
arguments: data_label: D9
type_infered: dword[4]
name: parse_data
arguments: data_label: D19
type_infered: string
name: parse_data
arguments: data_label: D26
type_infered: f64

tool: {"D9": [0x7fffffff,0x00000000,0x00000000,0x00000000]}

tool: D19: Error in routine POLINT

tool: D26: 0

assistant: ```c
void polint(xa, ya, n, x, y, dy) float xa[], ya[], x, *y, *dy;
int n;
{
  int i, m, ns = 1;
  float den, dif, dift, ho, hp, w;
  float *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);
}
```