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);
}
}
}
}
- for(std::size_t i {1}; i + 1 < SecondaryStructuresStatusMap.size(); ++i){
- if (SecondaryStructuresStatusMap[i].getStatus(secondaryStructureTypes::Loop)){
- 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)){
+// 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) ;
+// }
+// }
- if (isTurn){
- SecondaryStructuresStatusMap[i].setStatus(secondaryStructureTypes::Turn);
- }
- else if (SecondaryStructuresStatusMap[i].getStatus(secondaryStructureTypes::Bend)){
- SecondaryStructuresStatusMap[i].setStatus(secondaryStructureTypes::Bend);
- }
- }
- }
+// if (isTurn){
+// SecondaryStructuresStatusMap[i].setStatus(secondaryStructureTypes::Turn);
+// }
+// else if (SecondaryStructuresStatusMap[i].getStatus(secondaryStructureTypes::Bend)){
+// SecondaryStructuresStatusMap[i].setStatus(secondaryStructureTypes::Bend);
+// }
+// }
+// }
}