[spaces.bochner] allow to evaluate out of bounds

dune/gdt/discretefunction/bochner.hh

@@ -12,6 +12,7 @@
 #define DUNE_GDT_DISCRETEFUNCTION_BOCHNER_HH
 
 #include <dune/xt/common/memory.hh>
+#include <dune/xt/common/math.hh>
 #include <dune/xt/la/container/vector-array/list.hh>
 #include <dune/xt/grid/search.hh>
 
@@ -78,14 +79,16 @@ public:
 
   DiscreteFunction<V, GV, r, rC, R> evaluate(const double& time) const
   {
-    const auto search_result = XT::Grid::make_entity_in_level_search(bochner_space_.temporal_space().grid_view())(
-        std::vector<double>(1, time));
+    const double t =
+        XT::Common::clamp(time, bochner_space_.time_interval().first, bochner_space_.time_interval().second);
+    const auto search_result =
+        XT::Grid::make_entity_in_level_search(bochner_space_.temporal_space().grid_view())(std::vector<double>(1, t));
     DUNE_THROW_IF(!search_result.at(0),
                   Exceptions::finite_element_error,
-                  "when evaluating timedependent function: could not find time_interval for time = " << time);
+                  "when evaluating timedependent function: could not find time_interval for time = " << t);
     const auto& time_interval = *search_result.at(0);
     const auto temporal_basis = bochner_space_.temporal_space().basis().localize(time_interval);
-    const auto time_in_reference_element = time_interval.geometry().local(time);
+    const auto time_in_reference_element = time_interval.geometry().local(t);
     const auto temporal_basis_values = temporal_basis->evaluate_set(time_in_reference_element);
     V result(bochner_space_.spatial_space().mapper().size(), 0.);
     const auto global_dof_indices = bochner_space_.temporal_space().mapper().global_indices(time_interval);

dune/gdt/spaces/bochner.hh

@@ -12,7 +12,9 @@
 #define DUNE_GDT_SPACES_BOCHNER_HH
 
 #include <dune/grid/onedgrid.hh>
+#include <dune/grid/common/rangegenerators.hh>
 
+#include <dune/xt/common/ranges.hh>
 #include <dune/gdt/spaces/h1/continuous-lagrange.hh>
 
 #include "interface.hh"
@@ -30,7 +32,16 @@ public:
     : spatial_space_(spatial_space)
     , temporal_grid_(std::forward<TemporalGridArgs>(temporal_grid_args)...)
     , temporal_space_(temporal_grid_.leafGridView(), /*order=*/1)
-  {}
+    , time_interval_(std::make_pair(std::numeric_limits<double>::max(), std::numeric_limits<double>::min()))
+  {
+    for (auto&& time_interval : elements(temporal_space_.grid_view())) {
+      for (auto ii : XT::Common::value_range(time_interval.geometry().corners())) {
+        const auto time_point = time_interval.geometry().corner(ii);
+        time_interval_.first = std::min(time_interval_.first, time_point[0]);
+        time_interval_.second = std::max(time_interval_.second, time_point[0]);
+      }
+    }
+  }
 
   const ContinuousLagrangeSpace<typename OneDGrid::LeafGridView>& temporal_space() const
   {
@@ -42,10 +53,16 @@ public:
     return spatial_space_;
   }
 
+  const std::pair<double, double>& time_interval() const
+  {
+    return time_interval_;
+  }
+
 private:
   const SpaceInterface<GV, r, rC, R>& spatial_space_;
   const OneDGrid temporal_grid_;
   const ContinuousLagrangeSpace<typename OneDGrid::LeafGridView> temporal_space_;
+  std::pair<double, double> time_interval_;
 }; // class BochnerSpace