Add update_collision_model to script interface (not working yet)

src/core/lb/LBNone.hpp

@@ -68,6 +68,7 @@ struct LBNone {
   void on_timestep_change() const { throw NoLBActive{}; }
   void on_temperature_change() const { throw NoLBActive{}; }
   void on_lees_edwards_change() const { throw NoLBActive{}; }
+  void update_collision_model() const { throw NoLBActive{}; }
 };
 
 } // namespace LB

src/core/lb/Solver.cpp

@@ -146,6 +146,12 @@ void Solver::on_lees_edwards_change() {
   }
 }
 
+void Solver::update_collision_model() {
+  if (impl->solver) {
+    std::visit([](auto &ptr) { ptr->update_collision_model(); }, *impl->solver);
+  }
+}
+
 bool Solver::is_gpu() const {
   check_solver(impl);
   return std::visit([](auto &ptr) { return ptr->is_gpu(); }, *impl->solver);

src/core/lb/Solver.hpp

@@ -169,6 +169,7 @@ struct Solver : public System::Leaf<Solver> {
   void on_timestep_change();
   void on_temperature_change();
   void on_lees_edwards_change();
+  void update_collision_model();
   void veto_boxl_change() const;
 
 private:

src/script_interface/walberla/LBFluid.cpp

@@ -188,6 +188,8 @@ void LBFluid::do_construct(VariantMap const &params) {
       throw std::domain_error("Parameter 'kinematic_viscosity' must be >= 0");
     }
     make_instance(params);
+    ::System::get_system().lb.update_collision_model();
+    /*
     auto const &system = ::System::get_system();
     if (auto le_protocol = system.lees_edwards->get_protocol()) {
       if (lb_kT != 0.) {
@@ -210,6 +212,7 @@ void LBFluid::do_construct(VariantMap const &params) {
       m_instance->set_collision_model(lb_kT, m_seed);
     }
     m_instance->ghost_communication(); // synchronize ghost layers
+    */
     m_instance->set_external_force(lb_ext_f);
     for (auto &vtk : m_vtk_writers) {
       vtk->attach_to_lattice(m_instance, get_latice_to_md_units_conversion());