diff options
| -rw-r--r-- | log-fourier_integrator.cpp | 17 |
1 files 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; |