bool isNearPeptide(const std::vector< gmx::RVec > &inputFrame, const std::vector< size_t > &inputIndex,
const std::vector< std::pair< std::string, size_t > > &inputColor, const double cutOff) {
for (size_t i {0}; i < inputColor.size(); ++i) {
- while (pairSearch.findNextPair(&pair)) {
- for (size_t j {0}; j < inputIndex.size(); ++j) {
- if ((inputFrame[inputIndex[j]] - inputFrame[inputColor[i].second]).norm() <= cutOff) {
- return true;
- }
+ for (size_t j {0}; j < inputIndex.size(); ++j) {
+ if ((inputFrame[inputIndex[j]] - inputFrame[inputColor[i].second]).norm() <= cutOff) {
+ return true;
}
}
}
colorsToPeptide.resize(colorsNames.size(), false);
std::cout << "\t\tWhich colors are \"near\" the pepride mass." << std::endl;
for (size_t i = 0; i < colorsNames.size(); ++i) {
- colorsToPeptide[i] = isNearPeptide(fr, pbc, trajectoryFrame, /*nb_, */index, colorsNames[i]);
+ colorsToPeptide[i] = isNearPeptide(trajectoryFrame, index, colorsNames[i], effRad);
}
// расчёт угла и среднего угла с ближайшими бета листами
std::vector< std::vector< bool > > colorsToBeta;