From 3b966afb9c276f741f5d8cea3ea7f3d9f37bdd0c Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Thu, 22 Feb 2024 16:08:51 +0000
Subject: [PATCH] BUG: snappyHexMesh: parallel inconsistency. Fixes #3106

---
 .../snappyHexMeshDriver/snappySnapDriverFeature.C         | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

diff --git a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriverFeature.C b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriverFeature.C
index 2e3b373f4f..e5e951c30d 100644
--- a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriverFeature.C
+++ b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriverFeature.C
@@ -6,7 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2011-2015 OpenFOAM Foundation
-    Copyright (C) 2015-2022 OpenCFD Ltd.
+    Copyright (C) 2015-2022,2024 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -513,7 +513,9 @@ void Foam::snappySnapDriver::calcNearestFacePointProperties
         forAll(pFaces, i)
         {
             label facei = pFaces[i];
-            if (isMasterFace[facei] && faceSurfaceGlobalRegion[facei] != -1)
+            label globalRegioni = faceSurfaceGlobalRegion[facei];
+
+            if (isMasterFace[pp.addressing()[facei]] && globalRegioni != -1)
             {
                 nFaces++;
             }
@@ -535,7 +537,7 @@ void Foam::snappySnapDriver::calcNearestFacePointProperties
             label facei = pFaces[i];
             label globalRegioni = faceSurfaceGlobalRegion[facei];
 
-            if (isMasterFace[facei] && globalRegioni != -1)
+            if (isMasterFace[pp.addressing()[facei]] && globalRegioni != -1)
             {
                 pNormals[nFaces] = faceSurfaceNormal[facei];
                 pDisp[nFaces] = faceDisp[facei];