src/conversion/polyDualMesh/polyDualMesh.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 Class
  26     Foam::polyDualMesh
  27
  28 Description
  29     Creates dual of polyMesh.
  30
  31     Every polyMesh point becomes a cell on the dual mesh
  32
  33     Every polyMesh cell and patchFace becomes a point on the dual mesh.
  34
  35
  36 SourceFiles
  37     polyDualMesh.C
  38
  39 \*---------------------------------------------------------------------------*/
  40
  41 #ifndef polyDualMesh_H
  42 #define polyDualMesh_H
  43
  44 #include "polyMesh.H"
  45 #include "labelIOList.H"
  46 #include "typeInfo.H"
  47
  48 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  49
  50 namespace Foam
  51 {
  52
  53 /*---------------------------------------------------------------------------*\
  54                            Class polyDualMesh Declaration
  55 \*---------------------------------------------------------------------------*/
  56
  57 class polyDualMesh
  58 :
  59     public polyMesh
  60 {
  61     // Private data
  62
  63         //- From polyMesh cell to my point
  64         labelIOList cellPoint_;
  65
  66         //- From polyMesh boundary face (face-mesh.nInternalFaces()) to my point
  67         labelIOList boundaryFacePoint_;
  68
  69
  70     // Private Member Functions
  71
  72         static labelList getFaceOrder
  73         (
  74             const labelList& faceOwner,
  75             const labelList& faceNeighbour,
  76             const cellList& cells,
  77             label& nInternalFaces
  78         );
  79
  80         static void getPointEdges
  81         (
  82             const primitivePatch& patch,
  83             const label faceI,
  84             const label pointI,
  85             label& e0,
  86             label& e1
  87         );
  88
  89         static labelList collectPatchSideFace
  90         (
  91             const polyPatch& patch,
  92             const label patchToDualOffset,
  93             const labelList& edgeToDualPoint,
  94             const labelList& pointToDualPoint,
  95             const label pointI,
  96
  97             label& edgeI
  98         );
  99
 100         static void collectPatchInternalFace
 101         (
 102             const polyPatch& patch,
 103             const label patchToDualOffset,
 104             const labelList& edgeToDualPoint,
 105
 106             const label pointI,
 107             const label startEdgeI,
 108
 109             labelList& dualFace2,
 110             labelList& featEdgeIndices2
 111         );
 112
 113         static void splitFace
 114         (
 115             const polyPatch& patch,
 116             const labelList& pointToDualPoint,
 117
 118             const label pointI,
 119             const labelList& dualFace,
 120             const labelList& featEdgeIndices,
 121
 122             DynamicList<face>& dualFaces,
 123             DynamicList<label>& dualOwner,
 124             DynamicList<label>& dualNeighbour,
 125             DynamicList<label>& dualRegion
 126         );
 127
 128         static void dualPatch
 129         (
 130             const polyPatch& patch,
 131             const label patchToDualOffset,
 132             const labelList& edgeToDualPoint,
 133             const labelList& pointToDualPoint,
 134
 135             const pointField& dualPoints,
 136
 137             DynamicList<face>& dualFaces,
 138             DynamicList<label>& dualOwner,
 139             DynamicList<label>& dualNeighbour,
 140             DynamicList<label>& dualRegion
 141         );
 142
 143         void calcDual
 144         (
 145             const polyMesh& mesh,
 146             const labelList& featureEdges,
 147             const labelList& featurePoints
 148         );
 149
 150
 151         //- Disallow default bitwise copy construct
 152         polyDualMesh(const polyDualMesh&);
 153
 154         //- Disallow default bitwise assignment
 155         void operator=(const polyDualMesh&);
 156
 157
 158 public:
 159
 160     //- Runtime type information
 161     ClassName("polyDualMesh");
 162
 163
 164     // Constructors
 165
 166         //- Construct from IOobject
 167         polyDualMesh(const IOobject&);
 168
 169         //- Construct from polyMesh and list of edges and points to represent.
 170         //  Feature edge and point labels are in local addressing of a patch
 171         //  over all boundary faces.
 172         polyDualMesh
 173         (
 174             const polyMesh&,
 175             const labelList& featureEdges,
 176             const labelList& featurePoints
 177         );
 178
 179         //- Construct from polyMesh and feature edge angle. Uses calcFeatures
 180         //  below to determine feature edges and points.
 181         polyDualMesh
 182         (
 183             const polyMesh&,
 184             const scalar featureCos
 185         );
 186
 187         //- Helper function to create feature edges and points based on
 188         //  feature angle and patches.
 189         static void calcFeatures
 190         (
 191             const polyMesh&,
 192             const scalar featureCos,
 193             labelList& featureEdges,
 194             labelList& featurePoints
 195         );
 196
 197
 198     // Destructor
 199
 200         ~polyDualMesh();
 201
 202
 203     // Member Functions
 204
 205         // Access
 206
 207             //- From polyMesh cell to dual point
 208             const labelIOList& cellPoint() const
 209             {
 210                 return cellPoint_;
 211             }
 212
 213             //- From polyMesh patch face to dual point
 214             const labelIOList& boundaryFacePoint() const
 215             {
 216                 return boundaryFacePoint_;
 217             }
 218 };
 219
 220 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 221
 222 } // End namespace Foam
 223
 224 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 225
 226 #endif
 227
 228 // ************************************************************************* //