src/meshTools/searchableSurface/searchableSphere.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 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::searchableSphere
  27
  28 Description
  29     Searching on sphere
  30
  31 SourceFiles
  32     searchableSphere.C
  33
  34 \*---------------------------------------------------------------------------*/
  35
  36 #ifndef searchableSphere_H
  37 #define searchableSphere_H
  38
  39 #include "treeBoundBox.H"
  40 #include "searchableSurface.H"
  41
  42 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  43
  44 namespace Foam
  45 {
  46
  47 // Forward declaration of classes
  48
  49 /*---------------------------------------------------------------------------*\
  50                            Class searchableSphere Declaration
  51 \*---------------------------------------------------------------------------*/
  52
  53 class searchableSphere
  54 :
  55     public searchableSurface
  56 {
  57 private:
  58
  59     // Private Member Data
  60
  61         //- Centre point
  62         const point centre_;
  63
  64         //- Radius squared
  65         const scalar radius_;
  66
  67         //- Names of regions
  68         mutable wordList regions_;
  69
  70
  71     // Private Member Functions
  72
  73         //- Find nearest point on sphere.
  74         pointIndexHit findNearest
  75         (
  76             const point& sample,
  77             const scalar nearestDistSqr
  78         ) const;
  79
  80         //- Find intersection with sphere
  81         void findLineAll
  82         (
  83             const point& start,
  84             const point& end,
  85             pointIndexHit& near,
  86             pointIndexHit& far
  87         ) const;
  88
  89
  90         //- Disallow default bitwise copy construct
  91         searchableSphere(const searchableSphere&);
  92
  93         //- Disallow default bitwise assignment
  94         void operator=(const searchableSphere&);
  95
  96
  97 public:
  98
  99     //- Runtime type information
 100     TypeName("searchableSphere");
 101
 102
 103     // Constructors
 104
 105         //- Construct from components
 106         searchableSphere(const IOobject& io, const point&, const scalar radius);
 107
 108         //- Construct from dictionary (used by searchableSurface)
 109         searchableSphere
 110         (
 111             const IOobject& io,
 112             const dictionary& dict
 113         );
 114
 115     // Destructor
 116
 117         virtual ~searchableSphere();
 118
 119
 120     // Member Functions
 121
 122         virtual const wordList& regions() const;
 123
 124         //- Whether supports volume type below
 125         virtual bool hasVolumeType() const
 126         {
 127             return true;
 128         }
 129
 130         //- Range of local indices that can be returned.
 131         virtual label size() const
 132         {
 133             return 1;
 134         }
 135
 136
 137         // Multiple point queries.
 138
 139             virtual void findNearest
 140             (
 141                 const pointField& sample,
 142                 const scalarField& nearestDistSqr,
 143                 List<pointIndexHit>&
 144             ) const;
 145
 146             virtual void findLine
 147             (
 148                 const pointField& start,
 149                 const pointField& end,
 150                 List<pointIndexHit>&
 151             ) const;
 152
 153             virtual void findLineAny
 154             (
 155                 const pointField& start,
 156                 const pointField& end,
 157                 List<pointIndexHit>&
 158             ) const;
 159
 160             //- Get all intersections in order from start to end.
 161             virtual void findLineAll
 162             (
 163                 const pointField& start,
 164                 const pointField& end,
 165                 List<List<pointIndexHit> >&
 166             ) const;
 167
 168             //- From a set of points and indices get the region
 169             virtual void getRegion
 170             (
 171                 const List<pointIndexHit>&,
 172                 labelList& region
 173             ) const;
 174
 175             //- From a set of points and indices get the normal
 176             virtual void getNormal
 177             (
 178                 const List<pointIndexHit>&,
 179                 vectorField& normal
 180             ) const;
 181
 182             //- Determine type (inside/outside/mixed) for point. unknown if
 183             //  cannot be determined (e.g. non-manifold surface)
 184             virtual void getVolumeType
 185             (
 186                 const pointField&,
 187                 List<volumeType>&
 188             ) const;
 189
 190
 191         // regIOobject implementation
 192
 193             bool writeData(Ostream&) const
 194             {
 195                 notImplemented("searchableSphere::writeData(Ostream&) const");
 196                 return false;
 197             }
 198
 199 };
 200
 201
 202 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 203
 204 } // End namespace Foam
 205
 206 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 207
 208 #endif
 209
 210 // ************************************************************************* //