Warning: do not compile

src/language/modules/SchemeModule.cpp

@@ -30,6 +30,7 @@
 #include <scheme/DiscreteFunctionVectorInterpoler.hpp>
 #include <scheme/ExternalBoundaryConditionDescriptor.hpp>
 #include <scheme/FixedBoundaryConditionDescriptor.hpp>
+#include <scheme/FluxingAdvectionSolver.hpp>
 #include <scheme/FourierBoundaryConditionDescriptor.hpp>
 #include <scheme/FreeBoundaryConditionDescriptor.hpp>
 #include <scheme/HyperelasticSolver.hpp>
@@ -589,6 +590,13 @@ SchemeModule::SchemeModule()
                                                  }
 
                                                  ));
+  this->_addBuiltinFunction("fluxing_advection", std::function(
+
+                                                   [](const std::shared_ptr<const IMesh> new_mesh,
+                                                      const std::shared_ptr<const DiscreteFunctionVariant>& old_q)
+                                                     -> std::shared_ptr<const DiscreteFunctionVariant> {
+                                                     return FluxingAdvectionSolverHandler(new_mesh, old_q);
+                                                   }));
 
   MathFunctionRegisterForVh{*this};
 }

src/scheme/FluxingAdvectionSolver.cpp

@@ -85,7 +85,7 @@ class FluxingAdvectionSolver
   FluxingAdvectionSolver(const std::shared_ptr<const MeshType> old_mesh,
                          const std::shared_ptr<const MeshType> new_mesh,
                          const DiscreteFunctionP0<Dimension, const double>& old_q)
-    : m_old_mesh{old_mesh}, m_new_mesh{new_mesh}, m_old_q{old_q}
+    : m_old_mesh{old_mesh}, m_new_mesh{new_mesh}
   {}
 
   ~FluxingAdvectionSolver() = default;