From: David Wells <drwells@email.unc.edu>
Date: Wed, 18 Nov 2020 22:04:43 +0000 (-0500)
Subject: Add simplex elements to the MF permitted list.
X-Git-Tag: v9.3.0-rc1~885^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F11203%2Fhead;p=dealii.git

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 802607bfd6..81280e4f48 100644
--- 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<