Debut Codage Roe composite

src/scheme/CMakeLists.txt

@@ -24,6 +24,7 @@ add_library(
   RusanovEulerianCompositeSolverTools.cpp
   RusanovEulerianCompositeSolver.cpp
   RusanovEulerianCompositeSolver_v2.cpp
+  RoeViscousFormEulerianCompositeSolver_v2.cpp
 )
 
 target_link_libraries(

src/scheme/RoeViscousFormEulerianCompositeSolver_v2.hpp

+#ifndef ROE_VISCOUS_FORM_EULERIAN_COMPOSITE_SOLVER_V2_HPP
+#define ROE_VISCOUS_FORM_EULERIAN_COMPOSITE_SOLVER_V2_HPP
+
+#include <mesh/MeshVariant.hpp>
+#include <scheme/DiscreteFunctionVariant.hpp>
+#include <scheme/IBoundaryConditionDescriptor.hpp>
+#include <scheme/RusanovEulerianCompositeSolverTools.hpp>
+
+#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>>
+roeViscousFormEulerianCompositeSolver_v2(
+  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 std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>>& bc_descriptor_list,
+  const double& dt,
+  const bool check = false);
+
+#endif   // ROE_VISCOUS_FORM_EULERIAN_COMPOSITE_SOLVER_V2_HPP

src/scheme/RusanovEulerianCompositeSolverTools.hpp

@@ -17,6 +17,9 @@ double EvaluateMaxEigenValueTimesNormalLengthInGivenDirection(   // const double
   const double& c_mean,
   const Rd& normal);
 
+// template<class Rpxd>
+// Rpxd Flux(const double& Rho, const Rd& u, const R& E, )
+
 double compute_dt(const std::shared_ptr<const DiscreteFunctionVariant>& u_v,
                   const std::shared_ptr<const DiscreteFunctionVariant>& c_v);