user: 请从给定的汇编指令推导出C源代码：

```asm
<splint>:
  endbr64
  push   %r15
  mov    %rdx,%r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %fs:0x28,%rdx
  mov    %rdx,0x18(%rsp)
  xor    %edx,%edx
  cmpq   $0x0,0x60(%rsp)
  mov    0x68(%rsp),%r10d
  je     L81
  mov    %rdi,%rbp
  lea    -0x1(%rcx),%rdi
  mov    %rsi,%r12
  mov    %r8,%r13
  mov    %rdi,0x8(%rsp)
  mov    %r9,%r14
  xor    %ebx,%ebx
  xor    %esi,%esi
  jmp    L19
  nopl   0x0(%rax)
L137:
  cmp    $0x1,%r10d
  je     L131
  mov    0x10(%rsp),%rdi
  movsd  0x0(%rbp,%rax,8),%xmm3
  lea    0x0(,%rax,8),%rdx
  movsd  0x0(%rbp,%rdi,8),%xmm1
  lea    0x0(,%rdi,8),%rcx
L53:
  movapd %xmm3,%xmm4
  subsd  %xmm0,%xmm3
  subsd  %xmm1,%xmm0
  subsd  %xmm1,%xmm4
  divsd  %xmm4,%xmm3
  divsd  %xmm4,%xmm0
  movapd %xmm3,%xmm2
  mulsd  %xmm3,%xmm2
  mulsd  %xmm3,%xmm2
  subsd  %xmm3,%xmm2
  mulsd  (%r15,%rcx,1),%xmm2
  mulsd  (%r12,%rcx,1),%xmm3
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  subsd  %xmm0,%xmm1
  mulsd  (%r15,%rdx,1),%xmm1
  mulsd  (%r12,%rdx,1),%xmm0
  addsd  %xmm2,%xmm1
  addsd  %xmm0,%xmm3
  mulsd  %xmm4,%xmm1
  mulsd  %xmm4,%xmm1
  divsd  D8(%rip),%xmm1
  addsd  %xmm3,%xmm1
  movsd  %xmm1,(%r14,%rbx,8)
  add    $0x1,%rbx
  cmp    %rbx,0x60(%rsp)
  je     L81
L19:
  movsd  0x0(%r13,%rbx,8),%xmm0
  test   %r10d,%r10d
  jne    L137
  mov    0x8(%rsp),%rdx
  lea    0x10(%rsp),%rcx
  mov    %rbp,%rdi
  mov    %r10d,0x68(%rsp)
  movsd  %xmm0,(%rsp)
  call   <lin_search@plt>
  movsd  (%rsp),%xmm0
  mov    0x68(%rsp),%r10d
  test   %rax,%rax
  je     L99
  mov    0x10(%rsp),%rsi
  movsd  0x0(%r13,%rbx,8),%xmm0
  lea    0x1(%rsi),%rax
  lea    0x0(,%rax,8),%rdx
  movsd  0x0(%rbp,%rax,8),%xmm3
  movsd  -0x8(%rbp,%rdx,1),%xmm1
  lea    -0x8(%rdx),%rcx
  jmp    L53
L131:
  lea    0x1(%rsi),%rax
  movsd  0x0(%rbp,%rax,8),%xmm3
  lea    0x0(,%rax,8),%rdx
  comisd %xmm3,%xmm0
  jae    L35
  movsd  -0x8(%rbp,%rdx,1),%xmm1
  lea    -0x8(%rdx),%rcx
  comisd %xmm0,%xmm1
  ja     L73
L4:
  mov    %rsi,0x10(%rsp)
  jmp    L53
  nopw   0x0(%rax,%rax,1)
L81:
  mov    0x18(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L72
  add    $0x28,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L99:
  mov    0x8(%rsp),%rdx
  xor    %esi,%esi
  mov    %rbp,%rdi
  mov    %r10d,0x68(%rsp)
  call   <b_search@plt>
  movsd  0x0(%r13,%rbx,8),%xmm0
  mov    0x68(%rsp),%r10d
  mov    %rax,0x10(%rsp)
  mov    %rax,%rsi
  lea    0x1(%rax),%rax
  lea    0x0(,%rax,8),%rdx
  movsd  0x0(%rbp,%rax,8),%xmm3
  movsd  -0x8(%rbp,%rdx,1),%xmm1
  lea    -0x8(%rdx),%rcx
  jmp    L53
  nopl   0x0(%rax)
L35:
  mov    %r10d,0x68(%rsp)
  mov    0x8(%rsp),%rdx
L75:
  mov    %rbp,%rdi
  call   <b_search@plt>
  movsd  0x0(%r13,%rbx,8),%xmm0
  mov    0x68(%rsp),%r10d
  mov    %rax,%rsi
  lea    0x1(%rax),%rax
  lea    0x0(,%rax,8),%rdx
  movsd  0x0(%rbp,%rax,8),%xmm3
  movsd  -0x8(%rbp,%rdx,1),%xmm1
  lea    -0x8(%rdx),%rcx
  jmp    L4
  nopw   0x0(%rax,%rax,1)
L73:
  mov    %rsi,%rdx
  mov    %r10d,0x68(%rsp)
  xor    %esi,%esi
  jmp    L75
L72:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: label: D8
datatype: f64

tool: D8: 6

assistant: ```c
int splint(const double *pxa, const double *pya, const double *py2a, size_t n, const double *px, double *py, size_t nx, int blookup) {
  size_t klo;
  size_t khi;
  size_t prev_idx = 0;
  double h;
  double b;
  double a;
  double x;
  for (size_t i = 0; i < nx; i++) {
    x = px[i];
    if (blookup == 0) {
      if (lin_search(pxa, x, prev_idx, n - 1, &klo)) {
        khi = klo + 1;
        prev_idx = klo;
      } else {
        klo = b_search(pxa, x, 0, n - 1);
        khi = klo + 1;
        prev_idx = klo;
      }
    } else if (blookup == 1) {
      if (x >= pxa[prev_idx + 1])
        prev_idx = b_search(pxa, x, prev_idx, n - 1);
      else if (x < pxa[prev_idx])
        prev_idx = b_search(pxa, x, 0, prev_idx);
      klo = prev_idx;
      khi = klo + 1;
    }
    h = pxa[khi] - pxa[klo];
    a = (pxa[khi] - px[i]) / h;
    b = (px[i] - pxa[klo]) / h;
    py[i] = a * pya[klo] + b * pya[khi] + ((a * a * a - a) * py2a[klo] + (b * b * b - b) * py2a[khi]) * h * h / 6.0;
  }
  return 0;
}
```