tests: Add SerialMLSDC test.

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

Author: memmett

examples/advection_diffusion/serial_mlsdc.cpp

@@ -18,7 +18,7 @@
 using namespace pfasst;
 using namespace pfasst::encap;
 
-int main(int argc, char** argv)
+error_map run_serial_mlsdc()
 {
   MLSDC<> mlsdc;
 
@@ -74,4 +74,13 @@ int main(int argc, char** argv)
   mlsdc.run();
 
   fftw_cleanup();
+
+  return sweeper->get_errors();
+}
+
+#ifndef PFASST_UNIT_TESTING
+int main(int argc, char** argv)
+{
+  run_serial_mlsdc();
 }
+#endif

tests/test-advection-diffusion.cpp

@@ -12,6 +12,7 @@
 
 #define PFASST_UNIT_TESTING
 #include "../examples/advection_diffusion/vanilla_sdc.cpp"
+#include "../examples/advection_diffusion/serial_mlsdc.cpp"
 #undef PFASST_UNIT_TESTING
 
 using namespace std;
@@ -29,7 +30,7 @@ MATCHER(DoubleLess, "")
 
 TEST(ErrorTest, VanillaSDC)
 {
-  typedef pair<pair<size_t, size_t>, double> vtype;
+  typedef error_map::value_type vtype;
 
   auto errors = run_vanilla_sdc();
   auto get_iter  = [](const vtype x) { return get<1>(get<0>(x)); };
@@ -50,6 +51,29 @@ TEST(ErrorTest, VanillaSDC)
   EXPECT_THAT(err, testing::Pointwise(DoubleLess(), tol));
 }
 
+TEST(ErrorTest, SerialMLSDC)
+{
+  typedef error_map::value_type vtype;
+
+  auto errors = run_serial_mlsdc();
+  auto get_iter  = [](const vtype x) { return get<1>(get<0>(x)); };
+  auto get_step  = [](const vtype x) { return get<0>(get<0>(x)); };
+  auto get_error = [](const vtype x) { return get<1>(x); };
+
+  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 = { 8e-11, 8e-11, 8e-11, 8e-11 };
+  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));
+}
+
 int main(int argc, char** argv)
 {
   testing::InitGoogleTest(&argc, argv);