Allow read access to data in Interpolated*Function

diff --git a/include/deal.II/base/function_lib.h b/include/deal.II/base/function_lib.h
index e4b9a83c60e125facf64be185efef535a099ccb4..d6f48579e191fca57726313c3971206c543bf4ca 100644 (file)
--- a/include/deal.II/base/function_lib.h
+++ b/include/deal.II/base/function_lib.h
@@ -1441,6 +1441,12 @@ namespace Functions
     virtual std::size_t
     memory_consumption() const override;
 
+    /**
+     * Return a reference to the internally stored data.
+     */
+    const Table<dim, double> &
+    get_data() const;
+
   protected:
     /**
      * Find the index in the table of the rectangle containing an input point
@@ -1552,6 +1558,12 @@ namespace Functions
     virtual std::size_t
     memory_consumption() const override;
 
+    /**
+     * Return a reference to the internally stored data.
+     */
+    const Table<dim, double> &
+    get_data() const;
+
   private:
     /**
      * The set of interval endpoints in each of the coordinate directions.

diff --git a/source/base/function_lib.cc b/source/base/function_lib.cc
index 6f25f9df574c21a8d2c53be07cd445bebcc44fa3..4df376fc24773fa547dd0e65be2b8451096d4749 100644 (file)
--- a/source/base/function_lib.cc
+++ b/source/base/function_lib.cc
@@ -2551,6 +2551,15 @@ namespace Functions
 
 
 
+  template <int dim>
+  const Table<dim, double> &
+  InterpolatedTensorProductGridData<dim>::get_data() const
+  {
+    return data_values;
+  }
+
+
+
   template <int dim>
   double
   InterpolatedTensorProductGridData<dim>::value(
@@ -2745,6 +2754,15 @@ namespace Functions
 
 
 
+  template <int dim>
+  const Table<dim, double> &
+  InterpolatedUniformGridData<dim>::get_data() const
+  {
+    return data_values;
+  }
+
+
+
   /* ---------------------- Polynomial ----------------------- */