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){
}
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);
}
}
}
}
if (isTurn){
+ std::cout << "Resi " << i << " is Turn" << std::endl;
SecondaryStructuresStatusMap[i].setStatus(secondaryStructureTypes::Turn);
}
}