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 390e0e0319806e4db4bbd80e7397a876fdf07e17..6e88dab7b0dc1e4ae6bebf6c77795cb74f23a636 100644 (file)
--- a/deal.II/source/fe/fe_trace.cc
+++ b/deal.II/source/fe/fe_trace.cc
@@ -101,6 +101,9 @@ FE_TraceQ<dim,spacedim>::FE_TraceQ (const unsigned int degree)
                   FiniteElementData<dim>(get_dpo_vector(degree), 1, degree, FiniteElementData<dim>::L2),
                   std::vector<bool>(1,true))
 {  
+  Assert (degree > 0,
+          ExcMessage ("FE_Trace can only be used for polynomial degrees "
+                      "greater than zero"));
   std::vector<unsigned int> renumber (this->dofs_per_face);
   FETools::hierarchic_to_lexicographic_numbering<dim-1> (degree, renumber);
   this->poly_space.set_numbering(renumber);
@@ -232,6 +235,7 @@ template <int dim, int spacedim>
 std::vector<unsigned int>
 FE_TraceQ<dim,spacedim>::get_dpo_vector (const unsigned int deg)
 {
+  AssertThrow(deg>0,ExcMessage("FE_TraceQ needs to be of degree > 0."));
   std::vector<unsigned int> dpo(dim+1, 1U);
   dpo[dim]=0;
   dpo[0]=1;