From 036890b5a52b936e8ea2eacfbccd00eb3f2ffa4b Mon Sep 17 00:00:00 2001 From: wolf Date: Mon, 24 Feb 2003 15:51:38 +0000 Subject: [PATCH] Check for primitivity in one function. git-svn-id: https://svn.dealii.org/trunk@7229 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/fe_tools.cc | 6 ++++++ 1 file changed, 6 insertions(+) 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 &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 &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 &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(), -- 2.39.5