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

```asm
<md_kernel>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x18,%rsp
  mov    %rdi,-0x68(%rbp)
  mov    %rsi,-0x70(%rbp)
  mov    %rdx,-0x78(%rbp)
  mov    %rcx,-0x80(%rbp)
  mov    %r8,-0x88(%rbp)
  mov    %r9,-0x90(%rbp)
  mov    0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x48(%rbp)
  mov    0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x44(%rbp)
  mov    0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x40(%rbp)
  movl   $0x0,-0x54(%rbp)
  movl   $0x0,-0x50(%rbp)
  movl   $0x0,-0x4c(%rbp)
  movl   $0x0,-0x58(%rbp)
  jmp    L9
L11:
  mov    0x18(%rbp),%eax
  shl    $0x5,%eax
  mov    %eax,%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x3c(%rbp)
  mov    -0x3c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x38(%rbp)
  mov    -0x3c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x34(%rbp)
  mov    -0x3c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x30(%rbp)
  mov    -0x48(%rbp),%eax
  sub    -0x38(%rbp),%eax
  mov    %eax,-0x2c(%rbp)
  mov    -0x44(%rbp),%eax
  sub    -0x34(%rbp),%eax
  mov    %eax,-0x28(%rbp)
  mov    -0x40(%rbp),%eax
  sub    -0x30(%rbp),%eax
  mov    %eax,-0x24(%rbp)
  mov    -0x2c(%rbp),%eax
  imul   %eax,%eax
  mov    %eax,-0x20(%rbp)
  mov    -0x28(%rbp),%eax
  imul   %eax,%eax
  mov    %eax,-0x1c(%rbp)
  mov    -0x24(%rbp),%eax
  imul   %eax,%eax
  mov    %eax,-0x18(%rbp)
  mov    -0x20(%rbp),%edx
  mov    -0x1c(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x14(%rbp)
  mov    -0x14(%rbp),%edx
  mov    -0x18(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x10(%rbp)
  mov    -0x10(%rbp),%eax
  mov    %eax,-0xc(%rbp)
  mov    -0xc(%rbp),%eax
  imul   %eax,%eax
  mov    -0xc(%rbp),%edx
  imul   %edx,%eax
  mov    %eax,-0x8(%rbp)
  mov    -0xc(%rbp),%eax
  imul   -0x8(%rbp),%eax
  mov    %eax,%ecx
  mov    -0x8(%rbp),%edx
  mov    %edx,%eax
  shl    $0x4,%eax
  sub    %edx,%eax
  sub    $0x14,%eax
  imul   %ecx,%eax
  mov    %eax,-0x4(%rbp)
  mov    -0x2c(%rbp),%eax
  imul   -0x4(%rbp),%eax
  add    %eax,-0x54(%rbp)
  mov    -0x28(%rbp),%eax
  imul   -0x4(%rbp),%eax
  add    %eax,-0x50(%rbp)
  mov    -0x24(%rbp),%eax
  imul   -0x4(%rbp),%eax
  add    %eax,-0x4c(%rbp)
  addl   $0x1,-0x58(%rbp)
L9:
  cmpl   $0x1f,-0x58(%rbp)
  jle    L11
  mov    0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x68(%rbp),%rax
  add    %rax,%rdx
  mov    -0x54(%rbp),%eax
  mov    %eax,(%rdx)
  mov    0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x70(%rbp),%rax
  add    %rax,%rdx
  mov    -0x50(%rbp),%eax
  mov    %eax,(%rdx)
  mov    0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x78(%rbp),%rax
  add    %rax,%rdx
  mov    -0x4c(%rbp),%eax
  mov    %eax,(%rdx)
  nop
  leave
  ret
```

assistant: ```c
void md_kernel(int d_force_x[32], int d_force_y[32], int d_force_z[32], int position_x[32], int position_y[32], int position_z[32], int NL[32 * 32], int i) {
  int j, jidx;
  int delx, dely, delz, r2inv, r2invTEMP, r2invTEMP2, r2invTEMP3, t1, t2, t3;
  int i_x = position_x[i];
  int i_y = position_y[i];
  int i_z = position_z[i];
  int fx = 0;
  int fy = 0;
  int fz = 0;
loop_j:
  for (j = 0; j < 32; j++) {
    jidx = NL[i * 32 + j];
    int j_x = position_x[jidx];
    int j_y = position_y[jidx];
    int j_z = position_z[jidx];
    delx = i_x - j_x;
    dely = i_y - j_y;
    delz = i_z - j_z;
    r2invTEMP = delx * delx;
    r2invTEMP2 = dely * dely;
    r2invTEMP3 = delz * delz;
    t1 = r2invTEMP + r2invTEMP2;
    t2 = t1 + r2invTEMP3;
    r2inv = t2;
    int r6inv = r2inv * r2inv * r2inv;
    int force = r2inv * r6inv * (15 * r6inv - 20);
    fx += delx * force;
    fy += dely * force;
    fz += delz * force;
  }
  d_force_x[i] = fx;
  d_force_y[i] = fy;
  d_force_z[i] = fz;
}
```