Initial commit for version 2.0.x patch release
[OpenFOAM-2.0.x.git] / src / OpenFOAM / meshes / primitiveShapes / triangle / intersection.H
blob168011a6e13c88afc415a9b434ee45d0e87438e7
1 /*---------------------------------------------------------------------------*\
2   =========                 |
3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
4    \\    /   O peration     |
5     \\  /    A nd           | Copyright (C) 2004-2010 OpenCFD Ltd.
6      \\/     M anipulation  |
7 -------------------------------------------------------------------------------
8 License
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
19     for more details.
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/>.
24 Class
25     Foam::intersection
27 Description
28     Foam::intersection
30 SourceFiles
31    intersection.C
33 \*---------------------------------------------------------------------------*/
35 #ifndef intersection_H
36 #define intersection_H
38 #include "scalar.H"
39 #include "NamedEnum.H"
41 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
43 namespace Foam
46 /*---------------------------------------------------------------------------*\
47                         Class intersection Declaration
48 \*---------------------------------------------------------------------------*/
50 class intersection
52 private:
54     // Static data
56         //- Relative tolerance for point in a plane.
57         //  The tolerance is relative to the object size.
58         //  Used to determine if a half-ray from a point in plane of triangle
59         // intersects the triangle or not.
60         static scalar planarTol_;
62 public:
64     enum direction
65     {
66         VECTOR,
67         CONTACT_SPHERE
68     };
70     enum algorithm
71     {
72         FULL_RAY,        // Intersecting with full ray
73         HALF_RAY,        // Intersecting with half ray
74         VISIBLE          // Intersecting with the visible side
75     };
77     // Static Member Functions
79         //- Direction names
80         static const NamedEnum<direction, 2> directionNames_;
82         //- Projection algorithm names
83         static const NamedEnum<algorithm, 3> algorithmNames_;
85         //- Return planar tolerance
86         static scalar planarTol()
87         {
88             return planarTol_;
89         }
91         //- Set the planar tolerance, returning the previous value
92         static scalar setPlanarTol(const scalar t)
93         {
94             if (t < -VSMALL)
95             {
96                 FatalErrorIn
97                 (
98                     "scalar intersection::setPlanarTol(const scalar t)"
99                 )   << "Negative planar tolerance.  This is not allowed."
100                     << abort(FatalError);
101             }
103             scalar oldTol = planarTol_;
104             planarTol_ = t;
106             return oldTol;
107         }
111 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
113 } // End namespace Foam
115 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
117 #endif
119 // ************************************************************************* //