Save half the data to the disk

1 files changed, 5 insertions, 2 deletions

diff --git a/get_energy.cpp b/get_energy.cpp
index 6e14d9b..ae77595 100644
--- a/get_energy.cpp
+++ b/get_energy.cpp
@@ -63,11 +63,14 @@ int main(int argc, char* argv[]) {
   while (y += Δy, y <= yₘₐₓ) {
     std::ifstream cfile(fourierFile("C", p, s, λ, τ₀, y, log2n, τₘₐₓ), std::ios::binary);
     if (cfile.is_open()) {
-      cfile.read((char*)(C.data()), C.size() * sizeof(Real));
+      cfile.read((char*)(C.data()), (C.size() / 2) * sizeof(Real));
       cfile.close();
+      for (unsigned i = 1; i < n; i++) {
+        C[2 * n - i] = C[i];
+      }
 
       std::ifstream rfile(fourierFile("R", p, s, λ, τ₀, y, log2n, τₘₐₓ), std::ios::binary);
-      rfile.read((char*)(R.data()), R.size() * sizeof(Real));
+      rfile.read((char*)(R.data()), (R.size() / 2) * sizeof(Real));
       rfile.close();
 
       Real e = energy(C, R, p, s, λ, y, Δτ);