From: Martin Kronbichler <kronbichler@lnm.mw.tum.de> Date: Thu, 17 Jul 2014 21:36:31 +0000 (+0000) Subject: Compute unit face support points of FE_TraceQ X-Git-Tag: v8.2.0-rc1~280 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=b570b9f2ccd2e08a8b295638c65955d4492572ef;p=dealii.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 <deal.II/base/config.h> #include <deal.II/base/tensor_product_polynomials.h> +#include <deal.II/base/quadrature_lib.h> #include <deal.II/fe/fe_poly_face.h> +#include <deal.II/fe/fe_q.h> #include <deal.II/fe/fe_tools.h> #include <deal.II/fe/fe_trace.h> @@ -42,6 +44,10 @@ FE_TraceQ<dim,spacedim>::FE_TraceQ (const unsigned int degree) std::vector<unsigned int> renumber (this->dofs_per_face); FETools::hierarchic_to_lexicographic_numbering<dim-1> (degree, renumber); this->poly_space.set_numbering(renumber); + + // Initialize face support points + FE_Q<dim,spacedim> fe_q(degree); + this->unit_face_support_points = fe_q.get_unit_face_support_points(); } @@ -73,7 +79,7 @@ FE_TraceQ<dim,spacedim>::get_name () const template <int dim, int spacedim> bool FE_TraceQ<dim,spacedim>::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));