Begin try for second order scheme 4 using periodic bc in 1d

src/language/modules/KineticSchemeModule.cpp

@@ -12,6 +12,7 @@
 #include <scheme/EulerKineticSolverLagrangeMultiD.hpp>
 #include <scheme/EulerKineticSolverLagrangeMultiDLocal.hpp>
 #include <scheme/EulerKineticSolverLagrangeMultiDLocalOrder2.hpp>
+#include <scheme/EulerKineticSolverLagrangeMultiDLocalOrder2Periodic.hpp>
 #include <scheme/EulerKineticSolverMeanFluxMood.hpp>
 #include <scheme/EulerKineticSolverMoodFD.hpp>
 #include <scheme/EulerKineticSolverMoodFV.hpp>
@@ -904,6 +905,27 @@ KineticSchemeModule::KineticSchemeModule()
 
                               ));
 
+  this->_addBuiltinFunction("euler_kinetic_lagrange_multiD_local_order2_periodic",
+                            std::function(
+
+                              [](const double& dt, const size_t& L, const size_t& M, const size_t& k,
+                                 const double& gamma, const double& eps, const size_t& space_order,
+                                 const size_t& time_order, 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>& p,
+                                 const std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>>&
+                                   bc_descriptor_list) -> std::tuple<std::shared_ptr<const MeshVariant>,
+                                                                     std::shared_ptr<const DiscreteFunctionVariant>,
+                                                                     std::shared_ptr<const DiscreteFunctionVariant>,
+                                                                     std::shared_ptr<const DiscreteFunctionVariant>> {
+                                return eulerKineticSolverLagrangeMultiDLocalOrder2(dt, L, M, k, gamma, eps, space_order,
+                                                                                   time_order, rho, u, E, p,
+                                                                                   bc_descriptor_list);
+                              }
+
+                              ));
+
   this->_addBuiltinFunction("euler_kinetic_acoustic_lagrange_FV",
                             std::function(
 

src/scheme/CMakeLists.txt

@@ -20,6 +20,7 @@ add_library(
   EulerKineticSolverLagrangeMultiD.cpp
   EulerKineticSolverLagrangeMultiDLocal.cpp
   EulerKineticSolverLagrangeMultiDLocalOrder2.cpp
+  EulerKineticSolverLagrangeMultiDLocalOrder2Periodic.cpp
   EulerKineticSolverMeanFluxMood.cpp
   EulerKineticSolverOneFluxMood.cpp
   EulerKineticSolverMoodFD.cpp

src/scheme/EulerKineticSolverLagrangeMultiDLocalOrder2Periodic.hpp

+#ifndef EULER_KINETIC_SOLVER_LAGRANGE_MULTI_D_LOCAL_ORDER_2PERIODIC_HPP
+#define EULER_KINETIC_SOLVER_LAGRANGE_MULTI_D_LOCAL_ORDER_2PERIODIC_HPP
+
+#include <language/utils/FunctionSymbolId.hpp>
+#include <scheme/DiscreteFunctionVariant.hpp>
+
+class IBoundaryConditionDescriptor;
+
+std::tuple<std::shared_ptr<const MeshVariant>,
+           std::shared_ptr<const DiscreteFunctionVariant>,
+           std::shared_ptr<const DiscreteFunctionVariant>,
+           std::shared_ptr<const DiscreteFunctionVariant>>
+eulerKineticSolverLagrangeMultiDLocalOrder2Periodic(
+  const double& dt,
+  const size_t& L,
+  const size_t& M,
+  const size_t& k,
+  const double& gamma,
+  const double& eps,
+  const size_t& space_order,
+  const size_t& time_order,
+  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>& p,
+  const std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>>& bc_descriptor_list);
+
+#endif   // EULER_KINETIC_SOLVER_LAGRANGE_MULTI_D_LOCAL_ORDER_2PERIODIC_HPP