user: What is the c source code of the assembly code below?

```asm
<contr_grad_eri>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xa0,%rsp
  mov    %edi,-0x34(%rbp)
  mov    %rsi,-0x40(%rbp)
  mov    %rdx,-0x48(%rbp)
  mov    %rcx,-0x50(%rbp)
  movsd  %xmm0,-0x58(%rbp)
  movsd  %xmm1,-0x60(%rbp)
  movsd  %xmm2,-0x68(%rbp)
  mov    %r8d,-0x38(%rbp)
  mov    %r9d,-0x6c(%rbp)
  movsd  %xmm3,-0x78(%rbp)
  movsd  %xmm4,-0x80(%rbp)
  movsd  %xmm5,-0x88(%rbp)
  movsd  %xmm6,-0x90(%rbp)
  movsd  %xmm7,-0x98(%rbp)
  cmpl   $0x2,0xe0(%rbp)
  je     L247
  cmpl   $0x2,0xe0(%rbp)
  jg     L110
  cmpl   $0x0,0xe0(%rbp)
  je     L53
  cmpl   $0x1,0xe0(%rbp)
  je     L6
  jmp    L110
L53:
  movl   $0x1,-0x1c(%rbp)
  movl   $0x0,-0x18(%rbp)
  movl   $0x0,-0x14(%rbp)
  jmp    L127
L6:
  movl   $0x0,-0x1c(%rbp)
  movl   $0x1,-0x18(%rbp)
  movl   $0x0,-0x14(%rbp)
  jmp    L127
L247:
  movl   $0x0,-0x1c(%rbp)
  movl   $0x0,-0x18(%rbp)
  movl   $0x1,-0x14(%rbp)
  jmp    L127
L110:
  lea    D6(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
  mov    $0xffffffff,%edi
  call   <exit@plt>
L127:
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x1c(%rbp),%eax
  imul   -0x38(%rbp),%eax
  mov    %eax,%edx
  mov    -0x18(%rbp),%eax
  imul   -0x6c(%rbp),%eax
  add    %eax,%edx
  mov    -0x14(%rbp),%eax
  imul   0x10(%rbp),%eax
  add    %edx,%eax
  test   %eax,%eax
  jle    L26
  mov    -0x34(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x2c(%rbp)
  jmp    L233
L85:
  mov    0x18(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x28(%rbp)
  jmp    L211
L78:
  mov    0x50(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x24(%rbp)
  jmp    L235
L79:
  mov    0x90(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x20(%rbp)
  jmp    L108
L36:
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x20(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0xa0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,-0xa0(%rbp)
  mov    -0x20(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x98(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm5
  mov    -0x20(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm4
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm3
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm11
  mov    0x10(%rbp),%eax
  sub    -0x14(%rbp),%eax
  mov    %eax,%edx
  mov    -0x6c(%rbp),%eax
  sub    -0x18(%rbp),%eax
  mov    %eax,%esi
  mov    -0x38(%rbp),%eax
  sub    -0x1c(%rbp),%eax
  mov    %eax,%edi
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x50(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rax),%xmm10
  mov    0x48(%rbp),%r9d
  mov    0x40(%rbp),%r10d
  mov    0x38(%rbp),%ecx
  movsd  -0x88(%rbp),%xmm7
  movsd  -0x80(%rbp),%xmm6
  movsd  -0x78(%rbp),%xmm12
  movsd  -0x68(%rbp),%xmm9
  movsd  -0x60(%rbp),%xmm8
  mov    -0x58(%rbp),%rax
  lea    -0x8(%rsp),%rsp
  movsd  %xmm5,(%rsp)
  mov    0xd8(%rbp),%r8d
  push   %r8
  mov    0xd0(%rbp),%r8d
  push   %r8
  mov    0xc8(%rbp),%r8d
  push   %r8
  lea    -0x8(%rsp),%rsp
  movsd  %xmm4,(%rsp)
  push   0xc0(%rbp)
  push   0xb8(%rbp)
  push   0xb0(%rbp)
  lea    -0x8(%rsp),%rsp
  movsd  %xmm3,(%rsp)
  mov    0x88(%rbp),%r8d
  push   %r8
  mov    0x80(%rbp),%r8d
  push   %r8
  mov    0x78(%rbp),%r8d
  push   %r8
  lea    -0x8(%rsp),%rsp
  movsd  %xmm2,(%rsp)
  push   0x70(%rbp)
  push   -0x98(%rbp)
  push   -0x90(%rbp)
  lea    -0x8(%rsp),%rsp
  movsd  %xmm1,(%rsp)
  lea    -0x8(%rsp),%rsp
  movsd  %xmm0,(%rsp)
  mov    %r10d,%r8d
  movapd %xmm12,%xmm5
  movapd %xmm11,%xmm4
  movapd %xmm10,%xmm3
  movapd %xmm9,%xmm2
  movapd %xmm8,%xmm1
  movq   %rax,%xmm0
  call   <eri@plt>
  add    $0x90,%rsp
  mulsd  -0xa0(%rbp),%xmm0
  movsd  -0x10(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  subl   $0x1,-0x20(%rbp)
L108:
  cmpl   $0x0,-0x20(%rbp)
  jns    L36
  subl   $0x1,-0x24(%rbp)
L235:
  cmpl   $0x0,-0x24(%rbp)
  jns    L79
  subl   $0x1,-0x28(%rbp)
L211:
  cmpl   $0x0,-0x28(%rbp)
  jns    L78
  subl   $0x1,-0x2c(%rbp)
L233:
  cmpl   $0x0,-0x2c(%rbp)
  jns    L85
  mov    -0x1c(%rbp),%eax
  imul   -0x38(%rbp),%eax
  mov    %eax,%edx
  mov    -0x18(%rbp),%eax
  imul   -0x6c(%rbp),%eax
  add    %eax,%edx
  mov    -0x14(%rbp),%eax
  imul   0x10(%rbp),%eax
  add    %edx,%eax
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  -0x10(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
L26:
  movsd  -0x10(%rbp),%xmm0
  movq   D2(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x34(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x2c(%rbp)
  jmp    L95
L129:
  mov    0x18(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x28(%rbp)
  jmp    L237
L138:
  mov    0x50(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x24(%rbp)
  jmp    L161
L76:
  mov    0x90(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x20(%rbp)
  jmp    L146
L179:
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x20(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0xa0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,-0xa0(%rbp)
  mov    -0x20(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x98(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm5
  mov    -0x20(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm4
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm3
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm11
  mov    0x10(%rbp),%edx
  mov    -0x14(%rbp),%eax
  lea    (%rdx,%rax,1),%r10d
  mov    -0x6c(%rbp),%edx
  mov    -0x18(%rbp),%eax
  lea    (%rdx,%rax,1),%esi
  mov    -0x38(%rbp),%edx
  mov    -0x1c(%rbp),%eax
  lea    (%rdx,%rax,1),%edi
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm10
  mov    0x48(%rbp),%r9d
  mov    0x40(%rbp),%r8d
  mov    0x38(%rbp),%edx
  movsd  -0x88(%rbp),%xmm7
  movsd  -0x80(%rbp),%xmm6
  movsd  -0x78(%rbp),%xmm12
  movsd  -0x68(%rbp),%xmm9
  movsd  -0x60(%rbp),%xmm8
  mov    -0x58(%rbp),%rax
  lea    -0x8(%rsp),%rsp
  movsd  %xmm5,(%rsp)
  mov    0xd8(%rbp),%ecx
  push   %rcx
  mov    0xd0(%rbp),%ecx
  push   %rcx
  mov    0xc8(%rbp),%ecx
  push   %rcx
  lea    -0x8(%rsp),%rsp
  movsd  %xmm4,(%rsp)
  push   0xc0(%rbp)
  push   0xb8(%rbp)
  push   0xb0(%rbp)
  lea    -0x8(%rsp),%rsp
  movsd  %xmm3,(%rsp)
  mov    0x88(%rbp),%ecx
  push   %rcx
  mov    0x80(%rbp),%ecx
  push   %rcx
  mov    0x78(%rbp),%ecx
  push   %rcx
  lea    -0x8(%rsp),%rsp
  movsd  %xmm2,(%rsp)
  push   0x70(%rbp)
  push   -0x98(%rbp)
  push   -0x90(%rbp)
  lea    -0x8(%rsp),%rsp
  movsd  %xmm1,(%rsp)
  lea    -0x8(%rsp),%rsp
  movsd  %xmm0,(%rsp)
  mov    %edx,%ecx
  movapd %xmm12,%xmm5
  movapd %xmm11,%xmm4
  mov    %r10d,%edx
  movapd %xmm10,%xmm3
  movapd %xmm9,%xmm2
  movapd %xmm8,%xmm1
  movq   %rax,%xmm0
  call   <eri@plt>
  add    $0x90,%rsp
  mulsd  -0xa0(%rbp),%xmm0
  movsd  -0x10(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  subl   $0x1,-0x20(%rbp)
L146:
  cmpl   $0x0,-0x20(%rbp)
  jns    L179
  subl   $0x1,-0x24(%rbp)
L161:
  cmpl   $0x0,-0x24(%rbp)
  jns    L76
  subl   $0x1,-0x28(%rbp)
L237:
  cmpl   $0x0,-0x28(%rbp)
  jns    L138
  subl   $0x1,-0x2c(%rbp)
L95:
  cmpl   $0x0,-0x2c(%rbp)
  jns    L129
  movsd  -0x8(%rbp),%xmm0
  addsd  -0x10(%rbp),%xmm0
  movsd  -0x10(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: read_data
arguments: data_tag: D6
data_type: string
name: read_data
arguments: data_tag: D2
data_type: f32[2]

tool: D6: contr_grad_eri: error unrecognized direction flag


tool: {"D2": [0.0,-0.0]}

assistant: ```c
double contr_grad_eri(int lena, double *aexps, double *acoefs, double *anorms, double xa, double ya, double za, int la, int ma, int na, int lenb, double *bexps, double *bcoefs, double *bnorms, double xb, double yb, double zb, int lb, int mb, int nb, int lenc, double *cexps, double *ccoefs, double *cnorms, double xc, double yc, double zc, int lc, int mc, int nc, int lend, double *dexps, double *dcoefs, double *dnorms, double xd, double yd, double zd, int ld, int md, int nd, int direction) {
  double val;
  int p, q, r, s;
  int l, m, n;
  double EE;
  switch (direction) {
  case 0:
    l = 1;
    m = 0;
    n = 0;
    break;
  case 1:
    l = 0;
    m = 1;
    n = 0;
    break;
  case 2:
    l = 0;
    m = 0;
    n = 1;
    break;
  default:
    printf("contr_grad_eri: error unrecognized direction flag\n");
    exit(-1);
  }
  val = 0.0;
  if ((l * la + m * ma + n * na) > 0) {
    for (p = lena - 1; p >= 0; p--)
      for (q = lenb - 1; q >= 0; q--)
        for (r = lenc - 1; r >= 0; r--)
          for (s = lend - 1; s >= 0; s--)
            val += acoefs[p] * bcoefs[q] * ccoefs[r] * dcoefs[s] * eri(xa, ya, za, anorms[p], la - l, ma - m, na - n, aexps[p], xb, yb, zb, bnorms[q], lb, mb, nb, bexps[q], xc, yc, zc, cnorms[r], lc, mc, nc, cexps[r], xd, yd, zd, dnorms[s], ld, md, nd, dexps[s]);
    ;
    val = val * (l * la + m * ma + n * na);
  }
  EE = -val;
  val = 0.0;
  for (p = lena - 1; p >= 0; p--)
    for (q = lenb - 1; q >= 0; q--)
      for (r = lenc - 1; r >= 0; r--)
        for (s = lend - 1; s >= 0; s--)
          val += aexps[p] * acoefs[p] * bcoefs[q] * ccoefs[r] * dcoefs[s] * eri(xa, ya, za, anorms[p], la + l, ma + m, na + n, aexps[p], xb, yb, zb, bnorms[q], lb, mb, nb, bexps[q], xc, yc, zc, cnorms[r], lc, mc, nc, cexps[r], xd, yd, zd, dnorms[s], ld, md, nd, dexps[s]);
  ;
  EE = EE + val + val;
  return EE;
}
```