config: Add abs/rel residual tolerances.

Signed-off-by: Matthew Emmett <[email protected]>

Author: memmett

include/pfasst/config.hpp

@@ -171,6 +171,9 @@ namespace pfasst
       options::add_option<double>("Duration", "tend", "final time of simulation");
       options::add_option<size_t>("Duration", "num_iters", "number of iterations");
 
+      options::add_option<double>("Tolerances", "abs_res_tol", "absolute residual tolerance");
+      options::add_option<double>("Tolerances", "rel_res_tol", "relative residual tolerance");
+
       options::get_instance().init();
     }
 

include/pfasst/encap/encap_sweeper.hpp

@@ -92,6 +92,13 @@ namespace pfasst
         }
         //! @}
 
+        //! @{
+        virtual void set_options() override
+        {
+          this->abs_residual_tol = time(config::get_value<double>("abs_res_tol", this->abs_residual_tol));
+          this->rel_residual_tol = time(config::get_value<double>("rel_res_tol", this->rel_residual_tol));
+        }
+
         virtual void setup(bool coarse) override
         {
           auto const nodes = this->quadrature->get_nodes();
@@ -113,6 +120,7 @@ namespace pfasst
             }
           }
         }
+        //! @}
 
         //! @{
         virtual void spread() override
@@ -121,7 +129,6 @@ namespace pfasst
             this->state[m]->copy(this->state[0]);
           }
         }
-        //! @}
 
         /**
          * Save current solution states.
@@ -137,6 +144,7 @@ namespace pfasst
             }
           }
         }
+        //! @}
 
         //! @{
         void set_quadrature(shared_ptr<IQuadrature<time>> quadrature)