Add simplex elements to the MF permitted list.

diff --git a/include/deal.II/matrix_free/shape_info.templates.h b/include/deal.II/matrix_free/shape_info.templates.h
index 802607bfd6022c4566a57bf39a530ae5586d92ac..81280e4f480239b4dd9c6019e9080d31d182f49d 100644 (file)
--- a/include/deal.II/matrix_free/shape_info.templates.h
+++ b/include/deal.II/matrix_free/shape_info.templates.h
@@ -106,6 +106,14 @@ namespace internal
             {
               const FE_Poly<dim, spacedim> *fe_poly_ptr =
                 dynamic_cast<const FE_Poly<dim, spacedim> *>(fe_ptr);
+#ifdef DEAL_II_WITH_SIMPLEX_SUPPORT
+              // Simplices are a special case since the polynomial family is not
+              // indicative of their support
+              if (dynamic_cast<const Simplex::FE_P<dim> *>(fe_poly_ptr) ||
+                  dynamic_cast<const Simplex::FE_DGP<dim> *>(fe_poly_ptr))
+                return true;
+#endif
+
               if (dynamic_cast<const TensorProductPolynomials<dim> *>(
                     &fe_poly_ptr->get_poly_space()) == nullptr &&
                   dynamic_cast<const TensorProductPolynomials<