std::vector<int> genQmmmIndices(const gmx_mtop_t &mtop)
{
- std::vector<int> output;
- int global_at = 0;
- unsigned char *grpnr = mtop.groups.grpnr[egcQMMM];
+ std::vector<int> output;
+ int global_at = 0;
+ const unsigned char *grpnr = mtop.groups.groupNumbers[SimulationAtomGroupType::QuantumMechanics].data();
for (const gmx_molblock_t &molb : mtop.molblock)
{
for (int mol = 0; mol < molb.nmol; ++mol)