WIP: Sort out boundaries with pull scheme

espressomd · Aug 22, 2023 · fc1903e · fc1903e
1 parent 50bfd3a
commit fc1903e
Show file tree

Hide file tree

Showing 62 changed files with 5,701 additions and 14,533 deletions.
diff --git a/maintainer/walberla_kernels/generate_lb_kernels.py b/maintainer/walberla_kernels/generate_lb_kernels.py
@@ -33,7 +33,7 @@
 import lbmpy.enums
 
 import lbmpy_walberla
-import lbmpy_espresso
+import lbmpy_walberla.additional_data_handler
 
 import lees_edwards
 import relaxation_rates
@@ -184,17 +184,17 @@ def paramlist(parameters, keys):
         )
 
     # boundary conditions
-    ubb_dynamic = lbmpy_espresso.UBB(
-        lambda *args: None, dim=3, data_type=config.data_type.default_factory())
-    ubb_data_handler = lbmpy_espresso.BounceBackSlipVelocityUBB(
-        method.stencil, ubb_dynamic)
-
+    ubb = lbmpy.boundaries.boundaryconditions.UBB(
+        lambda *args: None, dim=3, data_type=config.data_type.default_factory(), adapt_velocity_to_force=True)
+    ubb_data_handler = lbmpy_walberla.additional_data_handler.default_additional_data_handler(
+        ubb, method, fields["pdfs"], target=target) 
     for _, target_suffix in paramlist(parameters, ("GPU", "CPU")):
         lbmpy_walberla.generate_boundary(
-            ctx, f"Dynamic_UBB_{precision_suffix}{target_suffix}", ubb_dynamic,
+            ctx, f"Dynamic_UBB_{precision_suffix}{target_suffix}", ubb,
             method, additional_data_handler=ubb_data_handler,
             streaming_pattern="pull", target=target)
 
+        # replace Walberla-specific real_t by float/double.
         with open(f"Dynamic_UBB_{precision_suffix}{target_suffix}.h", "r+") as f:
             content = f.read()
             f.seek(0)

diff --git a/maintainer/walberla_kernels/lbmpy_espresso.py b/maintainer/walberla_kernels/lbmpy_espresso.py
diff --git a/...rla_bridge/src/electrokinetics/generated_kernels/AdvectiveFluxKernel_double_precision.cpp b/...rla_bridge/src/electrokinetics/generated_kernels/AdvectiveFluxKernel_double_precision.cpp