const TopologyInformation & /*top*/)
{
std::vector< int > a;
- a.resize(0);
index_all.resize(0);
index.resize(0);
for (int i = 0; i < sel_.size(); i++) {
- index_all.push_back(a);
+ a.resize(0);
ConstArrayRef< int > atomind = sel_[i].atomIndices();
for (ConstArrayRef< int >::iterator ai = atomind.begin(); (ai < atomind.end()); ai++) {
- index_all.back().push_back(*ai);
+ a.push_back(*ai);
index.push_back(*ai);
}
+ index_all.push_back(a);
}
std::cout << "First part finished";
}