Make offset visible with an impulsion

apps/showcases/LeesEdwards/LeesEdwards.cpp

@@ -161,6 +161,18 @@ int main(int argc, char** argv)
    for (auto& block : *blocks)
       setterSweep(&block);
 
+   // create impulsion
+   {
+      auto &block = *(blocks->begin());
+      auto pdf_field = block.template getData<PdfField_T>(pdfFieldID);
+      auto const dim_x = blocks->getNumberOfXCells(block);
+      auto const dim_y = blocks->getNumberOfYCells(block);
+      Cell const cell1{dim_x / 2, dim_y - 4, 0u};
+      Cell const cell2{dim_x / 2, dim_y - 5, 0u};
+      pdf_field->get(cell1, cell_idx_t(2)) = real_t(1e8);
+      pdf_field->get(cell2, cell_idx_t(2)) = real_t(1e8);
+   }
+
    // create time loop
    SweepTimeloop timeloop(blocks->getBlockStorage(), timesteps);
 

apps/showcases/LeesEdwards/LeesEdwards.prm

 
 Parameters 
 {
-	timesteps       10001;
-	vtkWriteFrequency 500;
+	timesteps       51;
+	vtkWriteFrequency 1;
 	remainingTimeLoggerFrequency 3; // in seconds
-	offset 10;
+	offset 16;
 }
 
 DomainSetup