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

```asm
<lagrange_interp_1Dfo>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x78,%rsp
  movsd  %xmm0,-0x68(%rbp)
  mov    %edi,-0x6c(%rbp)
  mov    %rsi,-0x78(%rbp)
  mov    %edx,-0x70(%rbp)
  mov    %rcx,-0x80(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  mov    %rsp,%rax
  mov    %rax,%rbx
  mov    -0x6c(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  sete   %al
  movzbl %al,%eax
  mov    %eax,%edi
  call   <assert@plt>
  mov    -0x70(%rbp),%eax
  cmp    -0x6c(%rbp),%eax
  setg   %al
  movzbl %al,%eax
  mov    %eax,%edi
  call   <assert@plt>
  mov    -0x6c(%rbp),%eax
  add    $0x1,%eax
  movslq %eax,%rdx
  sub    $0x1,%rdx
  mov    %rdx,-0x58(%rbp)
  movslq %eax,%rdx
  mov    %rdx,-0x90(%rbp)
  movq   $0x0,-0x88(%rbp)
  movslq %eax,%rdx
  mov    %rdx,-0xa0(%rbp)
  movq   $0x0,-0x98(%rbp)
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%edi
  mov    $0x0,%edx
  div    %rdi
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L17:
  cmp    %rdx,%rsp
  je     L80
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L17
L80:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L74
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L74:
  mov    %rsp,%rax
  add    $0x7,%rax
  shr    $0x3,%rax
  shl    $0x3,%rax
  mov    %rax,-0x50(%rbp)
  mov    -0x6c(%rbp),%eax
  add    $0x1,%eax
  movslq %eax,%rdx
  sub    $0x1,%rdx
  mov    %rdx,-0x48(%rbp)
  movslq %eax,%rdx
  mov    %rdx,%r14
  mov    $0x0,%r15d
  movslq %eax,%rdx
  mov    %rdx,%r12
  mov    $0x0,%r13d
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%esi
  mov    $0x0,%edx
  div    %rsi
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L14:
  cmp    %rdx,%rsp
  je     L85
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L14
L85:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L57
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L57:
  mov    %rsp,%rax
  add    $0x7,%rax
  shr    $0x3,%rax
  shl    $0x3,%rax
  mov    %rax,-0x40(%rbp)
  mov    -0x6c(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  mov    %eax,%edx
  mov    -0x70(%rbp),%esi
  mov    -0x78(%rbp),%rcx
  mov    -0x68(%rbp),%rax
  mov    %rcx,%rdi
  movq   %rax,%xmm0
  call   <find_xb@plt>
  mov    %eax,-0x5c(%rbp)
  movl   $0x0,-0x60(%rbp)
  jmp    L1
L26:
  mov    -0x5c(%rbp),%edx
  mov    -0x60(%rbp),%eax
  add    %eax,%edx
  mov    -0x6c(%rbp),%eax
  mov    %eax,%ecx
  shr    $0x1f,%ecx
  add    %ecx,%eax
  sar    %eax
  neg    %eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x50(%rbp),%rax
  mov    -0x60(%rbp),%edx
  movslq %edx,%rdx
  movsd  %xmm0,(%rax,%rdx,8)
  mov    -0x5c(%rbp),%edx
  mov    -0x60(%rbp),%eax
  add    %eax,%edx
  mov    -0x6c(%rbp),%eax
  mov    %eax,%ecx
  shr    $0x1f,%ecx
  add    %ecx,%eax
  sar    %eax
  neg    %eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x40(%rbp),%rax
  mov    -0x60(%rbp),%edx
  movslq %edx,%rdx
  movsd  %xmm0,(%rax,%rdx,8)
  addl   $0x1,-0x60(%rbp)
L1:
  mov    -0x6c(%rbp),%eax
  cmp    -0x60(%rbp),%eax
  jge    L26
  mov    -0x6c(%rbp),%eax
  lea    0x1(%rax),%esi
  mov    -0x40(%rbp),%rdx
  mov    -0x50(%rbp),%rcx
  mov    -0x68(%rbp),%rax
  mov    %rcx,%rdi
  movq   %rax,%xmm0
  call   <lagrange_interp_1D@plt>
  movq   %xmm0,%rax
  mov    %rbx,%rsp
  mov    -0x38(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L79
  call   <__stack_chk_fail@plt>
L79:
  movq   %rax,%xmm0
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
```

assistant: ```c
double lagrange_interp_1Dfo(double x, int order, const double x_values[], int num_x, const double y_values[]) {
  assert(order % 2 == 0);
  assert(num_x > order);
  double mx[order + 1];
  double my[order + 1];
  int i;
  int ix = find_xb(x, x_values, num_x, order / 2);
  for (i = 0; i < order + 1; i++) {
    mx[i] = x_values[ix + i - order / 2];
    my[i] = y_values[ix + i - order / 2];
  }
  return lagrange_interp_1D(x, mx, order + 1, my);
}
```