[Algorit] Terminate time stepping properly

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_,