From: Martin Kronbichler <kronbichler@lnm.mw.tum.de>
Date: Sat, 15 May 2021 08:18:34 +0000 (+0200)
Subject: Let TensorProductPolynomials return 1D polynomial space
X-Git-Tag: v9.3.0-rc1~69^2~2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d2abc8f04615400b03ceb67f4716f4d2a921b4d0;p=dealii.git

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 da0ce0e73a..a8666cc80f 100644
--- 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 da6ecac666..8bb03217cf 100644
--- 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 -------------- */