summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAdrian Kummerlaender2019-01-31 14:26:35 +0100
committerAdrian Kummerlaender2019-06-24 15:17:09 +0200
commit39e22863be1e50bd8e6b36f3d3ecbd1bec09e6a0 (patch)
tree5a304e6a3c4603eeb3b057773f0f15cedde999b0
parent5e888fc13f38c94777963d72fc9e391cd4fa477a (diff)
downloadgrid_refinement_openlb-39e22863be1e50bd8e6b36f3d3ecbd1bec09e6a0.tar
grid_refinement_openlb-39e22863be1e50bd8e6b36f3d3ecbd1bec09e6a0.tar.gz
grid_refinement_openlb-39e22863be1e50bd8e6b36f3d3ecbd1bec09e6a0.tar.bz2
grid_refinement_openlb-39e22863be1e50bd8e6b36f3d3ecbd1bec09e6a0.tar.xz
grid_refinement_openlb-39e22863be1e50bd8e6b36f3d3ecbd1bec09e6a0.zip
Move geometry print to prepareLattice
This way statistics reflect the cells disabled by refinement.
-rw-r--r--apps/adrian/cylinder2d/cylinder2d.cpp12
1 files changed, 6 insertions, 6 deletions
diff --git a/apps/adrian/cylinder2d/cylinder2d.cpp b/apps/adrian/cylinder2d/cylinder2d.cpp
index 0893928..56866c4 100644
--- a/apps/adrian/cylinder2d/cylinder2d.cpp
+++ b/apps/adrian/cylinder2d/cylinder2d.cpp
@@ -93,7 +93,6 @@ void prepareGeometry(Grid2D<T,DESCRIPTOR>& grid)
sGeometry.clean();
sGeometry.innerClean();
sGeometry.checkForErrors();
- sGeometry.print();
clout << "Prepare Geometry ... OK" << std::endl;
}
@@ -149,6 +148,7 @@ void prepareLattice(Grid2D<T,DESCRIPTOR>& grid)
sLattice.initialize();
clout << "Prepare lattice ... OK" << std::endl;
+ sGeometry.print();
}
void setBoundaryValues(Grid2D<T,DESCRIPTOR>& grid, int iT)
@@ -227,10 +227,12 @@ int main(int argc, char* argv[])
auto& fineGrid = coarseGrid.refine(fineOrigin, fineExtend);
prepareGeometry(fineGrid);
- const Vector<T,2> fineExtendB {10*coarseDeltaX, ly};
- const Vector<T,2> fineOriginB {lx-10*coarseDeltaX, 0};
+ disableRefinedArea(coarseGrid, fineGrid);
+
+ const Vector<T,2> fineOutflowExtend {5*coarseDeltaX, ly};
+ const Vector<T,2> fineOutflowOrigin {lx-5*coarseDeltaX, 0};
- auto& fineOutflowGrid = coarseGrid.refine(fineOriginB, fineExtendB, false);
+ auto& fineOutflowGrid = coarseGrid.refine(fineOutflowOrigin, fineOutflowExtend, false);
prepareGeometry(fineOutflowGrid);
{
@@ -249,8 +251,6 @@ int main(int argc, char* argv[])
coarseGrid.getSuperGeometry().rename(4,0,refined);
}
- disableRefinedArea(coarseGrid, fineGrid);
-
const Vector<T,2> fineExtend2 {4*cylinderD, 2*cylinderD};
const Vector<T,2> fineOrigin2 {1*cylinderD, 2*cylinderD-fineExtend2[1]/2};