FreeFOAM The Cross-Platform CFD Toolkit
faceZoneToFaceZone.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 "faceZoneToFaceZone.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(faceZoneToFaceZone, 0);
38 
39 addToRunTimeSelectionTable(topoSetSource, faceZoneToFaceZone, word);
40 
41 addToRunTimeSelectionTable(topoSetSource, faceZoneToFaceZone, istream);
42 
43 }
44 
45 
46 Foam::topoSetSource::addToUsageTable Foam::faceZoneToFaceZone::usage_
47 (
48  faceZoneToFaceZone::typeName,
49  "\n Usage: faceZoneToFaceZone <faceZone>\n\n"
50  " Select all faces in the faceZone\n\n"
51 );
52 
53 
54 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
55 
56 // Construct from components
58 (
59  const polyMesh& mesh,
60  const word& setName
61 )
62 :
63  topoSetSource(mesh),
64  setName_(setName)
65 {}
66 
67 
68 // Construct from dictionary
70 (
71  const polyMesh& mesh,
72  const dictionary& dict
73 )
74 :
75  topoSetSource(mesh),
76  setName_(dict.lookup("zone"))
77 {}
78 
79 
80 // Construct from Istream
82 (
83  const polyMesh& mesh,
84  Istream& is
85 )
86 :
87  topoSetSource(mesh),
88  setName_(checkIs(is))
89 {}
90 
91 
92 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
93 
95 {}
96 
97 
98 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
99 
101 (
102  const topoSetSource::setAction action,
103  topoSet& set
104 ) const
105 {
106  if (!isA<faceZoneSet>(set))
107  {
108  WarningIn
109  (
110  "faceZoneToFaceZone::applyToSet(const topoSetSource::setAction"
111  ", topoSet"
112  ) << "Operation only allowed on a faceZoneSet." << endl;
113  }
114  else
115  {
116  faceZoneSet& fSet = refCast<faceZoneSet>(set);
117 
118  if ((action == topoSetSource::NEW) || (action == topoSetSource::ADD))
119  {
120  Info<< " Adding all faces from faceZone " << setName_ << " ..."
121  << endl;
122 
123  // Load the set
124  faceZoneSet loadedSet(mesh_, setName_);
125 
126  DynamicList<label> newAddressing(fSet.addressing());
127  DynamicList<bool> newFlipMap(fSet.flipMap());
128 
129  forAll(loadedSet.addressing(), i)
130  {
131  if (!fSet.found(loadedSet.addressing()[i]))
132  {
133  newAddressing.append(loadedSet.addressing()[i]);
134  newFlipMap.append(loadedSet.flipMap()[i]);
135  }
136  }
137  fSet.addressing().transfer(newAddressing);
138  fSet.flipMap().transfer(newFlipMap);
139  fSet.updateSet();
140  }
141  else if (action == topoSetSource::DELETE)
142  {
143  Info<< " Removing all faces from faceZone " << setName_ << " ..."
144  << endl;
145 
146  // Load the set
147  faceZoneSet loadedSet(mesh_, setName_);
148 
149  DynamicList<label> newAddressing(fSet.addressing().size());
150  DynamicList<bool> newFlipMap(fSet.flipMap().size());
151 
152  forAll(fSet.addressing(), i)
153  {
154  if (!loadedSet.found(fSet.addressing()[i]))
155  {
156  newAddressing.append(fSet.addressing()[i]);
157  newFlipMap.append(fSet.flipMap()[i]);
158  }
159  }
160  fSet.addressing().transfer(newAddressing);
161  fSet.flipMap().transfer(newFlipMap);
162  fSet.updateSet();
163  }
164  }
165 }
166 
167 
168 // ************************ vim: set sw=4 sts=4 et: ************************ //