WIP: Adapt UBB to pull scheme

Author: RudolfWeeber

maintainer/walberla_kernels/generate_lb_kernels.py

@@ -193,7 +193,7 @@ def paramlist(parameters, keys):
         lbmpy_walberla.generate_boundary(
             ctx, f"Dynamic_UBB_{precision_suffix}{target_suffix}", ubb_dynamic,
             method, additional_data_handler=ubb_data_handler,
-            streaming_pattern="push", target=target)
+            streaming_pattern="pull", target=target)
 
         with open(f"Dynamic_UBB_{precision_suffix}{target_suffix}.h", "r+") as f:
             content = f.read()

src/walberla_bridge/src/electrokinetics/generated_kernels/AdvectiveFluxKernel_double_precision.cpp

@@ -1625,9 +1625,9 @@ static FUNC_PREFIX void advectivefluxkernel_double_precision_advectivefluxkernel
 } // namespace internal_5255e1c780a944d646f270232511968b
 
 void AdvectiveFluxKernel_double_precision::run(IBlock *block) {
+  auto j = block->getData<field::GhostLayerField<double, 13>>(jID);
   auto rho = block->getData<field::GhostLayerField<double, 1>>(rhoID);
   auto u = block->getData<field::GhostLayerField<double, 3>>(uID);
-  auto j = block->getData<field::GhostLayerField<double, 13>>(jID);
 
   WALBERLA_ASSERT_GREATER_EQUAL(-1, -int_c(j->nrOfGhostLayers()));
   double *RESTRICT const _data_j = j->dataAt(-1, -1, -1, 0);
@@ -1664,9 +1664,9 @@ void AdvectiveFluxKernel_double_precision::runOnCellInterval(const shared_ptr<St
   if (ci.empty())
     return;
 
+  auto j = block->getData<field::GhostLayerField<double, 13>>(jID);
   auto rho = block->getData<field::GhostLayerField<double, 1>>(rhoID);
   auto u = block->getData<field::GhostLayerField<double, 3>>(uID);
-  auto j = block->getData<field::GhostLayerField<double, 13>>(jID);
 
   WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin() - 1, -int_c(j->nrOfGhostLayers()));
   WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin() - 1, -int_c(j->nrOfGhostLayers()));

src/walberla_bridge/src/electrokinetics/generated_kernels/AdvectiveFluxKernel_single_precision.cpp

@@ -1625,9 +1625,9 @@ static FUNC_PREFIX void advectivefluxkernel_single_precision_advectivefluxkernel
 } // namespace internal_47df4b171f276b8c3a55fc08d45e245e
 
 void AdvectiveFluxKernel_single_precision::run(IBlock *block) {
-  auto j = block->getData<field::GhostLayerField<float, 13>>(jID);
   auto rho = block->getData<field::GhostLayerField<float, 1>>(rhoID);
   auto u = block->getData<field::GhostLayerField<float, 3>>(uID);
+  auto j = block->getData<field::GhostLayerField<float, 13>>(jID);
 
   WALBERLA_ASSERT_GREATER_EQUAL(-1, -int_c(j->nrOfGhostLayers()));
   float *RESTRICT const _data_j = j->dataAt(-1, -1, -1, 0);
@@ -1664,9 +1664,9 @@ void AdvectiveFluxKernel_single_precision::runOnCellInterval(const shared_ptr<St
   if (ci.empty())
     return;
 
-  auto j = block->getData<field::GhostLayerField<float, 13>>(jID);
   auto rho = block->getData<field::GhostLayerField<float, 1>>(rhoID);
   auto u = block->getData<field::GhostLayerField<float, 3>>(uID);
+  auto j = block->getData<field::GhostLayerField<float, 13>>(jID);
 
   WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin() - 1, -int_c(j->nrOfGhostLayers()));
   WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin() - 1, -int_c(j->nrOfGhostLayers()));

src/walberla_bridge/src/electrokinetics/generated_kernels/ContinuityKernel_double_precision.cpp

@@ -108,8 +108,8 @@ static FUNC_PREFIX void continuitykernel_double_precision_continuitykernel_doubl
 } // namespace internal_5c5e903f8ea7925cf790d7c2318b2c56
 
 void ContinuityKernel_double_precision::run(IBlock *block) {
-  auto rho = block->getData<field::GhostLayerField<double, 1>>(rhoID);
   auto j = block->getData<field::GhostLayerField<double, 13>>(jID);
+  auto rho = block->getData<field::GhostLayerField<double, 1>>(rhoID);
 
   WALBERLA_ASSERT_GREATER_EQUAL(-1, -int_c(j->nrOfGhostLayers()));
   double *RESTRICT const _data_j = j->dataAt(-1, -1, -1, 0);
@@ -140,8 +140,8 @@ void ContinuityKernel_double_precision::runOnCellInterval(const shared_ptr<Struc
   if (ci.empty())
     return;
 
-  auto rho = block->getData<field::GhostLayerField<double, 1>>(rhoID);
   auto j = block->getData<field::GhostLayerField<double, 13>>(jID);
+  auto rho = block->getData<field::GhostLayerField<double, 1>>(rhoID);
 
   WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin() - 1, -int_c(j->nrOfGhostLayers()));
   WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin() - 1, -int_c(j->nrOfGhostLayers()));

src/walberla_bridge/src/electrokinetics/generated_kernels/ContinuityKernel_single_precision.cpp

@@ -108,8 +108,8 @@ static FUNC_PREFIX void continuitykernel_single_precision_continuitykernel_singl
 } // namespace internal_990034b4e4dd57d2802b4bcb5f716e46
 
 void ContinuityKernel_single_precision::run(IBlock *block) {
-  auto j = block->getData<field::GhostLayerField<float, 13>>(jID);
   auto rho = block->getData<field::GhostLayerField<float, 1>>(rhoID);
+  auto j = block->getData<field::GhostLayerField<float, 13>>(jID);
 
   WALBERLA_ASSERT_GREATER_EQUAL(-1, -int_c(j->nrOfGhostLayers()));
   float *RESTRICT const _data_j = j->dataAt(-1, -1, -1, 0);
@@ -140,8 +140,8 @@ void ContinuityKernel_single_precision::runOnCellInterval(const shared_ptr<Struc
   if (ci.empty())
     return;
 
-  auto j = block->getData<field::GhostLayerField<float, 13>>(jID);
   auto rho = block->getData<field::GhostLayerField<float, 1>>(rhoID);
+  auto j = block->getData<field::GhostLayerField<float, 13>>(jID);
 
   WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin() - 1, -int_c(j->nrOfGhostLayers()));
   WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin() - 1, -int_c(j->nrOfGhostLayers()));