Berk introduced this in
e82d62ed0. It won't have a purpose in the
Verlet scheme, because no charge groups exist. Even if it does still
have a purpose, tpxio.c is not the place to do the transformation.
Change-Id: I1a2d79d004f64b6fdb772f712ac16616516741c3
``GMX_NO_PULLVIR``
when set, do not add virial contribution to COM pull forces.
-``GMX_NOCHARGEGROUPS``
- disables multi-atom charge groups, i.e. each atom
- in all non-solvent molecules is assigned its own charge group.
-
``GMX_NOPREDICT``
shell positions are not predicted.
{
gmx_mtop_finalize(mtop);
}
-
- if (file_version >= 57)
- {
- char *env;
- int ienv;
- env = getenv("GMX_NOCHARGEGROUPS");
- if (env != NULL)
- {
- sscanf(env, "%d", &ienv);
- fprintf(stderr, "\nFound env.var. GMX_NOCHARGEGROUPS = %d\n",
- ienv);
- if (ienv > 0)
- {
- fprintf(stderr,
- "Will make single atomic charge groups in non-solvent%s\n",
- ienv > 1 ? " and solvent" : "");
- gmx_mtop_make_atomic_charge_groups(mtop, ienv == 1);
- }
- fprintf(stderr, "\n");
- }
- }
}
return ePBC;
return atoms;
}
-void gmx_mtop_make_atomic_charge_groups(gmx_mtop_t *mtop,
- gmx_bool bKeepSingleMolCG)
-{
- int mb, cg;
- t_block *cgs_mol;
-
- for (mb = 0; mb < mtop->nmolblock; mb++)
- {
- cgs_mol = &mtop->moltype[mtop->molblock[mb].type].cgs;
- if (!(bKeepSingleMolCG && cgs_mol->nr == 1))
- {
- cgs_mol->nr = mtop->molblock[mb].natoms_mol;
- cgs_mol->nalloc_index = cgs_mol->nr + 1;
- srenew(cgs_mol->index, cgs_mol->nalloc_index);
- for (cg = 0; cg < cgs_mol->nr+1; cg++)
- {
- cgs_mol->index[cg] = cg;
- }
- }
- }
-}
-
/*
* The cat routines below are old code from src/kernel/topcat.c
*/
gmx_mtop_global_atoms(const struct gmx_mtop_t *mtop);
-/* Make all charge groups the size of one atom.
- * When bKeepSingleMolCG==TRUE keep charge groups for molecules
- * that consist of a single charge group.
- */
-void
-gmx_mtop_make_atomic_charge_groups(struct gmx_mtop_t *mtop, gmx_bool bKeepSingleMolCG);
-
-
/* Generate a 'local' topology for the whole system.
* When ir!=NULL the free energy interactions will be sorted to the end.
*/