From 03ab577adb93e91c8a1868cbb7dfbfccc81c61c6 Mon Sep 17 00:00:00 2001
From: Philipp Suffa <philipp.suffa@fau.de>
Date: Fri, 28 Apr 2023 15:38:13 +0200
Subject: [PATCH] Setup empty channel Benchmark

---
 .../ListLBM/HybridBenchmark/HybridBenchmark.cpp          | 6 ++++--
 apps/benchmarks/ListLBM/HybridBenchmark/input.py         | 9 ++++++++-
 2 files changed, 12 insertions(+), 3 deletions(-)

diff --git a/apps/benchmarks/ListLBM/HybridBenchmark/HybridBenchmark.cpp b/apps/benchmarks/ListLBM/HybridBenchmark/HybridBenchmark.cpp
index c0283918c..3fd0b1cb5 100644
--- a/apps/benchmarks/ListLBM/HybridBenchmark/HybridBenchmark.cpp
+++ b/apps/benchmarks/ListLBM/HybridBenchmark/HybridBenchmark.cpp
@@ -535,7 +535,6 @@ int main(int argc, char **argv)
             vtk::writeDomainDecomposition(blocks, "domain_decomposition", "vtk_out", "write_call", true, true, 0);
          }
 
-         lbm::PerformanceEvaluation< FlagField_T > performance(blocks, flagFieldId, fluidFlagUID);
 
          WALBERLA_LOG_INFO_ON_ROOT("Simulating ListLBM:"
                                    "\n timesteps:                  "
@@ -570,6 +569,8 @@ int main(int argc, char **argv)
             WALBERLA_LOG_INFO_ON_ROOT("Simulation finished")
             real_t time = simTimer.max();
             WALBERLA_MPI_SECTION() { walberla::mpi::reduceInplace(time, walberla::mpi::MAX); }
+
+            lbm::PerformanceEvaluation< FlagField_T > performance(blocks, flagFieldId, fluidFlagUID);
             performance.logResultOnRoot(timesteps, time);
 
             const auto reducedTimeloopTiming = timeloopTiming.getReduced();
@@ -753,7 +754,6 @@ int main(int argc, char **argv)
 
          }
 
-         lbm::PerformanceEvaluation< FlagField_T > performance(blocks, flagFieldId, fluidFlagUID);
 
          WALBERLA_LOG_INFO_ON_ROOT("Simulating ListLBM:"
                                       "\n timesteps:                  " << timesteps
@@ -788,6 +788,8 @@ int main(int argc, char **argv)
             WALBERLA_LOG_INFO_ON_ROOT("Simulation finished")
             real_t time = simTimer.max();
             WALBERLA_MPI_SECTION() { walberla::mpi::reduceInplace(time, walberla::mpi::MAX); }
+
+            lbm::PerformanceEvaluation< FlagField_T > performance(blocks, flagFieldId, fluidFlagUID);
             performance.logResultOnRoot(timesteps, time);
 
             const auto reducedTimeloopTiming = timeloopTiming.getReduced();
diff --git a/apps/benchmarks/ListLBM/HybridBenchmark/input.py b/apps/benchmarks/ListLBM/HybridBenchmark/input.py
index 935faae91..ca4678b7f 100644
--- a/apps/benchmarks/ListLBM/HybridBenchmark/input.py
+++ b/apps/benchmarks/ListLBM/HybridBenchmark/input.py
@@ -94,6 +94,12 @@ def porosity_benchmark():
         scenario = Scenario(porosity=porosity, geometry_setup="randomNoslip", inflow_velocity=(0,0,0), run_boundaries=False)
         scenarios.add(scenario)
 
+
+def emptyChannel():
+    scenarios = wlb.ScenarioManager()
+    scenario = Scenario(porosity=1.0, vtk_write_frequency=0, geometry_setup="randomNoslip", cells_per_block=(64, 64, 64), time_step_strategy="kernelOnly", run_boundaries=False, porositySwitch=1.1)
+    scenarios.add(scenario)
+
 def randomNoslip():
     scenarios = wlb.ScenarioManager()
     scenario = Scenario(porosity=0.9, vtk_write_frequency=50, geometry_setup="randomNoslip", inflow_velocity=(0,0,0))
@@ -122,5 +128,6 @@ def particleBed():
 
 #randomNoslip()
 #spheres()
-Artery()
+#Artery()
 #particleBed()
+emptyChannel()
\ No newline at end of file
-- 
GitLab