From: Max <Infinity2573@gmail.com>
Date: Sun, 25 Sep 2022 13:45:43 +0000 (+0300)
Subject: Нашёл в чём вся хуйня. Исправляю
X-Git-Url: http://biod.pnpi.spb.ru/gitweb/?a=commitdiff_plain;h=7defe9120cad776caa253f5d1a691b98bd44165b;p=alexxy%2Fgromacs-dssp.git

Нашёл в чём вся хуйня. Исправляю
---

diff --git a/src/dssptools.cpp b/src/dssptools.cpp
index 0dc84e8..8d63288 100644
--- a/src/dssptools.cpp
+++ b/src/dssptools.cpp
@@ -661,7 +661,8 @@ void DsspTool::calculateHBondEnergy(ResInfo& Donor,
                    float prevCODist {CalculateAtomicDistances(Donor.prevResi->getIndex(backboneAtomTypes::AtomC), Donor.prevResi->getIndex(backboneAtomTypes::AtomO), fr, pbc)};
                    for (int i{XX}; i <= ZZ; ++i){
                        float prevCO = fr.x[Donor.prevResi->getIndex(backboneAtomTypes::AtomC)][i] - fr.x[Donor.prevResi->getIndex(backboneAtomTypes::AtomO)][i];
-                       atomH[i] = prevCO / prevCODist;
+                       atomH[i] = fr.x[Donor.getIndex(backboneAtomTypes::AtomH)][i]; // Но на самом деле берутся координаты N
+                       atomH[i] += prevCO / prevCODist;
                    }
                    distanceHO = CalculateAtomicDistances(atomH, Acceptor.getIndex(backboneAtomTypes::AtomO), fr, pbc);
                    distanceHC = CalculateAtomicDistances(atomH, Acceptor.getIndex(backboneAtomTypes::AtomC), fr, pbc);