From cc70aa51c3f65f3d771cd027595ef2f586bb3d95 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Thu, 3 Apr 2025 17:57:00 -0300
Subject: Remove allocation from loop to speed up energy colculation

---
 get_energy.cpp | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/get_energy.cpp b/get_energy.cpp
index b1527dc..b6e835f 100644
--- a/get_energy.cpp
+++ b/get_energy.cpp
@@ -57,10 +57,10 @@ int main(int argc, char* argv[]) {
 
   FourierTransform fft(n, Δω, Δτ, FFTW_ESTIMATE);
 
-  while (y += Δy, y <= yₘₐₓ) {
-    std::vector<Real> C(2 * n);
-    std::vector<Real> R(2 * n);
+  std::vector<Real> C(2 * n);
+  std::vector<Real> R(2 * n);
 
+  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));
-- 
cgit v1.2.3-70-g09d2