From 3bc4210c9a3c10ffe61f0237bba101f9700b7668 Mon Sep 17 00:00:00 2001
From: Adrian Kummerlaender
Date: Thu, 31 Jan 2019 15:14:03 +0100
Subject: Report total number of active cells in refinement hierarchy

---
 src/refinement/grid2D.h  |  2 ++
 src/refinement/grid2D.hh | 10 ++++++++++
 2 files changed, 12 insertions(+)

(limited to 'src')

diff --git a/src/refinement/grid2D.h b/src/refinement/grid2D.h
index 37e764a..95324f7 100644
--- a/src/refinement/grid2D.h
+++ b/src/refinement/grid2D.h
@@ -110,6 +110,8 @@ public:
 
   RefiningGrid2D<T,DESCRIPTOR>& refine(Vector<T,2> origin, Vector<T,2> extend, bool addCouplers=true);
 
+  std::size_t getActiveVoxelN() const;
+
 };
 
 template <typename T, template<typename> class DESCRIPTOR>
diff --git a/src/refinement/grid2D.hh b/src/refinement/grid2D.hh
index 941bff7..504808a 100644
--- a/src/refinement/grid2D.hh
+++ b/src/refinement/grid2D.hh
@@ -217,6 +217,16 @@ Vector<T,2> Grid2D<T,DESCRIPTOR>::alignExtendToGrid(Vector<T,2> extend) const
   return util::floor(extend / deltaX) * deltaX;
 }
 
+template <typename T, template<typename> class DESCRIPTOR>
+std::size_t Grid2D<T,DESCRIPTOR>::getActiveVoxelN() const
+{
+  std::size_t n = _geometry->getStatistics().getNvoxel();
+  for (const auto& grid : _fineGrids) {
+    n += grid->getActiveVoxelN();
+  }
+  return n;
+}
+
 template <typename T, template<typename> class DESCRIPTOR>
 RefiningGrid2D<T,DESCRIPTOR>& Grid2D<T,DESCRIPTOR>::refine(
   Vector<T,2> wantedOrigin, Vector<T,2> wantedExtend, bool addCouplers)
-- 
cgit v1.2.3