edgeFaceCirculator.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | www.openfoam.com
     \\/     M anipulation  |
-------------------------------------------------------------------------------
    Copyright (C) 2011-2017 OpenFOAM Foundation
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.
 
    OpenFOAM is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
 
    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.
 
    You should have received a copy of the GNU General Public License
    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
Class
    Foam::edgeFaceCirculator
 
Description
    Walks from starting face around edge.
 
    Implicit description of edge:
    - face
    - index in face. edge is always between f[index] and f[index+1]
    - direction (cell to walk into)
 
    -# Use in-place: \n
        \code
            edgeFaceCirculator circ(..);
            // Optionally rotate to beginning: circ.setCanonical();
 
            // Walk
            do
            {
                Info<< "face:" << circ.face() << endl;
                ++circ;
            }
            while (circ != circ.end());
        \endcode
 
    -# Use like STL iterator: \n
        \code
            edgeFaceCirculator circ(..);
            for
            (
                edgeFaceCirculator iter = circ.begin();
                iter != circ.end();
                ++iter
            )
            {
                Info<< "face:" << iter.face() << endl;
            }
        \endcode
 
 
SourceFiles
    edgeFaceCirculator.C
 
\*---------------------------------------------------------------------------*/
 
#ifndef edgeFaceCirculator_H
#define edgeFaceCirculator_H
 
#include "face.H"
#include "ListOps.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
 
// Forward declaration of classes
class primitiveMesh;
 
/*---------------------------------------------------------------------------*\
                           Class edgeFaceCirculator Declaration
\*---------------------------------------------------------------------------*/
 
class edgeFaceCirculator
{
    // Private Member Data
 
        //- The underlying mesh pointer
        const primitiveMesh* meshPtr_;
 
        //- Is current side of face the owner side?
        bool ownerSide_;
 
        //- Is boundary edge?
        bool isBoundaryEdge_;
 
        //- Current face
        label faceLabel_;
 
        //- Edge (between index and index+1 on faces[faceLabel_]
        label index_;
 
        //- Starting face so we know when to stop. Used when circulating over
        //  internal edges.
        label startFaceLabel_;
 
 
    // Private Member Functions
 
        //- The underlying mesh
        inline const primitiveMesh& mesh() const;
 
        //- Set to end() iterator
        inline void setEnd();
 
        //- Check and set faceLabel_ and ownerSide_
        inline void setFace(const label facei, const label celli);
 
        //- Set faceLabel_ to be the other face on the cell that uses the edge.
        inline void otherFace(const label celli);
 
        //- Default construct - this is also an end iterator
        inline edgeFaceCirculator();
 
 
public:
 
    // Constructors
 
        //- Construct from components
        inline edgeFaceCirculator
        (
            const primitiveMesh& mesh,
            const label faceLabel,
            const bool ownerSide,
            const label index,
            const bool isBoundaryEdge
        );
 
        //- Construct as copy
        inline edgeFaceCirculator(const edgeFaceCirculator& circ);
 
 
    // Member Functions
 
        //- Helper: find index in face of edge or -1. Index is such that edge is
        //  between f[index] and f[index+1]
        inline static label getMinIndex
        (
            const face& f,
            const label v0,
            const label v1
        );
 
        //- Return the face label, -1 for end iterator
        inline label faceLabel() const;
 
        //- Return true if the face label corresponds to an internal face
        inline bool isInternalFace() const;
 
        inline bool ownerSide() const;
 
        inline label index() const;
 
        //- Helper: get the neighbouring cell according to the ownerSide.
        //  Returns -1 if on neighbourside of boundary face.
        inline label cellLabel() const;
 
        //- Helper: return true if normal of generated face points along
        //  edge from v0 to v1. (v0 and v1 have to be on edge)
        inline bool sameOrder(const label v0, const label v1) const;
 
        //- Set edge to a unique state so different ones can be compared.
        //  Internal edge: minimum face index.
        //  Boundary edge: walk back until boundary face.
        inline void setCanonical();
 
 
    // Member Operators
 
        inline void operator=(const edgeFaceCirculator& iter);
 
        inline bool operator==(const edgeFaceCirculator& iter) const;
 
        inline bool operator!=(const edgeFaceCirculator& iter) const;
 
        //- Step to next face. Uses no edge addressing!
        inline edgeFaceCirculator& operator++();
 
        //- Iterator set to the beginning face. For internal edges this is
        //  the current face. For boundary edges this is the first boundary face
        //  reached from walking back (i.e. in opposite direction to ++)
        inline edgeFaceCirculator begin() const;
        inline edgeFaceCirculator cbegin() const;
 
        //- Iterator set to beyond the end of the walk.
        inline const edgeFaceCirculator end() const;
        inline const edgeFaceCirculator cend() const;
};
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
} // End namespace Foam
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#include "edgeFaceCirculatorI.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif
 
// ************************************************************************* //