user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<rys_contr_coulomb>:
  endbr64
  push   %r15
  mov    %edi,%eax
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xe8,%rsp
  mov    0x1d8(%rsp),%r14
  movsd  %xmm0,0xa8(%rsp)
  movsd  %xmm6,0xb0(%rsp)
  movsd  %xmm7,0xb8(%rsp)
  test   %edi,%edi
  jle    L77
  mov    %rdx,%rdi
  sub    $0x1,%eax
  mov    %rsi,0x80(%rsp)
  mov    0x120(%rsp),%rsi
  lea    0x8(%rdi,%rax,8),%rax
  mov    %rdi,%r15
  mov    0x138(%rsp),%rdi
  mov    %rcx,%rdx
  mov    %rax,0xd8(%rsp)
  movapd %xmm3,%xmm0
  movapd %xmm5,%xmm8
  mov    0x128(%rsp),%eax
  mov    %rsi,0x88(%rsp)
  pxor   %xmm3,%xmm3
  mov    %r15,%r12
  mov    0x170(%rsp),%rsi
  sub    $0x1,%eax
  movsd  %xmm3,(%rsp)
  movapd %xmm0,%xmm5
  movapd %xmm4,%xmm6
  lea    0x8(%rdi,%rax,8),%rax
  movapd %xmm8,%xmm7
  mov    %r14,%r15
  mov    %r9,0x90(%rsp)
  mov    %rax,0xd0(%rsp)
  mov    0x160(%rsp),%eax
  mov    %r8,0x98(%rsp)
  sub    $0x1,%eax
  mov    %rdx,0xa0(%rsp)
  lea    0x8(%rsi,%rax,8),%rax
  mov    %rax,0xc8(%rsp)
  movslq 0x1a0(%rsp),%rax
  mov    %rax,0xc0(%rsp)
L50:
  mov    0x128(%rsp),%eax
  test   %eax,%eax
  jle    L23
  mov    0x130(%rsp),%rax
  mov    %rax,0x78(%rsp)
  mov    0x158(%rsp),%rax
  mov    %rax,0x70(%rsp)
  mov    0x150(%rsp),%rax
  mov    %rax,0x68(%rsp)
  mov    0x148(%rsp),%rax
  mov    %rax,0x60(%rsp)
  mov    0x140(%rsp),%rax
  mov    %rax,0x58(%rsp)
  mov    0x138(%rsp),%rax
  mov    %rax,0x50(%rsp)
  nopl   0x0(%rax,%rax,1)
L35:
  mov    0x160(%rsp),%edx
  test   %edx,%edx
  jle    L71
  mov    0x168(%rsp),%rax
  mov    0x198(%rsp),%r14
  mov    0x190(%rsp),%r13
  mov    0x188(%rsp),%rbp
  mov    %rax,0x48(%rsp)
  mov    0x178(%rsp),%rax
  mov    %rax,0x40(%rsp)
  mov    0x170(%rsp),%rax
  mov    %rax,0x38(%rsp)
  xchg   %ax,%ax
L21:
  mov    0x1a0(%rsp),%ecx
  xor    %ebx,%ebx
  test   %ecx,%ecx
  jle    L32
  nopl   0x0(%rax)
L83:
  mov    0x50(%rsp),%rax
  movsd  (%r12),%xmm0
  mulsd  (%rax),%xmm0
  mov    0x38(%rsp),%rax
  mulsd  (%rax),%xmm0
  mov    0x1b0(%rsp),%rax
  mulsd  (%rax,%rbx,8),%xmm0
  mov    0x1a8(%rsp),%rax
  movsd  %xmm0,0x8(%rsp)
  push   (%rax,%rbx,8)
  mov    0x1f0(%rsp),%rax
  mov    (%rax,%rbx,4),%eax
  push   %rax
  mov    0x1f0(%rsp),%rax
  mov    (%rax,%rbx,4),%eax
  push   %rax
  mov    (%r15,%rbx,4),%eax
  push   %rax
  mov    0x1d8(%rsp),%rax
  push   (%rax,%rbx,8)
  add    $0x1,%rbx
  push   0x1f8(%rsp)
  push   0x1f8(%rsp)
  push   0x1f8(%rsp)
  mov    0x88(%rsp),%rax
  push   (%rax)
  mov    (%r14),%eax
  push   %rax
  mov    0x0(%r13),%eax
  push   %rax
  mov    0x0(%rbp),%eax
  push   %rax
  mov    0xa0(%rsp),%rax
  push   (%rax)
  push   0x1e8(%rsp)
  push   0x128(%rsp)
  push   0x128(%rsp)
  mov    0xf8(%rsp),%rax
  push   (%rax)
  mov    0xe0(%rsp),%rax
  push   (%rax)
  mov    0x100(%rsp),%rax
  mov    (%rax),%r9d
  mov    0xf8(%rsp),%rax
  movsd  %xmm7,0xc0(%rsp)
  movsd  %xmm6,0xb8(%rsp)
  mov    (%rax),%r8d
  mov    0xf0(%rsp),%rax
  movsd  %xmm5,0xb0(%rsp)
  mov    (%rax),%ecx
  mov    0x110(%rsp),%rax
  movsd  (%rax),%xmm4
  mov    0x118(%rsp),%rax
  mov    (%rax),%edx
  mov    0x120(%rsp),%rax
  movsd  %xmm2,0xa8(%rsp)
  movsd  0x138(%rsp),%xmm0
  movsd  %xmm1,0xa0(%rsp)
  mov    (%rax),%esi
  mov    0x128(%rsp),%rax
  mov    (%rax),%edi
  mov    0x130(%rsp),%rax
  movsd  (%rax),%xmm3
  call   <rys_coulomb_repulsion@plt>
  mulsd  0x98(%rsp),%xmm0
  addsd  0x90(%rsp),%xmm0
  movsd  %xmm0,0x90(%rsp)
  add    $0x90,%rsp
  movsd  0x10(%rsp),%xmm1
  movsd  0x18(%rsp),%xmm2
  cmp    %rbx,0xc0(%rsp)
  movsd  0x20(%rsp),%xmm5
  movsd  0x28(%rsp),%xmm6
  movsd  0x30(%rsp),%xmm7
  jne    L83
L32:
  addq   $0x8,0x38(%rsp)
  add    $0x4,%r14
  add    $0x4,%r13
  add    $0x4,%rbp
  addq   $0x8,0x48(%rsp)
  mov    0x38(%rsp),%rax
  addq   $0x8,0x40(%rsp)
  cmp    %rax,0xc8(%rsp)
  jne    L21
L71:
  addq   $0x8,0x50(%rsp)
  mov    0x50(%rsp),%rax
  addq   $0x8,0x78(%rsp)
  addq   $0x4,0x70(%rsp)
  addq   $0x4,0x68(%rsp)
  addq   $0x4,0x60(%rsp)
  addq   $0x8,0x58(%rsp)
  cmp    0xd0(%rsp),%rax
  jne    L35
L23:
  addq   $0x8,0x80(%rsp)
  add    $0x8,%r12
  addq   $0x4,0x88(%rsp)
  addq   $0x4,0x90(%rsp)
  addq   $0x4,0x98(%rsp)
  addq   $0x8,0xa0(%rsp)
  cmp    0xd8(%rsp),%r12
  jne    L50
L55:
  movsd  (%rsp),%xmm0
  add    $0xe8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L77:
  pxor   %xmm4,%xmm4
  movsd  %xmm4,(%rsp)
  jmp    L55
```

assistant: ```c
double rys_contr_coulomb(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) {
  double val = 0.;
  int i, j, k, l;
  for (i = 0; i < lena; i++) {
    for (j = 0; j < lenb; j++) {
      for (k = 0; k < lenc; k++) {
        for (l = 0; l < lend; l++) {
          val += acoefs[i] * bcoefs[j] * ccoefs[k] * dcoefs[l] * rys_coulomb_repulsion(xa, ya, za, anorms[i], la[i], ma[i], na[i], aexps[i], xb, yb, zb, bnorms[j], lb[j], mb[j], nb[j], bexps[j], xc, yc, zc, cnorms[k], lc[k], mc[k], nc[k], cexps[k], xd, yd, zd, dnorms[l], ld[l], md[l], nd[l], dexps[l]);
        }
      }
    }
  }
  return val;
}
```