Add a bunch of tests for quadratic PolynomialReconstruction

src/scheme/test_reconstruction.cpp

@@ -7,10 +7,36 @@ void
 test_reconstruction(const std::vector<std::shared_ptr<const DiscreteFunctionVariant>>& discrete_function_variant_list,
                     uint64_t degree)
 {
-  PolynomialReconstructionDescriptor descriptor{PolynomialReconstructionDescriptor::IntegrationMethod::element, degree};
-  {
-    const size_t nb_loops = 1;
-    std::cout << "** variable list at once (" << nb_loops << " times)\n";
+  std::vector descriptor_list =
+    {PolynomialReconstructionDescriptor{PolynomialReconstructionDescriptor::IntegrationMethod::cell_center, degree},
+     PolynomialReconstructionDescriptor{PolynomialReconstructionDescriptor::IntegrationMethod::element, degree},
+     PolynomialReconstructionDescriptor{PolynomialReconstructionDescriptor::IntegrationMethod::boundary, degree}};
+
+  [[maybe_unused]] auto x =
+    PolynomialReconstruction{
+      PolynomialReconstructionDescriptor{PolynomialReconstructionDescriptor::IntegrationMethod::cell_center, degree}}
+      .build(discrete_function_variant_list);
+
+  for (auto&& descriptor : descriptor_list) {
+    std::string method_name;
+    switch (descriptor.integrationMethod()) {
+    case PolynomialReconstructionDescriptor::IntegrationMethod::element: {
+      method_name = "element";
+      break;
+    }
+    case PolynomialReconstructionDescriptor::IntegrationMethod::boundary: {
+      method_name = "boundary";
+      break;
+    }
+    case PolynomialReconstructionDescriptor::IntegrationMethod::cell_center: {
+      method_name = "cell_center";
+      break;
+    }
+    }
+
+    const size_t nb_loops = 50;
+    std::cout << "** variable list at once (" << nb_loops << " times) using " << rang::fgB::yellow << method_name
+              << rang::fg::reset << "\n";
 
     PolynomialReconstruction reconstruction{descriptor};
     Timer t;