added some usefull comments
authorAnatoly Titov <toluk@omrb.pnpi.spb.ru>
Fri, 8 Sep 2017 09:46:51 +0000 (12:46 +0300)
committerAnatoly Titov <toluk@omrb.pnpi.spb.ru>
Fri, 8 Sep 2017 09:46:51 +0000 (12:46 +0300)
src/spacetimecorr.cpp

index f754ce32e459f2d56f4e676401e293469f201e58..b50b74207bff0a23723b82cf12ecb1fd1467ee92 100644 (file)
@@ -232,7 +232,7 @@ void
 Domains::initAfterFirstFrame(const TrajectoryAnalysisSettings       &settings,
                              const t_trxframe                       &fr)
 {
-    std::cout << "Init After First Frame\n\n";
+    std::cout << "Init After First Frame - start\n\n";
     ConstArrayRef< int > atomind  = sel_[0].atomIndices();
     index.resize(0);
     for (ConstArrayRef<int>::iterator ai = atomind.begin(); (ai < atomind.end()); ai++) {
@@ -314,7 +314,7 @@ Domains::initAfterFirstFrame(const TrajectoryAnalysisSettings       &settings,
     frankenstein_trajectory.back() = trajectory.back();
 
     frames++;
-    std::cout << "Init Analysis Finish\n\n";
+    std::cout << "Init After First Frame - end\n\n";
 }
 
 void
@@ -381,7 +381,7 @@ Domains::analyzeFrame(int frnr, const t_trxframe &fr, t_pbc *pbc,
 void
 Domains::finishAnalysis(int nframes)
 {
-    std::cout << "Finish Analysis\n\n";
+    std::cout << "Finish Analysis - start\n\n";
 
     //
     //make_pdb_trajectory(frankenstein_trajectory, "/home/toluk/Develop/samples/reca_rd/10k_frames.pdb");
@@ -400,13 +400,11 @@ Domains::finishAnalysis(int nframes)
         }
     }
 
-    std::cout << "\n" ;
+    std::cout << "Correlation's evaluation - start\n\n";
     #pragma omp parallel
     {
         #pragma omp for schedule(dynamic)
         for (int i = m; i <= k; i += 1) {
-            std::cout << " k = " << i << " ";
-
             rvec *medx, *medy, temp_zero;
             snew(medx, index.size());
             snew(medy, index.size());
@@ -476,16 +474,13 @@ Domains::finishAnalysis(int nframes)
             sfree(medy);
         }
     }
-    std::cout << "\n" ;
-
-    std::cout << "printing correlation's file\n" ;
+    std::cout << "Correlation's evaluation - end\n\n";
 
     //
     //make_correlation_file(crltns, "corrs_vectors.txt", m);
     //
 
-    std::cout << "done\n" ;
-
+    std::cout << "Correlation's overview - start\n\n";
     int imax = 0;
     for (int i1 = 0; i1 < 100; i1++) {
         int i5 = 0;
@@ -522,12 +517,14 @@ Domains::finishAnalysis(int nframes)
     for (std::set< std::pair< int, int > >::iterator it = crr_prs.begin(); it != crr_prs.end(); it++) {
         corr_prs.push_back(*it);
     }
+    std::cout << "Correlation's overview - end\n\n";
     //make_best_corrs_graphics(crltns, corr_prs, index, "best_graphics.txt");
 
     std::vector< std::vector< int > > graph;
     graph.resize(index.size());
     rvec temp;
 
+    std::cout << "Path's finding - start\n\n";
     for (int i = 1; i <= k; i++) {
         for (int j = 0; j < index.size(); j++) {
             for (int f = 0; f < index.size(); f++) {
@@ -629,13 +626,13 @@ Domains::finishAnalysis(int nframes)
             rout_new.push_back(rout_old[i]);
         }
     }
-    std::cout << "\nfinished routs\n\n" ;
+    std::cout << "Path's finding - end\n\n";
 
     //
     //make_rout_file(crl_border, index, rout_new, "routs.txt");
     //
 
-    std::cout << "\n\n\n job done \n\n\n";
+    std::cout << "Finish Analysis - end\n\n";
 }
 
 void