[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