Let TensorProductPolynomials return 1D polynomial space

diff --git a/include/deal.II/base/tensor_product_polynomials.h b/include/deal.II/base/tensor_product_polynomials.h
index da0ce0e73a44deaded44fa80cb2a9bfff89d7a45..a8666cc80f4f359591d8cd0a8b8a63075cd275db 100644 (file)
--- a/include/deal.II/base/tensor_product_polynomials.h
+++ b/include/deal.II/base/tensor_product_polynomials.h
@@ -84,9 +84,9 @@ public:
 
   /**
    * Constructor. <tt>pols</tt> is a vector of objects that should be derived
-   * or otherwise convertible to one-dimensional polynomial objects of type @p
-   * PolynomialType (template argument of class). It will be copied element by
-   * element into a private variable.
+   * or otherwise convertible to one-dimensional polynomial objects of type
+   * `PolynomialType` (template argument of class). It will be copied element
+   * by element into a protected member variable.
    */
   template <class Pol>
   TensorProductPolynomials(const std::vector<Pol> &pols);
@@ -245,6 +245,13 @@ public:
   virtual std::size_t
   memory_consumption() const override;
 
+  /**
+   * Return a copy of the underlying one-dimensional polynomials given to the
+   * constructor of this class.
+   */
+  std::vector<PolynomialType>
+  get_underlying_polynomials() const;
+
 protected:
   /**
    * Copy of the vector <tt>pols</tt> of polynomials given to the constructor.

diff --git a/source/base/tensor_product_polynomials.cc b/source/base/tensor_product_polynomials.cc
index da6ecac6668de94548969053b758587de5a3a8e7..8bb03217cf84bd02f0b239238ceee6c8f9796bd0 100644 (file)
--- a/source/base/tensor_product_polynomials.cc
+++ b/source/base/tensor_product_polynomials.cc
@@ -504,6 +504,16 @@ TensorProductPolynomials<dim, PolynomialType>::memory_consumption() const
 
 
 
+template <int dim, typename PolynomialType>
+std::vector<PolynomialType>
+TensorProductPolynomials<dim, PolynomialType>::get_underlying_polynomials()
+  const
+{
+  return polynomials;
+}
+
+
+
 /* ------------------- AnisotropicPolynomials -------------- */