- debug
authorAnatoly <Titov_AI@pnpi.nrcki.ru>
Wed, 21 Oct 2020 15:32:35 +0000 (18:32 +0300)
committerAnatoly <Titov_AI@pnpi.nrcki.ru>
Wed, 21 Oct 2020 15:32:35 +0000 (18:32 +0300)
src/colorvec.cpp

index afda82bd66c60938ee4b046d084c368f5810a913..8d60cd52937261d70b12712496bb6d8641ad2979 100644 (file)
@@ -176,11 +176,6 @@ inline void betaListsRVecsEvaluation(const std::vector< gmx::RVec > &frame, cons
     temp.resize(0);
     gmx::RVec tempA;
     for (const auto &i : inputBetaLists) {
-        std::cout << i.size() << std::endl;
-        std::cout << frame[inputCA[i[i.size() - 1]]][0] << std::endl;
-        std::cout << frame[inputCA[i[i.size() - 2]]][0] << std::endl;
-        std::cout << frame[inputCA[i[0]]][0] << std::endl;
-        std::cout << frame[inputCA[i[1]]][0] << std::endl;
         tempA = frame[inputCA[i[i.size() - 1]]] + frame[inputCA[i[i.size() - 2]]] - frame[inputCA[i[1]]] - frame[inputCA[i[0]]];
         tempA /= tempA.norm();
         temp.push_back(tempA);
@@ -495,6 +490,10 @@ colorVec::analyzeFrame(int frnr, const t_trxframe &fr, t_pbc *pbc, gmx::Trajecto
     std::cout << "\n\ntest001\n";
     betaListsRVecsEvaluation(trajectoryFrame, betaLists[frnr], betaListsRVecs, indexCA);
     std::cout << "\n\ntest002\n";
+
+    std::cout << aminoacidsIndex.size();
+
+
     for (size_t i = 0; i < colorsToPeptide.size(); ++i) {
         if (colorsToPeptide[i]) {
             searchNearBetaLists(trajectoryFrame, betaLists[frnr], colorsNames[i], colorsToBeta[i], aminoacidsIndex, effRad * 9000);