for (std::size_t k {1}; k < stride; ++k){
if( SecondaryStructuresStatusMap[j + k].getStatus(i) == HelixPositions::None ){
SecondaryStructuresStatusMap[j + k].setStatus(HelixPositions::Middle, i);
- SecondaryStructuresStatusMap[j + k].setStatus(secondaryStructureTypes::Turn);
+// SecondaryStructuresStatusMap[j + k].setStatus(secondaryStructureTypes::Turn);
}
}
}
}
-// for(std::size_t i {1}; i + 1 < SecondaryStructuresStatusMap.size(); ++i){
-// if (SecondaryStructuresStatusMap[i].getStatus() == secondaryStructureTypes::Loop || SecondaryStructuresStatusMap[i].getStatus(secondaryStructureTypes::Bend)){
-// bool isTurn = false;
-// for(const turnsTypes &j : {turnsTypes::Turn_3, turnsTypes::Turn_4, turnsTypes::Turn_5}){
-// std::size_t stride {static_cast<std::size_t>(i) + 3};
-// for(std::size_t k {1}; k < stride; ++k){
-// isTurn = (i >= k) && (SecondaryStructuresStatusMap[i - k].getStatus(j) == HelixPositions::Start || SecondaryStructuresStatusMap[i - k].getStatus(j) == HelixPositions::Start_AND_End) ;
-// }
-// }
+ for(std::size_t i {1}; i + 1 < SecondaryStructuresStatusMap.size(); ++i){
+ if (SecondaryStructuresStatusMap[i].getStatus() == secondaryStructureTypes::Loop || SecondaryStructuresStatusMap[i].getStatus(secondaryStructureTypes::Bend)){
+ bool isTurn = false;
+ for(const turnsTypes &j : {turnsTypes::Turn_3, turnsTypes::Turn_4, turnsTypes::Turn_5}){
+ std::size_t stride {static_cast<std::size_t>(j) + 3};
+ for(std::size_t k {1}; k < stride; ++k){
+ isTurn = (i >= k) && (SecondaryStructuresStatusMap[i - k].getStatus(j) == HelixPositions::Start || SecondaryStructuresStatusMap[i - k].getStatus(j) == HelixPositions::Start_AND_End) ;
+ }
+ }
-// if (isTurn){
-// SecondaryStructuresStatusMap[i].setStatus(secondaryStructureTypes::Turn);
-// }
-// }
-// }
+ if (isTurn){
+ SecondaryStructuresStatusMap[i].setStatus(secondaryStructureTypes::Turn);
+ }
+ }
+ }
}