From: buerg Date: Sun, 29 Jan 2012 10:52:21 +0000 (+0000) Subject: Add asserts for FESystem. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=4267a5713549ea1c6e3aa7397c556cb96b8a3e26;p=dealii-svn.git Add asserts for FESystem. git-svn-id: https://svn.dealii.org/trunk@24950 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/include/deal.II/numerics/vectors.templates.h b/deal.II/include/deal.II/numerics/vectors.templates.h index 578c21b52b..0bb9187283 100644 --- a/deal.II/include/deal.II/numerics/vectors.templates.h +++ b/deal.II/include/deal.II/numerics/vectors.templates.h @@ -3307,13 +3307,19 @@ namespace VectorTools if (dynamic_cast*> (&cell->get_fe ()) != 0) return; - // this is only + // This is only // implemented, if the // FE is a Nedelec - // element - typedef FiniteElement FEL; - AssertThrow (dynamic_cast*> (&cell->get_fe ()) != 0, - typename FEL::ExcInterpolationNotImplemented ()); + // element. If the FE is + // a FESystem we cannot + // check this. + if (dynamic_cast*> (&cell->get_fe ()) == 0) + { + typedef FiniteElement FEL; + + AssertThrow (dynamic_cast*> (&cell->get_fe ()) != 0, + typename FEL::ExcInterpolationNotImplemented ()); + } for (unsigned int dof = 0; dof < dofs_per_face; ++dof) dof_values[dof] = 0.0;