looking for segmentation error

diff --git a/src/spacetimecorr.cpp b/src/spacetimecorr.cpp
index 96997915e3d85bb762ef737db237fa201cb5a8db..e58b8fda9748070070292b71be56d55be8233d33 100644 (file)
--- a/src/spacetimecorr.cpp
+++ b/src/spacetimecorr.cpp
@@ -110,7 +110,7 @@ void make_best_corrs_graphics(std::vector< std::vector< std::vector< double > >
     file = std::fopen(file_name, "w+");
     for (int i = 0; i < corr_pairs.size(); i++) {
         std::fprintf(file, "%3d %3d", indx[corr_pairs[i].first] + 1, indx[corr_pairs[i].second] + 1);
-        for (int j = 0; j < correlations.size(); j++) {
+        for (int j = 1; j < correlations.size(); j++) {
             std::fprintf(file, "%3.2f ", correlations[j][corr_pairs[i].first][corr_pairs[i].second]);
         }
         std::fprintf(file, "\n");
@@ -225,12 +225,14 @@ Domains::initAnalysis(const TrajectoryAnalysisSettings &settings,
                       const TopologyInformation        &top)
 {
     domains_ePBC = top.ePBC();
+    std::cout << "Init Analysis\n\n";
 }
 
 void
 Domains::initAfterFirstFrame(const TrajectoryAnalysisSettings       &settings,
                              const t_trxframe                       &fr)
 {
+    std::cout << "Init After First Frame\n\n";
     ConstArrayRef< int > atomind  = sel_[0].atomIndices();
     index.resize(0);
     for (ConstArrayRef<int>::iterator ai = atomind.begin(); (ai < atomind.end()); ai++) {
@@ -312,12 +314,14 @@ Domains::initAfterFirstFrame(const TrajectoryAnalysisSettings       &settings,
     frankenstein_trajectory.back() = trajectory.back();
 
     frames++;
+    std::cout << "Init Analysis Finish\n\n";
 }
 
 void
 Domains::analyzeFrame(int frnr, const t_trxframe &fr, t_pbc *pbc,
                       TrajectoryAnalysisModuleData *pdata)
 {
+    std::cout << "Analyze Frame\n\n";
     trajectory.resize(trajectory.size() + 1);
     trajectory.back().resize(index.size());
 
@@ -373,11 +377,13 @@ Domains::analyzeFrame(int frnr, const t_trxframe &fr, t_pbc *pbc,
     }
 
     frames++;
+    std::cout << "Finishe Analyze Frame\n\n";
 }
 
 void
 Domains::finishAnalysis(int nframes)
 {
+    std::cout << "Finish Analysis\n\n";
     make_pdb_trajectory(frankenstein_trajectory, "/home/toluk/Develop/samples/reca_rd/10k_frames.pdb");
 
     std::vector< std::vector< std::vector< double > > > crltns;