From: Max Date: Tue, 27 Sep 2022 12:32:35 +0000 (+0300) Subject: Добавил сиауты, чтобы найти возвможные ошибки. X-Git-Url: http://biod.pnpi.spb.ru/gitweb/?a=commitdiff_plain;h=f7233d2742e531552a2348621d7100bc039b7b79;p=alexxy%2Fgromacs-dssp.git Добавил сиауты, чтобы найти возвможные ошибки. --- diff --git a/src/dssptools.cpp b/src/dssptools.cpp index ac28492..f900da9 100644 --- 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(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); /* Пропустил проверку того, что заменяемая ак - петля */