Implement Simplex::FE_Poly::get_constant_modes

diff --git a/include/deal.II/simplex/fe_lib.h b/include/deal.II/simplex/fe_lib.h
index e7c45a8ed97a858a3ac128d307114f509ea93706..201620d5437588ea148a63c00d140936c9b79f75 100644 (file)
--- a/include/deal.II/simplex/fe_lib.h
+++ b/include/deal.II/simplex/fe_lib.h
@@ -43,6 +43,13 @@ namespace Simplex
     FE_Poly(const unsigned int               degree,
             const std::vector<unsigned int> &dpo_vector);
 
+    /**
+     * Return a list of constant modes of the element. For this element, the
+     * list consists of true arguments for all components.
+     */
+    std::pair<Table<2, bool>, std::vector<unsigned int>>
+    get_constant_modes() const override;
+
   private:
     /**
      * @copydoc dealii::FiniteElement::convert_generalized_support_point_values_to_dof_values()

diff --git a/source/simplex/fe_lib.cc b/source/simplex/fe_lib.cc
index abd6e5904290d3147afa265f8ec017deb70211e4..93fbed41cbc690fe64472b8d17c9f158221e0cd0 100644 (file)
--- a/source/simplex/fe_lib.cc
+++ b/source/simplex/fe_lib.cc
@@ -228,6 +228,18 @@ namespace Simplex
 
 
 
+  template <int dim, int spacedim>
+  std::pair<Table<2, bool>, std::vector<unsigned int>>
+  FE_Poly<dim, spacedim>::get_constant_modes() const
+  {
+    Table<2, bool> constant_modes(1, this->n_dofs_per_cell());
+    constant_modes.fill(true);
+    return std::pair<Table<2, bool>, std::vector<unsigned int>>(
+      constant_modes, std::vector<unsigned int>(1, 0));
+  }
+
+
+
   template <int dim, int spacedim>
   void
   FE_Poly<dim, spacedim>::