Check for primitivity in one function.

git-svn-id: https://svn.dealii.org/trunk@7229 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/deal.II/source/fe/fe_tools.cc b/deal.II/deal.II/source/fe/fe_tools.cc
index 52e24ec3b7103743c2bf93ef08cc8568944dcf87..b53f65c97393d0c30ac337ac1efdcf2da99c9940 100644 (file)
--- a/deal.II/deal.II/source/fe/fe_tools.cc
+++ b/deal.II/deal.II/source/fe/fe_tools.cc
@@ -35,6 +35,8 @@ void FETools::get_interpolation_matrix(const FiniteElement<dim> &fe1,
 {
   Assert (fe1.n_components() == fe2.n_components(),
          ExcDimensionMismatch(fe1.n_components(), fe2.n_components()));
+  Assert (fe1.is_primitive() == true, ExcFEMustBePrimitive());
+  Assert (fe2.is_primitive() == true, ExcFEMustBePrimitive());  
   Assert(interpolation_matrix.m()==fe2.dofs_per_cell &&
         interpolation_matrix.n()==fe1.dofs_per_cell,
         ExcMatrixDimensionMismatch(interpolation_matrix.m(),
@@ -74,6 +76,8 @@ void FETools::get_back_interpolation_matrix(const FiniteElement<dim> &fe1,
 {
   Assert (fe1.n_components() == fe2.n_components(),
          ExcDimensionMismatch(fe1.n_components(), fe2.n_components()));
+  Assert (fe1.is_primitive() == true, ExcFEMustBePrimitive());
+  Assert (fe2.is_primitive() == true, ExcFEMustBePrimitive());  
   Assert(interpolation_matrix.m()==fe1.dofs_per_cell &&
         interpolation_matrix.n()==fe1.dofs_per_cell, 
         ExcMatrixDimensionMismatch(interpolation_matrix.m(),
@@ -100,6 +104,8 @@ void FETools::get_interpolation_difference_matrix(const FiniteElement<dim> &fe1,
 {
   Assert (fe1.n_components() == fe2.n_components(),
          ExcDimensionMismatch(fe1.n_components(), fe2.n_components()));
+  Assert (fe1.is_primitive() == true, ExcFEMustBePrimitive());
+  Assert (fe2.is_primitive() == true, ExcFEMustBePrimitive());  
   Assert(difference_matrix.m()==fe1.dofs_per_cell &&
         difference_matrix.n()==fe1.dofs_per_cell, 
         ExcMatrixDimensionMismatch(difference_matrix.m(),