std::cout << "Testing Helix_" << static_cast<std::size_t>(i) + 3 << std::endl;
std::size_t stride {static_cast<std::size_t>(i) + 3};
for(std::size_t j {0}; j + static_cast<std::size_t>(i) < SecondaryStructuresStatusMap.size(); ++j){
std::cout << "Testing Helix_" << static_cast<std::size_t>(i) + 3 << std::endl;
std::size_t stride {static_cast<std::size_t>(i) + 3};
for(std::size_t j {0}; j + static_cast<std::size_t>(i) < SecondaryStructuresStatusMap.size(); ++j){
if ( hasHBondBetween(j, j + (static_cast<std::size_t>(i))) && NoChainBreaksBetween(j, j + stride) ){
if ( hasHBondBetween(j, j + (static_cast<std::size_t>(i))) && NoChainBreaksBetween(j, j + stride) ){
SecondaryStructuresStatusMap[j + static_cast<std::size_t>(i)].setStatus(HelixPositions::End, i);
for (std::size_t k {1}; k < (static_cast<std::size_t>(i)); ++k){
SecondaryStructuresStatusMap[j + static_cast<std::size_t>(i)].setStatus(HelixPositions::End, i);
for (std::size_t k {1}; k < (static_cast<std::size_t>(i)); ++k){
- std::cout << j << " is HELIX" << std::endl;
- if ( empty ){
- for(std::size_t k {0}; k < (static_cast<std::size_t>(i)); ++k ){
+ if ( empty || Helix == secondaryStructureTypes::Helix_4 ){
+ for(std::size_t k {0}; k < stride ; ++k ){