From 800a3a73161af1ceb58feefc6e6808cbcfaf0c9a Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 14 May 2025 20:02:41 -0300
Subject: Make sure β doesn't drift
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 get_energy.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

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));
-- 
cgit v1.2.3-70-g09d2