Synchronize for MPI

src/scheme/EulerKineticSolver2D.cpp

@@ -8,12 +8,15 @@
 #include <mesh/MeshFlatNodeBoundary.hpp>
 #include <mesh/MeshNodeBoundary.hpp>
 #include <mesh/MeshVariant.hpp>
+#include <scheme/DiscreteFunctionDPkVariant.hpp>
 #include <scheme/DiscreteFunctionDescriptorP0Vector.hpp>
 #include <scheme/DiscreteFunctionP0Vector.hpp>
 #include <scheme/DiscreteFunctionUtils.hpp>
 #include <scheme/IBoundaryConditionDescriptor.hpp>
 #include <scheme/IDiscreteFunctionDescriptor.hpp>
 #include <scheme/InflowListBoundaryConditionDescriptor.hpp>
+#include <scheme/PolynomialReconstruction.hpp>
+#include <scheme/PolynomialReconstructionDescriptor.hpp>
 #include <scheme/SymmetryBoundaryConditionDescriptor.hpp>
 
 #include <analysis/GaussLegendreQuadratureDescriptor.hpp>
@@ -827,6 +830,21 @@ class EulerKineticSolver2D
     const DiscreteFunctionP0<const double> rho_u2 = compute_PFn(Fn_rho_u2);
     const DiscreteFunctionP0<const double> rho_E  = compute_PFn(Fn_rho_E);
 
+    std::vector<std::shared_ptr<const IBoundaryDescriptor>> symmetry_boundary_descriptor_list;
+    // Reconstruction uniquement à l'intérieur
+    // for (auto&& bc_descriptor : bc_descriptor_list) {
+    //   if (bc_descriptor->type() == IBoundaryConditionDescriptor::Type::symmetry) {
+    //     symmetry_boundary_descriptor_list.push_back(bc_descriptor->boundaryDescriptor_shared());
+    //   }
+    // }
+
+    PolynomialReconstructionDescriptor reconstruction_descriptor(IntegrationMethodType::cell_center, 1,
+                                                                 symmetry_boundary_descriptor_list);
+
+    auto reconstruction =
+      PolynomialReconstruction{reconstruction_descriptor}.build(Fn_rho, Fn_rho_u1, Fn_rho_u2, Fn_rho_E);
+    auto DPk_Fn_rho = reconstruction[0]->get<DiscreteFunctionDPkVector<Dimension, const double>>();
+
     // NodeArray<double> Fr_rho    = compute_Flux1(Fn_rho);
     // NodeArray<double> Fr_rho_u1 = compute_Flux1(Fn_rho_u1);
     // NodeArray<double> Fr_rho_u2 = compute_Flux1(Fn_rho_u2);
@@ -844,6 +862,11 @@ class EulerKineticSolver2D
 
     apply_scalar_bc(Fr_rho, Fr_rho_u1, Fr_rho_u2, Fr_rho_E, rho, rho_u1, rho_u2, rho_E, bc_list);
 
+    synchronize(Fr_rho);
+    synchronize(Fr_rho_u1);
+    synchronize(Fr_rho_u2);
+    synchronize(Fr_rho_E);
+
     // DiscreteFunctionP0Vector<double> deltaLFn_rho    = compute_deltaLFn(Fr_rho);
     // DiscreteFunctionP0Vector<double> deltaLFn_rho_u1 = compute_deltaLFn(Fr_rho_u1);
     // DiscreteFunctionP0Vector<double> deltaLFn_rho_u2 = compute_deltaLFn(Fr_rho_u2);