Fix, test recursive grid refinement

1 files changed, 38 insertions, 3 deletions

diff --git a/apps/adrian/poiseuille2d/poiseuille2d.cpp b/apps/adrian/poiseuille2d/poiseuille2d.cpp
index 20d2585..a6d7bd1 100644
--- a/apps/adrian/poiseuille2d/poiseuille2d.cpp
+++ b/apps/adrian/poiseuille2d/poiseuille2d.cpp
@@ -161,7 +161,7 @@ void getResults(const std::string& prefix,
     vtmWriter.createMasterFile();
   }
 
-  if (iT%10==0) {
+  if (iT%20==0) {
     vtmWriter.write(iT);
   }
 
@@ -220,6 +220,19 @@ int main(int argc, char* argv[])
   sOnLatticeBoundaryCondition2D<T, DESCRIPTOR> fineSBoundaryCondition(fineGrid->getSuperLattice());
   createLocalBoundaryCondition2D<T, DESCRIPTOR>(fineSBoundaryCondition);
 
+  const Vector<T,2> wantedFineExtend2 {0.4, 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;
+
+  auto fineGrid2 = &fineGrid->refine(fineOrigin2, fineExtend2);
+  prepareGeometry(fineGrid2->getConverter(), fineGrid2->getSuperGeometry());
+
+  const T fine2DeltaX = fineGrid->getConverter().getPhysDeltaX();
+
+  IndicatorCuboid2D<T> refinedOverlap2(fineExtend2 - 4*fine2DeltaX, fineOrigin2 + 2*fine2DeltaX);
+  fineGrid->getSuperGeometry().rename(1,0,refinedOverlap2);
+  fineGrid->getSuperGeometry().rename(2,0,refinedOverlap2);
+
   prepareLattice(
     fineGrid->getConverter(),
     fineGrid->getSuperLattice(),
@@ -227,12 +240,26 @@ int main(int argc, char* argv[])
     fineSBoundaryCondition,
     fineGrid->getSuperGeometry());
 
+  BGKdynamics<T, DESCRIPTOR> fineBulkDynamics2(
+    fineGrid2->getConverter().getLatticeRelaxationFrequency(),
+    instances::getBulkMomenta<T, DESCRIPTOR>());
+
+  sOnLatticeBoundaryCondition2D<T, DESCRIPTOR> fineSBoundaryCondition2(fineGrid2->getSuperLattice());
+  createLocalBoundaryCondition2D<T, DESCRIPTOR>(fineSBoundaryCondition2);
+
+  prepareLattice(
+    fineGrid2->getConverter(),
+    fineGrid2->getSuperLattice(),
+    fineBulkDynamics2,
+    fineSBoundaryCondition2,
+    fineGrid2->getSuperGeometry());
+
   clout << "starting simulation..." << endl;
   Timer<T> timer(
     coarseGrid->getConverter().getLatticeTime(maxPhysT),
     coarseGrid->getSuperGeometry().getStatistics().getNvoxel());
   util::ValueTracer<T> converge(
-    fineGrid->getConverter().getLatticeTime(physInterval),
+    fineGrid2->getConverter().getLatticeTime(physInterval),
     residuum);
   timer.start();
 
@@ -260,8 +287,16 @@ int main(int argc, char* argv[])
       fineGrid->getSuperGeometry(),
       timer,
       converge.hasConverged());
+    getResults(
+      "fine2_",
+      fineGrid2->getSuperLattice(),
+      fineGrid2->getConverter(),
+      iT,
+      fineGrid2->getSuperGeometry(),
+      timer,
+      converge.hasConverged());
 
-    converge.takeValue(fineGrid->getSuperLattice().getStatistics().getAverageEnergy(), true);
+    converge.takeValue(fineGrid2->getSuperLattice().getStatistics().getAverageEnergy(), true);
   }
 
   timer.stop();