]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Convert FE_TraceQ to Gauss-Lobatto points for collaboration with FE_Q.
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Fri, 8 Apr 2016 15:20:24 +0000 (17:20 +0200)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Sat, 9 Apr 2016 13:49:09 +0000 (15:49 +0200)
source/fe/fe_trace.cc

index 87ed272a22b12abc0689a7fa8bab5326e56f3f1a..1e0509ee291db0ee38f5f74f3539cd8d04ec5a30 100644 (file)
@@ -35,7 +35,7 @@ template <int dim, int spacedim>
 FE_TraceQ<dim,spacedim>::FE_TraceQ (const unsigned int degree)
   :
   FE_PolyFace<TensorProductPolynomials<dim-1>, dim, spacedim> (
-    TensorProductPolynomials<dim-1>(Polynomials::LagrangeEquidistant::generate_complete_basis(degree)),
+    TensorProductPolynomials<dim-1>(Polynomials::generate_complete_Lagrange_basis(QGaussLobatto<1>(degree+1).get_points())),
     FiniteElementData<dim>(get_dpo_vector(degree), 1, degree, FiniteElementData<dim>::L2),
     std::vector<bool>(1,true)),
   fe_q (degree)

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.