From: angela.klewinghaus Date: Mon, 11 Nov 2013 13:10:56 +0000 (+0000) Subject: add Asserts for polynomial degree of FE_TraceQ X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=4afbbd93c00663ac29f57edaed9bbd0a74c6bde8;p=dealii-svn.git add Asserts for polynomial degree of FE_TraceQ git-svn-id: https://svn.dealii.org/branches/branch_face_elements@31623 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/source/fe/fe_trace.cc b/deal.II/source/fe/fe_trace.cc index 390e0e0319..6e88dab7b0 100644 --- a/deal.II/source/fe/fe_trace.cc +++ b/deal.II/source/fe/fe_trace.cc @@ -101,6 +101,9 @@ FE_TraceQ::FE_TraceQ (const unsigned int degree) FiniteElementData(get_dpo_vector(degree), 1, degree, FiniteElementData::L2), std::vector(1,true)) { + Assert (degree > 0, + ExcMessage ("FE_Trace can only be used for polynomial degrees " + "greater than zero")); std::vector renumber (this->dofs_per_face); FETools::hierarchic_to_lexicographic_numbering (degree, renumber); this->poly_space.set_numbering(renumber); @@ -232,6 +235,7 @@ template std::vector FE_TraceQ::get_dpo_vector (const unsigned int deg) { + AssertThrow(deg>0,ExcMessage("FE_TraceQ needs to be of degree > 0.")); std::vector dpo(dim+1, 1U); dpo[dim]=0; dpo[0]=1;