From 9fdeb69b004e8bb538d3e6a51b78e1827213b718 Mon Sep 17 00:00:00 2001
From: Philipp Suffa <philipp.suffa@fau.de>
Date: Tue, 2 Apr 2024 16:49:55 +0200
Subject: [PATCH] Small changes for cluster[skip_ci]
---
.../ListLBM/HybridBenchmark/HybridBenchmark.cpp | 9 +++++----
apps/benchmarks/ListLBM/HybridBenchmark/input.py | 2 +-
2 files changed, 6 insertions(+), 5 deletions(-)
diff --git a/apps/benchmarks/ListLBM/HybridBenchmark/HybridBenchmark.cpp b/apps/benchmarks/ListLBM/HybridBenchmark/HybridBenchmark.cpp
index 890b9e44e..d06124bda 100644
--- a/apps/benchmarks/ListLBM/HybridBenchmark/HybridBenchmark.cpp
+++ b/apps/benchmarks/ListLBM/HybridBenchmark/HybridBenchmark.cpp
@@ -220,6 +220,9 @@ void gatherAndPrintPorosityStats(weak_ptr<StructuredBlockForest> forest, BlockDa
std::vector<real_t> processPorositiesVec;
processPorositiesVec = mpi::gather( averageProcessPorosity);
WALBERLA_ROOT_SECTION() {
+ for (auto p : processPorositiesVec) {
+ WALBERLA_LOG_INFO(p)
+ }
real_t sum = std::accumulate(processPorositiesVec.begin(), processPorositiesVec.end(), 0.0);
real_t mean = real_c(sum) / real_c(processPorositiesVec.size());
@@ -398,14 +401,12 @@ int main(int argc, char **argv)
bfc.setWorkloadMemorySUIDAssignmentFunction( meshWorkloadMemory );
auto setupForest = bfc.createSetupBlockForest( cellsPerBlock );
+ WALBERLA_LOG_INFO_ON_ROOT("Number of blocks is <" << setupForest->getXSize() << "," << setupForest->getYSize() << "," << setupForest->getZSize() << ">, without excluded blocks its " << setupForest->getNumberOfBlocks())
+ WALBERLA_LOG_INFO_ON_ROOT("Number of cells per Block is " << cellsPerBlock)
if(writeDomainDecompositionAndReturn) {
- WALBERLA_LOG_INFO_ON_ROOT("Number of blocks is <" << setupForest->getXSize() << "," << setupForest->getYSize() << "," << setupForest->getZSize() << ">, without excluded blocks its " << setupForest->getNumberOfBlocks())
- WALBERLA_LOG_INFO_ON_ROOT("Number of cells per Block is " << cellsPerBlock)
WALBERLA_ROOT_SECTION() { setupForest->writeVTKOutput("SetupBlockForest"); }
return EXIT_SUCCESS;
}
-
-
blocks = bfc.createStructuredBlockForest(cellsPerBlock);
flagFieldId = field::addFlagFieldToStorage< FlagField_T >(blocks, "flag field");
diff --git a/apps/benchmarks/ListLBM/HybridBenchmark/input.py b/apps/benchmarks/ListLBM/HybridBenchmark/input.py
index 78daea054..18fc6382f 100644
--- a/apps/benchmarks/ListLBM/HybridBenchmark/input.py
+++ b/apps/benchmarks/ListLBM/HybridBenchmark/input.py
@@ -119,7 +119,7 @@ def Artery():
scenarios = wlb.ScenarioManager()
#mesh_file = "Artery.obj"
mesh_file = "coronary_colored_medium.obj"
- scenario = Scenario(dx=0.0257, cells_per_block=(128, 128, 128), vtk_write_frequency=0, geometry_setup="artery", mesh_file=mesh_file, timesteps=20, omega=1.7, porosity_switch=1.0, run_hybrid=False, time_step_strategy="Overlap", run_boundaries=True, balance_load=True)
+ scenario = Scenario(dx=0.257, cells_per_block=(32, 32, 32), vtk_write_frequency=0, geometry_setup="artery", mesh_file=mesh_file, timesteps=20, omega=1.7, porosity_switch=1.0, run_hybrid=False, time_step_strategy="Overlap", run_boundaries=True, balance_load=True)
scenarios.add(scenario)
def ArterySparseVsDense():
--
GitLab