Report total number of active cells in refinement hierarchy

2 files changed, 4 insertions, 2 deletions

diff --git a/apps/adrian/cylinder2d/cylinder2d.cpp b/apps/adrian/cylinder2d/cylinder2d.cpp
index 56866c4..a0732e4 100644
--- a/apps/adrian/cylinder2d/cylinder2d.cpp
+++ b/apps/adrian/cylinder2d/cylinder2d.cpp
@@ -251,7 +251,7 @@ int main(int argc, char* argv[])
     coarseGrid.getSuperGeometry().rename(4,0,refined);
   }
 
-  const Vector<T,2> fineExtend2 {4*cylinderD, 2*cylinderD};
+  const Vector<T,2> fineExtend2 {3.6*cylinderD, 2*cylinderD};
   const Vector<T,2> fineOrigin2 {1*cylinderD, 2*cylinderD-fineExtend2[1]/2};
 
   auto& fineGrid2 = fineGrid.refine(fineOrigin2, fineExtend2);
@@ -264,12 +264,13 @@ int main(int argc, char* argv[])
   prepareLattice(fineOutflowGrid);
   prepareLattice(fineGrid2);
 
+  clout << "Total number of active cells: " << coarseGrid.getActiveVoxelN() << endl;
   clout << "Starting simulation..." << endl;
   Timer<T> timer(
     coarseGrid.getConverter().getLatticeTime(maxPhysT),
     coarseGrid.getSuperGeometry().getStatistics().getNvoxel());
   timer.start();
-  const int statIter = coarseGrid.getConverter().getLatticeTime(0.5);
+  const int statIter = coarseGrid.getConverter().getLatticeTime(1);
 
   for (int iT = 0; iT < coarseGrid.getConverter().getLatticeTime(maxPhysT); ++iT) {
     setBoundaryValues(coarseGrid, iT);
diff --git a/apps/adrian/poiseuille2d/poiseuille2d.cpp b/apps/adrian/poiseuille2d/poiseuille2d.cpp
index 4091905..5a47f96 100644
--- a/apps/adrian/poiseuille2d/poiseuille2d.cpp
+++ b/apps/adrian/poiseuille2d/poiseuille2d.cpp
@@ -248,6 +248,7 @@ int main(int argc, char* argv[])
 
   prepareLattice(fineGrid);
 
+  clout << "Total number of active cells: " << coarseGrid.getActiveVoxelN() << endl;
   clout << "starting simulation..." << endl;
   Timer<T> timer(
     coarseGrid.getConverter().getLatticeTime(maxPhysT),