user: Converte this assembly snippet into its equivalent C code:

```asm
<three_center_1D>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xa0,%rsp
  movsd  %xmm0,-0x58(%rbp)
  mov    %edi,-0x5c(%rbp)
  movsd  %xmm1,-0x68(%rbp)
  movsd  %xmm2,-0x70(%rbp)
  mov    %esi,-0x60(%rbp)
  movsd  %xmm3,-0x78(%rbp)
  movsd  %xmm4,-0x80(%rbp)
  mov    %edx,-0x84(%rbp)
  movsd  %xmm5,-0x90(%rbp)
  movsd  -0x68(%rbp),%xmm0
  addsd  -0x78(%rbp),%xmm0
  movsd  -0x90(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x68(%rbp),%xmm0
  movq   D14(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mulsd  -0x78(%rbp),%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0x58(%rbp),%xmm0
  subsd  -0x70(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    $0x2,%esi
  mov    %eax,%edi
  call   <pow@plt>
  mulsd  -0x98(%rbp),%xmm0
  divsd  -0x30(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0x68(%rbp),%xmm0
  movq   D14(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0x90(%rbp),%xmm2
  movsd  %xmm2,-0xa0(%rbp)
  movsd  -0x58(%rbp),%xmm0
  subsd  -0x80(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    $0x2,%esi
  mov    %eax,%edi
  call   <pow@plt>
  mulsd  -0xa0(%rbp),%xmm0
  movapd %xmm0,%xmm3
  divsd  -0x30(%rbp),%xmm3
  movq   %xmm3,%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm4
  mulsd  -0x98(%rbp),%xmm4
  movsd  %xmm4,-0x98(%rbp)
  movsd  -0x78(%rbp),%xmm0
  movq   D14(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mulsd  -0x90(%rbp),%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  movsd  -0x70(%rbp),%xmm0
  subsd  -0x80(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    $0x2,%esi
  mov    %eax,%edi
  call   <pow@plt>
  mulsd  -0xa0(%rbp),%xmm0
  divsd  -0x30(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  mulsd  -0x98(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x68(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x58(%rbp),%xmm1
  movsd  -0x78(%rbp),%xmm0
  mulsd  -0x70(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x90(%rbp),%xmm0
  mulsd  -0x80(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  divsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x20(%rbp),%xmm0
  subsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x20(%rbp),%xmm0
  subsd  -0x70(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x20(%rbp),%xmm0
  subsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movl   $0x0,-0x48(%rbp)
  jmp    L20
L2:
  movl   $0x0,-0x44(%rbp)
  jmp    L63
L38:
  movl   $0x0,-0x40(%rbp)
  jmp    L39
L54:
  mov    -0x48(%rbp),%edx
  mov    -0x44(%rbp),%eax
  add    %eax,%edx
  mov    -0x40(%rbp),%eax
  add    %edx,%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  mov    %eax,-0x3c(%rbp)
  mov    -0x48(%rbp),%edx
  mov    -0x44(%rbp),%eax
  add    %eax,%edx
  mov    -0x40(%rbp),%eax
  add    %edx,%eax
  and    $0x1,%eax
  test   %eax,%eax
  jne    L59
  mov    -0x48(%rbp),%edx
  mov    -0x5c(%rbp),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <binomial@plt>
  movsd  %xmm0,-0x98(%rbp)
  mov    -0x44(%rbp),%edx
  mov    -0x60(%rbp),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <binomial@plt>
  movapd %xmm0,%xmm7
  mulsd  -0x98(%rbp),%xmm7
  movsd  %xmm7,-0x98(%rbp)
  mov    -0x40(%rbp),%edx
  mov    -0x84(%rbp),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <binomial@plt>
  mulsd  -0x98(%rbp),%xmm0
  movsd  %xmm0,-0x98(%rbp)
  mov    -0x5c(%rbp),%eax
  sub    -0x48(%rbp),%eax
  mov    %eax,%edx
  movsd  -0x18(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <pow@plt>
  movapd %xmm0,%xmm3
  mulsd  -0x98(%rbp),%xmm3
  movsd  %xmm3,-0x98(%rbp)
  mov    -0x60(%rbp),%eax
  sub    -0x44(%rbp),%eax
  mov    %eax,%edx
  movsd  -0x10(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <pow@plt>
  mulsd  -0x98(%rbp),%xmm0
  movsd  %xmm0,-0x98(%rbp)
  mov    -0x84(%rbp),%eax
  sub    -0x40(%rbp),%eax
  mov    %eax,%edx
  movsd  -0x8(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <pow@plt>
  movapd %xmm0,%xmm6
  mulsd  -0x98(%rbp),%xmm6
  movsd  %xmm6,-0x98(%rbp)
  mov    -0x3c(%rbp),%eax
  add    %eax,%eax
  sub    $0x1,%eax
  mov    %eax,%edi
  call   <fact2@plt>
  mulsd  -0x98(%rbp),%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0x30(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  cvttsd2si %xmm0,%eax
  mov    -0x3c(%rbp),%edx
  mov    %edx,%esi
  mov    %eax,%edi
  call   <pow@plt>
  movsd  -0x98(%rbp),%xmm6
  divsd  %xmm0,%xmm6
  movsd  %xmm6,-0x98(%rbp)
  movsd  D17(%rip),%xmm0
  divsd  -0x30(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  mulsd  -0x98(%rbp),%xmm0
  movsd  -0x38(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
L59:
  addl   $0x1,-0x40(%rbp)
L39:
  mov    -0x84(%rbp),%eax
  cmp    -0x40(%rbp),%eax
  jge    L54
  addl   $0x1,-0x44(%rbp)
L63:
  mov    -0x60(%rbp),%eax
  cmp    -0x44(%rbp),%eax
  jge    L38
  addl   $0x1,-0x48(%rbp)
L20:
  mov    -0x5c(%rbp),%eax
  cmp    -0x48(%rbp),%eax
  jge    L2
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x38(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: parse_data
arguments: label: D14
data_type: byte[8]
name: parse_data
arguments: label: D17
data_type: f64

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

tool: D17: 3.141592653589793

assistant: ```c
double three_center_1D(double xi, int ai, double alphai, double xj, int aj, double alphaj, double xk, int ak, double alphak) {
  double gamma, dx, px, xpi, xpj, xpk, intgl;
  int q, r, s, n;
  gamma = alphai + alphaj + alphak;
  dx = exp(-alphai * alphaj * pow(xi - xj, 2) / gamma) * exp(-alphai * alphak * pow(xi - xk, 2) / gamma) * exp(-alphaj * alphak * pow(xj - xk, 2) / gamma);
  px = (alphai * xi + alphaj * xj + alphak * xk) / gamma;
  xpi = px - xi;
  xpj = px - xj;
  xpk = px - xk;
  intgl = 0;
  for (q = 0; q < ai + 1; q++) {
    for (r = 0; r < aj + 1; r++) {
      for (s = 0; s < ak + 1; s++) {
        n = (q + r + s) / 2;
        if ((q + r + s) % 2 == 0) {
          intgl += binomial(ai, q) * binomial(aj, r) * binomial(ak, s) * pow(xpi, ai - q) * pow(xpj, aj - r) * pow(xpk, ak - s) * fact2(2 * n - 1) / pow(2 * gamma, n) * sqrt(3.14159265358979323846 / gamma);
        }
      }
    }
  }
  return dx * intgl;
}
```