Fix grid alignment of refinement area

1 files changed, 7 insertions, 7 deletions

diff --git a/apps/adrian/poiseuille2d/poiseuille2d.cpp b/apps/adrian/poiseuille2d/poiseuille2d.cpp
index df166e1..b16a19b 100644
--- a/apps/adrian/poiseuille2d/poiseuille2d.cpp
+++ b/apps/adrian/poiseuille2d/poiseuille2d.cpp
@@ -37,9 +37,9 @@ typedef double T;
 
 const T lx  = 8.0;            // length of the channel
 const T ly  = 2.0;            // height of the channel
-const int N = 64;             // resolution of the model
-const T Re  = 500.;           // Reynolds number
-const T baseTau = 0.57;       // Relaxation time of coarsest grid
+const int N = 50;             // resolution of the model
+const T Re  = 200.;           // Reynolds number
+const T baseTau = 0.8;        // Relaxation time of coarsest grid
 const T maxPhysT = 60.;       // max. simulation time in s, SI unit
 const T physInterval = 0.25;  // interval for the convergence check in s
 const T residuum = 1e-5;      // residuum for the convergence check
@@ -188,8 +188,8 @@ int main(int argc, char* argv[])
   prepareGeometry(coarseGrid->getConverter(), coarseGrid->getSuperGeometry());
 
   const Vector<T,2> wantedFineExtend {3.0, 1.5};
-  const Vector<T,2> fineOrigin = coarseGrid->alignToGrid({0.8, (ly-wantedFineExtend[1])/2});
-  const Vector<T,2> fineExtend = coarseGrid->alignToGrid(fineOrigin + wantedFineExtend) - fineOrigin;
+  const Vector<T,2> fineOrigin = coarseGrid->alignLocationToGrid({0.8, (ly-wantedFineExtend[1])/2});
+  const Vector<T,2> fineExtend = coarseGrid->alignExtendToGrid(wantedFineExtend);
 
   auto fineGrid = &coarseGrid->refine(fineOrigin, fineExtend);
   prepareGeometry(fineGrid->getConverter(), fineGrid->getSuperGeometry());
@@ -222,8 +222,8 @@ int main(int argc, char* argv[])
   createLocalBoundaryCondition2D<T, DESCRIPTOR>(fineSBoundaryCondition);
 
   const Vector<T,2> wantedFineExtend2 {0.6, 0.4};
-  const Vector<T,2> fineOrigin2 = fineGrid->alignToGrid({1.05, (ly-wantedFineExtend2[1])/2});
-  const Vector<T,2> fineExtend2 = fineGrid->alignToGrid(fineOrigin2 + wantedFineExtend2) - fineOrigin2;
+  const Vector<T,2> fineOrigin2 = fineGrid->alignLocationToGrid({1.05, (ly-wantedFineExtend2[1])/2});
+  const Vector<T,2> fineExtend2 = fineGrid->alignExtendToGrid(wantedFineExtend2);
 
   auto fineGrid2 = &fineGrid->refine(fineOrigin2, fineExtend2);
   prepareGeometry(fineGrid2->getConverter(), fineGrid2->getSuperGeometry());