From: David Wells <drwells@email.unc.edu>
Date: Fri, 15 Jan 2021 00:21:41 +0000 (-0500)
Subject: Switch ScalarPolynomial<dim> -> BarycentricPolynomials<dim>
X-Git-Tag: v9.3.0-rc1~587^2~1
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d031fef8c62400ed46079ce79c2a6799f9137a4e;p=dealii.git

Switch ScalarPolynomial<dim> -> BarycentricPolynomials<dim>
---

diff --git a/include/deal.II/simplex/polynomials.h b/include/deal.II/simplex/polynomials.h
index b3cd474e95..31acf50362 100644
--- a/include/deal.II/simplex/polynomials.h
+++ b/include/deal.II/simplex/polynomials.h
@@ -21,6 +21,8 @@
 
 #include <deal.II/base/scalar_polynomials_base.h>
 
+#include <deal.II/simplex/barycentric_polynomials.h>
+
 DEAL_II_NAMESPACE_OPEN
 
 /**
@@ -264,14 +266,14 @@ namespace Simplex
 
   private:
     /**
-     * Scalar polynomial defined on a triangle.
+     * Scalar polynomials defined on a triangle.
      */
-    const ScalarPolynomial<2> poly_tri;
+    const BarycentricPolynomials<2> poly_tri;
 
     /**
-     * Scalar polynomial defined on a line.
+     * Scalar polynomials defined on a line.
      */
-    const ScalarPolynomial<1> poly_line;
+    const BarycentricPolynomials<1> poly_line;
   };
 
 
diff --git a/source/base/qprojector.cc b/source/base/qprojector.cc
index 85a997f514..858a03c8f6 100644
--- a/source/base/qprojector.cc
+++ b/source/base/qprojector.cc
@@ -18,7 +18,7 @@
 #include <deal.II/base/qprojector.h>
 #include <deal.II/base/tensor_product_polynomials.h>
 
-#include <deal.II/simplex/polynomials.h>
+#include <deal.II/simplex/barycentric_polynomials.h>
 
 DEAL_II_NAMESPACE_OPEN
 
@@ -529,7 +529,7 @@ QProjector<2>::project_to_all_faces(
 
       // linear polynomial to map the reference quadrature points correctly
       // on faces
-      const Simplex::ScalarPolynomial<1> poly(1);
+      const auto poly = Simplex::BarycentricPolynomials<1>::get_fe_p_basis(1);
 
       // new (projected) quadrature points and weights
       std::vector<Point<2>> points;
@@ -668,8 +668,8 @@ QProjector<3>::project_to_all_faces(
     std::vector<Point<3>> points;
     std::vector<double>   weights;
 
-    const Simplex::ScalarPolynomial<2> poly_tri(1);
-    const TensorProductPolynomials<2>  poly_quad(
+    const auto poly_tri = Simplex::BarycentricPolynomials<2>::get_fe_p_basis(1);
+    const TensorProductPolynomials<2> poly_quad(
       Polynomials::generate_complete_Lagrange_basis(
         {Point<1>(0.0), Point<1>(1.0)}));
 
diff --git a/source/simplex/fe_lib.cc b/source/simplex/fe_lib.cc
index d06188f983..61c14b5a43 100644
--- a/source/simplex/fe_lib.cc
+++ b/source/simplex/fe_lib.cc
@@ -283,7 +283,7 @@ namespace Simplex
     const std::vector<unsigned int> &                 dpo_vector,
     const typename FiniteElementData<dim>::Conformity conformity)
     : dealii::FE_Poly<dim, spacedim>(
-        Simplex::ScalarPolynomial<dim>(degree),
+        BarycentricPolynomials<dim>::get_fe_p_basis(degree),
         FiniteElementData<dim>(dpo_vector,
                                dim == 2 ? ReferenceCell::Type::Tri :
                                           ReferenceCell::Type::Tet,
diff --git a/source/simplex/polynomials.cc b/source/simplex/polynomials.cc
index 4d8dff665a..58910a9a87 100644
--- a/source/simplex/polynomials.cc
+++ b/source/simplex/polynomials.cc
@@ -14,6 +14,7 @@
 // ---------------------------------------------------------------------
 
 
+#include <deal.II/simplex/barycentric_polynomials.h>
 #include <deal.II/simplex/polynomials.h>
 
 DEAL_II_NAMESPACE_OPEN
@@ -636,8 +637,8 @@ namespace Simplex
   ScalarWedgePolynomial<dim>::ScalarWedgePolynomial(const unsigned int degree)
     : ScalarPolynomialsBase<dim>(degree,
                                  compute_n_polynomials_wedge(dim, degree))
-    , poly_tri(degree)
-    , poly_line(degree)
+    , poly_tri(BarycentricPolynomials<2>::get_fe_p_basis(degree))
+    , poly_line(BarycentricPolynomials<1>::get_fe_p_basis(degree))
   {}
 
 
diff --git a/tests/simplex/polynomials_01.cc b/tests/simplex/polynomials_01.cc
index eb8c82cfcc..c501094a14 100644
--- a/tests/simplex/polynomials_01.cc
+++ b/tests/simplex/polynomials_01.cc
@@ -14,8 +14,8 @@
 // ---------------------------------------------------------------------
 
 
-// Test Simplex::ScalarPolynomial on an the points of an arbitrary quadrature
-// rule.
+// Test Simplex::BarycentricPolynomials on an the points of an arbitrary
+// quadrature rule.
 
 
 #include <deal.II/base/quadrature_lib.h>
@@ -30,8 +30,9 @@ template <int dim>
 void
 test(const unsigned int degree)
 {
-  Simplex::ScalarPolynomial<dim> poly(degree);
-  QSimplex<dim>                  quad(QGauss<dim>(degree + 1));
+  const auto poly =
+    Simplex::BarycentricPolynomials<dim>::get_fe_p_basis(degree);
+  QSimplex<dim> quad(QGauss<dim>(degree + 1));
 
   std::vector<double>         values(poly.n());
   std::vector<Tensor<1, dim>> grads(poly.n());