localfunction of discsretefunction could not be stored without the discretefunction object

src/amdis/gridfunctions/DiscreteFunction.inc.hpp

@@ -23,9 +23,9 @@ private:
 
 public:
   LocalFunction(DiscreteFunction const& globalFunction)
-    : globalFunction_(&globalFunction)
-    , localView_(globalFunction_->basis().localView())
-    , subTree_(&child(localView_.tree(), globalFunction_->treePath()))
+    : globalFunction_(globalFunction)
+    , localView_(globalFunction_.basis().localView())
+    , subTree_(&child(localView_.tree(), globalFunction_.treePath()))
   {}
 
   /// \brief Bind the LocalView to the element
@@ -48,7 +48,7 @@ public:
   /// \brief Create a LocalFunction representing the gradient. \relates GradientLocalFunction
   friend GradientLocalFunction derivative(LocalFunction const& localFunction)
   {
-    return GradientLocalFunction{*localFunction.globalFunction_};
+    return GradientLocalFunction{localFunction.globalFunction_};
   }
 
   /// \brief The \ref polynomialDegree() of the LocalFunctions
@@ -66,7 +66,7 @@ public:
   }
 
 private:
-  DiscreteFunction const* globalFunction_;
+  DiscreteFunction globalFunction_;
   LocalView localView_;
   SubTree const* subTree_;
   bool bound_ = false;
@@ -94,9 +94,9 @@ private:
 
 public:
   GradientLocalFunction(DiscreteFunction const& globalFunction)
-    : globalFunction_(&globalFunction)
-    , localView_(globalFunction_->basis().localView())
-    , subTree_(&child(localView_.tree(), globalFunction_->treePath()))
+    : globalFunction_(globalFunction)
+    , localView_(globalFunction_.basis().localView())
+    , subTree_(&child(localView_.tree(), globalFunction_.treePath()))
   {}
 
   void bind(Element const& element)
@@ -130,7 +130,7 @@ public:
   }
 
 private:
-  DiscreteFunction const* globalFunction_;
+  DiscreteFunction globalFunction_;
   LocalView localView_;
   SubTree const* subTree_;
   Dune::Std::optional<Geometry> geometry_;
@@ -145,8 +145,8 @@ LocalFunction::operator()(Domain const& x) const
   assert( bound_ );
   auto y = Range(0);
 
-  auto&& coefficients = *globalFunction_->dofVector_;
-  auto&& nodeToRangeEntry = globalFunction_->nodeToRangeEntry_;
+  auto&& coefficients = *globalFunction_.dofVector_;
+  auto&& nodeToRangeEntry = globalFunction_.nodeToRangeEntry_;
   forEachLeafNode_(*subTree_, [&,this](auto const& node, auto const& tp)
   {
     auto&& fe = node.finiteElement();
@@ -191,8 +191,8 @@ GradientLocalFunction::operator()(Domain const& x) const
   for (std::size_t j = 0; j < dy.size(); ++j)
     dy[j] = 0;
 
-  auto&& coefficients = *globalFunction_->dofVector_;
-  auto&& nodeToRangeEntry = globalFunction_->nodeToRangeEntry_;
+  auto&& coefficients = *globalFunction_.dofVector_;
+  auto&& nodeToRangeEntry = globalFunction_.nodeToRangeEntry_;
   forEachLeafNode_(*subTree_, [&,this](auto const& node, auto const& tp)
   {
     // TODO: may DOFVectorView::Range to FieldVector type if necessary