user: Demonstrate the C code that generates the assembly output below:

```asm
<lagrange_interp_2Dfo>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x88,%rsp
  movsd  %xmm0,-0x98(%rbp)
  movsd  %xmm1,-0xa0(%rbp)
  mov    %edi,%ebx
  mov    %edi,-0x74(%rbp)
  mov    %esi,%r15d
  mov    %esi,-0x44(%rbp)
  mov    %rdx,%r12
  mov    %ecx,-0x54(%rbp)
  mov    %r8,%r13
  mov    %r9d,-0x58(%rbp)
  mov    0x10(%rbp),%r14
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  xor    $0x1,%edi
  and    $0x1,%edi
  call   <assert@plt>
  mov    %r15d,%edi
  xor    $0x1,%edi
  and    $0x1,%edi
  call   <assert@plt>
  mov    -0x54(%rbp),%eax
  cmp    %eax,%ebx
  setl   %dil
  movzbl %dil,%edi
  call   <assert@plt>
  mov    -0x58(%rbp),%eax
  cmp    %eax,%r15d
  setl   %dil
  movzbl %dil,%edi
  call   <assert@plt>
  lea    0x1(%rbx),%r15d
  movslq %r15d,%rax
  lea    0xf(,%rax,8),%rax
  mov    %rax,%rcx
  and    $0xfffffffffffffff0,%rcx
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rdx
  sub    %rax,%rdx
L96:
  cmp    %rdx,%rsp
  je     L87
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L96
L87:
  mov    %rcx,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L5
  orq    $0x0,-0x8(%rsp,%rax,1)
L5:
  mov    %rsp,-0x80(%rbp)
  mov    -0x44(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x60(%rbp)
  cltq
  lea    0xf(,%rax,8),%rax
  mov    %rax,%rdx
  and    $0xfffffffffffffff0,%rdx
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rcx
  sub    %rax,%rcx
L72:
  cmp    %rcx,%rsp
  je     L0
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L72
L0:
  mov    %rdx,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L22
  orq    $0x0,-0x8(%rsp,%rax,1)
L22:
  mov    %rsp,-0x68(%rbp)
  mov    -0x60(%rbp),%eax
  imul   %r15d,%eax
  cltq
  lea    0xf(,%rax,8),%rdx
  mov    %rdx,%rax
  and    $0xfffffffffffffff0,%rax
  and    $0xfffffffffffff000,%rdx
  mov    %rsp,%rcx
  sub    %rdx,%rcx
L61:
  cmp    %rcx,%rsp
  je     L77
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L61
L77:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  test   %rdx,%rdx
  je     L68
  orq    $0x0,-0x8(%rsp,%rdx,1)
L68:
  mov    %rsp,-0xa8(%rbp)
  mov    -0x74(%rbp),%edi
  mov    %edi,%eax
  shr    $0x1f,%eax
  add    %edi,%eax
  sar    %eax
  mov    %eax,-0x48(%rbp)
  mov    %eax,%edx
  mov    -0x54(%rbp),%esi
  mov    %r12,%rdi
  movsd  -0x98(%rbp),%xmm0
  call   <find_xb@plt>
  mov    %eax,%ebx
  mov    -0x44(%rbp),%ecx
  mov    %ecx,%eax
  shr    $0x1f,%eax
  add    %ecx,%eax
  sar    %eax
  mov    %eax,-0x8c(%rbp)
  mov    %eax,%edx
  mov    -0x58(%rbp),%esi
  mov    %r13,%rdi
  movsd  -0xa0(%rbp),%xmm0
  call   <find_xb@plt>
  mov    %eax,%ecx
  mov    %eax,-0x90(%rbp)
  cmpl   $0x0,-0x74(%rbp)
  js     L102
  mov    -0x48(%rbp),%edi
  movslq %edi,%rdx
  movslq %ebx,%rax
  sub    %rdx,%rax
  lea    (%r12,%rax,8),%rsi
  mov    %rsi,-0x88(%rbp)
  movslq -0x8c(%rbp),%rdx
  movslq %ecx,%rax
  sub    %rdx,%rax
  lea    0x0(%r13,%rax,8),%r13
  movq   $0x0,-0x70(%rbp)
  sub    %edi,%ebx
  mov    %ebx,-0x78(%rbp)
  mov    %r15d,-0x5c(%rbp)
  mov    -0xa8(%rbp),%r15
L35:
  mov    -0x70(%rbp),%rax
  mov    %eax,-0x48(%rbp)
  mov    -0x78(%rbp),%esi
  add    %eax,%esi
  mov    %esi,-0x4c(%rbp)
  mov    -0x88(%rbp),%rdi
  movsd  (%rdi,%rax,8),%xmm0
  mov    -0x80(%rbp),%rcx
  movsd  %xmm0,(%rcx,%rax,8)
  cmpl   $0x0,-0x44(%rbp)
  js     L138
  mov    $0x0,%ebx
  mov    -0x90(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  mov    %eax,-0x50(%rbp)
L95:
  movsd  0x0(%r13,%rbx,8),%xmm0
  mov    -0x68(%rbp),%rax
  movsd  %xmm0,(%rax,%rbx,8)
  mov    -0x60(%rbp),%ecx
  mov    %ebx,%edx
  mov    -0x5c(%rbp),%esi
  mov    -0x48(%rbp),%edi
  call   <index_2D_to_1D@plt>
  mov    %eax,%r12d
  mov    -0x50(%rbp),%eax
  lea    (%rax,%rbx,1),%edx
  mov    -0x58(%rbp),%ecx
  mov    -0x54(%rbp),%esi
  mov    -0x4c(%rbp),%edi
  call   <index_2D_to_1D@plt>
  movslq %r12d,%r12
  cltq
  movsd  (%r14,%rax,8),%xmm0
  movsd  %xmm0,(%r15,%r12,8)
  add    $0x1,%rbx
  cmp    %ebx,-0x44(%rbp)
  jge    L95
L138:
  addq   $0x1,-0x70(%rbp)
  mov    -0x70(%rbp),%rax
  cmp    %eax,-0x74(%rbp)
  jge    L35
  mov    -0x5c(%rbp),%r15d
L102:
  mov    -0xa8(%rbp),%r8
  mov    -0x60(%rbp),%ecx
  mov    -0x68(%rbp),%rdx
  mov    %r15d,%esi
  mov    -0x80(%rbp),%rdi
  movsd  -0xa0(%rbp),%xmm1
  movsd  -0x98(%rbp),%xmm0
  call   <lagrange_interp_2D@plt>
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L99
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
L99:
  call   <__stack_chk_fail@plt>
```

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);
}
```