Продолжнаю выяснять причину расхождение через coutы
authorГорелов Сергей Васильевич <gorelov_sv@pnpi.nrcki.ru>
Mon, 3 Oct 2022 09:46:30 +0000 (12:46 +0300)
committerГорелов Сергей Васильевич <gorelov_sv@pnpi.nrcki.ru>
Mon, 3 Oct 2022 09:46:30 +0000 (12:46 +0300)
src/dssptools.cpp

index 77b7e32d4f1f4d2b95405ceacf48e62a7fa8e590..301624b4f80cf566833568378e06efdf0832a413 100644 (file)
@@ -766,6 +766,13 @@ float DsspTool::CalculateDihedralAngle(const int &A, const int &B, const int &C,
         y[i] = (vBD[j] * x[k]) - (vBD[k] * x[j]);
     }
 
+    std::cout << "v12 = " << vBA[0] << ", " << vBA[1] << ", " << vBA[2] << std::endl;
+    std::cout << "v43 = " << vDC[0] << ", " << vDC[1] << ", " << vDC[2] << std::endl;
+    std::cout << "z = " << vBD[0] << ", " << vBD[1] << ", " << vBD[2] << std::endl;
+    std::cout << "p = " << p[0] << ", " << p[1] << ", " << p[2] << std::endl;
+    std::cout << "x = " << x[0] << ", " << x[1] << ", " << x[2] << std::endl;
+    std::cout << "y = " << y[0] << ", " << y[1] << ", " << y[2] << std::endl;
+
     u = (x[XX] * x[XX]) + (x[YY] * x[YY]) + (x[ZZ] * x[ZZ]);
     v = (y[XX] * y[XX]) + (y[YY] * y[YY]) + (y[ZZ] * y[ZZ]);