gmx::AnalysisNeighborhoodSearch nbsearch = nbhood.initSearch(pbc, gmx::AnalysisNeighborhoodPositions(fr.x, fr.natoms));
gmx::AnalysisNeighborhoodPair pair;
for (const auto &i : inputColor) {
gmx::AnalysisNeighborhoodSearch nbsearch = nbhood.initSearch(pbc, gmx::AnalysisNeighborhoodPositions(fr.x, fr.natoms));
gmx::AnalysisNeighborhoodPair pair;
for (const auto &i : inputColor) {
gmx::AnalysisNeighborhoodPairSearch pairSearch = nbsearch.startPairSearch(inputFrame[i.second].as_vec());
while (pairSearch.findNextPair(&pair)) {
for (const auto &j : inputIndex) {
gmx::AnalysisNeighborhoodPairSearch pairSearch = nbsearch.startPairSearch(inputFrame[i.second].as_vec());
while (pairSearch.findNextPair(&pair)) {
for (const auto &j : inputIndex) {