FreeFOAM The Cross-Platform CFD Toolkit
setToFaceZone.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) 1991-2010 OpenCFD Ltd.
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 "setToFaceZone.H"
27 #include <OpenFOAM/polyMesh.H>
28 #include <meshTools/faceZoneSet.H>
29 
31 
32 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
33 
34 namespace Foam
35 {
36 
37 defineTypeNameAndDebug(setToFaceZone, 0);
38 
39 addToRunTimeSelectionTable(topoSetSource, setToFaceZone, word);
40 
41 addToRunTimeSelectionTable(topoSetSource, setToFaceZone, istream);
42 
43 }
44 
45 
46 Foam::topoSetSource::addToUsageTable Foam::setToFaceZone::usage_
47 (
48  setToFaceZone::typeName,
49  "\n Usage: setToFaceZone <faceSet>\n\n"
50  " Select all faces in the faceSet."
51  " Sets flipMap.\n\n"
52 );
53 
54 
55 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
56 
57 // Construct from components
59 (
60  const polyMesh& mesh,
61  const word& setName
62 )
63 :
64  topoSetSource(mesh),
65  setName_(setName)
66 {}
67 
68 
69 // Construct from dictionary
71 (
72  const polyMesh& mesh,
73  const dictionary& dict
74 )
75 :
76  topoSetSource(mesh),
77  setName_(dict.lookup("faceSet"))
78 {}
79 
80 
81 // Construct from Istream
83 (
84  const polyMesh& mesh,
85  Istream& is
86 )
87 :
88  topoSetSource(mesh),
89  setName_(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  {
109  WarningIn
110  (
111  "setToFaceZone::applyToSet(const topoSetSource::setAction"
112  ", topoSet"
113  ) << "Operation only allowed on a faceZoneSet." << endl;
114  }
115  else
116  {
117  faceZoneSet& fzSet = refCast<faceZoneSet>(set);
118 
119  if ((action == topoSetSource::NEW) || (action == topoSetSource::ADD))
120  {
121  Info<< " Adding all faces from faceSet " << setName_
122  << " ..." << endl;
123 
124  // Load the sets
125  faceSet fSet(mesh_, setName_);
126 
127  // Start off from copy
128  DynamicList<label> newAddressing(fzSet.addressing());
129  DynamicList<bool> newFlipMap(fzSet.flipMap());
130 
131  forAllConstIter(faceSet, fSet, iter)
132  {
133  label faceI = iter.key();
134 
135  if (!fzSet.found(faceI))
136  {
137  newAddressing.append(faceI);
138  newFlipMap.append(false);
139  }
140  }
141 
142  fzSet.addressing().transfer(newAddressing);
143  fzSet.flipMap().transfer(newFlipMap);
144  fzSet.updateSet();
145  }
146  else if (action == topoSetSource::DELETE)
147  {
148  Info<< " Removing all faces from faceSet " << setName_
149  << " ..." << endl;
150 
151  // Load the set
152  faceSet loadedSet(mesh_, setName_);
153 
154  // Start off empty
155  DynamicList<label> newAddressing(fzSet.addressing().size());
156  DynamicList<bool> newFlipMap(fzSet.flipMap().size());
157 
158  forAll(fzSet.addressing(), i)
159  {
160  if (!loadedSet.found(fzSet.addressing()[i]))
161  {
162  newAddressing.append(fzSet.addressing()[i]);
163  newFlipMap.append(fzSet.flipMap()[i]);
164  }
165  }
166  fzSet.addressing().transfer(newAddressing);
167  fzSet.flipMap().transfer(newFlipMap);
168  fzSet.updateSet();
169  }
170  }
171 }
172 
173 
174 // ************************ vim: set sw=4 sts=4 et: ************************ //