1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
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.
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
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/>.
30 Constructs feature edges and feature points, which are edges/points with
31 and angle > given specification.
32 Can be asked for 'superFaces' which can be used to see if a cell is a
38 \*---------------------------------------------------------------------------*/
40 #ifndef cellFeatures_H
41 #define cellFeatures_H
44 #include "labelList.H"
48 #include "DynamicList.H"
50 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
55 // Forward declaration of classes
58 /*---------------------------------------------------------------------------*\
59 Class cellFeatures Declaration
60 \*---------------------------------------------------------------------------*/
66 const primitiveMesh& mesh_;
68 //- cos of angle between two connected faces or two connected edges on
69 // same face before edge/point is 'feature'.
75 labelHashSet featureEdge_;
77 //- (demand driven) Faces after removing internal points&edges
78 mutable faceList* facesPtr_;
80 //- New to old face mapping
81 mutable List<DynamicList<label> > faceMap_;
84 // Private Member Functions
86 bool faceAlignedEdge(const label, const label) const;
90 const Map<label>& toSuperFace,
91 const label superFaceI,
92 const label thisEdgeI,
96 bool isCellFeatureEdge(const scalar, const label) const;
101 const label superFaceI,
102 Map<label>& toSuperFace
105 void calcSuperFaces() const;
108 //- Disallow default bitwise copy construct
109 cellFeatures(const cellFeatures&);
111 //- Disallow default bitwise assignment
112 void operator=(const cellFeatures&);
118 //- Construct from cell in mesh
121 const primitiveMesh&,
122 const scalar minCos, // angle to use for feature recognition.
135 const labelHashSet& featureEdge() const
140 const faceList& faces() const
149 //- New to old faceMap. Guaranteed to be shrunk.
150 const List<DynamicList<label> >& faceMap() const
162 //- Is edge a feature edge (uniquely determined since on cell
163 // only two faces sharing edge)
164 bool isFeatureEdge(const label edgeI) const
166 return featureEdge().found(edgeI);
169 //- Are two edges connected at feature point?
170 // Is local to face since point might be seen as feature point
171 // from one face but not from another.
172 bool isFeaturePoint(const label edge0, const label edge1) const;
174 //- Is vertexI on faceI used by two edges that form feature
176 bool isFeatureVertex(const label faceI, const label vertI) const;
181 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
183 } // End namespace Foam
185 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
189 // ************************************************************************* //