From 75dc30ef643a62c1864435c8655b505da05a4bcb Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 10 May 2025 18:14:03 -0300
Subject: Make sure things stop on nan

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

diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp
index 12a3516..627f49e 100644
--- a/log-fourier_integrator.cpp
+++ b/log-fourier_integrator.cpp
@@ -147,6 +147,8 @@ int main(int argc, char* argv[]) {
 
       Real C₀ = Cₜ₊₁[0];
 
+      if (!std::isnan(Cₜ₊₁[0])) {
+
       bool trigger0 = false;
       bool trigger1 = false;
       for (unsigned i = 0; i < N; i++) {
@@ -189,6 +191,7 @@ int main(int argc, char* argv[]) {
         if (Cₜ₊₁[N-1-i] > Cmax) Cmax = Cₜ₊₁[N-1-i];
         Cₜ₊₁[N-1-i] = Cmax;
       }
+      }
 
       μₜ *= pow(tanh(C₀-1)+1, x);
 
-- 
cgit v1.2.3-70-g09d2