if (bMol)
gmx_rmpbc(gpbc,natoms,box,x[cur]);
+ /* calculate the molecules' centers of masses and put them into xa */
+ if (bMol)
+ calc_mol_com(gnx[0],index[0],&top->mols,&top->atoms, x[cur],xa[cur]);
+
/* first remove the periodic boundary condition crossings */
for(i=0;i<curr->ngrp;i++)
{
&top->atoms, com);
}
- /* calculate the molecules' centers of masses and put them into xa */
- if (bMol)
- calc_mol_com(gnx[0],index[0],&top->mols,&top->atoms, x[cur],xa[cur]);
-
/* loop over all groups in index file */
for(i=0; (i<curr->ngrp); i++)
{