src/dynamicMesh/dynamicFvMesh/dynamicTopoFvMesh/convexSetAlgorithm/tetIntersection.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright held by original author
   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     tetIntersection
  27
  28 Description
  29     Implementation of the tetrahedron intersection algorithm given in:
  30
  31     D.H. Eberly, 3D Game Engine Design: A Practical Approach to Real-time
  32     Computer Graphics, Morgan Kaufmann, 2001.
  33
  34     Geometric Tools, LLC
  35     Distributed under the Boost Software License, Version 1.0.
  36     http://www.boost.org/LICENSE_1_0.txt
  37
  38 Implemented by
  39     Sandeep Menon
  40     University of Massachusetts Amherst
  41
  42 SourceFiles
  43     tetIntersectionI.H
  44
  45 \*---------------------------------------------------------------------------*/
  46
  47 #ifndef tetIntersection_H
  48 #define tetIntersection_H
  49
  50 #include "point.H"
  51 #include "label.H"
  52 #include "Tuple2.H"
  53 #include "FixedList.H"
  54 #include "DynamicList.H"
  55
  56 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  57
  58 namespace Foam
  59 {
  60
  61 /*---------------------------------------------------------------------------*\
  62                        Class tetIntersection Declaration
  63 \*---------------------------------------------------------------------------*/
  64
  65 class tetIntersection
  66 {
  67     // Private data
  68
  69         //- Const reference to clipping tetrahedron
  70         const FixedList<point, 4>& clipTet_;
  71
  72         //- Hessian-normal plane definition
  73         typedef Tuple2<vector, scalar> hPlane;
  74
  75         FixedList<hPlane, 4> clipPlanes_;
  76
  77         //- Magnitude of clipping tetrahedron
  78         scalar clipTetMag_;
  79
  80         //- Tetrahedra used as temporaries
  81         DynamicList<FixedList<point, 4> > inside_;
  82
  83         //- All intersection tets
  84         DynamicList<FixedList<point, 4> > allTets_;
  85
  86     // Private Member Functions
  87
  88         //- Disallow default bitwise copy construct
  89         tetIntersection(const tetIntersection&);
  90
  91         //- Disallow default bitwise assignment
  92         void operator=(const tetIntersection&);
  93
  94         //- Compute clip-planes
  95         inline void computeClipPlanes();
  96
  97         //- Split and decompose
  98         inline void splitAndDecompose
  99         (
 100             const label tetPlaneIndex,
 101             FixedList<point, 4>& tetra,
 102             DynamicList<FixedList<point, 4> >& decompTets
 103         ) const;
 104
 105 public:
 106
 107     // Constructors
 108
 109         //- Construct from components
 110         inline tetIntersection(const FixedList<point, 4>& clipTet);
 111
 112
 113     // Destructor
 114
 115         inline ~tetIntersection();
 116
 117
 118     // Member Functions
 119
 120         //- Return magnitude of clipping tetrahedron
 121         inline scalar clipTetMag() const;
 122
 123         //- Evaluate for intersections against input tetrahedron
 124         inline bool evaluate(const FixedList<point, 4>& subjectTet);
 125
 126         //- Return intersections
 127         inline const DynamicList<FixedList<point, 4> >& getIntersection() const;
 128
 129         //- Evaluate and return volume / centroid
 130         inline void getVolumeAndCentre(scalar& volume, vector& centre) const;
 131 };
 132
 133 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 134
 135 } // End namespace Foam
 136
 137 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 138
 139 #include "tetIntersectionI.H"
 140
 141 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 142
 143 #endif
 144
 145 // ************************************************************************* //