[gmx-developers] double sum over neighbor list

[mashayak at votca.org]

Hi

 

I am tyring to compute total virial, i.e. sum_i sum_j (rij*Fij). To achieve this, I am making modifications in the gmx_nb_generic_kernel() function in nb_generic.c. In it, I see that atoms I and j are accessed from neighbor lists to determine forces on atom I . My question is that in this loop do I have to take care of double counting pair forces, i.e. will the pair of atom I and j be accessed twice or just once?

 

Thanks

Sikandar
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-developers/attachments/20130531/d27a19ac/attachment.html>