changed window output

diff --git a/src/spirals.cpp b/src/spirals.cpp
index 4e603a35a32d122d8cdc87b17d39e1a99c98e5b1..03177ad78b50249652bf6888f73f70a076b47cfc 100644 (file)
--- a/src/spirals.cpp
+++ b/src/spirals.cpp
@@ -498,25 +498,25 @@ Spirals::finishAnalysis(int /*nframes*/)
 
     file = std::fopen("LocalSteps_Rspiral_Nmonomers.txt", "w+");
     for (int m = 0; m < window_circles.front().size(); m++) {
-        for (int i = 0; i < window_circles[m].size(); i++) {
-            std::fprintf(file, "Frame # %6.2d | window # %3.2d\n", i, m);
+        for (int i = 0; i < window_circles.size(); i++) {
+            std::fprintf(file, "Frame # %6.2d | window # %3.2d\n", m, i);
             std::fprintf(file, "Spiral steps:\n");
-            for (int j = 0; j < window_circles[m][i].size(); j++) {
-                std::fprintf(file, "%3.2f ", std::sqrt(distance2(window_kernel[m][i].krnl[window_circles[m][i][j].front() - 1], window_kernel[m][i].krnl[window_circles[m][i][j].back() - 1])));
+            for (int j = 0; j < window_circles[i][m].size(); j++) {
+                std::fprintf(file, "%3.2f ", std::sqrt(distance2(kernel[m].krnl[window_circles[i][m][j].front() - 1], kernel[m].krnl[window_circles[i][m][j].back() - 1])));
             }
             std::fprintf(file, "\n");
             std::fprintf(file, "Spiral radii\n");
-            for (int j = 0; j < window_circles[m][i].size(); j++) {
+            for (int j = 0; j < window_circles[i][m].size(); j++) {
                 float temp = 0;
-                for (int k = 0; k < window_circles[m][i][j].size(); k++) {
-                    temp += std::sqrt(distance2(window_kernel[m][i].krnl[window_circles[m][i][j][k] - 1], monomers[i][window_circles[m][i][j][k] - 1 + m]));
-                    std::fprintf(file, "%3.2f ", std::sqrt(distance2(window_kernel[m][i].krnl[window_circles[m][i][j][k] - 1], monomers[i][window_circles[m][i][j][k] - 1 + m])));
+                for (int k = 0; k < window_circles[i][m][j].size(); k++) {
+                    temp += std::sqrt(distance2(kernel[m].krnl[window_circles[i][m][j][k] - 1], monomers[i][window_circles[i][m][j][k] - 1 + i]));
+                    std::fprintf(file, "%3.2f ", std::sqrt(distance2(kernel[m].krnl[window_circles[i][m][j][k] - 1], monomers[i][window_circles[i][m][j][k] - 1 + i])));
                 }
-                std::fprintf(file, "average: %3.2f\n", temp / window_circles[m][i][j].size());
+                std::fprintf(file, "average: %3.2f\n", temp / window_circles[i][m][j].size());
             }
             std::fprintf(file, "# of monomers per coil:\n");
-            for (int j = 0; j < window_circles[m][i].size(); j++) {
-                std::fprintf(file, "%3.2d ", window_circles[m][i][j].size());
+            for (int j = 0; j < window_circles[i][m].size(); j++) {
+                std::fprintf(file, "%3.2d ", window_circles[i][m][j].size());
             }
             std::fprintf(file, "\n\n");
         }