Add boundary condition support for fluxing advection

Allowed bc tyes are
- inflow
- outflow
- symmetry

Inflow is not yet implemented for variables of type P0Vector

src/language/modules/SchemeModule.cpp

@@ -677,13 +677,14 @@ SchemeModule::SchemeModule()
                                                  }
 
                                                  ));
-  this->_addBuiltinFunction("fluxing_advection",
-                            std::function(
+
+  this->_addBuiltinFunction("fluxing_advection", std::function(
 
                                                    [](const std::shared_ptr<const IMesh> new_mesh,
-                                 const std::vector<std::shared_ptr<const DiscreteFunctionVariant>>& remapped_variables)
+                                                      const std::vector<std::shared_ptr<const VariableBCDescriptor>>&
+                                                        remapped_quantity_with_bc)
                                                      -> std::vector<std::shared_ptr<const DiscreteFunctionVariant>> {
-                                return advectByFluxing(new_mesh, remapped_variables);
+                                                     return advectByFluxing(new_mesh, remapped_quantity_with_bc);
                                                    }
 
                                                    ));

src/scheme/FluxingAdvectionSolver.hpp

@@ -3,6 +3,7 @@
 
 #include <language/utils/FunctionSymbolId.hpp>
 #include <scheme/DiscreteFunctionVariant.hpp>
+#include <scheme/VariableBCDescriptor.hpp>
 
 #include <vector>
 
@@ -10,4 +11,8 @@ std::vector<std::shared_ptr<const DiscreteFunctionVariant>> advectByFluxing(
   const std::shared_ptr<const IMesh> new_mesh,
   const std::vector<std::shared_ptr<const DiscreteFunctionVariant>>& remapped_variables);
 
+std::vector<std::shared_ptr<const DiscreteFunctionVariant>> advectByFluxing(
+  const std::shared_ptr<const IMesh> new_mesh,
+  const std::vector<std::shared_ptr<const VariableBCDescriptor>>& remapped_variables_with_bc);
+
 #endif   // FLUXING_ADVECION_SOLVER_HPP