From dcf8c1ab2008b65ff8a15685b7f94e787c127195 Mon Sep 17 00:00:00 2001
From: HOCH PHILIPPE <philippe.hoch@gmail.com>
Date: Sat, 16 Nov 2024 20:41:15 +0100
Subject: [PATCH] Fix parallelism
---
src/scheme/RusanovEulerianCompositeSolver_v2.cpp | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp
index ea8ca8d01..c78fe14ac 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();
--
GitLab