phase 1 seems to be done

diff --git a/src/spirals.cpp b/src/spirals.cpp
index d3767d4a202a5121985d43e320cd382b6710ea5b..8f8bdb3c99be7d356dfd4b81fcf9125057381e31 100644 (file)
--- a/src/spirals.cpp
+++ b/src/spirals.cpp
@@ -155,7 +155,7 @@ double fp3 (double x0, double y0, double z0, double p1, double p2, double p3, st
 void linear_kernel_search (double &x0, double &y0, double &z0, double &p1, double &p2, double &p3, std::vector< RVec > x, double epsi) {
     double F = 1, FX = 0, FX0 = 0, FY0 = 0, FZ0 = 0, FP1 = 0, FP2 = 0, FP3 = 0;
     double L1, L2, L3, L4, L5, L6;
-    int count = 0;
+    //int count = 0;
     while (std::abs(F - FX) >= epsi) {
         F = FX;
 
@@ -197,9 +197,9 @@ void linear_kernel_search (double &x0, double &y0, double &z0, double &p1, doubl
         p1 -= L4 * FP1;
         p2 -= L5 * FP2;
         p3 -= L6 * FP3;
-        count++;
+        //count++;
     }
-    std::cout << count << "\n";
+    //std::cout << count << "\n";
 }
 
 RVec kernel_pro (double x0, double y0, double z0, double p1, double p2, double p3, RVec x) {
@@ -338,7 +338,7 @@ Spirals::analyzeFrame(int frnr, const t_trxframe &fr, t_pbc *pbc,
     //linear_kernel_search(mid[0], mid[1], mid[2], arrow[0], arrow[1], arrow[2], temp, epsi); //изменить формат функции для изменения значений переменных
     linear_kernel_search(t1, t2, t3, t4, t5, t6, temp, epsi);
 
-    std::cout << t1 << " " << t2 << " " << t3 << " " << t4 << " " << t5 << " " << t6 << "\n";
+    //std::cout << t1 << " " << t2 << " " << t3 << " " << t4 << " " << t5 << " " << t6 << "\n";
 
     mid[0] = t1;
     mid[1] = t2;