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

```asm
<scalars_lub_poly_full>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x320,%rsp
  mov    %edi,-0x2e4(%rbp)
  movsd  %xmm0,-0x2f0(%rbp)
  movsd  %xmm1,-0x2f8(%rbp)
  movsd  %xmm2,-0x300(%rbp)
  mov    %rsi,-0x308(%rbp)
  mov    %rdx,-0x310(%rbp)
  mov    %ecx,-0x2e8(%rbp)
  mov    %r8d,-0x314(%rbp)
  mov    %r9,-0x320(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movl   $0x0,-0x2d4(%rbp)
  jmp    L32
L49:
  mov    -0x2d4(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x320(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x2d4(%rbp)
L32:
  cmpl   $0x2b,-0x2d4(%rbp)
  jle    L49
  lea    -0x2d0(%rbp),%r8
  mov    -0x314(%rbp),%ecx
  mov    -0x2e8(%rbp),%edx
  mov    -0x308(%rbp),%rsi
  movsd  -0x300(%rbp),%xmm1
  movsd  -0x2f8(%rbp),%xmm0
  mov    -0x2f0(%rbp),%rdi
  mov    -0x2e4(%rbp),%eax
  mov    %r8,%r9
  mov    $0x1,%r8d
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rdi,%xmm0
  mov    %eax,%edi
  mov    $0x3,%eax
  call   <twobody_scalars_res@plt>
  lea    -0x220(%rbp),%r8
  mov    -0x314(%rbp),%ecx
  mov    -0x2e8(%rbp),%edx
  mov    -0x310(%rbp),%rsi
  movsd  -0x2f8(%rbp),%xmm1
  movsd  -0x300(%rbp),%xmm0
  mov    -0x2f0(%rbp),%rdi
  mov    -0x2e4(%rbp),%eax
  mov    %r8,%r9
  mov    $0x1,%r8d
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rdi,%xmm0
  mov    %eax,%edi
  mov    $0x3,%eax
  call   <twobody_scalars_res@plt>
  cmpl   $0x0,-0x2e4(%rbp)
  jne    L28
  lea    -0x170(%rbp),%rdx
  movsd  -0x300(%rbp),%xmm1
  movsd  -0x2f8(%rbp),%xmm0
  mov    -0x2f0(%rbp),%rax
  mov    %rdx,%rdi
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  mov    $0x3,%eax
  call   <scalars_minv_f_poly@plt>
  jmp    L78
L28:
  cmpl   $0x1,-0x2e4(%rbp)
  jne    L88
  lea    -0x170(%rbp),%rdx
  movsd  -0x300(%rbp),%xmm1
  movsd  -0x2f8(%rbp),%xmm0
  mov    -0x2f0(%rbp),%rax
  mov    %rdx,%rdi
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  mov    $0x3,%eax
  call   <scalars_minv_ft_poly@plt>
  jmp    L78
L88:
  lea    -0x170(%rbp),%rdx
  movsd  -0x300(%rbp),%xmm1
  movsd  -0x2f8(%rbp),%xmm0
  mov    -0x2f0(%rbp),%rax
  mov    %rdx,%rdi
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  mov    $0x3,%eax
  call   <scalars_minv_fts_poly@plt>
L78:
  movsd  -0x2d0(%rbp),%xmm0
  movsd  -0x170(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x320(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  -0x2c8(%rbp),%xmm0
  movsd  -0x168(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x8,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x218(%rbp),%xmm0
  movsd  -0x160(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x10,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x220(%rbp),%xmm0
  movsd  -0x158(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x18,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x2c0(%rbp),%xmm0
  movsd  -0x150(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x20,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x2b8(%rbp),%xmm0
  movsd  -0x148(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x28,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x208(%rbp),%xmm0
  movsd  -0x140(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x30,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x210(%rbp),%xmm0
  movsd  -0x138(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x38,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  cmpl   $0x0,-0x2e4(%rbp)
  je     L89
  movsd  -0x2b0(%rbp),%xmm0
  movsd  -0x130(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x40,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x2a8(%rbp),%xmm0
  movsd  -0x128(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x48,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1f8(%rbp),%xmm0
  movsd  -0x120(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x50,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x200(%rbp),%xmm0
  movsd  -0x118(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x58,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x2a0(%rbp),%xmm0
  movsd  -0x110(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x60,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x298(%rbp),%xmm0
  movsd  -0x108(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x68,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1e8(%rbp),%xmm0
  movsd  -0x100(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x70,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1f0(%rbp),%xmm0
  movsd  -0xf8(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x78,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x290(%rbp),%xmm0
  movsd  -0xf0(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  sub    $0xffffffffffffff80,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x288(%rbp),%xmm0
  movsd  -0xe8(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x88,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1d8(%rbp),%xmm0
  movsd  -0xe0(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x90,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1e0(%rbp),%xmm0
  movsd  -0xd8(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x98,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  cmpl   $0x1,-0x2e4(%rbp)
  je     L39
  movsd  -0x280(%rbp),%xmm0
  movsd  -0xd0(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xa0,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x278(%rbp),%xmm0
  movsd  -0xc8(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xa8,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1c8(%rbp),%xmm0
  movsd  -0xc0(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xb0,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1d0(%rbp),%xmm0
  movsd  -0xb8(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xb8,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x270(%rbp),%xmm0
  movsd  -0xb0(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xc0,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x268(%rbp),%xmm0
  movsd  -0xa8(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xc8,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1b8(%rbp),%xmm0
  movsd  -0xa0(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xd0,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1c0(%rbp),%xmm0
  movsd  -0x98(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xd8,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x260(%rbp),%xmm0
  movsd  -0x90(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xe0,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x258(%rbp),%xmm0
  movsd  -0x88(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xe8,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1a8(%rbp),%xmm0
  movsd  -0x80(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xf0,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1b0(%rbp),%xmm0
  movsd  -0x78(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0xf8,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x250(%rbp),%xmm0
  movsd  -0x70(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x100,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x248(%rbp),%xmm0
  movsd  -0x68(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x108,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x198(%rbp),%xmm0
  movsd  -0x60(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x110,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1a0(%rbp),%xmm0
  movsd  -0x58(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x118,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x240(%rbp),%xmm0
  movsd  -0x50(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x120,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x238(%rbp),%xmm0
  movsd  -0x48(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x128,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x188(%rbp),%xmm0
  movsd  -0x40(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x130,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x190(%rbp),%xmm0
  movsd  -0x38(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x138,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x230(%rbp),%xmm0
  movsd  -0x30(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x140,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x228(%rbp),%xmm0
  movsd  -0x28(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x148,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x178(%rbp),%xmm0
  movsd  -0x20(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x150,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x180(%rbp),%xmm0
  movsd  -0x18(%rbp),%xmm1
  mov    -0x320(%rbp),%rax
  add    $0x158,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  jmp    L71
L89:
  nop
  jmp    L71
L39:
  nop
L71:
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L85
  call   <__stack_chk_fail@plt>
L85:
  leave
  ret
```

assistant: ```c
void scalars_lub_poly_full(int version, double r, double a1, double a2, struct twobody_f *f12, struct twobody_f *f21, int n, int flag_lub, double *lub) {
  int i;
  for (i = 0; i < 44; i++) {
    lub[i] = 0.0;
  }
  double res12[22];
  twobody_scalars_res(version, r, a1, a2, f12, n, flag_lub, 1, res12);
  double res21[22];
  twobody_scalars_res(version, r, a2, a1, f21, n, flag_lub, 1, res21);
  double minv[44];
  if (version == 0) {
    scalars_minv_f_poly(r, a1, a2, minv);
  } else if (version == 1) {
    scalars_minv_ft_poly(r, a1, a2, minv);
  } else {
    scalars_minv_fts_poly(r, a1, a2, minv);
  }
  lub[0] = res12[0] - minv[0];
  lub[1] = res12[1] - minv[1];
  lub[2] = res21[1] - minv[2];
  lub[3] = res21[0] - minv[3];
  lub[4] = res12[2] - minv[4];
  lub[5] = res12[3] - minv[5];
  lub[6] = res21[3] - minv[6];
  lub[7] = res21[2] - minv[7];
  if (version == 0)
    return;
  lub[8] = res12[4] - minv[8];
  lub[9] = res12[5] - minv[9];
  lub[10] = res21[5] - minv[10];
  lub[11] = res21[4] - minv[11];
  lub[12] = res12[6] - minv[12];
  lub[13] = res12[7] - minv[13];
  lub[14] = res21[7] - minv[14];
  lub[15] = res21[6] - minv[15];
  lub[16] = res12[8] - minv[16];
  lub[17] = res12[9] - minv[17];
  lub[18] = res21[9] - minv[18];
  lub[19] = res21[8] - minv[19];
  if (version == 1)
    return;
  lub[20] = res12[10] - minv[20];
  lub[21] = res12[11] - minv[21];
  lub[22] = res21[11] - minv[22];
  lub[23] = res21[10] - minv[23];
  lub[24] = res12[12] - minv[24];
  lub[25] = res12[13] - minv[25];
  lub[26] = res21[13] - minv[26];
  lub[27] = res21[12] - minv[27];
  lub[28] = res12[14] - minv[28];
  lub[29] = res12[15] - minv[29];
  lub[30] = res21[15] - minv[30];
  lub[31] = res21[14] - minv[31];
  lub[32] = res12[16] - minv[32];
  lub[33] = res12[17] - minv[33];
  lub[34] = res21[17] - minv[34];
  lub[35] = res21[16] - minv[35];
  lub[36] = res12[18] - minv[36];
  lub[37] = res12[19] - minv[37];
  lub[38] = res21[19] - minv[38];
  lub[39] = res21[18] - minv[39];
  lub[40] = res12[20] - minv[40];
  lub[41] = res12[21] - minv[41];
  lub[42] = res21[21] - minv[42];
  lub[43] = res21[20] - minv[43];
}
```