applications/utilities/mesh/generation/blockMesh/curvedEdges/polyLine.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 Description
  26     polyLineEdge class : defines a curvedEdge in terms of a series of
  27     straight line segments
  28
  29 \*---------------------------------------------------------------------------*/
  30
  31 #include "error.H"
  32
  33 #include "polyLine.H"
  34
  35 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  36
  37 namespace Foam
  38 {
  39
  40 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  41
  42 // calcDistances generates the distances_ lookup table (cumulative
  43 // distance along the line) from the individual vectors to the points
  44
  45 void polyLine::calcDistances()
  46 {
  47     distances_[0] = 0.0;
  48
  49     for (label i=1; i<distances_.size(); i++)
  50     {
  51         distances_[i] =
  52             mag(controlPoints_[i] - controlPoints_[i-1])
  53           + distances_[i-1];
  54     }
  55
  56     lineLength_ = distances_[distances_.size()-1];
  57
  58     for (label i=1; i<distances_.size(); i++)
  59     {
  60         distances_[i] /= lineLength_;
  61     }
  62 }
  63
  64
  65 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  66
  67 polyLine::polyLine(const pointField& ps)
  68 :
  69     controlPoints_(ps),
  70     distances_(ps.size())
  71 {
  72     if (ps.size())
  73     {
  74         calcDistances();
  75     }
  76 }
  77
  78
  79 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  80
  81 vector polyLine::position(const scalar lambda) const
  82 {
  83     // check range of lambda
  84
  85     if (lambda < 0 || lambda > 1)
  86     {
  87         FatalErrorIn("polyLine::position(const scalar)")
  88             << "Parameter out of range, "
  89             << "lambda = " << lambda
  90             << abort(FatalError);
  91     }
  92
  93     // Quick calc of endpoints
  94
  95     if (lambda < SMALL)
  96     {
  97         return controlPoints_[0];
  98     }
  99     else if (lambda > 1 - SMALL)
 100     {
 101         return controlPoints_[controlPoints_.size()-1];
 102     }
 103
 104
 105     // search table of cumulative distance to find which linesegment we
 106     // are on
 107
 108     label i(0);
 109     do
 110     {
 111         i++;
 112     } while (distances_[i] < lambda);
 113
 114     i--;               // we overshot!
 115
 116     // construct position vector
 117     scalar offsetDist =
 118         (lambda - distances_[i])
 119        /(distances_[i+1] - distances_[i]);
 120
 121     vector offsetV = controlPoints_[i+1] - controlPoints_[i];
 122
 123     return controlPoints_[i] + offsetDist*offsetV;
 124 }
 125
 126
 127 scalar polyLine::length() const
 128 {
 129     return lineLength_;
 130 }
 131
 132
 133 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 134
 135 } // End namespace Foam
 136
 137 // ************************************************************************* //