src/mesh/cfMesh/utilities/tetrahedra/tetCreatorOctree/tetCreatorOctree.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | cfMesh: A library for mesh generation
   4    \\    /   O peration     |
   5     \\  /    A nd           | Author: Franjo Juretic (franjo.juretic@c-fields.com)
   6      \\/     M anipulation  | Copyright (C) Creative Fields, Ltd.
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of cfMesh.
  10
  11     cfMesh 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 3 of the License, or (at your
  14     option) any later version.
  15
  16     cfMesh 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 cfMesh.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Class
  25     tetCreatorOctree
  26
  27 Description
  28     A class which creates tets from the octree structure
  29
  30 SourceFiles
  31     tetCreatorOctree.C
  32
  33 \*---------------------------------------------------------------------------*/
  34
  35 #ifndef tetCreatorOctree_H
  36 #define tetCreatorOctree_H
  37
  38 #include "polyMeshGenModifier.H"
  39 #include "partTet.H"
  40 #include "meshOctreeAddressing.H"
  41 #include "meshOctreeCubeCoordinates.H"
  42 #include "VRWGraph.H"
  43 #include "labelLongList.H"
  44
  45 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  46
  47 namespace Foam
  48 {
  49
  50 /*---------------------------------------------------------------------------*\
  51                     Class tetCreatorOctree Declaration
  52 \*---------------------------------------------------------------------------*/
  53
  54 class tetCreatorOctree
  55 {
  56     // Private data
  57         //- helper for searching coordinates of cubes around an edge
  58         static const meshOctreeCubeCoordinates edgeCoordinates_[12][4];
  59
  60         //- helper for finding face centres of cubes sharing an edge
  61         static const label faceCentreHelper_[3][4];
  62
  63         //- reference to the octree
  64         meshOctreeAddressing octreeCheck_;
  65
  66         //- points of the tetrahedrisation
  67         LongList<point> tetPoints_;
  68
  69         //- tetrahedra
  70         LongList<partTet> tets_;
  71
  72         //- octree leaves sorted according to their level
  73         List<labelLongList> sortedLeaves_;
  74
  75         //- node labels of vertices created inside split-hex boxes
  76         VRWGraph* subNodeLabelsPtr_;
  77
  78         //- cube centre label
  79         labelList* cubeLabelPtr_;
  80
  81         //- cube face label
  82         VRWGraph* faceCentreLabelPtr_;
  83
  84         //- are tets created or not
  85         bool created_;
  86
  87     // Private member functions
  88         //- deletes all pointer data
  89         void clearOut();
  90
  91         //- find elements which will be used as mesh cells
  92         void selectElements();
  93
  94     // Private member functions for creating tetrahedra
  95         //- create tetPoints_ and necessary addressing
  96         void createPointsAndAddressing();
  97
  98         //- create faceCentreLabelPtr_
  99         void createFaceCentreLabels();
 100
 101         //- create tetrahedra from faces, owner and neighbour
 102         void createTetsAroundEdges();
 103
 104         //- create tetrahedra from split faces
 105         void createTetsFromSplitFaces();
 106
 107         //- create tetrahedra from faces which share at least one edge
 108         //- with a refined cube
 109         void createTetsFromFacesWithCentreNode();
 110
 111         //- create tetrahedra from faces with split edges
 112         void createTetsAroundSplitEdges();
 113
 114         //- helper funtion which checks validity of a created tet and appends
 115         //- it if valid
 116         inline void checkAndAppendTet(const partTet);
 117
 118         //- function containing the workflow
 119         void createTets();
 120
 121     // Private copy constructor
 122         //- Disallow default bitwise copy construct
 123         tetCreatorOctree(const tetCreatorOctree&);
 124
 125         //- Disallow default bitwise assignment
 126         void operator=(const tetCreatorOctree&);
 127
 128 public:
 129
 130     // Constructors
 131
 132         //- Construct from octree and IOdictionary
 133         tetCreatorOctree
 134         (
 135             const meshOctree& octree,
 136             const IOdictionary& meshDict
 137         );
 138
 139     // Destructor
 140
 141         ~tetCreatorOctree();
 142
 143
 144     // Member Functions
 145         const LongList<point>& tetPoints() const
 146         {
 147             if( !created_ )
 148                 FatalErrorIn
 149                 (
 150                     "const LongList<point>& tetPoints() const"
 151                 ) << "Tets are not created!" << exit(FatalError);
 152
 153             return tetPoints_;
 154         }
 155
 156         const LongList<partTet>& tets() const
 157         {
 158             if( !created_ )
 159                 FatalErrorIn
 160                 (
 161                     "const LongList<point>& tets() const"
 162                 ) << "Tets are not created!" << exit(FatalError);
 163
 164             return tets_;
 165         }
 166 };
 167
 168
 169 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 170
 171 } // End namespace Foam
 172
 173 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 174
 175 #include "tetCreatorOctreeI.H"
 176
 177 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 178
 179 #endif
 180
 181 // ************************************************************************* //