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

```asm
<overlap>:
  endbr64
  push   %r15
  movapd %xmm3,%xmm8
  movapd %xmm2,%xmm9
  mov    %edi,%r15d
  push   %r14
  movapd %xmm1,%xmm11
  movapd %xmm5,%xmm10
  mov    %r8d,%r14d
  push   %r13
  movapd %xmm7,%xmm5
  movapd %xmm8,%xmm2
  movapd %xmm9,%xmm1
  push   %r12
  movapd %xmm10,%xmm3
  mov    %r9d,%r12d
  mov    %esi,%r13d
  push   %rbp
  mov    %edx,%ebp
  push   %rbx
  mov    %ecx,%ebx
  sub    $0x88,%rsp
  movsd  %xmm0,0x10(%rsp)
  movapd %xmm11,%xmm0
  movsd  %xmm4,0x18(%rsp)
  movapd %xmm6,%xmm4
  movsd  %xmm7,0x78(%rsp)
  movsd  %xmm8,0x38(%rsp)
  movsd  %xmm9,0x30(%rsp)
  movsd  %xmm6,0x48(%rsp)
  movsd  %xmm10,0x40(%rsp)
  movsd  %xmm11,0x20(%rsp)
  call   <dist2@plt>
  movsd  0x10(%rsp),%xmm5
  movsd  0x18(%rsp),%xmm6
  movsd  0x40(%rsp),%xmm10
  movsd  0x20(%rsp),%xmm11
  movsd  %xmm0,0x28(%rsp)
  movapd %xmm5,%xmm2
  movapd %xmm5,%xmm0
  addsd  %xmm6,%xmm2
  movapd %xmm10,%xmm3
  movapd %xmm11,%xmm1
  movsd  %xmm10,0x70(%rsp)
  movsd  %xmm11,0x68(%rsp)
  movsd  %xmm2,0x8(%rsp)
  movapd %xmm6,%xmm2
  call   <product_center_1D@plt>
  movsd  0x48(%rsp),%xmm6
  movsd  0x30(%rsp),%xmm9
  movsd  0x18(%rsp),%xmm2
  movsd  %xmm0,0x60(%rsp)
  movsd  0x10(%rsp),%xmm0
  movapd %xmm6,%xmm3
  movapd %xmm9,%xmm1
  movsd  %xmm6,0x58(%rsp)
  movsd  %xmm9,0x50(%rsp)
  call   <product_center_1D@plt>
  movsd  0x38(%rsp),%xmm8
  movsd  0x78(%rsp),%xmm7
  movsd  0x18(%rsp),%xmm2
  movsd  %xmm0,0x48(%rsp)
  movsd  0x10(%rsp),%xmm0
  movapd %xmm8,%xmm1
  movapd %xmm7,%xmm3
  movsd  %xmm7,0x40(%rsp)
  call   <product_center_1D@plt>
  movsd  D20(%rip),%xmm1
  movsd  %xmm0,0x30(%rsp)
  movsd  D2(%rip),%xmm0
  divsd  0x8(%rsp),%xmm0
  call   <pow@plt>
  movsd  0x10(%rsp),%xmm5
  xorpd  D17(%rip),%xmm5
  movsd  %xmm0,0x20(%rsp)
  movsd  0x18(%rsp),%xmm0
  mulsd  %xmm5,%xmm0
  mulsd  0x28(%rsp),%xmm0
  divsd  0x8(%rsp),%xmm0
  call   <exp@plt>
  mov    %ebx,%esi
  mov    %r15d,%edi
  movsd  0x20(%rsp),%xmm7
  movsd  0x60(%rsp),%xmm5
  movsd  0x68(%rsp),%xmm11
  movsd  0x70(%rsp),%xmm10
  mulsd  %xmm0,%xmm7
  movapd %xmm5,%xmm1
  movapd %xmm5,%xmm0
  movsd  0x8(%rsp),%xmm2
  subsd  %xmm10,%xmm1
  subsd  %xmm11,%xmm0
  movsd  %xmm7,0x10(%rsp)
  call   <overlap_1D@plt>
  movsd  0x48(%rsp),%xmm4
  mov    %r14d,%esi
  movsd  0x50(%rsp),%xmm9
  movsd  0x58(%rsp),%xmm6
  movsd  0x8(%rsp),%xmm2
  movsd  %xmm0,0x18(%rsp)
  mov    %r13d,%edi
  movapd %xmm4,%xmm3
  movapd %xmm4,%xmm1
  subsd  %xmm9,%xmm3
  subsd  %xmm6,%xmm1
  movapd %xmm3,%xmm0
  call   <overlap_1D@plt>
  movsd  0x30(%rsp),%xmm3
  mov    %r12d,%esi
  mov    %ebp,%edi
  movsd  0x40(%rsp),%xmm7
  movsd  0x38(%rsp),%xmm8
  movsd  %xmm0,0x20(%rsp)
  movapd %xmm3,%xmm5
  movsd  0x8(%rsp),%xmm2
  subsd  %xmm7,%xmm3
  subsd  %xmm8,%xmm5
  movapd %xmm3,%xmm1
  movapd %xmm5,%xmm0
  call   <overlap_1D@plt>
  movsd  0x10(%rsp),%xmm7
  mulsd  0x18(%rsp),%xmm7
  movapd %xmm0,%xmm1
  movsd  0x20(%rsp),%xmm0
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  mulsd  %xmm7,%xmm0
  mulsd  %xmm1,%xmm0
  ret
```

tool_calls: name: extract_symbolic_data
arguments: label: D20
type_infered: f64
name: extract_symbolic_data
arguments: label: D2
type_infered: f64
name: extract_symbolic_data
arguments: label: D17
type_infered: byte[16]

tool: D20: 1.5

tool: D2: 3.141592653589793

tool: {"D17": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

assistant: ```c
double overlap(double alpha1, int l1, int m1, int n1, double xa, double ya, double za, double alpha2, int l2, int m2, int n2, double xb, double yb, double zb) {
  double rab2, gamma, xp, yp, zp, pre, wx, wy, wz;
  rab2 = dist2(xa, ya, za, xb, yb, zb);
  gamma = alpha1 + alpha2;
  xp = product_center_1D(alpha1, xa, alpha2, xb);
  yp = product_center_1D(alpha1, ya, alpha2, yb);
  zp = product_center_1D(alpha1, za, alpha2, zb);
  pre = pow(3.14159265358979323846 / gamma, 1.5) * exp(-alpha1 * alpha2 * rab2 / gamma);
  wx = overlap_1D(l1, l2, xp - xa, xp - xb, gamma);
  wy = overlap_1D(m1, m2, yp - ya, yp - yb, gamma);
  wz = overlap_1D(n1, n2, zp - za, zp - zb, gamma);
  return pre * wx * wy * wz;
}
```