summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2025-05-10 17:00:00 -0300
committerJaron Kent-Dobias <jaron@kent-dobias.com>2025-05-10 17:00:00 -0300
commit6878bab16270855178d1b2b0fb22285894f62827 (patch)
tree21631d03ac7292ba4de72d5e21a1e2aa3efe8664
parented91704b302cbe3500378810aac6acad342e0cc9 (diff)
downloadcode-6878bab16270855178d1b2b0fb22285894f62827.tar.gz
code-6878bab16270855178d1b2b0fb22285894f62827.tar.bz2
code-6878bab16270855178d1b2b0fb22285894f62827.zip
Tweak
-rw-r--r--log-fourier_integrator.cpp17
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;