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 --- apps/adrian/cylinder2d/cylinder2d.cpp | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) (limited to 'apps/adrian/cylinder2d/cylinder2d.cpp') 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 fineExtend2 {4*cylinderD, 2*cylinderD}; + const Vector fineExtend2 {3.6*cylinderD, 2*cylinderD}; const Vector 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 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); -- cgit v1.2.3