From: kronbichler Date: Thu, 17 Jul 2014 21:36:31 +0000 (+0000) Subject: Compute unit face support points of FE_TraceQ X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=44f1aeee663c4f0b4da3a6a68302f976daa23849;p=dealii-svn.git Compute unit face support points of FE_TraceQ git-svn-id: https://svn.dealii.org/trunk@33194 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/source/fe/fe_trace.cc b/deal.II/source/fe/fe_trace.cc index af79e48461..4906f0d54e 100644 --- a/deal.II/source/fe/fe_trace.cc +++ b/deal.II/source/fe/fe_trace.cc @@ -16,7 +16,9 @@ #include #include +#include #include +#include #include #include @@ -42,6 +44,10 @@ FE_TraceQ::FE_TraceQ (const unsigned int degree) std::vector renumber (this->dofs_per_face); FETools::hierarchic_to_lexicographic_numbering (degree, renumber); this->poly_space.set_numbering(renumber); + + // Initialize face support points + FE_Q fe_q(degree); + this->unit_face_support_points = fe_q.get_unit_face_support_points(); } @@ -73,7 +79,7 @@ FE_TraceQ::get_name () const template bool FE_TraceQ::has_support_on_face (const unsigned int shape_index, - const unsigned int face_index) const + const unsigned int face_index) const { Assert (shape_index < this->dofs_per_cell, ExcIndexRange (shape_index, 0, this->dofs_per_cell));