From 52e50b5decf27e8618539d8682e1d05c0dceae4d Mon Sep 17 00:00:00 2001 From: mattijs Date: Wed, 12 Sep 2012 10:55:13 +0100 Subject: [PATCH] ENH: polyTopoChange: check on deleting cells without deleting/changing surrounding faces --- .../polyTopoChange/polyTopoChange.C | 15 +++++++++++++++ 1 file changed, 15 insertions(+) diff --git a/src/dynamicMesh/polyTopoChange/polyTopoChange/polyTopoChange.C b/src/dynamicMesh/polyTopoChange/polyTopoChange/polyTopoChange.C index b9b7bde14b..07618de976 100644 --- a/src/dynamicMesh/polyTopoChange/polyTopoChange/polyTopoChange.C +++ b/src/dynamicMesh/polyTopoChange/polyTopoChange/polyTopoChange.C @@ -495,6 +495,21 @@ void Foam::polyTopoChange::makeCells for (label faceI = 0; faceI < nActiveFaces; faceI++) { + if (faceOwner_[faceI] < 0) + { + FatalErrorIn + ( + "polyTopoChange::makeCells\n" + "(\n" + " const label,\n" + " labelList&,\n" + " labelList&\n" + ") const\n" + ) << "Face " << faceI << " is active but its owner has" + << " been deleted. This is usually due to deleting cells" + << " without modifying exposed faces to be boundary faces." + << exit(FatalError); + } nNbrs[faceOwner_[faceI]]++; } for (label faceI = 0; faceI < nActiveFaces; faceI++)