src/meshTools/sets/cellSources/nbrToCell/nbrToCell.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright held by original author
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "nbrToCell.H"
  28 #include "polyMesh.H"
  29
  30 #include "addToRunTimeSelectionTable.H"
  31
  32 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36
  37 defineTypeNameAndDebug(nbrToCell, 0);
  38
  39 addToRunTimeSelectionTable(topoSetSource, nbrToCell, word);
  40
  41 addToRunTimeSelectionTable(topoSetSource, nbrToCell, istream);
  42
  43 }
  44
  45
  46 Foam::topoSetSource::addToUsageTable Foam::nbrToCell::usage_
  47 (
  48     nbrToCell::typeName,
  49     "\n    Usage: nbrToCell <nNeighbours>\n\n"
  50     "    Select all cells with <= nNeighbours neighbouring cells\n\n"
  51 );
  52
  53
  54 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  55
  56 void Foam::nbrToCell::combine(topoSet& set, const bool add) const
  57 {
  58     const cellList& cells = mesh().cells();
  59     const polyBoundaryMesh& patches = mesh_.boundaryMesh();
  60
  61     boolList isCoupled(mesh_.nFaces()-mesh_.nInternalFaces(), false);
  62
  63     forAll(patches, patchI)
  64     {
  65         const polyPatch& pp = patches[patchI];
  66
  67         if (pp.coupled())
  68         {
  69             label faceI = pp.start();
  70             forAll(pp, i)
  71             {
  72                 isCoupled[faceI-mesh_.nInternalFaces()] = true;
  73                 faceI++;
  74             }
  75         }
  76     }
  77
  78     forAll(cells, cellI)
  79     {
  80         const cell& cFaces = cells[cellI];
  81
  82         label nNbrCells = 0;
  83
  84         forAll(cFaces, i)
  85         {
  86             label faceI = cFaces[i];
  87
  88             if (mesh_.isInternalFace(faceI))
  89             {
  90                 nNbrCells++;
  91             }
  92             else if (isCoupled[faceI-mesh_.nInternalFaces()])
  93             {
  94                 nNbrCells++;
  95             }
  96         }
  97
  98         if (nNbrCells <= minNbrs_)
  99         {
 100             addOrDelete(set, cellI, add);
 101         }
 102     }
 103 }
 104
 105
 106 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 107
 108 // Construct from components
 109 Foam::nbrToCell::nbrToCell
 110 (
 111     const polyMesh& mesh,
 112     const label minNbrs
 113 )
 114 :
 115     topoSetSource(mesh),
 116     minNbrs_(minNbrs)
 117 {}
 118
 119
 120 // Construct from dictionary
 121 Foam::nbrToCell::nbrToCell
 122 (
 123     const polyMesh& mesh,
 124     const dictionary& dict
 125 )
 126 :
 127     topoSetSource(mesh),
 128     minNbrs_(readLabel(dict.lookup("neighbours")))
 129 {}
 130
 131
 132 // Construct from Istream
 133 Foam::nbrToCell::nbrToCell
 134 (
 135     const polyMesh& mesh,
 136     Istream& is
 137 )
 138 :
 139     topoSetSource(mesh),
 140     minNbrs_(readLabel(checkIs(is)))
 141 {}
 142
 143
 144 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 145
 146 Foam::nbrToCell::~nbrToCell()
 147 {}
 148
 149
 150 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 151
 152 void Foam::nbrToCell::applyToSet
 153 (
 154     const topoSetSource::setAction action,
 155     topoSet& set
 156 ) const
 157 {
 158     if ((action == topoSetSource::NEW) || (action == topoSetSource::ADD))
 159     {
 160         Info<< "    Adding cells with only " << minNbrs_ << " or less"
 161                 " neighbouring cells" << " ..." << endl;
 162
 163         combine(set, true);
 164     }
 165     else if (action == topoSetSource::DELETE)
 166     {
 167         Info<< "    Removing cells with only " << minNbrs_ << " or less"
 168                 " neighbouring cells" << " ..." << endl;
 169
 170         combine(set, false);
 171     }
 172 }
 173
 174
 175 // ************************************************************************* //