[gmx-developers] remove periodicity
Gaurav Goel
gauravgoeluta at gmail.com
Thu Jul 1 22:22:41 CEST 2010
If I know the coordinates at time 't-del_t' and 't', I can remove the
periodicity as follows:
/* Remove periodicity */
...
...
while(xcur[ind][m]-xprev[ind][m] <= -hbox[m])
rvec_inc(xcur[ind],box[m]);
while(xcur[ind][m]-xprev[ind][m] > hbox[m])
rvec_dec(xcur[ind],box[m]);
...
...
-> If I only know the coordinates at 't', is it still possible to
remove periodicity?
In the update() subroutine, I want to be able to make a list of all
atoms within a thin slab at a fixed distance from the edge of the box.
E.g., all particles in a slab of width L/10, with the x-coordinate of
the center of slab at a distance L/4 from the yz-face of box.
thanks,
Gaurav
More information about the gromacs.org_gmx-developers
mailing list