diff --git a/apps/benchmarks/ListLBM/HybridBenchmark/HybridBenchmark.cpp b/apps/benchmarks/ListLBM/HybridBenchmark/HybridBenchmark.cpp index c0283918c74a47c87baab7611f9f8b591bf617c8..3fd0b1cb5e0df7068e1ea973a58e2698e6d40f67 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 935faae91cf6f9f318b763b1c8e3f1db26060dac..ca4678b7fcbaa01d83f2a11ee0c6670942f5d940 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