Add stencil descriptor to polynomial reconstruction descriptor

This is optional by now. It may
- become mandatory
- evolve to take into account better default value it is deduced from
reconstruction degree

src/scheme/PolynomialReconstructionDescriptor.hpp

@@ -2,6 +2,7 @@
 #define POLYNOMIAL_RECONSTRUCTION_DESCRIPTOR_HPP
 
 #include <mesh/IBoundaryDescriptor.hpp>
+#include <mesh/StencilDescriptor.hpp>
 #include <scheme/IntegrationMethodType.hpp>
 #include <utils/PugsMacros.hpp>
 
@@ -17,6 +18,7 @@ class PolynomialReconstructionDescriptor
  private:
   IntegrationMethodType m_integration_method;
   size_t m_degree;
+  StencilDescriptor m_stencil_descriptor;
 
   BoundaryDescriptorList m_symmetry_boundary_descriptor_list;
 
@@ -79,11 +81,34 @@ class PolynomialReconstructionDescriptor
     m_row_weighting = row_weighting;
   }
 
+  PolynomialReconstructionDescriptor(const IntegrationMethodType integration_method, const size_t degree)
+    : m_integration_method{integration_method},
+      m_degree{degree},
+      m_stencil_descriptor(degree, StencilDescriptor::Type::by_nodes)
+  {}
+
+  PolynomialReconstructionDescriptor(const IntegrationMethodType integration_method,
+                                     const size_t degree,
+                                     const BoundaryDescriptorList& symmetry_boundary_descriptor_list)
+    : m_integration_method{integration_method},
+      m_degree{degree},
+      m_stencil_descriptor(degree, StencilDescriptor::Type::by_nodes),
+      m_symmetry_boundary_descriptor_list(symmetry_boundary_descriptor_list)
+  {}
+
+  PolynomialReconstructionDescriptor(const IntegrationMethodType integration_method,
+                                     const size_t degree,
+                                     const StencilDescriptor& stencil_descriptor)
+    : m_integration_method{integration_method}, m_degree{degree}, m_stencil_descriptor{stencil_descriptor}
+  {}
+
   PolynomialReconstructionDescriptor(const IntegrationMethodType integration_method,
                                      const size_t degree,
-                                     const BoundaryDescriptorList& symmetry_boundary_descriptor_list = {})
+                                     const StencilDescriptor& stencil_descriptor,
+                                     const BoundaryDescriptorList& symmetry_boundary_descriptor_list)
     : m_integration_method{integration_method},
       m_degree{degree},
+      m_stencil_descriptor{stencil_descriptor},
       m_symmetry_boundary_descriptor_list(symmetry_boundary_descriptor_list)
   {}