examples: boris: enable residual control for MLSDC

Author: torbjoernk

examples/boris/boris_mlsdc.cpp

@@ -22,7 +22,8 @@ namespace pfasst
     {
       template<typename scalar>
       error_map<scalar> run_boris_sdc(const size_t nsteps, const scalar dt, const size_t nnodes,
-                                      const size_t nparticles, const size_t niters)
+                                      const size_t nparticles, const size_t niters,
+                                      const double abs_res_tol, const double rel_res_tol)
       {
         MLSDC<> controller;
 
@@ -33,6 +34,7 @@ namespace pfasst
           make_shared<bindings::WrapperSimplePhysicsSolver<double, double>>();
         bindings::setup(dynamic_pointer_cast<bindings::WrapperSimplePhysicsSolver<double, double>>(impl_solver));
 
+        // fine level
         auto quad1        = quadrature::quadrature_factory<double>(nnodes,
                                                                    quadrature::QuadratureType::GaussLobatto);
         auto factory1     = make_shared<ParticleCloudFactory<double>>(nparticles, 3, mass, charge);
@@ -43,8 +45,10 @@ namespace pfasst
         auto transfer1    = make_shared<InjectiveTransfer<double, double>>();
         sweeper1->set_quadrature(quad1);
         sweeper1->set_factory(factory1);
+        sweeper1->set_residual_tolerances(abs_res_tol, rel_res_tol);
         controller.add_level(sweeper1, transfer1);
 
+        // coarse level
         auto quad2        = quadrature::quadrature_factory<double>(nnodes,
                                                                    quadrature::QuadratureType::GaussLobatto);
         auto factory2     = make_shared<ParticleCloudFactory<double>>(nparticles, 3, mass, charge);
@@ -95,13 +99,17 @@ int main(int argc, char** argv)
   const size_t nnodes     = pfasst::config::get_value<size_t>("num_nodes", 5);
   const size_t nparticles = pfasst::config::get_value<size_t>("num_particles", 1);
   const size_t niters     = pfasst::config::get_value<size_t>("num_iter", 2);
+  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);
 
   LOG(INFO) << "nsteps=" << nsteps << ", "
             << "dt=" << dt << ", "
             << "nnodes=" << nnodes << ", "
             << "nparticles=" << nparticles << ", "
-            << "niter=" << niters;
+            << "niter=" << niters << ", "
+            << "abs res=" << abs_res_tol << ", "
+            << "rel res=" << rel_res_tol;
 
-  pfasst::examples::boris::run_boris_sdc<double>(nsteps, dt, nnodes, nparticles, niters);
+  pfasst::examples::boris::run_boris_sdc<double>(nsteps, dt, nnodes, nparticles, niters, abs_res_tol, rel_res_tol);
 }
 #endif