asdsad

diff --git a/src/dssptools.cpp b/src/dssptools.cpp
index 6dcadbf417bf36ae1e3971f1561c55b6a2bb6a68..86fde8f691d1b63ac67a1cf763647701037ac040 100644 (file)
--- a/src/dssptools.cpp
+++ b/src/dssptools.cpp
@@ -615,10 +615,10 @@ void DsspTool::calculateHBondEnergy(ResInfo& Donor,
 
            if (initParams.addHydrogens){
                if (Donor.prevResi != nullptr && Donor.prevResi->getIndex(backboneAtomTypes::AtomC) && Donor.prevResi->getIndex(backboneAtomTypes::AtomO)){
-                   std::cout << "On donor " << Donor.info->nr << *(Donor.info->name) << std::endl;
-                   std::cout <<"Prev donor is " << Donor.prevResi->info->nr << *(Donor.prevResi->info->name) << std::endl;
-                   std::cout <<"Prev C index is " << Donor.prevResi->getIndex(backboneAtomTypes::AtomC) << std::endl;
-                   std::cout <<"Prev O index is " << Donor.prevResi->getIndex(backboneAtomTypes::AtomO) << std::endl;
+//                   std::cout << "On donor " << Donor.info->nr << *(Donor.info->name) << std::endl;
+//                   std::cout <<"Prev donor is " << Donor.prevResi->info->nr << *(Donor.prevResi->info->name) << std::endl;
+//                   std::cout <<"Prev C index is " << Donor.prevResi->getIndex(backboneAtomTypes::AtomC) << std::endl;
+//                   std::cout <<"Prev O index is " << Donor.prevResi->getIndex(backboneAtomTypes::AtomO) << std::endl;
                    rvec atomH{};
                    float prevCODist {CalculateAtomicDistances(Donor.prevResi->getIndex(backboneAtomTypes::AtomC), Donor.prevResi->getIndex(backboneAtomTypes::AtomO), fr, pbc)};
                    for (int i{XX}; i <= ZZ; ++i){
@@ -747,6 +747,10 @@ void DsspTool::initAnalysis(/*const TrajectoryAnalysisSettings &settings,*/const
            IndexMap[i - 1].nextResi = &IndexMap[i];
 
            IndexMap[i].prevResi = &IndexMap[i - 1];
+
+           std::cout << "Resi " << IndexMap[i].info->nr << *(IndexMap[i].info->name) << std::endl;
+           std::cout <<"Prev resi is " << IndexMap[i-1].prevResi->info->nr << *(IndexMap[i-1].prevResi->info->name) << std::endl;
+
        }
        std::string atomname(*(top.atoms()->atomname[static_cast<std::size_t>(*ai)]));
        if (atomname == backboneAtomTypeNames[backboneAtomTypes::AtomCA])