diff --git a/applications/utilities/mesh/manipulation/checkMesh/checkTopology.C b/applications/utilities/mesh/manipulation/checkMesh/checkTopology.C index 7e395be879..8af38b9394 100644 --- a/applications/utilities/mesh/manipulation/checkMesh/checkTopology.C +++ b/applications/utilities/mesh/manipulation/checkMesh/checkTopology.C @@ -1,3 +1,28 @@ +/*---------------------------------------------------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | + \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation + \\/ M anipulation | +------------------------------------------------------------------------------- +License + This file is part of OpenFOAM. + + OpenFOAM is free software: you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 of the License, or + (at your option) any later version. + + OpenFOAM is distributed in the hope that it will be useful, but WITHOUT + ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or + FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License + for more details. + + You should have received a copy of the GNU General Public License + along with OpenFOAM. If not, see . + +\*---------------------------------------------------------------------------*/ + #include "checkTopology.H" #include "polyMesh.H" #include "Time.H" @@ -9,6 +34,8 @@ #include "emptyPolyPatch.H" #include "processorPolyPatch.H" +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + Foam::label Foam::checkTopology ( const polyMesh& mesh, @@ -287,6 +314,29 @@ Foam::label Foam::checkTopology rs ); ctr.write(); + + for (label i = 0; i < rs.nRegions(); i++) + { + cellSet rCells + ( + mesh, + "cellToRegion:" + Foam::name(i), + mesh.nCells()/100 + ); + + forAll(rs, j) + { + if (rs[j] == i) + { + rCells.insert(j); + } + } + + Info<< " <