diff options
| -rw-r--r-- | log-fourier_integrator.cpp | 3 |
1 files changed, 3 insertions, 0 deletions
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); |