Adding generic Rusanov solver at arbitrary degree reconstruction (not yet fully implemented)

src/scheme/CMakeLists.txt

@@ -29,6 +29,7 @@ add_library(
   RoeViscousFormEulerianCompositeSolver_v2.cpp
   RusanovEulerianCompositeSolver_o2.cpp
   RusanovEulerianCompositeSolver_v2_o2.cpp
+  RusanovEulerianCompositeSolver_v2_order_n.cpp
   RoeViscousFormEulerianCompositeSolver_v2_o2.cpp
 )
 

src/scheme/RusanovEulerianCompositeSolver_v2_order_n.hpp

+#ifndef RUSANOV_EULERIAN_COMPOSITE_SOLVER_V2_ORDER_N_HPP
+#define RUSANOV_EULERIAN_COMPOSITE_SOLVER_V2_ORDER_N_HPP
+
+#include <memory>
+#include <mesh/MeshVariant.hpp>
+#include <scheme/CellbyCellLimitation.hpp>
+#include <scheme/DiscreteFunctionVariant.hpp>
+#include <scheme/IBoundaryConditionDescriptor.hpp>
+#include <scheme/RusanovEulerianCompositeSolverTools.hpp>
+#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>>
+rusanovEulerianCompositeSolver_v2_order_n(
+  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 bool check = false);
+
+#endif   // RUSANOV_EULERIAN_COMPOSITE_SOLVER_V2_ORDER_N_HPP