applications/utilities/surface/surfaceFind/surfaceFind.C

   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 Description
  26     Finds nearest triangle and vertex.
  27
  28 \*---------------------------------------------------------------------------*/
  29
  30 #include "triSurface.H"
  31 #include "argList.H"
  32 #include "OFstream.H"
  33
  34 #ifndef namespaceFoam
  35 #define namespaceFoam
  36     using namespace Foam;
  37 #endif
  38
  39
  40 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  41 // Main program:
  42
  43 int main(int argc, char *argv[])
  44 {
  45     argList::noParallel();
  46     argList::validArgs.clear();
  47     argList::validOptions.insert("x", "X");
  48     argList::validOptions.insert("y", "Y");
  49     argList::validOptions.insert("z", "Z");
  50
  51     argList::validArgs.append("surface file");
  52
  53     argList args(argc, argv);
  54
  55     point samplePt
  56     (
  57         args.optionRead<scalar>("x"),
  58         args.optionRead<scalar>("y"),
  59         args.optionRead<scalar>("z")
  60     );
  61     Info<< "Looking for nearest face/vertex to " << samplePt << endl;
  62
  63
  64     Info<< "Reading surf1 ..." << endl;
  65     triSurface surf1(args.additionalArgs()[0]);
  66
  67     //
  68     // Nearest vertex
  69     //
  70
  71     const pointField& localPoints = surf1.localPoints();
  72
  73     label minIndex = -1;
  74     scalar minDist = GREAT;
  75
  76     forAll(localPoints, pointI)
  77     {
  78         const scalar dist = mag(localPoints[pointI] - samplePt);
  79         if (dist < minDist)
  80         {
  81             minDist = dist;
  82             minIndex = pointI;
  83         }
  84     }
  85
  86     Info<< "Nearest vertex:" << endl
  87         << "    index      :" << minIndex << " (in localPoints)" << endl
  88         << "    index      :" << surf1.meshPoints()[minIndex]
  89         << " (in points)" << endl
  90         << "    coordinates:" << localPoints[minIndex] << endl
  91         << endl;
  92
  93     //
  94     // Nearest face
  95     //
  96
  97     const pointField& points = surf1.points();
  98
  99     minIndex = -1;
 100     minDist = GREAT;
 101
 102     forAll(surf1, faceI)
 103     {
 104         const labelledTri& f = surf1[faceI];
 105         const point centre = f.centre(points);
 106
 107         const scalar dist = mag(centre - samplePt);
 108         if (dist < minDist)
 109         {
 110             minDist = dist;
 111             minIndex = faceI;
 112         }
 113     }
 114
 115     const labelledTri& f = surf1[minIndex];
 116
 117     Info<< "Face with nearest centre:" << endl
 118         << "    index        :" << minIndex << endl
 119         << "    centre       :" << f.centre(points) << endl
 120         << "    face         :" << f << endl
 121         << "    vertex coords:" << points[f[0]] << " "
 122         << points[f[1]] << " " << points[f[2]] << endl
 123         << endl;
 124
 125
 126     Info << "End\n" << endl;
 127
 128     return 0;
 129 }
 130
 131
 132 // ************************************************************************* //