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 (C) 2011 OpenFOAM Foundation
   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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "nbrToCell.H"
  27 #include "polyMesh.H"
  28
  29 #include "addToRunTimeSelectionTable.H"
  30
  31 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  32
  33 namespace Foam
  34 {
  35
  36 defineTypeNameAndDebug(nbrToCell, 0);
  37
  38 addToRunTimeSelectionTable(topoSetSource, nbrToCell, word);
  39
  40 addToRunTimeSelectionTable(topoSetSource, nbrToCell, istream);
  41
  42 }
  43
  44
  45 Foam::topoSetSource::addToUsageTable Foam::nbrToCell::usage_
  46 (
  47     nbrToCell::typeName,
  48     "\n    Usage: nbrToCell <nNeighbours>\n\n"
  49     "    Select all cells with <= nNeighbours neighbouring cells\n\n"
  50 );
  51
  52
  53 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  54
  55 void Foam::nbrToCell::combine(topoSet& set, const bool add) const
  56 {
  57     const cellList& cells = mesh().cells();
  58     const polyBoundaryMesh& patches = mesh_.boundaryMesh();
  59
  60     boolList isCoupled(mesh_.nFaces()-mesh_.nInternalFaces(), false);
  61
  62     forAll(patches, patchI)
  63     {
  64         const polyPatch& pp = patches[patchI];
  65
  66         if (pp.coupled())
  67         {
  68             label faceI = pp.start();
  69             forAll(pp, i)
  70             {
  71                 isCoupled[faceI-mesh_.nInternalFaces()] = true;
  72                 faceI++;
  73             }
  74         }
  75     }
  76
  77     forAll(cells, cellI)
  78     {
  79         const cell& cFaces = cells[cellI];
  80
  81         label nNbrCells = 0;
  82
  83         forAll(cFaces, i)
  84         {
  85             label faceI = cFaces[i];
  86
  87             if (mesh_.isInternalFace(faceI))
  88             {
  89                 nNbrCells++;
  90             }
  91             else if (isCoupled[faceI-mesh_.nInternalFaces()])
  92             {
  93                 nNbrCells++;
  94             }
  95         }
  96
  97         if (nNbrCells <= minNbrs_)
  98         {
  99             addOrDelete(set, cellI, add);
 100         }
 101     }
 102 }
 103
 104
 105 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 106
 107 // Construct from components
 108 Foam::nbrToCell::nbrToCell
 109 (
 110     const polyMesh& mesh,
 111     const label minNbrs
 112 )
 113 :
 114     topoSetSource(mesh),
 115     minNbrs_(minNbrs)
 116 {}
 117
 118
 119 // Construct from dictionary
 120 Foam::nbrToCell::nbrToCell
 121 (
 122     const polyMesh& mesh,
 123     const dictionary& dict
 124 )
 125 :
 126     topoSetSource(mesh),
 127     minNbrs_(readLabel(dict.lookup("neighbours")))
 128 {}
 129
 130
 131 // Construct from Istream
 132 Foam::nbrToCell::nbrToCell
 133 (
 134     const polyMesh& mesh,
 135     Istream& is
 136 )
 137 :
 138     topoSetSource(mesh),
 139     minNbrs_(readLabel(checkIs(is)))
 140 {}
 141
 142
 143 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 144
 145 Foam::nbrToCell::~nbrToCell()
 146 {}
 147
 148
 149 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 150
 151 void Foam::nbrToCell::applyToSet
 152 (
 153     const topoSetSource::setAction action,
 154     topoSet& set
 155 ) const
 156 {
 157     if ((action == topoSetSource::NEW) || (action == topoSetSource::ADD))
 158     {
 159         Info<< "    Adding cells with only " << minNbrs_ << " or less"
 160                 " neighbouring cells" << " ..." << endl;
 161
 162         combine(set, true);
 163     }
 164     else if (action == topoSetSource::DELETE)
 165     {
 166         Info<< "    Removing cells with only " << minNbrs_ << " or less"
 167                 " neighbouring cells" << " ..." << endl;
 168
 169         combine(set, false);
 170     }
 171 }
 172
 173
 174 // ************************************************************************* //