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 " << j << " and " << j + stride << 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 " << j << " and " << j + stride << std::endl;
if( SecondaryStructuresStatusMap[j + k].getStatus(i) == HelixPositions::None ){
SecondaryStructuresStatusMap[j + k].setStatus(HelixPositions::Middle, i);
}
if( SecondaryStructuresStatusMap[j + k].getStatus(i) == HelixPositions::None ){
SecondaryStructuresStatusMap[j + k].setStatus(HelixPositions::Middle, i);
}
{
gmx::RVec r{ 0, 0, 0 };
pbc_dx(pbc, fr.x[A], fr.x[B], r.as_vec());
{
gmx::RVec r{ 0, 0, 0 };
pbc_dx(pbc, fr.x[A], fr.x[B], r.as_vec());
{
gmx::RVec r{ 0, 0, 0 };
pbc_dx(pbc, A, fr.x[B], r.as_vec());
{
gmx::RVec r{ 0, 0, 0 };
pbc_dx(pbc, A, fr.x[B], r.as_vec());