Вернул, как было master
authorMax <Infinity2573@gmail.com>
Tue, 11 Oct 2022 21:24:58 +0000 (00:24 +0300)
committerMax <Infinity2573@gmail.com>
Tue, 11 Oct 2022 21:24:58 +0000 (00:24 +0300)
src/dssptools.cpp

index 93063ea75d70310faf6ec0062692ffc146b0d9f0..55156463533ab32a609c2ff829d567774a975047 100644 (file)
@@ -964,12 +964,7 @@ void DsspTool::calculateHBondEnergy(ResInfo& Donor,
                 * ((1 / distanceNO) + (1 / distanceHC) - (1 / distanceHO) - (1 / distanceNC));
         }
 
-               std::cout << "CA-CA distance: " << CalculateAtomicDistances(
-                                Donor.getIndex(backboneAtomTypes::AtomCA), Acceptor.getIndex(backboneAtomTypes::AtomCA), fr, pbc) << std::endl;
-               std::cout << "N-O distance: " << distanceNO << std::endl;
-               std::cout << "N-C distance: " << distanceNC << std::endl;
-               std::cout << "H-O distance: " << distanceHO << std::endl;
-               std::cout << "H-C distance: " << distanceHC << std::endl;
+
 
         HbondEnergy = std::round(HbondEnergy * 1000) / 1000;
 
@@ -977,10 +972,18 @@ void DsspTool::calculateHBondEnergy(ResInfo& Donor,
             HbondEnergy = minEnergy;
         }
 
-        std::cout << "Calculated energy = " << HbondEnergy << std::endl;
+
 
     } // TODO
 
+    std::cout << "CA-CA distance: " << CalculateAtomicDistances(
+        Donor.getIndex(backboneAtomTypes::AtomCA), Acceptor.getIndex(backboneAtomTypes::AtomCA), fr, pbc) << std::endl;
+    std::cout << "N-O distance: " << distanceNO << std::endl;
+    std::cout << "N-C distance: " << distanceNC << std::endl;
+    std::cout << "H-O distance: " << distanceHO << std::endl;
+    std::cout << "H-C distance: " << distanceHC << std::endl;
+    std::cout << "Calculated energy = " << HbondEnergy << std::endl;
+
 
 
 //    else{
@@ -1017,7 +1020,7 @@ float DsspTool::CalculateAtomicDistances(const int &A, const int &B, const t_trx
 {
    gmx::RVec r{ 0, 0, 0 };
    pbc_dx(pbc, fr.x[A], fr.x[B], r.as_vec());
-   return r.norm() * gmx::c_nm2A * gmx::c_nm2A; // НЕ ТРОГАТЬ
+   return r.norm() * gmx::c_nm2A; // НЕ ТРОГАТЬ
 }
 
 /* Calculate Distance From B to A, where A is only fake coordinates */
@@ -1025,7 +1028,7 @@ float DsspTool::CalculateAtomicDistances(const rvec &A, const int &B, const t_tr
 {
    gmx::RVec r{ 0, 0, 0 };
    pbc_dx(pbc, A, fr.x[B], r.as_vec());
-   return r.norm() * gmx::c_nm2A * gmx::c_nm2A; // НЕ ТРОГАТЬ
+   return r.norm() * gmx::c_nm2A; // НЕ ТРОГАТЬ
 }
 
 void DsspTool::initAnalysis(/*const TrajectoryAnalysisSettings &settings,*/const TopologyInformation& top, const initParameters &initParamz)