src/meshTools/octree/octreeDataEdges.C

   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 \*---------------------------------------------------------------------------*/
  26
  27 #include "octreeDataEdges.H"
  28
  29 #include "line.H"
  30 #include "labelList.H"
  31 #include "octree.H"
  32 #include "linePointRef.H"
  33 #include "pointHit.H"
  34
  35 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  36
  37 defineTypeNameAndDebug(Foam::octreeDataEdges, 0);
  38
  39 Foam::scalar Foam::octreeDataEdges::tol = 1E-6;
  40
  41
  42 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  43
  44
  45 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  46
  47 // Construct from selected edges. Bounding box calculated.
  48 Foam::octreeDataEdges::octreeDataEdges
  49 (
  50     const edgeList& edges,
  51     const pointField& points,
  52     const labelList& edgeLabels
  53 )
  54 :
  55     edges_(edges),
  56     points_(points),
  57     edgeLabels_(edgeLabels),
  58     allBb_(edgeLabels_.size())
  59 {
  60     // Generate tight fitting bounding box
  61     forAll(edgeLabels_, i)
  62     {
  63         label edgeI = edgeLabels_[i];
  64
  65         const edge& e = edges_[edgeI];
  66
  67         const point& a = points_[e.start()];
  68         const point& b = points_[e.end()];
  69
  70         allBb_[i].min() = min(a, b);
  71         allBb_[i].max() = max(a, b);
  72     }
  73 }
  74
  75
  76 // Construct as copy
  77 Foam::octreeDataEdges::octreeDataEdges(const octreeDataEdges& shapes)
  78 :
  79     edges_(shapes.edges()),
  80     points_(shapes.points()),
  81     edgeLabels_(shapes.edgeLabels()),
  82     allBb_(shapes.allBb())
  83 {}
  84
  85
  86 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  87
  88 Foam::octreeDataEdges::~octreeDataEdges()
  89 {}
  90
  91
  92 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  93
  94 Foam::label Foam::octreeDataEdges::getSampleType
  95 (
  96     const octree<octreeDataEdges>&,
  97     const point&
  98 ) const
  99 {
 100     return octree<octreeDataEdges>::UNKNOWN;
 101 }
 102
 103
 104 bool Foam::octreeDataEdges::overlaps
 105 (
 106     const label index,
 107     const treeBoundBox& sampleBb
 108 ) const
 109 {
 110     return sampleBb.intersects(allBb_[index]);
 111 }
 112
 113
 114 bool Foam::octreeDataEdges::contains
 115 (
 116     const label,
 117     const point&
 118 ) const
 119 {
 120     notImplemented
 121     (
 122         "octreeDataEdges::contains(const label, const point&)"
 123     );
 124     return false;
 125 }
 126
 127
 128 bool Foam::octreeDataEdges::intersects
 129 (
 130     const label,
 131     const point&,
 132     const point&,
 133     point&
 134 ) const
 135 {
 136     notImplemented
 137     (
 138         "octreeDataEdges::intersects(const label, const point&"
 139         ", const point&, point&)"
 140     );
 141     return false;
 142 }
 143
 144
 145 bool Foam::octreeDataEdges::findTightest
 146 (
 147     const label index,
 148     const point& sample,
 149     treeBoundBox& tightest
 150 ) const
 151 {
 152     // Get nearest and furthest away vertex
 153     point myNear, myFar;
 154     allBb_[index].calcExtremities(sample, myNear, myFar);
 155
 156     const point dist = myFar - sample;
 157     scalar myFarDist = mag(dist);
 158
 159     point tightestNear, tightestFar;
 160     tightest.calcExtremities(sample, tightestNear, tightestFar);
 161
 162     scalar tightestFarDist = mag(tightestFar - sample);
 163
 164     if (tightestFarDist < myFarDist)
 165     {
 166         // Keep current tightest.
 167         return false;
 168     }
 169     else
 170     {
 171         // Construct bb around sample and myFar
 172         const point dist2(fabs(dist.x()), fabs(dist.y()), fabs(dist.z()));
 173
 174         tightest.min() = sample - dist2;
 175         tightest.max() = sample + dist2;
 176
 177         return true;
 178     }
 179 }
 180
 181
 182 // Determine numerical value of sign of sample compared to shape at index
 183 Foam::scalar Foam::octreeDataEdges::calcSign
 184 (
 185     const label,
 186     const point&,
 187     point& n
 188 ) const
 189 {
 190     n = vector::zero;
 191
 192     return 1;
 193 }
 194
 195
 196 // Calculate nearest point on/in shapei
 197 Foam::scalar Foam::octreeDataEdges::calcNearest
 198 (
 199     const label index,
 200     const point& sample,
 201     point& nearest
 202 ) const
 203 {
 204     const edge& e = edges_[edgeLabels_[index]];
 205
 206     pointHit nearHit = e.line(points_).nearestDist(sample);
 207
 208     nearest = nearHit.rawPoint();
 209
 210     return nearHit.distance();
 211 }
 212
 213
 214 // Calculate nearest point on/in shapei
 215 Foam::scalar Foam::octreeDataEdges::calcNearest
 216 (
 217     const label index,
 218     const linePointRef& sampleLine,
 219     point& sampleLinePt,
 220     point& shapePt
 221 ) const
 222 {
 223     const edge& e = edges_[edgeLabels_[index]];
 224
 225     linePointRef edgeLine(e.line(points_));
 226
 227     return edgeLine.nearestDist(sampleLine, shapePt, sampleLinePt);
 228 }
 229
 230
 231 void Foam::octreeDataEdges::write
 232 (
 233     Ostream& os,
 234     const label index
 235 ) const
 236 {
 237     os << edgeLabels_[index] << " " << allBb_[index];
 238 }
 239
 240
 241 // ************************************************************************* //