setAndNormalToFaceZone.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2013-2018 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 "setAndNormalToFaceZone.H"
27 #include "polyMesh.H"
28 #include "faceZoneSet.H"
29 
31 
32 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
33 
34 namespace Foam
35 {
36 
37 defineTypeNameAndDebug(setAndNormalToFaceZone, 0);
38 
39 addToRunTimeSelectionTable(topoSetSource, setAndNormalToFaceZone, word);
40 
41 addToRunTimeSelectionTable(topoSetSource, setAndNormalToFaceZone, istream);
42 
43 }
44 
45 
46 Foam::topoSetSource::addToUsageTable Foam::setAndNormalToFaceZone::usage_
47 (
48  setAndNormalToFaceZone::typeName,
49  "\n Usage: setAndNormalToFaceZone <faceSet> <normal>\n\n"
50  " Select all faces in the faceSet and orient using normal.\n\n"
51 );
52 
53 
54 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
55 
57 (
58  const polyMesh& mesh,
59  const word& setName,
60  const vector& normal
61 )
62 :
63  topoSetSource(mesh),
64  setName_(setName),
65  normal_(normal)
66 {}
67 
68 
70 (
71  const polyMesh& mesh,
72  const dictionary& dict
73 )
74 :
75  topoSetSource(mesh),
76  setName_(dict.lookup("faceSet")),
77  normal_(dict.lookup("normal"))
78 {}
79 
80 
82 (
83  const polyMesh& mesh,
84  Istream& is
85 )
86 :
87  topoSetSource(mesh),
88  setName_(checkIs(is)),
89  normal_(checkIs(is))
90 {}
91 
92 
93 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
94 
96 {}
97 
98 
99 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
100 
102 (
103  const topoSetSource::setAction action,
104  topoSet& set
105 ) const
106 {
107  if (!isA<faceZoneSet>(set))
108  {
110  << "Operation only allowed on a faceZoneSet." << endl;
111  }
112  else
113  {
114  faceZoneSet& fzSet = refCast<faceZoneSet>(set);
115 
116  if ((action == topoSetSource::NEW) || (action == topoSetSource::ADD))
117  {
118  Info<< " Adding all faces from faceSet " << setName_
119  << " ..." << endl;
120 
121  // Load the sets
122  faceSet fSet(mesh_, setName_);
123 
124  // Start off from copy
125  DynamicList<label> newAddressing(fzSet.addressing());
126  DynamicList<bool> newFlipMap(fzSet.flipMap());
127 
128  const faceList& faces = mesh_.faces();
129  const pointField& points = mesh_.points();
130 
131  forAllConstIter(faceSet, fSet, iter)
132  {
133  label facei = iter.key();
134 
135  if (!fzSet.found(facei))
136  {
137  newAddressing.append(facei);
138 
139  const vector a = faces[facei].area(points);
140  if ((a & normal_) > 0)
141  {
142  newFlipMap.append(false);
143  }
144  else
145  {
146  newFlipMap.append(true);
147  }
148  }
149  }
150 
151  fzSet.addressing().transfer(newAddressing);
152  fzSet.flipMap().transfer(newFlipMap);
153  fzSet.updateSet();
154  }
155  else if (action == topoSetSource::DELETE)
156  {
157  Info<< " Removing all faces from faceSet " << setName_
158  << " ..." << endl;
159 
160  // Load the set
161  faceSet loadedSet(mesh_, setName_);
162 
163  // Start off empty
164  DynamicList<label> newAddressing(fzSet.addressing().size());
165  DynamicList<bool> newFlipMap(fzSet.flipMap().size());
166 
167  forAll(fzSet.addressing(), i)
168  {
169  if (!loadedSet.found(fzSet.addressing()[i]))
170  {
171  newAddressing.append(fzSet.addressing()[i]);
172  newFlipMap.append(fzSet.flipMap()[i]);
173  }
174  }
175  fzSet.addressing().transfer(newAddressing);
176  fzSet.flipMap().transfer(newFlipMap);
177  fzSet.updateSet();
178  }
179  }
180 }
181 
182 
183 // ************************************************************************* //
setAndNormalToFaceZone(const polyMesh &mesh, const word &setName, const vector &normal)
Construct from components.
const boolList & flipMap() const
Definition: faceZoneSet.H:118
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:434
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
A list of face labels.
Definition: faceSet.H:48
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:158
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
Definition: Istream.H:57
void size(const label)
Override size to be inconsistent with allocated storage.
Definition: ListI.H:163
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:256
Macros for easy insertion into run-time selection tables.
Base class of a source for a topoSet.
Definition: topoSetSource.H:63
const pointField & points
bool found(const Key &) const
Return true if hashedEntry is found in table.
Definition: HashTable.C:113
A class for handling words, derived from string.
Definition: word.H:59
setAction
Enumeration defining the valid actions.
Definition: topoSetSource.H:82
DynamicList< T, SizeInc, SizeMult, SizeDiv > & append(const T &)
Append an element at the end of the list.
Definition: DynamicListI.H:296
forAllConstIter(PtrDictionary< phaseModel >, mixture.phases(), phase)
Definition: pEqn.H:29
virtual void applyToSet(const topoSetSource::setAction action, topoSet &) const
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
defineTypeNameAndDebug(combustionModel, 0)
void updateSet()
Sort addressing and make faceSet part consistent with addressing.
Definition: faceZoneSet.C:51
Like faceSet but -reads data from faceZone -updates faceZone when writing.
Definition: faceZoneSet.H:49
General set of labels of mesh quantity (points, cells, faces).
Definition: topoSet.H:61
#define WarningInFunction
Report a warning using Foam::Warning.
Class with constructor to add usage string to table.
messageStream Info
Mesh consisting of general polyhedral cells.
Definition: polyMesh.H:74
const labelList & addressing() const
Definition: faceZoneSet.H:107
virtual ~setAndNormalToFaceZone()
Destructor.
void transfer(List< T > &)
Transfer the contents of the argument List into this list.
Definition: List.C:342
Namespace for OpenFOAM.
ITstream & lookup(const word &, bool recursive=false, bool patternMatch=true) const
Find and return an entry data stream.
Definition: dictionary.C:583