diff --git a/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C b/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C
index df95c968187d5cabfefd7d407a793e525e65b807..7bd5bd0d4b17c640e2de955336b8af91e25b92de 100644
--- a/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C
+++ b/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C
@@ -210,6 +210,7 @@ Foam::labelList Foam::polyMeshTetDecomposition::findFaceBasePts
 )
 {
     const labelList& pOwner = mesh.faceOwner();
+    const vectorField& pC = mesh.cellCentres();
 
     // Find a suitable base point for each face, considering both
     // cells for interface faces or those on coupled patches
@@ -228,7 +229,7 @@ Foam::labelList Foam::polyMeshTetDecomposition::findFaceBasePts
     for(label faceI = nInternalFaces; faceI < mesh.nFaces(); faceI++)
     {
         neighbourCellCentres[faceI - nInternalFaces] =
-            mesh.cellCentres()[pOwner[faceI]];
+            pC[pOwner[faceI]];
     }
 
     syncTools::swapBoundaryFacePositions(mesh, neighbourCellCentres);