aminoacidsIndex.resize(0);
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)));
+ 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]);
}
}