od->nref = 0;
for (int i = 0; i < mtop->natoms; i++)
{
- if (getGroupType(mtop->groups, egcORFIT, i) == 0)
+ if (getGroupType(mtop->groups, SimulationAtomGroupType::OrientationRestraintsFit, i) == 0)
{
od->nref++;
}
{
const t_atom &local = atomP.atom();
int i = atomP.globalAtomNumber();
- if (mtop->groups.grpnr[egcORFIT] == nullptr ||
- mtop->groups.grpnr[egcORFIT][i] == 0)
+ if (mtop->groups.groupNumbers[SimulationAtomGroupType::OrientationRestraintsFit].empty() ||
+ mtop->groups.groupNumbers[SimulationAtomGroupType::OrientationRestraintsFit][i] == 0)
{
/* Not correct for free-energy with changing masses */
od->mref[j] = local.m;