From 654a713ebbaaff8958a026e8fd83bd180b5b611d Mon Sep 17 00:00:00 2001 From: Max Date: Tue, 27 Sep 2022 16:03:02 +0300 Subject: [PATCH] =?utf8?q?=D0=92=D1=80=D0=BE=D0=B4=D0=B5=20=D0=BA=D0=B0?= =?utf8?q?=D0=BA=20=D1=80=D0=B0=D0=B1=D0=BE=D1=82=D0=B0=D0=B5=D1=82.=20?= =?utf8?q?=D0=A2=D0=B5=D0=BF=D0=B5=D1=80=D1=8C=20=D0=BD=D0=B0=D0=B4=D0=BE?= =?utf8?q?=20=D0=B1=D1=83=D0=B4=D0=B5=D1=82=20=D0=BD=D0=B0=D1=81=D1=82?= =?utf8?q?=D1=80=D0=B0=D0=B8=D0=B2=D0=B0=D1=82=D1=8C=20=D0=B5=D0=B3=D0=BE?= =?utf8?q?=20=D1=80=D0=B0=D0=B1=D0=BE=D1=82=D1=83.?= MIME-Version: 1.0 Content-Type: text/plain; charset=utf8 Content-Transfer-Encoding: 8bit --- src/dssptools.cpp | 12 +++++------- 1 file changed, 5 insertions(+), 7 deletions(-) diff --git a/src/dssptools.cpp b/src/dssptools.cpp index f900da9..62f8695 100644 --- a/src/dssptools.cpp +++ b/src/dssptools.cpp @@ -781,10 +781,10 @@ float DsspTool::CalculateDihedralAngle(const int &A, const int &B, const int &C, void DsspTool::calculateDihedrals(const t_trxframe &fr, const t_pbc *pbc){ const float epsilon = 29; - const float phi_min = -75 - epsilon; - const float phi_max = -75 + epsilon; - const float psi_min = 145 - epsilon; - const float psi_max = 145 + epsilon; + const float phi_min = -75 - epsilon; // -104 + const float phi_max = -75 + epsilon; // -46 + const float psi_min = 145 - epsilon; // 116 + const float psi_max = 145 + epsilon; // 176 std::vector phi, psi; phi.resize(0); psi.resize(0); @@ -804,7 +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; +// std::cout << "For " << i << " phi = " << phi[i] << ", psi = " << psi[i] << std::endl; } for (std::size_t i = 1; i + 3 < IndexMap.size(); ++i){ @@ -865,8 +865,6 @@ 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); /* Пропустил проверку того, что заменяемая ак - петля */ -- 2.22.0