pointToCell.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2011-2016 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 "pointToCell.H"
27 #include "polyMesh.H"
28 #include "pointSet.H"
29 
31 
32 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
33 
34 namespace Foam
35 {
36  defineTypeNameAndDebug(pointToCell, 0);
37  addToRunTimeSelectionTable(topoSetSource, pointToCell, word);
38  addToRunTimeSelectionTable(topoSetSource, pointToCell, istream);
39 
40  template<>
41  const char* Foam::NamedEnum
42  <
44  2
45  >::names[] =
46  {
47  "any",
48  "edge"
49  };
50 }
51 
52 
53 Foam::topoSetSource::addToUsageTable Foam::pointToCell::usage_
54 (
55  pointToCell::typeName,
56  "\n Usage: pointToCell <pointSet> any|edge\n\n"
57  " Select all cells with any point ('any') or any edge ('edge')"
58  " in the pointSet\n\n"
59 );
60 
62  Foam::pointToCell::pointActionNames_;
63 
64 
65 // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
66 
67 void Foam::pointToCell::combine(topoSet& set, const bool add) const
68 {
69  // Load the set
70  pointSet loadedSet(mesh_, setName_);
71 
72 
73  // Handle any selection
74  if (option_ == ANY)
75  {
76  forAllConstIter(pointSet, loadedSet, iter)
77  {
78  const label pointi = iter.key();
79  const labelList& pCells = mesh_.pointCells()[pointi];
80 
81  forAll(pCells, pCelli)
82  {
83  addOrDelete(set, pCells[pCelli], add);
84  }
85  }
86  }
87  else if (option_ == EDGE)
88  {
89  const faceList& faces = mesh_.faces();
90  forAll(faces, facei)
91  {
92  const face& f = faces[facei];
93 
94  forAll(f, fp)
95  {
96  if (loadedSet.found(f[fp]) && loadedSet.found(f.nextLabel(fp)))
97  {
98  addOrDelete(set, mesh_.faceOwner()[facei], add);
99  if (mesh_.isInternalFace(facei))
100  {
101  addOrDelete(set, mesh_.faceNeighbour()[facei], add);
102  }
103  }
104  }
105  }
106  }
107 }
108 
109 
110 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
111 
112 // Construct from components
114 (
115  const polyMesh& mesh,
116  const word& setName,
117  const pointAction option
118 )
119 :
120  topoSetSource(mesh),
121  setName_(setName),
122  option_(option)
123 {}
124 
125 
126 // Construct from dictionary
128 (
129  const polyMesh& mesh,
130  const dictionary& dict
131 )
132 :
133  topoSetSource(mesh),
134  setName_(dict.lookup("set")),
135  option_(pointActionNames_.read(dict.lookup("option")))
136 {}
137 
138 
139 // Construct from Istream
141 (
142  const polyMesh& mesh,
143  Istream& is
144 )
145 :
146  topoSetSource(mesh),
147  setName_(checkIs(is)),
148  option_(pointActionNames_.read(checkIs(is)))
149 {}
150 
151 
152 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
153 
155 {}
156 
157 
158 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
159 
161 (
162  const topoSetSource::setAction action,
163  topoSet& set
164 ) const
165 {
166  if ((action == topoSetSource::NEW) || (action == topoSetSource::ADD))
167  {
168  Info<< " Adding cells according to pointSet " << setName_
169  << " ..." << endl;
170 
171  combine(set, true);
172  }
173  else if (action == topoSetSource::DELETE)
174  {
175  Info<< " Removing cells according to pointSet " << setName_
176  << " ..." << endl;
177 
178  combine(set, false);
179  }
180 }
181 
182 
183 // ************************************************************************* //
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:428
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
virtual ~pointToCell()
Destructor.
Definition: pointToCell.C:154
A face is a list of labels corresponding to mesh vertices.
Definition: face.H:75
A set of point labels.
Definition: pointSet.H:48
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:137
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
Definition: Istream.H:57
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:253
Initialise the NamedEnum HashTable from the static list of names.
Definition: NamedEnum.H:52
AccessType combine(const List< T > &, AccessOp aop=accessOp< T >())
Combines sublists into one big list.
Definition: ListListOps.C:34
Macros for easy insertion into run-time selection tables.
Base class of a source for a topoSet.
Definition: topoSetSource.H:63
virtual void applyToSet(const topoSetSource::setAction action, topoSet &) const
Definition: pointToCell.C:161
A class for handling words, derived from string.
Definition: word.H:59
setAction
Enumeration defining the valid actions.
Definition: topoSetSource.H:82
pointAction
Enumeration defining the valid options.
Definition: pointToCell.H:55
forAllConstIter(PtrDictionary< phaseModel >, mixture.phases(), phase)
Definition: pEqn.H:29
bool found(const Key &) const
Return true if hashedEntry is found in table.
Definition: HashTable.C:109
void add(FieldField< Field1, typename typeOfSum< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
defineTypeNameAndDebug(combustionModel, 0)
labelList f(nPoints)
General set of labels of mesh quantity (points, cells, faces).
Definition: topoSet.H:61
Class with constructor to add usage string to table.
messageStream Info
label nextLabel(const label i) const
Next vertex on face.
Definition: faceI.H:117
Mesh consisting of general polyhedral cells.
Definition: polyMesh.H:74
pointToCell(const polyMesh &mesh, const word &setName, const pointAction option)
Construct from components.
Definition: pointToCell.C:114
Namespace for OpenFOAM.
ITstream & lookup(const word &, bool recursive=false, bool patternMatch=true) const
Find and return an entry data stream.
Definition: dictionary.C:451