ajout d'une fonction sources pour flux rusanov

src/language/modules/CompositeSchemeSourcesModule.cpp

@@ -18,7 +18,33 @@
 #include <scheme/RusanovEulerianCompositeSolver_v2_o2.hpp>
 #include <scheme/RusanovEulerianCompositeSolver_v2_order_n.hpp>
 
-CompositeSchemeSourcesModule::CompositeSchemeSourcesModule() {}
+CompositeSchemeSourcesModule::CompositeSchemeSourcesModule(){
+                              
+   this->_addBuiltinFunction("rusanov_eulerian_composite_solver_v2_with_source",
+                            std::function(
+                              [](const std::shared_ptr<const DiscreteFunctionVariant>& rho,
+                                 const std::shared_ptr<const DiscreteFunctionVariant>& u,
+                                 const std::shared_ptr<const DiscreteFunctionVariant>& E, const double& gamma,
+                                 const std::shared_ptr<const DiscreteFunctionVariant>& c,
+                                 const std::shared_ptr<const DiscreteFunctionVariant>& p,   // const size_t& degree,
+                                 const std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>>&
+                                   bc_descriptor_list,
+
+                                 const double& dt,
+                                 const std::shared_ptr<const DiscreteFunctionVariant>& function_id1,
+                                 const std::shared_ptr<const DiscreteFunctionVariant>& function_id2,
+                                 const std::shared_ptr<const DiscreteFunctionVariant>& function_id3)
+                                 
+                                 			-> std::tuple<std::shared_ptr<const DiscreteFunctionVariant>,
+                                                           std::shared_ptr<const DiscreteFunctionVariant>,
+                                                           std::shared_ptr<const DiscreteFunctionVariant>> {
+                                return rusanovEulerianCompositeSolver_v2_with_source(rho, u, E, gamma, c, p,   // degree,
+                                                                         bc_descriptor_list, dt, function_id1, function_id2,function_id3);
+                              }
+
+                              ));
+
+}
 
 void
 CompositeSchemeSourcesModule::registerOperators() const

src/scheme/RusanovEulerianCompositeSolver_v2.hpp

@@ -27,4 +27,22 @@ rusanovEulerianCompositeSolver_v2(
   const double& dt,
   const bool check = false);
   
+std::tuple<std::shared_ptr<const DiscreteFunctionVariant>,
+           std::shared_ptr<const DiscreteFunctionVariant>,
+           std::shared_ptr<const DiscreteFunctionVariant>>
+rusanovEulerianCompositeSolver_v2_with_source(
+  const std::shared_ptr<const DiscreteFunctionVariant>& rho,
+  const std::shared_ptr<const DiscreteFunctionVariant>& u,
+  const std::shared_ptr<const DiscreteFunctionVariant>& E,
+  const double& gamma,
+  const std::shared_ptr<const DiscreteFunctionVariant>& c,
+  const std::shared_ptr<const DiscreteFunctionVariant>& p,
+  // const size_t& degree,
+  const std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>>& bc_descriptor_list,
+  const double& dt,
+  const std::shared_ptr<const DiscreteFunctionVariant>& function_id1,
+  const std::shared_ptr<const DiscreteFunctionVariant>& function_id2,
+  const std::shared_ptr<const DiscreteFunctionVariant>& function_id3,
+  const bool check = false);
+
 #endif   // RUSANOV_EULERIAN_COMPOSITE_SOLVER_V2_HPP