Add support for calling a function for each subgrid of Grid2D

1 files changed, 13 insertions, 0 deletions

diff --git a/src/refinement/grid2D.hh b/src/refinement/grid2D.hh
index 6e9405e..6d16648 100644
--- a/src/refinement/grid2D.hh
+++ b/src/refinement/grid2D.hh
@@ -259,6 +259,19 @@ Vector<T,2> Grid2D<T,DESCRIPTOR>::alignExtendToGrid(Vector<T,2> extend) const
 }
 
 template <typename T, template<typename> class DESCRIPTOR>
+void Grid2D<T,DESCRIPTOR>::forEachGrid(
+  const std::string& id,
+  std::function<void(Grid2D<T,DESCRIPTOR>&,const std::string&)>&& f)
+{
+  f(*this, id);
+  for (std::size_t i = 0; i < _fineGrids.size(); ++i) {
+    _fineGrids[i]->forEachGrid(id + "_" + std::to_string(i),
+                               std::forward<decltype(f)>(f));
+
+  }
+}
+
+template <typename T, template<typename> class DESCRIPTOR>
 std::size_t Grid2D<T,DESCRIPTOR>::getActiveVoxelN() const
 {
   std::size_t n = _geometry->getStatistics().getNvoxel();