std::vector<int> qmmmAtoms;
for (int i = 0; i < numQmmmGroups; i++)
{
- for (const AtomProxy &atomP : AtomRange(mtop))
+ for (const AtomProxy atomP : AtomRange(mtop))
{
int index = atomP.globalAtomNumber();
if (getGroupType(groups, egcQMMM, index) == i)