for (size_t i = 0; i < inputIndex.size(); ++i) {
std::cout << top.atoms()->atom[inputIndex[i]].resind << std::endl;
aminoacidsIndex.resize(std::max(aminoacidsIndex.size(), static_cast< size_t >(top.atoms()->atom[inputIndex[i]].resind + 1)));
- aminoacidsIndex[top.atoms()->atom[inputIndex[i]].resind - 1].push_back(inputIndex[i]);
+ aminoacidsIndex[top.atoms()->atom[inputIndex[i]].resind].push_back(inputIndex[i]);
}
}