From ed91704b302cbe3500378810aac6acad342e0cc9 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 10 May 2025 16:54:21 -0300
Subject: Enforce monotony of R

---
 log-fourier_integrator.cpp | 7 +++++++
 1 file changed, 7 insertions(+)

(limited to 'log-fourier_integrator.cpp')

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