tests: Add RelativeAdaptiveErrorTest to test_advection_diffusion.cpp.

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

Author: memmett

examples/advection_diffusion/vanilla_sdc.cpp

@@ -23,7 +23,7 @@ namespace pfasst
   {
     namespace advection_diffusion
     {
-      error_map run_vanilla_sdc(double abs_residual_tol)
+      error_map run_vanilla_sdc(double abs_residual_tol, double rel_residual_tol=0.0)
       {
         SDC<> sdc;
 
@@ -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, rel_residual_tol);
 
         sdc.add_level(sweeper);
         sdc.set_duration(0.0, nsteps*dt, dt, niters);

tests/examples/advection_diffusion/test_advection_diffusion.cpp

@@ -70,6 +70,29 @@ TEST(ErrorTest, VanillaSDC)
   }
 }
 
+TEST(RelativeAdaptiveErrorTest, VanillaSDC)
+{
+  typedef error_map::value_type vtype;
+
+  auto get_iter  = [](const vtype x) { return get<1>(get<0>(x)); };
+  auto get_error = [](const vtype x) { return get<1>(x); };
+
+  auto errors = run_vanilla_sdc(0.0, 1.e-6);
+  auto max_iter = get_iter(*std::max_element(errors.begin(), errors.end(),
+                                             [get_iter](const vtype p1, const vtype p2) { return get_iter(p1) < get_iter(p2); }));
+
+  vector<double> tol = { 5e-8, 5e-8, 5e-8, 5e-8 };
+  vector<double> err;
+  for (auto& x: errors) {
+    if (get_iter(x) == max_iter) {
+      err.push_back(get_error(x));
+    }
+  }
+
+  EXPECT_THAT(err, testing::Pointwise(DoubleLess(), tol));
+  ASSERT_EQ(max_iter, (size_t) 2);
+}
+
 TEST(ErrorTest, SerialMLSDC)
 {
   typedef error_map::value_type vtype;