user: Please show me the C implementation that matches the assembly code listed below:

```asm
<twobody_XC_slip_16>:
  endbr64
  lea    0x8(%rdx),%rdi
  mov    %rdx,%rcx
  push   %rbp
  xor    %eax,%eax
  push   %rbx
  and    $0xfffffffffffffff8,%rdi
  mov    %rsi,%rbx
  mov    $0x5,%esi
  sub    %rdi,%rcx
  mov    %rdx,%rbp
  add    $0x88,%ecx
  sub    $0x28,%rsp
  shr    $0x3,%ecx
  movq   $0x0,(%rdx)
  movq   $0x0,0x80(%rdx)
  movsd  %xmm0,0x18(%rsp)
  rep stos %rax,%es:(%rdi)
  mov    $0xfffffffe,%edi
  call   <SL_G1@plt>
  xor    %edi,%edi
  mov    $0x3,%esi
  movsd  %xmm0,(%rsp)
  call   <SL_G1@plt>
  movsd  (%rsp),%xmm2
  xor    %edi,%edi
  mov    $0x3,%esi
  mulsd  D52(%rip),%xmm2
  movapd %xmm0,%xmm3
  movapd %xmm2,%xmm0
  mulsd  %xmm3,%xmm0
  movsd  %xmm0,(%rsp)
  call   <SL_G1@plt>
  movsd  (%rsp),%xmm5
  xor    %edi,%edi
  mov    $0x3,%esi
  mulsd  %xmm0,%xmm5
  movsd  %xmm5,(%rsp)
  call   <SL_G2@plt>
  mulsd  (%rsp),%xmm0
  xor    %edi,%edi
  mov    $0x3,%esi
  movsd  %xmm0,(%rsp)
  call   <SL_G2@plt>
  mulsd  (%rsp),%xmm0
  mov    $0x4,%esi
  mov    $0xffffffff,%edi
  movsd  %xmm0,0x30(%rbp)
  call   <SL_G1@plt>
  mov    $0x4,%esi
  mov    $0xffffffff,%edi
  movsd  %xmm0,(%rsp)
  call   <SL_G2@plt>
  movsd  (%rsp),%xmm2
  xor    %edi,%edi
  mov    $0x3,%esi
  mulsd  D25(%rip),%xmm2
  movapd %xmm0,%xmm3
  movapd %xmm2,%xmm0
  mulsd  %xmm3,%xmm0
  movsd  %xmm0,(%rsp)
  call   <SL_G1@plt>
  movsd  (%rsp),%xmm4
  xor    %edi,%edi
  mov    $0x3,%esi
  mulsd  %xmm0,%xmm4
  movsd  %xmm4,(%rsp)
  call   <SL_G1@plt>
  mulsd  (%rsp),%xmm0
  xor    %edi,%edi
  mov    $0x3,%esi
  movsd  %xmm0,(%rsp)
  call   <SL_G2@plt>
  mulsd  (%rsp),%xmm0
  xor    %edi,%edi
  mov    $0x3,%esi
  movsd  %xmm0,0x40(%rbp)
  call   <SL_G1@plt>
  xor    %edi,%edi
  mov    $0x3,%esi
  movsd  %xmm0,(%rsp)
  call   <SL_G1@plt>
  movsd  (%rsp),%xmm2
  xor    %edi,%edi
  mov    $0x3,%esi
  mulsd  D7(%rip),%xmm2
  movapd %xmm0,%xmm3
  movapd %xmm2,%xmm0
  mulsd  %xmm3,%xmm0
  movsd  %xmm0,(%rsp)
  call   <SL_G1@plt>
  movsd  (%rsp),%xmm6
  mov    $0x5,%esi
  mov    $0xfffffffe,%edi
  mulsd  %xmm0,%xmm6
  movsd  %xmm6,(%rsp)
  call   <SL_G2@plt>
  xor    %edi,%edi
  mov    $0x3,%esi
  movsd  %xmm0,0x8(%rsp)
  call   <SL_G2@plt>
  movsd  0x8(%rsp),%xmm2
  mov    $0x4,%esi
  mov    $0xffffffff,%edi
  mulsd  D1(%rip),%xmm2
  movapd %xmm0,%xmm3
  movapd %xmm2,%xmm0
  mulsd  %xmm3,%xmm0
  movsd  %xmm0,0x8(%rsp)
  call   <SL_G2@plt>
  mov    $0x4,%esi
  mov    $0xffffffff,%edi
  movsd  %xmm0,0x10(%rsp)
  call   <SL_G2@plt>
  movsd  0x10(%rsp),%xmm2
  mov    $0x8,%esi
  mov    $0xfffffffb,%edi
  mulsd  D17(%rip),%xmm2
  movapd %xmm0,%xmm3
  movapd %xmm2,%xmm0
  mulsd  %xmm3,%xmm0
  addsd  0x8(%rsp),%xmm0
  mulsd  (%rsp),%xmm0
  movsd  %xmm0,0x50(%rbp)
  call   <SL_G2@plt>
  mov    $0x3,%esi
  xor    %edi,%edi
  movsd  %xmm0,(%rsp)
  call   <SL_G1@plt>
  movsd  (%rsp),%xmm2
  mov    $0x3,%esi
  xor    %edi,%edi
  mulsd  D24(%rip),%xmm2
  movapd %xmm0,%xmm3
  movapd %xmm2,%xmm0
  mulsd  %xmm3,%xmm0
  movsd  %xmm0,(%rsp)
  call   <SL_G1@plt>
  mulsd  (%rsp),%xmm0
  movsd  0x18(%rsp),%xmm1
  movsd  %xmm0,0x68(%rbp)
  pxor   %xmm0,%xmm0
  movq   $0x0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x80(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x78(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x70(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x68(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x60(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x58(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x50(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x48(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x40(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x38(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x30(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x28(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x20(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x18(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x10(%rbp),%xmm0
  movsd  %xmm0,0x80(%rbx)
  mulsd  %xmm1,%xmm0
  addsd  0x8(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  %xmm0,0x80(%rbx)
  addsd  0x0(%rbp),%xmm1
  movsd  %xmm1,0x80(%rbx)
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  ret
```

tool_calls: name: disassemble_data
arguments: tag: D52
data_type: f64
name: disassemble_data
arguments: tag: D25
data_type: f64
name: disassemble_data
arguments: tag: D7
data_type: f64
name: disassemble_data
arguments: tag: D1
data_type: f64
name: disassemble_data
arguments: tag: D17
data_type: f64
name: disassemble_data
arguments: tag: D24
data_type: f64

tool: D52: 393216

tool: D25: 1572864

tool: D7: 196608

tool: D1: 4

tool: D17: 3

tool: D24: 1376256

assistant: ```c
void twobody_XC_slip_16(int n, double l, double *f, double *fk) {
  int q;
  for (q = 0; q <= 16; q++)
    fk[q] = 0.0;
  fk[6] = 393216.0 * SL_G1(-2, 5) * SL_G1(0, 3) * SL_G1(0, 3) * SL_G2(0, 3) * SL_G2(0, 3);
  fk[8] = 1572864.0 * SL_G1(-1, 4) * SL_G2(-1, 4) * SL_G1(0, 3) * SL_G1(0, 3) * SL_G2(0, 3);
  fk[10] = 196608.0 * SL_G1(0, 3) * SL_G1(0, 3) * SL_G1(0, 3) * (4.0 * SL_G2(-2, 5) * SL_G2(0, 3) + 3.0 * SL_G2(-1, 4) * SL_G2(-1, 4));
  fk[13] = 1376256.0 * SL_G2(-5, 8) * SL_G1(0, 3) * SL_G1(0, 3);
  f[16] = 0.0;
  for (q = 0; q <= 16; q++)
    f[16] = f[16] * l + fk[16 - q];
}
```