diff --git a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp index ea8ca8d017e8264aad64266a6251ab90e6e69633..c78fe14ace7c48b33a7819de4582b8b3422e8329 100644 --- a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp +++ b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp @@ -12,6 +12,7 @@ #include <mesh/MeshNodeBoundary.hpp> #include <mesh/MeshTraits.hpp> #include <mesh/MeshVariant.hpp> +#include <mesh/SubItemValuePerItemUtils.hpp> #include <scheme/DiscreteFunctionUtils.hpp> #include <scheme/InflowListBoundaryConditionDescriptor.hpp> #include <variant> @@ -1056,6 +1057,7 @@ class RusanovEulerianCompositeSolver_v2 } }); + synchronize(Gjr); if (checkLocalConservation) { auto is_boundary_node = p_mesh->connectivity().isBoundaryNode(); @@ -1128,6 +1130,7 @@ class RusanovEulerianCompositeSolver_v2 Gjf[j][l] *= 1. / face_to_cell.size(); } }); + synchronize(Gjf); if (checkLocalConservation) { auto is_boundary_face = p_mesh->connectivity().isBoundaryFace(); @@ -1210,6 +1213,7 @@ class RusanovEulerianCompositeSolver_v2 Gje[j][l] *= 1. / edge_to_cell.size(); } }); + synchronize(Gje); if (checkLocalConservation) { auto is_boundary_edge = p_mesh->connectivity().isBoundaryEdge();