Save every correlation function.

1 files changed, 16 insertions, 10 deletions

diff --git a/fourier_integrator.cpp b/fourier_integrator.cpp
index dd1b2d3..db5e07e 100644
--- a/fourier_integrator.cpp
+++ b/fourier_integrator.cpp
@@ -198,23 +198,29 @@ int main(int argc, char* argv[]) {
 
       z *= Cnew[0];
 
-      Real energy = 0;
-
-      for (unsigned i = 0; i < n; i++) {
-        energy += y * R[i] * df(λ, p, s, C[i]) * M_PI * Δτ;
-      }
-
-      std::cerr << "y " << y << " iteration " << it << ": " << sqrt(ΔC / C.size()) << " " << pow(y, 2) / z << " " << energy << std::endl;
 
       if (it > maxIterations) {
         it = 0;
         γ /= 2;
       }
     }
-  }
 
-  for (unsigned i = 0; i < C.size(); i++) {
-    std::cout << i * Δτ * M_PI << " " << C[i] << std::endl;
+    Real energy = 0;
+
+    for (unsigned i = 0; i < n; i++) {
+      energy += y * R[i] * df(λ, p, s, C[i]) * M_PI * Δτ;
+    }
+
+    std::cerr << "y " << y << " " << energy << std::endl;
+
+    for (unsigned i = 0; i < C.size(); i++) {
+      std::cout << i * Δτ * M_PI << " " << C[i] << " ";
+    }
+    std::cout << std::endl;
+    for (unsigned i = 0; i < R.size(); i++) {
+      std::cout << i * Δτ * M_PI << " " << R[i] << " ";
+    }
+    std::cout << std::endl;
   }
 
   return 0;