]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Compute unit face support points of FE_TraceQ
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Thu, 17 Jul 2014 21:36:31 +0000 (21:36 +0000)
committerMatthias Maier <tamiko@kyomu.43-1.org>
Sun, 20 Jul 2014 22:16:56 +0000 (00:16 +0200)
git-svn-id: https://svn.dealii.org/trunk@33194 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/source/fe/fe_trace.cc

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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.