From: wolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Mon, 24 Feb 2003 15:51:38 +0000 (+0000)
Subject: Check for primitivity in one function.
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=036890b5a52b936e8ea2eacfbccd00eb3f2ffa4b;p=dealii-svn.git

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 52e24ec3b7..b53f65c973 100644
--- 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(),