added: stopping criterion in time stepping

src/ASM/TimeDomain.h

@@ -25,15 +25,16 @@ struct TimeDomain
   double  dt; //!< Current timestep (or load parameter) increment
   double dtn; //!< Previous timestep (or load parameter) increment
   double CFL; //!< Current CFL number (used by CFD simulators)
+  double incNorm; //!< Norm of change between timesteps
   int     it; //!< Current iteration within current time/load step
   char first; //!< If \e true, this is the first load/time step
 
   //! \brief Default constructor.
   explicit TimeDomain(int i = 0, bool f = true) : it(i), first(f)
-  { t = dt = dtn = CFL = 0.0; }
+  { t = dt = dtn = CFL = incNorm = 0.0; }
   //! \brief Constructor for linear problems with fixed (non-zero) time.
   explicit TimeDomain(double t0) : t(t0), it(0), first(true)
-  { dt = dtn = CFL = 0.0; }
+  { dt = dtn = CFL = incNorm = 0.0; }
 };
 
 #endif

src/SIM/TimeStep.C

@@ -32,6 +32,7 @@ TimeStep::TimeStep () : step(0), iter(time.it), lstep(0)
   f1 = 1.5;
   f2 = 0.25;
   maxStep = niter = 0;
+  stop_tol = 0.0;
   stepIt = mySteps.end();
 }
 
@@ -51,6 +52,7 @@ TimeStep& TimeStep::operator= (const TimeStep& ts)
   maxStep = ts.maxStep;
   niter = ts.niter;
   iter = ts.iter;
+  stop_tol = ts.stop_tol;
 
   time = ts.time;
   mySteps = ts.mySteps;
@@ -110,6 +112,7 @@ bool TimeStep::parse (const TiXmlElement* elem)
   utl::getAttribute(elem,"dtMax",dtMax);
   utl::getAttribute(elem,"maxCFL",maxCFL);
   utl::getAttribute(elem,"nInitStep",nInitStep);
+  utl::getAttribute(elem,"stop_tolerance",stop_tol);
   utl::getAttribute(elem,"maxStep",maxStep);
   utl::getAttribute(elem,"f1",f1);
   utl::getAttribute(elem,"f2",f2);
@@ -266,6 +269,12 @@ bool TimeStep::increment ()
                << std::endl;
     return false;
   }
+  else if (time.incNorm > 0.0 && time.incNorm < stop_tol)
+  {
+    IFEM::cout <<"\n  ** Terminating, solution increment norm "
+               << time.incNorm << " less than tolerance " << stop_tol << std::endl;
+    return false;
+  }
 
   niter = iter;
   time.t += time.dt;

src/SIM/TimeStep.h

@@ -62,6 +62,9 @@ class TimeStep
   //! \brief Returns \e true if the end of the simulation has been reached.
   bool finished() const;
 
+  //! \brief Returns configured stopping tolerance.
+  double stopTolerance() const { return stop_tol; }
+
   //! \brief Serialize internal state for restarting purposes.
   //! \param data Container for serialized data
   bool serialize(std::map<std::string,std::string>& data) const;
@@ -87,6 +90,8 @@ class TimeStep
   double f1;     //!< Scale factor for increased time step size
   double f2;     //!< Scale factor for reduced time step size
 
+  double stop_tol; //!< Stop time step loop if incNorm < stop_tol
+
   typedef std::pair<std::vector<double>,double> Step; //!< Time step definition
   typedef std::vector<Step> TimeSteps;                //!< Time step container