faceSet.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) 2011-2016 OpenFOAM Foundation
9  Copyright (C) 2016-2018 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 "faceSet.H"
30 #include "polyMesh.H"
31 #include "mapPolyMesh.H"
32 #include "syncTools.H"
33 #include "mapDistributePolyMesh.H"
35 
36 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
37 
38 namespace Foam
39 {
40  defineTypeNameAndDebug(faceSet, 0);
41  addToRunTimeSelectionTable(topoSet, faceSet, word);
42  addToRunTimeSelectionTable(topoSet, faceSet, size);
43  addToRunTimeSelectionTable(topoSet, faceSet, set);
44 }
45 
46 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
47 
49 :
50  topoSet(obj, typeName)
51 {}
52 
53 
55 (
56  const polyMesh& mesh,
57  const word& name,
58  readOption r,
59  writeOption w
60 )
61 :
62  topoSet(mesh, typeName, name, r, w)
63 {
64  check(mesh.nFaces());
65 }
66 
67 
69 (
70  const polyMesh& mesh,
71  const word& name,
72  const label size,
73  writeOption w
74 )
75 :
76  topoSet(mesh, name, size, w)
77 {}
78 
79 
81 (
82  const polyMesh& mesh,
83  const word& name,
84  const topoSet& set,
85  writeOption w
86 )
87 :
88  topoSet(mesh, name, set, w)
89 {}
90 
91 
93 (
94  const polyMesh& mesh,
95  const word& name,
96  const labelHashSet& labels,
97  writeOption w
98 )
99 :
100  topoSet(mesh, name, labels, w)
101 {}
102 
103 
105 (
106  const polyMesh& mesh,
107  const word& name,
108  labelHashSet&& labels,
109  writeOption w
110 )
111 :
112  topoSet(mesh, name, std::move(labels), w)
113 {}
114 
115 
117 (
118  const polyMesh& mesh,
119  const word& name,
120  const labelUList& labels,
121  writeOption w
122 )
123 :
124  topoSet(mesh, name, labels, w)
125 {}
126 
127 
128 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
129 
131 {
132  labelHashSet& labels = *this;
133 
134  // Convert to boolList
135  // TBD: could change to using bitSet for the synchronization
136 
137  const label len = mesh.nFaces();
138 
139  boolList contents(len, false);
140 
141  for (const label facei : labels)
142  {
143  contents.set(facei);
144  }
145 
147 
148 
149  // Update labelHashSet
150 
151  for (label i=0; i < len; ++i)
152  {
153  if (contents.test(i))
154  {
155  labels.set(i);
156  }
157  }
158 }
159 
160 
161 Foam::label Foam::faceSet::maxSize(const polyMesh& mesh) const
162 {
163  return mesh.nFaces();
164 }
165 
166 
168 {
169  updateLabels(morphMap.reverseFaceMap());
170 }
171 
172 
174 {
175  labelHashSet& labels = *this;
176 
177  boolList contents(map.nOldFaces(), false);
178 
179  for (const label facei : labels)
180  {
181  contents.set(facei);
182  }
183 
184  map.distributeFaceData(contents);
185 
186  // The new length
187  const label len = contents.size();
188 
189  // Count - as per BitOps::count(contents)
190  label n = 0;
191  for (label i=0; i < len; ++i)
192  {
193  if (contents.test(i))
194  {
195  ++n;
196  }
197  }
198 
199  // Update labelHashSet
200 
201  labels.clear();
202  labels.resize(2*n);
203 
204  for (label i=0; i < len; ++i)
205  {
206  if (contents.test(i))
207  {
208  labels.set(i);
209  }
210  }
211 }
212 
213 
215 (
216  Ostream& os,
217  const primitiveMesh& mesh,
218  const label maxLen
219 ) const
220 {
221  topoSet::writeDebug(os, mesh.faceCentres(), maxLen);
222 }
223 
224 
225 // ************************************************************************* //
Foam::topoSet::writeDebug
void writeDebug(Ostream &os, const label maxElem, topoSet::const_iterator &iter, label &elemI) const
Write part of contents nicely formatted. Prints labels only.
Definition: topoSet.C:221
Foam::addToRunTimeSelectionTable
addToRunTimeSelectionTable(decompositionMethod, kahipDecomp, dictionary)
Foam::IOobject
Defines the attributes of an object for which implicit objectRegistry management is supported,...
Definition: IOobject.H:104
Foam::word
A class for handling words, derived from Foam::string.
Definition: word.H:62
Foam::primitiveMesh::nFaces
label nFaces() const
Number of mesh faces.
Definition: primitiveMeshI.H:90
mapPolyMesh.H
Foam::faceSet::distribute
virtual void distribute(const mapDistributePolyMesh &map)
Update any stored data for mesh redistribution.
Definition: faceSet.C:173
Foam::orEqOp
Definition: ops.H:86
Foam::faceSet::maxSize
virtual label maxSize(const polyMesh &mesh) const
Return max index+1.
Definition: faceSet.C:161
polyMesh.H
Foam::HashSet< label, Hash< label > >
syncTools.H
Foam::polyMesh
Mesh consisting of general polyhedral cells.
Definition: polyMesh.H:77
Foam::faceSet::sync
virtual void sync(const polyMesh &mesh)
Sync faceSet across coupled patches.
Definition: faceSet.C:130
Foam::IOobject::writeOption
writeOption
Enumeration defining the write options.
Definition: IOobject.H:127
mapDistributePolyMesh.H
n
label n
Definition: TABSMDCalcMethod2.H:31
Foam::syncTools::syncFaceList
static void syncFaceList(const polyMesh &mesh, UList< T > &faceValues, const CombineOp &cop)
Synchronize values on all mesh faces.
Definition: syncTools.H:390
Foam::faceSet::writeDebug
virtual void writeDebug(Ostream &os, const primitiveMesh &, const label maxLen) const
Write maxLen items with label and coordinates.
Definition: faceSet.C:215
Foam::name
word name(const complex &c)
Return string representation of complex.
Definition: complex.C:76
faceSet.H
Foam::topoSet
General set of labels of mesh quantity (points, cells, faces).
Definition: topoSet.H:63
Foam::faceSet::faceSet
faceSet(const IOobject &obj)
Construct from IOobject.
Definition: faceSet.C:48
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:6
addToRunTimeSelectionTable.H
Macros for easy insertion into run-time selection tables.
Foam::mapDistributePolyMesh::distributeFaceData
void distributeFaceData(List< T > &lst) const
Distribute list of face data.
Definition: mapDistributePolyMesh.H:249
Foam
Namespace for OpenFOAM.
Definition: atmBoundaryLayer.C:33
Foam::mapPolyMesh::reverseFaceMap
const labelList & reverseFaceMap() const
Reverse face map.
Definition: mapPolyMesh.H:500
Foam::List< bool >
Foam::primitiveMesh::faceCentres
const vectorField & faceCentres() const
Definition: primitiveMeshFaceCentresAndAreas.C:144
Foam::UList< label >
Foam::List::set
std::enable_if< std::is_same< bool, TypeT >::value, bool >::type set(const label i, bool val=true)
A bitSet::set() method for a list of bool.
Definition: List.H:325
Foam::mapPolyMesh
Class containing mesh-to-mesh mapping information after a change in polyMesh topology.
Definition: mapPolyMesh.H:160
Foam::mapDistributePolyMesh::nOldFaces
label nOldFaces() const
Number of faces in mesh before distribution.
Definition: mapDistributePolyMesh.H:181
Foam::faceSet::updateMesh
virtual void updateMesh(const mapPolyMesh &morphMap)
Update any stored data for new labels.
Definition: faceSet.C:167
Foam::IOobject::readOption
readOption
Enumeration defining the read options.
Definition: IOobject.H:118
Foam::Ostream
An Ostream is an abstract base class for all output systems (streams, files, token lists,...
Definition: Ostream.H:56
Foam::mapDistributePolyMesh
Class containing mesh-to-mesh mapping information after a mesh distribution where we send parts of me...
Definition: mapDistributePolyMesh.H:66
Foam::defineTypeNameAndDebug
defineTypeNameAndDebug(combustionModel, 0)
Foam::primitiveMesh
Cell-face mesh analysis engine.
Definition: primitiveMesh.H:78