Make sure β doesn't drift

1 files changed, 1 insertions, 1 deletions

diff --git a/get_energy.cpp b/get_energy.cpp
index 70b8eb8..deee26b 100644
--- a/get_energy.cpp
+++ b/get_energy.cpp
@@ -63,7 +63,7 @@ int main(int argc, char* argv[]) {
 
   std::cout << std::setprecision(15);
 
-  while (y += Δy, y <= yₘₐₓ) {
+  while (y = std::round(1e6 * (y + Δy)) / 1e6, 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() / 2) * sizeof(Real));