From 6878bab16270855178d1b2b0fb22285894f62827 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 10 May 2025 17:00:00 -0300
Subject: Tweak

---
 log-fourier_integrator.cpp | 17 +++++++++++------
 1 file changed, 11 insertions(+), 6 deletions(-)

diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp
index e127fc1..1d64d69 100644
--- a/log-fourier_integrator.cpp
+++ b/log-fourier_integrator.cpp
@@ -147,6 +147,15 @@ int main(int argc, char* argv[]) {
 
       Real C₀ = Cₜ₊₁[0];
 
+      bool trigger0 = false;
+      bool trigger1 = false;
+      for (unsigned i = 0; i < N; i++) {
+        if (Rₜ₊₁[i] < ε || trigger0) {
+          Rₜ₊₁[i] = 0;
+          trigger0 = true;
+        }
+      }
+
 
       Real Rmax = 0;
       for (unsigned i = 0; i < N; i++) {
@@ -154,8 +163,8 @@ int main(int argc, char* argv[]) {
         Rₜ₊₁[N-1-i] = Rmax;
       }
 
-      bool trigger0 = false;
-      bool trigger1 = false;
+      trigger0 = false;
+      trigger1 = false;
       for (unsigned i = 0; i < N; i++) {
         if (Cₜ₊₁[i] < ε || trigger0) {
           Cₜ₊₁[i] = 0;
@@ -169,10 +178,6 @@ int main(int argc, char* argv[]) {
       trigger0 = false;
       trigger1 = false;
       for (unsigned i = 0; i < N; i++) {
-        if (Rₜ₊₁[i] < ε || trigger0) {
-          Rₜ₊₁[i] = 0;
-          trigger0 = true;
-        }
         if (Rₜ₊₁[N-1-i] > 1 - ε || trigger1) {
           Rₜ₊₁[N-1-i] = 1;
           trigger1 = true;
-- 
cgit v1.2.3-70-g09d2