Codage Rusanov 2D 3D composite version 1

src/language/modules/SchemeModule.cpp

@@ -520,6 +520,24 @@ SchemeModule::SchemeModule()
 
                               ));
 
+  this->_addBuiltinFunction("rusanov_eulerian_composite_solver_version1_with_checks",
+                            std::function(
+
+                              [](const std::shared_ptr<const DiscreteFunctionVariant>& rho,
+                                 const std::shared_ptr<const DiscreteFunctionVariant>& u,
+                                 const std::shared_ptr<const DiscreteFunctionVariant>& E,
+                                 const std::shared_ptr<const DiscreteFunctionVariant>& c,
+                                 const std::shared_ptr<const DiscreteFunctionVariant>& p,
+                                 const std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>>&
+                                   bc_descriptor_list,
+                                 const double& dt) -> std::tuple<std::shared_ptr<const DiscreteFunctionVariant>,
+                                                                 std::shared_ptr<const DiscreteFunctionVariant>,
+                                                                 std::shared_ptr<const DiscreteFunctionVariant>> {
+                                return rusanovEulerianCompositeSolver(rho, u, E, c, p, bc_descriptor_list, dt, true);
+                              }
+
+                              ));
+
   this->_addBuiltinFunction("eucclhyd_fluxes",
                             std::function(
 
@@ -775,6 +793,13 @@ SchemeModule::SchemeModule()
 
                                              ));
 
+  this->_addBuiltinFunction("compute_dt", std::function(
+
+                                            [](const std::shared_ptr<const DiscreteFunctionVariant>& u,
+                                               const std::shared_ptr<const DiscreteFunctionVariant>& c) -> double {
+                                              return compute_dt(u, c);
+                                            }));
+
   this->_addBuiltinFunction("cell_volume",
                             std::function(
 

src/scheme/CMakeLists.txt

@@ -4,6 +4,16 @@ add_library(
   PugsScheme
   AcousticSolver.cpp
   AcousticCompositeSolver.cpp
+#  JacobianAndStructuralInfoForSystemsofEquations.cpp
+#  ApproximateRiemannCompositeSolver.cpp
+#  ApproximateRiemannCompositeSolver_FluxForm.cpp
+#  ApproximateRiemannCompositeSolver_ViscousForm.cpp
+#  Roe_FluxForm_v2_CompositeSolver.cpp
+#  Roe_ViscousForm_v2_CompositeSolver.cpp
+#  VFFC_FluxForm_v1_CompositeSolver.cpp
+#  VFFC_ViscousForm_v1_CompositeSolver.cpp
+#  VFFC_FluxForm_v2_CompositeSolver.cpp
+#  VFFC_ViscousForm_v2_CompositeSolver.cpp
   HyperelasticSolver.cpp
   DiscreteFunctionIntegrator.cpp
   DiscreteFunctionInterpoler.cpp

src/scheme/RusanovEulerianCompositeSolver.hpp

@@ -8,6 +8,9 @@
 #include <memory>
 #include <vector>
 
+double compute_dt(const std::shared_ptr<const DiscreteFunctionVariant>& u_v,
+                  const std::shared_ptr<const DiscreteFunctionVariant>& c_v);
+
 std::tuple<std::shared_ptr<const DiscreteFunctionVariant>,
            std::shared_ptr<const DiscreteFunctionVariant>,
            std::shared_ptr<const DiscreteFunctionVariant>>
@@ -18,6 +21,7 @@ rusanovEulerianCompositeSolver(
   const std::shared_ptr<const DiscreteFunctionVariant>& c,
   const std::shared_ptr<const DiscreteFunctionVariant>& p,
   const std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>>& bc_descriptor_list,
-  const double& dt);
+  const double& dt,
+  const bool check = false);
 
 #endif   // RUSANOV_EULERIAN_COMPOSITE_SOLVER_HPP