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;
rvec atomH{};
float prevCODist {CalculateAtomicDistances(Donor.prevResi->getIndex(backboneAtomTypes::AtomC), Donor.prevResi->getIndex(backboneAtomTypes::AtomO), fr, pbc)};
for (int i{XX}; i <= ZZ; ++i){