FreeFOAM The Cross-Platform CFD Toolkit
octreeDataTriSurfaceTreeLeaf.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 1991-2010 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
13  the Free Software Foundation, either version 3 of the License, or
14  (at your 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, see <http://www.gnu.org/licenses/>.
23 
24 Description
25 
26 \*---------------------------------------------------------------------------*/
27 
30 
31 // * * * * * * * * * * * * * Template Specialisations * * * * * * * * * * * //
32 
33 template<>
35 (
36  const octreeDataTriSurface& shapes,
37  const point& sample,
38  treeBoundBox& tightest,
39  label& tightestI,
40  scalar& tightestDist
41 ) const
42 {
43  // Some aliases
44  const treeBoundBoxList& allBb = shapes.allBb();
45  point& min = tightest.min();
46  point& max = tightest.max();
47 
48  point nearest;
49 
50  bool changed = false;
51  forAll(indices_, i)
52  {
53  label faceI = indices_[i];
54 
55  // Quick rejection test.
56  if (tightest.overlaps(allBb[faceI]))
57  {
58  // Full calculation
59  scalar dist = shapes.calcNearest(faceI, sample, nearest);
60 
61  if (dist < tightestDist)
62  {
63  // Update bb (centered around sample, span is dist)
64  min.x() = sample.x() - dist;
65  min.y() = sample.y() - dist;
66  min.z() = sample.z() - dist;
67 
68  max.x() = sample.x() + dist;
69  max.y() = sample.y() + dist;
70  max.z() = sample.z() + dist;
71 
72  tightestI = faceI;
73  tightestDist = dist;
74 
75  changed = true;
76  }
77  }
78  }
79  return changed;
80 }
81 
82 
83 // ************************ vim: set sw=4 sts=4 et: ************************ //