user: 下面的汇编指令对应的C源代码是什么:

```asm
<lagrange_interp_2Dfo>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0xd8,%rsp
  movsd  %xmm0,-0x88(%rbp)
  movsd  %xmm1,-0x90(%rbp)
  mov    %edi,-0x94(%rbp)
  mov    %esi,-0x98(%rbp)
  mov    %rdx,-0xa0(%rbp)
  mov    %ecx,-0xa4(%rbp)
  mov    %r8,-0xb0(%rbp)
  mov    %r9d,-0xa8(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0xb8(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  mov    %rsp,%rax
  mov    %rax,%rbx
  mov    -0x94(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  sete   %al
  movzbl %al,%eax
  mov    %eax,%edi
  call   <assert@plt>
  mov    -0x98(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  sete   %al
  movzbl %al,%eax
  mov    %eax,%edi
  call   <assert@plt>
  mov    -0xa4(%rbp),%eax
  cmp    -0x94(%rbp),%eax
  setg   %al
  movzbl %al,%eax
  mov    %eax,%edi
  call   <assert@plt>
  mov    -0xa8(%rbp),%eax
  cmp    -0x98(%rbp),%eax
  setg   %al
  movzbl %al,%eax
  mov    %eax,%edi
  call   <assert@plt>
  mov    -0x94(%rbp),%eax
  lea    0x1(%rax),%edx
  movslq %edx,%rax
  sub    $0x1,%rax
  mov    %rax,-0x68(%rbp)
  movslq %edx,%rax
  mov    %rax,-0xd0(%rbp)
  movq   $0x0,-0xc8(%rbp)
  movslq %edx,%rax
  mov    %rax,-0xe0(%rbp)
  movq   $0x0,-0xd8(%rbp)
  movslq %edx,%rax
  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
L82:
  cmp    %rdx,%rsp
  je     L85
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L82
L85:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L40
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L40:
  mov    %rsp,%rax
  add    $0x7,%rax
  shr    $0x3,%rax
  shl    $0x3,%rax
  mov    %rax,-0x60(%rbp)
  mov    -0x98(%rbp),%eax
  add    $0x1,%eax
  movslq %eax,%rdx
  sub    $0x1,%rdx
  mov    %rdx,-0x58(%rbp)
  movslq %eax,%rdx
  mov    %rdx,-0xf0(%rbp)
  movq   $0x0,-0xe8(%rbp)
  movslq %eax,%rdx
  mov    %rdx,-0x100(%rbp)
  movq   $0x0,-0xf8(%rbp)
  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
L36:
  cmp    %rdx,%rsp
  je     L46
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L36
L46:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L118
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L118:
  mov    %rsp,%rax
  add    $0x7,%rax
  shr    $0x3,%rax
  shl    $0x3,%rax
  mov    %rax,-0x50(%rbp)
  mov    -0x94(%rbp),%eax
  lea    0x1(%rax),%edx
  mov    -0x98(%rbp),%eax
  add    $0x1,%eax
  imul   %edx,%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,%edi
  mov    $0x0,%edx
  div    %rdi
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L61:
  cmp    %rdx,%rsp
  je     L5
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L61
L5:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L49
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L49:
  mov    %rsp,%rax
  add    $0x7,%rax
  shr    $0x3,%rax
  shl    $0x3,%rax
  mov    %rax,-0x40(%rbp)
  mov    -0x94(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  mov    %eax,%edx
  mov    -0xa4(%rbp),%esi
  mov    -0xa0(%rbp),%rcx
  mov    -0x88(%rbp),%rax
  mov    %rcx,%rdi
  movq   %rax,%xmm0
  call   <find_xb@plt>
  mov    %eax,-0x78(%rbp)
  mov    -0x98(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  mov    %eax,%edx
  mov    -0xa8(%rbp),%esi
  mov    -0xb0(%rbp),%rcx
  mov    -0x90(%rbp),%rax
  mov    %rcx,%rdi
  movq   %rax,%xmm0
  call   <find_xb@plt>
  mov    %eax,-0x74(%rbp)
  movl   $0x0,-0x80(%rbp)
  jmp    L32
L8:
  mov    -0x78(%rbp),%edx
  mov    -0x80(%rbp),%eax
  add    %eax,%edx
  mov    -0x94(%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    -0xa0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x60(%rbp),%rax
  mov    -0x80(%rbp),%edx
  movslq %edx,%rdx
  movsd  %xmm0,(%rax,%rdx,8)
  movl   $0x0,-0x7c(%rbp)
  jmp    L127
L34:
  mov    -0x74(%rbp),%edx
  mov    -0x7c(%rbp),%eax
  add    %eax,%edx
  mov    -0x98(%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    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x50(%rbp),%rax
  mov    -0x7c(%rbp),%edx
  movslq %edx,%rdx
  movsd  %xmm0,(%rax,%rdx,8)
  mov    -0x98(%rbp),%eax
  lea    0x1(%rax),%ecx
  mov    -0x94(%rbp),%eax
  lea    0x1(%rax),%esi
  mov    -0x7c(%rbp),%edx
  mov    -0x80(%rbp),%eax
  mov    %eax,%edi
  call   <index_2D_to_1D@plt>
  mov    %eax,-0x70(%rbp)
  mov    -0x74(%rbp),%edx
  mov    -0x7c(%rbp),%eax
  add    %eax,%edx
  mov    -0x98(%rbp),%eax
  mov    %eax,%ecx
  shr    $0x1f,%ecx
  add    %ecx,%eax
  sar    %eax
  neg    %eax
  lea    (%rdx,%rax,1),%esi
  mov    -0x78(%rbp),%edx
  mov    -0x80(%rbp),%eax
  add    %eax,%edx
  mov    -0x94(%rbp),%eax
  mov    %eax,%ecx
  shr    $0x1f,%ecx
  add    %ecx,%eax
  sar    %eax
  neg    %eax
  lea    (%rdx,%rax,1),%edi
  mov    -0xa8(%rbp),%edx
  mov    -0xa4(%rbp),%eax
  mov    %edx,%ecx
  mov    %esi,%edx
  mov    %eax,%esi
  call   <index_2D_to_1D@plt>
  mov    %eax,-0x6c(%rbp)
  mov    -0x6c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x40(%rbp),%rax
  mov    -0x70(%rbp),%edx
  movslq %edx,%rdx
  movsd  %xmm0,(%rax,%rdx,8)
  addl   $0x1,-0x7c(%rbp)
L127:
  mov    -0x98(%rbp),%eax
  cmp    -0x7c(%rbp),%eax
  jge    L34
  addl   $0x1,-0x80(%rbp)
L32:
  mov    -0x94(%rbp),%eax
  cmp    -0x80(%rbp),%eax
  jge    L8
  mov    -0x98(%rbp),%eax
  lea    0x1(%rax),%ecx
  mov    -0x94(%rbp),%eax
  lea    0x1(%rax),%esi
  mov    -0x40(%rbp),%r8
  mov    -0x50(%rbp),%rdx
  mov    -0x60(%rbp),%rdi
  movsd  -0x90(%rbp),%xmm0
  mov    -0x88(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <lagrange_interp_2D@plt>
  movq   %xmm0,%rax
  mov    %rbx,%rsp
  mov    -0x38(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L48
  call   <__stack_chk_fail@plt>
L48:
  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_2Dfo(double x, double y, int order_x, int order_y, const double x_values[], int num_x, const double y_values[], int num_y, const double z_values[]) {
  assert(order_x % 2 == 0);
  assert(order_y % 2 == 0);
  assert(num_x > order_x);
  assert(num_y > order_y);
  double mx[order_x + 1];
  double my[order_y + 1];
  double mz[(order_x + 1) * (order_y + 1)];
  int ix = find_xb(x, x_values, num_x, order_x / 2);
  int iy = find_xb(y, y_values, num_y, order_y / 2);
  int iz_small, iz_big;
  for (int i = 0; i < order_x + 1; i++) {
    mx[i] = x_values[ix + i - order_x / 2];
    for (int j = 0; j < order_y + 1; j++) {
      my[j] = y_values[iy + j - order_y / 2];
      iz_small = index_2D_to_1D(i, order_x + 1, j, order_y + 1);
      iz_big = index_2D_to_1D(ix + i - order_x / 2, num_x, iy + j - order_y / 2, num_y);
      mz[iz_small] = z_values[iz_big];
    }
  }
  return lagrange_interp_2D(x, y, mx, order_x + 1, my, order_y + 1, mz);
}
```