Добавил сиауты, чтобы найти возвможные ошибки.

diff --git a/src/dssptools.cpp b/src/dssptools.cpp
index ac28492238d6905e5f8ff1c68bd6ab268b20ca3e..f900da930a1a563679e7d8390e53f7928137b817 100644 (file)
--- a/src/dssptools.cpp
+++ b/src/dssptools.cpp
@@ -804,6 +804,7 @@ void DsspTool::calculateDihedrals(const t_trxframe &fr, const t_pbc *pbc){
                                         static_cast<int>(IndexMap[i + 1].getIndex(backboneAtomTypes::AtomN)),
                                         fr,
                                         pbc);
+        std::cout << "For " << i << " phi = " << phi[i] << ", psi = " << psi[i] << std::endl;
     }
 
     for (std::size_t i = 1; i + 3 < IndexMap.size(); ++i){
@@ -836,7 +837,7 @@ void DsspTool::calculateDihedrals(const t_trxframe &fr, const t_pbc *pbc){
             /* Пропустил проверку того, что заменяемая ак - петля */
             PatternSearch.SecondaryStructuresStatusMap[i].setStatus(secondaryStructureTypes::Helix_PP);
             PatternSearch.SecondaryStructuresStatusMap[i + 1].setStatus(secondaryStructureTypes::Helix_PP);
-        break;
+            break;
         }
         case 3:{
             if (phi_min > phi[i] or phi[i] > phi_max or
@@ -864,6 +865,8 @@ void DsspTool::calculateDihedrals(const t_trxframe &fr, const t_pbc *pbc){
                     break;
             }
 
+            std::cout << "ПОЛИПРОЛИИИН" << std::endl;
+
             PatternSearch.SecondaryStructuresStatusMap[i + 1].setStatus(HelixPositions::Middle, turnsTypes::Turn_PP);
             PatternSearch.SecondaryStructuresStatusMap[i + 2].setStatus(HelixPositions::End, turnsTypes::Turn_PP);
             /* Пропустил проверку того, что заменяемая ак - петля */