examples: advec: enable cmd arg for residual tolerances

Author: torbjoernk

examples/advection_diffusion/serial_mlsdc.cpp

@@ -39,6 +39,9 @@ namespace pfasst
         size_t nnodes = config::get_value<size_t>("num_nodes", 5);
         size_t ndofs  = config::get_value<size_t>("spatial_dofs", 128);
 
+        const double abs_res_tol = pfasst::config::get_value<double>("abs_res_tol", 0.0);
+        const double rel_res_tol = pfasst::config::get_value<double>("rel_res_tol", 0.0);
+
         /*
          * build space/time discretisation levels and add them to mlsdc
          * controller.  this loop adds the finest level first, and
@@ -53,6 +56,7 @@ namespace pfasst
 
           sweeper->set_quadrature(quad);
           sweeper->set_factory(factory);
+          sweeper->set_residual_tolerances(abs_res_tol, rel_res_tol);
 
           mlsdc.add_level(sweeper, transfer);
 

examples/advection_diffusion/vanilla_sdc.cpp

@@ -41,7 +41,7 @@ namespace pfasst
 
         sweeper->set_quadrature(quad);
         sweeper->set_factory(factory);
-        sweeper->set_residual_tolerances(abs_residual_tol, 0.0);
+        sweeper->set_residual_tolerances(abs_residual_tol, pfasst::config::get_value<double>("rel_res_tol", 0.0));
 
         sdc.add_level(sweeper);
         sdc.set_duration(0.0, nsteps*dt, dt, niters);
@@ -68,6 +68,8 @@ int main(int argc, char** argv)
   pfasst::init(argc, argv);
   pfasst::log::add_custom_logger("Advec");
 
-  pfasst::examples::advection_diffusion::run_vanilla_sdc(0.0);
+  const double abs_res_tol = pfasst::config::get_value<double>("abs_res_tol", 0.0);
+
+  pfasst::examples::advection_diffusion::run_vanilla_sdc(abs_res_tol);
 }
 #endif