added output files

diff --git a/src/spirals.cpp b/src/spirals.cpp
index 15debf0b5391e1ff3184b470df42adbe55367cda..de214054d534e2dc9fe581b8f0ec04224777cb01 100644 (file)
--- a/src/spirals.cpp
+++ b/src/spirals.cpp
@@ -269,7 +269,43 @@ Spirals::analyzeFrame(int frnr, const t_trxframe &fr, t_pbc *pbc,
 void
 Spirals::finishAnalysis(int /*nframes*/)
 {
+    /*
+     *  kernel.resize(frames + 1);
+        circles.resize(frames + 1);
+        spiral_dist.resize(frames + 1);
+     */
+    FILE *file;
+    file = std::fopen("kernel.txt", "w+");
+    for (int i = 0; i < kernel.size(); i++) {
+        for (int j = 0; j < kernel[i].size(); j++) {
+            std::fprintf(file, "%3.2f 3.2f 3.2f 3.2f\n", kernel[i][j].x, kernel[i][j].y, kernel[i][j].z, kernel[i][j].r);
+        }
+        std::fprintf(file, "\n");
+    }
+    std::fclose(file);
+
+    FILE *file;
+    file = std::fopen("circles.txt", "w+");
+    for (int i = 0; i < circles.size(); i++) {
+        for (int j = 0; j < circles[i].size(); j++) {
+            for (int k = 0; k < circles[i][j].size(); k++) {
+                std::fprintf(file, "%3.2f ", circles[i][j][k]);
+            }
+            std::fprintf(file, "\n");
+        }
+        std::fprintf(file, "\n");
+    }
+    std::fclose(file);
 
+    FILE *file;
+    file = std::fopen("spiral_dist.txt", "w+");
+    for (int i = 0; i < spiral_dist.size(); i++) {
+        for (int j = 0; j < spiral_dist[i].size(); j++) {
+            std::fprintf(file, "%3.2f\n", spiral_dist[i][j]);
+        }
+        std::fprintf(file, "\n");
+    }
+    std::fclose(file);
 }
 
 void