From: Daniel Arndt <arndtd@ornl.gov>
Date: Thu, 21 May 2020 16:19:07 +0000 (-0400)
Subject: Fix PolynomialsP
X-Git-Tag: v9.3.0-rc1~1594^2~5
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=aa1f909c02d9476291d53f2faefbba28d9470ce3;p=dealii.git

Fix PolynomialsP
---

diff --git a/include/deal.II/base/polynomials_p.h b/include/deal.II/base/polynomials_p.h
index cc445d5791..5f2dc4b488 100644
--- a/include/deal.II/base/polynomials_p.h
+++ b/include/deal.II/base/polynomials_p.h
@@ -23,6 +23,7 @@
 #include <deal.II/base/point.h>
 #include <deal.II/base/polynomial.h>
 #include <deal.II/base/polynomial_space.h>
+#include <deal.II/base/std_cxx14/memory.h>
 #include <deal.II/base/tensor.h>
 
 #include <vector>
@@ -66,8 +67,8 @@ public:
    * Note, that this number is <tt>PolynomialSpace::degree()-1</tt>, compare
    * definition in PolynomialSpace.
    */
-  unsigned int
-  degree() const;
+  virtual unsigned int
+  degree() const override;
 
   /**
    * For the <tt>n</tt>th polynomial $p_n(x,y,z)=x^i y^j z^k$ this function
@@ -78,6 +79,12 @@ public:
   std::array<unsigned int, dim>
   directional_degrees(unsigned int n) const;
 
+  std::unique_ptr<ScalarPolynomialsBase<dim>>
+  clone() const override
+  {
+    return std_cxx14::make_unique<PolynomialsP<dim>>(*this);
+  }
+
 private:
   /**
    * Fills the <tt>index_map</tt>.
diff --git a/include/deal.II/base/scalar_polynomials_base.h b/include/deal.II/base/scalar_polynomials_base.h
index c039e4b7cb..10a123ef1b 100644
--- a/include/deal.II/base/scalar_polynomials_base.h
+++ b/include/deal.II/base/scalar_polynomials_base.h
@@ -193,7 +193,7 @@ public:
    * class. A derived class may override this if its value is different from
    * @p my_degree.
    */
-  unsigned int
+  virtual unsigned int
   degree() const;
 
   /**