src/mesh/blockMesh/curvedEdges/curvedEdge.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-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 \*---------------------------------------------------------------------------*/
  25
  26 #include "error.H"
  27 #include "curvedEdge.H"
  28
  29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33     defineTypeNameAndDebug(curvedEdge, 0);
  34     defineRunTimeSelectionTable(curvedEdge, Istream);
  35 }
  36
  37
  38 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  39
  40 Foam::curvedEdge::curvedEdge
  41 (
  42     const pointField& points,
  43     const label start,
  44     const label end
  45 )
  46 :
  47     points_(points),
  48     start_(start),
  49     end_(end)
  50 {}
  51
  52
  53 Foam::curvedEdge::curvedEdge(const pointField& points, Istream& is)
  54 :
  55     points_(points),
  56     start_(readLabel(is)),
  57     end_(readLabel(is))
  58 {}
  59
  60
  61 Foam::curvedEdge::curvedEdge(const curvedEdge& c)
  62 :
  63     points_(c.points_),
  64     start_(c.start_),
  65     end_(c.end_)
  66 {}
  67
  68
  69 Foam::autoPtr<Foam::curvedEdge> Foam::curvedEdge::clone() const
  70 {
  71     notImplemented("curvedEdge::clone() const");
  72     return autoPtr<curvedEdge>(NULL);
  73 }
  74
  75
  76 Foam::autoPtr<Foam::curvedEdge> Foam::curvedEdge::New
  77 (
  78     const pointField& points,
  79     Istream& is
  80 )
  81 {
  82     if (debug)
  83     {
  84         Info<< "curvedEdge::New(const pointField&, Istream&) : "
  85             << "constructing curvedEdge"
  86             << endl;
  87     }
  88
  89     const word edgeType(is);
  90
  91     IstreamConstructorTable::iterator cstrIter =
  92         IstreamConstructorTablePtr_->find(edgeType);
  93
  94     if (cstrIter == IstreamConstructorTablePtr_->end())
  95     {
  96         FatalErrorIn("curvedEdge::New(const pointField&, Istream&)")
  97             << "Unknown curvedEdge type "
  98             << edgeType << nl << nl
  99             << "Valid curvedEdge types are" << endl
 100             << IstreamConstructorTablePtr_->sortedToc()
 101             << abort(FatalError);
 102     }
 103
 104     return autoPtr<curvedEdge>(cstrIter()(points, is));
 105 }
 106
 107
 108 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 109
 110 Foam::pointField Foam::curvedEdge::appendEndPoints
 111 (
 112     const pointField& points,
 113     const label start,
 114     const label end,
 115     const pointField& otherKnots
 116 )
 117 {
 118     pointField allKnots(otherKnots.size() + 2);
 119
 120     // start/end knots
 121     allKnots[0] = points[start];
 122     allKnots[otherKnots.size() + 1] = points[end];
 123
 124     // intermediate knots
 125     forAll(otherKnots, knotI)
 126     {
 127         allKnots[knotI+1] = otherKnots[knotI];
 128     }
 129
 130     return allKnots;
 131 }
 132
 133
 134 // * * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * //
 135
 136 void Foam::curvedEdge::operator=(const curvedEdge&)
 137 {
 138     notImplemented("void curvedEdge::operator=(const curvedEdge&)");
 139 }
 140
 141
 142 Foam::Ostream& Foam::operator<<(Ostream& os, const curvedEdge& p)
 143 {
 144     os << p.start_ << tab << p.end_ << endl;
 145
 146     return os;
 147 }
 148
 149
 150 // ************************************************************************* //