From 7d16074ea49e563bd4734c582b34b7da28d9b38c Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Sun, 17 Aug 2014 19:35:36 +0200
Subject: [PATCH] [Algorit] Terminate time stepping properly
---
src/adaptivetimestepper.hh | 7 ++-----
1 file changed, 2 insertions(+), 5 deletions(-)
diff --git a/src/adaptivetimestepper.hh b/src/adaptivetimestepper.hh
index 1b2a819d..0a034172 100644
--- a/src/adaptivetimestepper.hh
+++ b/src/adaptivetimestepper.hh
@@ -39,15 +39,12 @@ template <class Factory, class UpdaterPair> class AdaptiveTimeStepper {
iterationWriter_ << std::endl;
}
- // FIXME
- bool reachedEnd() { return relativeTime_ > 1.0 - 1e-10; }
+ bool reachedEnd() { return relativeTime_ >= 1.0; }
bool coarsen() {
bool didCoarsen = false;
- // FIXME: for a constant function, e.g., we will not only overstep but
- // diverge
- while (true) {
+ while (relativeTime_ + relativeTau_ < 1.0) {
R2_ = clonePair(R1_);
{
MyCoupledTimeStepper coupledTimeStepper(finalTime_, factory_, parset_,
--
GitLab