gdfgfdgfgd

diff --git a/src/dssptools.cpp b/src/dssptools.cpp
index 7c1d7c21366db52e2a22574199dc2048383db11c..1f09bfcda8b4f58e5a19ef7d1faa35b3236a1c66 100644 (file)
--- a/src/dssptools.cpp
+++ b/src/dssptools.cpp
@@ -284,7 +284,7 @@ void secondaryStructures::analyzeTurnsAndHelicesPatterns(){
         std::cout << "Testing Helix_" << stride << std::endl;
         for(std::size_t j {0}; j + stride < SecondaryStructuresStatusMap.size(); ++j){
             std::cout << "Testing " << j << " and " << j + stride << std::endl;
-            if ( hasHBondBetween(j, j + stride) && NoChainBreaksBetween(j, j + stride) ){
+            if ( hasHBondBetween(j + stride, j) && NoChainBreaksBetween(j, j + stride) ){
                 std::cout << j << " and " << j + stride << " has hbond!" << std::endl;
                 SecondaryStructuresStatusMap[j + stride].setStatus(HelixPositions::End, i);
 
@@ -713,9 +713,9 @@ void DsspTool::calculateHBondEnergy(ResInfo& Donor,
 
        HbondEnergy = std::round(HbondEnergy * 1000) / 1000;
 
-       if ( HbondEnergy < minEnergy ){
-            HbondEnergy = minEnergy;
-       }
+//       if ( HbondEnergy < minEnergy ){ // I don't think that this is correct
+//            HbondEnergy = minEnergy;
+//       }
 
 //       std::cout << "Calculated energy = " << HbondEnergy << std::endl;
     }