Implement BoxConstrainedQuadraticDirectionalRestriction::operator()

dune/tnnmg/functionals/boxconstrainedquadraticfunctional.hh

@@ -212,7 +212,9 @@ public:
 
   Range operator()(const Vector& v) const
   {
-    DUNE_THROW(Dune::NotImplemented, "Evaluation of QuadraticCoordinateRestriction not implemented");
+    if (checkInsideBox(v, lowerObstacle(), upperObstacle()))
+      return Base::operator()(v);
+    return std::numeric_limits<Range>::max();
   }
 
   const LowerObstacle& lowerObstacle() const