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 af79e48461b20783b729e7276fb054292f4a6910..4906f0d54ef49f67ea5b7f929f210cf5fa1713ca 100644 (file)
--- 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));