setAndNormalToFaceZone.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 | www.openfoam.com
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8  Copyright (C) 2013-2016 OpenFOAM Foundation
9  Copyright (C) 2018-2020 OpenCFD Ltd.
10 -------------------------------------------------------------------------------
11 License
12  This file is part of OpenFOAM.
13 
14  OpenFOAM is free software: you can redistribute it and/or modify it
15  under the terms of the GNU General Public License as published by
16  the Free Software Foundation, either version 3 of the License, or
17  (at your option) any later version.
18 
19  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
20  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
21  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
22  for more details.
23 
24  You should have received a copy of the GNU General Public License
25  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
26 
27 \*---------------------------------------------------------------------------*/
28 
29 #include "setAndNormalToFaceZone.H"
30 #include "polyMesh.H"
31 #include "faceZoneSet.H"
33 
34 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
35 
36 namespace Foam
37 {
38  defineTypeNameAndDebug(setAndNormalToFaceZone, 0);
39  addToRunTimeSelectionTable(topoSetSource, setAndNormalToFaceZone, word);
40  addToRunTimeSelectionTable(topoSetSource, setAndNormalToFaceZone, istream);
41 
43  (
44  topoSetFaceZoneSource,
45  setAndNormalToFaceZone,
46  word
47  );
49  (
50  topoSetFaceZoneSource,
51  setAndNormalToFaceZone,
52  istream
53  );
54 }
55 
56 
57 Foam::topoSetSource::addToUsageTable Foam::setAndNormalToFaceZone::usage_
58 (
59  setAndNormalToFaceZone::typeName,
60  "\n Usage: setAndNormalToFaceZone <faceSet> <normal>\n\n"
61  " Select all faces in the faceSet and orient using normal.\n\n"
62 );
63 
64 
65 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
66 
68 (
69  const polyMesh& mesh,
70  const word& setName,
71  const vector& normal
72 )
73 :
75  setName_(setName),
76  normal_(normal)
77 {}
78 
79 
81 (
82  const polyMesh& mesh,
83  const dictionary& dict
84 )
85 :
87  setName_(dict.get<word>("faceSet")),
88  normal_(dict.get<vector>("normal"))
89 {}
90 
91 
93 (
94  const polyMesh& mesh,
95  Istream& is
96 )
97 :
99  setName_(checkIs(is)),
100  normal_(checkIs(is))
101 {}
102 
103 
104 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
105 
107 (
108  const topoSetSource::setAction action,
109  topoSet& set
110 ) const
111 {
112  if (!isA<faceZoneSet>(set))
113  {
115  << "Operation only allowed on a faceZoneSet." << endl;
116  return;
117  }
118  else
119  {
120  faceZoneSet& zoneSet = refCast<faceZoneSet>(set);
121 
122  if (action == topoSetSource::ADD || action == topoSetSource::NEW)
123  {
124  if (verbose_)
125  {
126  Info<< " Adding all faces from faceSet " << setName_
127  << " ..." << endl;
128  }
129 
130  // Load the sets
131  faceSet loadedSet(mesh_, setName_);
132  labelHashSet& faceIds = loadedSet;
133 
134  // Start off from copy
135  DynamicList<label> newAddressing(zoneSet.addressing());
136  DynamicList<bool> newFlipMap(zoneSet.flipMap());
137 
138  const faceList& faces = mesh_.faces();
139  const pointField& points = mesh_.points();
140 
141  for (const label facei : faceIds)
142  {
143  if (!zoneSet.found(facei))
144  {
145  newAddressing.append(facei);
146 
147  const vector n = faces[facei].areaNormal(points);
148  if ((n & normal_) > 0)
149  {
150  newFlipMap.append(false);
151  }
152  else
153  {
154  newFlipMap.append(true);
155  }
156  }
157  }
158 
159  zoneSet.addressing().transfer(newAddressing);
160  zoneSet.flipMap().transfer(newFlipMap);
161  zoneSet.updateSet();
162  }
163  else if (action == topoSetSource::SUBTRACT)
164  {
165  if (verbose_)
166  {
167  Info<< " Removing all faces from faceSet " << setName_
168  << " ..." << endl;
169  }
170 
171  // Load the set
172  faceSet loadedSet(mesh_, setName_);
173 
174  // Start off empty
175  DynamicList<label> newAddressing(zoneSet.addressing().size());
176  DynamicList<bool> newFlipMap(zoneSet.flipMap().size());
177 
178  forAll(zoneSet.addressing(), i)
179  {
180  if (!loadedSet.found(zoneSet.addressing()[i]))
181  {
182  newAddressing.append(zoneSet.addressing()[i]);
183  newFlipMap.append(zoneSet.flipMap()[i]);
184  }
185  }
186  zoneSet.addressing().transfer(newAddressing);
187  zoneSet.flipMap().transfer(newFlipMap);
188  zoneSet.updateSet();
189  }
190  }
191 }
192 
193 
194 // ************************************************************************* //
Foam::addToRunTimeSelectionTable
addToRunTimeSelectionTable(decompositionMethod, kahipDecomp, dictionary)
Foam::faceZoneSet::addressing
const labelList & addressing() const
Definition: faceZoneSet.H:106
Foam::BitOps::set
void set(List< bool > &bools, const labelRange &range)
Set the specified range 'on' in a boolList.
Definition: BitOps.C:37
Foam::word
A class for handling words, derived from Foam::string.
Definition: word.H:65
Foam::DynamicList< label >
faceZoneSet.H
Foam::faceZoneSet::updateSet
void updateSet()
Sort addressing and make faceSet part consistent with addressing.
Definition: faceZoneSet.C:51
Foam::topoSetSource::addToUsageTable
Class with constructor to add usage string to table.
Definition: topoSetSource.H:129
Foam::topoSetFaceZoneSource
The topoSetFaceZoneSource is a intermediate class for handling topoSet sources for selecting face zon...
Definition: topoSetFaceZoneSource.H:57
Foam::faceSet
A list of face labels.
Definition: faceSet.H:51
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:369
Foam::topoSetSource::setAction
setAction
Enumeration defining the valid actions.
Definition: topoSetSource.H:100
polyMesh.H
Foam::HashSet< label, Hash< label > >
Foam::faceZoneSet::flipMap
const boolList & flipMap() const
Definition: faceZoneSet.H:117
Foam::polyMesh
Mesh consisting of general polyhedral cells.
Definition: polyMesh.H:77
forAll
#define forAll(list, i)
Loop across all elements in list.
Definition: stdFoam.H:296
n
label n
Definition: TABSMDCalcMethod2.H:31
Foam::topoSet::found
virtual bool found(const label id) const
Has the given index?
Definition: topoSet.C:508
Foam::Field< vector >
Foam::Istream
An Istream is an abstract base class for all input systems (streams, files, token lists etc)....
Definition: Istream.H:61
Foam::Info
messageStream Info
Information stream (stdout output on master, null elsewhere)
Foam::DynamicList::append
DynamicList< T, SizeMin > & append(const T &val)
Append an element to the end of this list.
Definition: DynamicListI.H:511
setAndNormalToFaceZone.H
Foam::List::transfer
void transfer(List< T > &list)
Definition: List.C:456
Foam::faceZoneSet
Like faceSet but -reads data from faceZone -updates faceZone when writing.
Definition: faceZoneSet.H:51
Foam::topoSet
General set of labels of mesh quantity (points, cells, faces).
Definition: topoSet.H:63
dict
dictionary dict
Definition: searchingEngine.H:14
Foam::dictionary
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Definition: dictionary.H:123
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:6
addToRunTimeSelectionTable.H
Macros for easy insertion into run-time selection tables.
Foam
Namespace for OpenFOAM.
Definition: atmBoundaryLayer.C:33
Foam::setAndNormalToFaceZone::applyToSet
virtual void applyToSet(const topoSetSource::setAction action, topoSet &set) const
Apply specified action to the topoSet.
Definition: setAndNormalToFaceZone.C:107
Foam::setAndNormalToFaceZone::setAndNormalToFaceZone
setAndNormalToFaceZone(const polyMesh &mesh, const word &setName, const vector &normal)
Construct from components.
Definition: setAndNormalToFaceZone.C:68
Foam::Vector< scalar >
Foam::List< face >
points
const pointField & points
Definition: gmvOutputHeader.H:1
Foam::defineTypeNameAndDebug
defineTypeNameAndDebug(combustionModel, 0)
WarningInFunction
#define WarningInFunction
Report a warning using Foam::Warning.
Definition: messageStream.H:328