lbboundaries reset after ``lbf.load_checkpoint()`` about espresso HOT 2 CLOSED

works like a charm

from espresso.

Thank you for this MWE!

When reloading LB populations, any LB boundary information is overwritten with zeros. There is functionality to regenerate the LB boundary information, but the callback call is missing. Here is the bugfix for 4.2.0 and 4.3-dev:

diff --git a/src/core/grid_based_algorithms/lb.cpp b/src/core/grid_based_algorithms/lb.cpp
index 504f3088f..a6e37cb58 100644
--- a/src/core/grid_based_algorithms/lb.cpp
+++ b/src/core/grid_based_algorithms/lb.cpp
@@ -33,4 +33,5 @@
 #include "communication.hpp"
 #include "errorhandling.hpp"
+#include "event.hpp"
 #include "grid.hpp"
 #include "grid_based_algorithms/lb_boundaries.hpp"
@@ -208,3 +209,4 @@ void lb_initialize_fields(std::vector<LB_FluidNode> &lb_fields,
 #endif // LB_BOUNDARIES
   }
+  on_lbboundary_change();
 }

Could you please confirm that your simulations are fixed with this change? From my side, both the testsuite and your MWE are successful.

This regression was introduced in 4.2.0, back when serialization of object lists was temporarily removed from ESPResSo. The checkpointing test worked around this issue by checking LB boundaries before the LB populations were reloaded. I'll rewrite the test to check LB boundary information both before and after the LB populations are reloaded.

from espresso.