ПП-хеликсы проверены. Проверяю теперь остальные хеликсы
authorГорелов Сергей Васильевич <gorelov_sv@pnpi.nrcki.ru>
Mon, 3 Oct 2022 10:44:55 +0000 (13:44 +0300)
committerГорелов Сергей Васильевич <gorelov_sv@pnpi.nrcki.ru>
Mon, 3 Oct 2022 10:44:55 +0000 (13:44 +0300)
src/dssptools.cpp

index 2e261eeceaed35dcc621790cab180d22af70c8c6..3520c6a8f319b5352d7380b1f9e064ad8f1cf8bd 100644 (file)
@@ -397,6 +397,7 @@ void secondaryStructures::analyzeTurnsAndHelicesPatterns(){
         std::size_t stride {static_cast<std::size_t>(i) + 3};
         for(std::size_t j {0}; j + stride < SecondaryStructuresStatusMap.size(); ++j){
             if ( hasHBondBetween(j + stride, j) && NoChainBreaksBetween(j, j + stride) ){
+                std::cout << "Resi " << j << " is Helix_" << stride << "start" << std::endl;
                 SecondaryStructuresStatusMap[j + stride].setStatus(HelixPositions::End, i);
 
                 for (std::size_t k {1}; k < stride; ++k){
@@ -442,6 +443,7 @@ void secondaryStructures::analyzeTurnsAndHelicesPatterns(){
                 }
                 if ( empty || Helix == secondaryStructureTypes::Helix_4 ){
                     for(std::size_t k {0}; k < stride; ++k ){
+                        std::cout << "Resi " << j << " is Helix_" << static_cast<std::size_t>(Helix) << std::endl;
                         SecondaryStructuresStatusMap[j + k].setStatus(Helix);
                     }
                 }
@@ -461,6 +463,7 @@ void secondaryStructures::analyzeTurnsAndHelicesPatterns(){
                 }
             }
             if (isTurn){
+                std::cout << "Resi " << i << " is Turn" << std::endl;
                 SecondaryStructuresStatusMap[i].setStatus(secondaryStructureTypes::Turn);
             }
         }