Add assert to MappingFE

diff --git a/include/deal.II/fe/mapping_fe.h b/include/deal.II/fe/mapping_fe.h
index 2c617b367d7335379a55c8b299b68734300b175c..aa96493fcb8e6161c51ce8c8ff7dd81697ebdb0d 100644 (file)
--- a/include/deal.II/fe/mapping_fe.h
+++ b/include/deal.II/fe/mapping_fe.h
@@ -49,7 +49,8 @@ DEAL_II_NAMESPACE_OPEN
  * equivalent to MappingQGeneric(degree). Please note that no optimizations
  * exploiting tensor-product structures of finite elements have been added here.
  *
- * @note Currently, only implemented for degree==1.
+ * @node Currently, only implemented for elements with tensor_degree==1 and
+ *   n_components==1.
  *
  * @ingroup simplex
  */

diff --git a/source/fe/mapping_fe.cc b/source/fe/mapping_fe.cc
index 746abb6785733e7157ca25adf210a0a983eefe61..16fa1546411b09afe2d65c0286b8b8518da88e0f 100644 (file)
--- a/source/fe/mapping_fe.cc
+++ b/source/fe/mapping_fe.cc
@@ -926,6 +926,8 @@ MappingFE<dim, spacedim>::MappingFE(const FiniteElement<dim, spacedim> &fe)
   Assert(polynomial_degree >= 1,
          ExcMessage("It only makes sense to create polynomial mappings "
                     "with a polynomial degree greater or equal to one."));
+  Assert(fe.tensor_degree() == 1, ExcNotImplemented());
+  Assert(fe.n_components() == 1, ExcNotImplemented());
 }