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));