ENH: featureEdgeMesh : merge from cvm

src/edgeMesh/featureEdgeMesh/featureEdgeMeshI.H

+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2010 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+inline Foam::label Foam::featureEdgeMesh::convexStart() const
+{
+    return convexStart_;
+}
+
+
+inline Foam::label Foam::featureEdgeMesh::concaveStart() const
+{
+	return concaveStart_;
+}
+
+
+inline Foam::label Foam::featureEdgeMesh::mixedStart() const
+{
+	return mixedStart_;
+}
+
+
+inline Foam::label Foam::featureEdgeMesh::nonFeatureStart() const
+{
+	return nonFeatureStart_;
+}
+
+
+inline Foam::label Foam::featureEdgeMesh::externalStart() const
+{
+	return externalStart_;
+}
+
+
+inline Foam::label Foam::featureEdgeMesh::internalStart() const
+{
+	return internalStart_;
+}
+
+
+inline Foam::label Foam::featureEdgeMesh::flatStart() const
+{
+	return flatStart_;
+}
+
+
+inline Foam::label Foam::featureEdgeMesh::openStart() const
+{
+	return openStart_;
+}
+
+
+inline Foam::label Foam::featureEdgeMesh::multipleStart() const
+{
+	return multipleStart_;
+}
+
+
+inline bool Foam::featureEdgeMesh::featurePoint(label ptI) const
+{
+    return ptI < nonFeatureStart_;
+}
+
+
+inline const Foam::vectorField& Foam::featureEdgeMesh::normals() const
+{
+    return normals_;
+}
+
+inline const Foam::vectorField& Foam::featureEdgeMesh::edgeDirections() const
+{
+    return edgeDirections_;
+}
+
+
+inline Foam::vector Foam::featureEdgeMesh::edgeDirection
+(
+    label edgeI,
+    label ptI
+) const
+{
+    const edge& e = edges()[edgeI];
+
+    if (ptI == e.start())
+    {
+        return edgeDirections()[edgeI];
+    }
+    else if (ptI == e.end())
+    {
+        return -edgeDirections()[edgeI];
+    }
+    else
+    {
+        FatalErrorIn("Foam::featureEdgeMesh::edgedirection")
+            << "Requested ptI " << ptI << " is not a point on the requested "
+            << "edgeI " << edgeI << ". edgeI start and end: "
+            << e.start() << " " << e.end()
+            << exit(FatalError);
+
+        return vector::zero;
+    }
+}
+
+
+inline const Foam::labelListList& Foam::featureEdgeMesh::edgeNormals() const
+{
+    return edgeNormals_;
+}
+
+
+inline Foam::vectorField Foam::featureEdgeMesh::edgeNormals
+(
+    const labelList& edgeNormIs
+) const
+{
+    vectorField norms(edgeNormIs.size());
+
+    forAll(edgeNormIs, i)
+    {
+        norms[i] = normals_[edgeNormIs[i]];
+    }
+
+    return norms;
+}
+
+
+inline Foam::vectorField Foam::featureEdgeMesh::edgeNormals(label edgeI) const
+{
+    return edgeNormals(edgeNormals_[edgeI]);
+}
+
+
+inline const Foam::labelListList&
+Foam::featureEdgeMesh::featurePointNormals() const
+{
+    return featurePointNormals_;
+}
+
+
+inline Foam::vectorField Foam::featureEdgeMesh::featurePointNormals
+(
+    label ptI
+) const
+{
+    if (!featurePoint(ptI))
+    {
+        WarningIn("vectorField featureEdgeMesh::featurePointNormals")
+            << "Requesting the normals of a non-feature point. "
+            << "Returned zero length vectorField."
+            << endl;
+
+        return vectorField(0);
+    }
+
+    labelList featPtNormIs(featurePointNormals_[ptI]);
+
+    vectorField norms(featPtNormIs.size());
+
+    forAll(featPtNormIs, i)
+    {
+        norms[i] = normals_[featPtNormIs[i]];
+    }
+
+    return norms;
+}
+
+
+inline const Foam::labelList& Foam::featureEdgeMesh::regionEdges() const
+{
+    return regionEdges_;
+}
+
+
+inline Foam::featureEdgeMesh::pointStatus Foam::featureEdgeMesh::getPointStatus
+(
+    label ptI
+) const
+{
+    if (ptI < concaveStart_)
+    {
+        return CONVEX;
+    }
+    else if (ptI < mixedStart_)
+    {
+        return CONCAVE;
+    }
+    else if (ptI < nonFeatureStart_)
+    {
+        return MIXED;
+    }
+    else
+    {
+        return NONFEATURE;
+    }
+}
+
+
+inline Foam::featureEdgeMesh::edgeStatus Foam::featureEdgeMesh::getEdgeStatus
+(
+    label edgeI
+) const
+{
+    if (edgeI < internalStart_)
+    {
+        return EXTERNAL;
+    }
+    else if (edgeI < flatStart_)
+    {
+        return INTERNAL;
+    }
+    else if (edgeI < openStart_)
+    {
+        return FLAT;
+    }
+    else if (edgeI < multipleStart_)
+    {
+        return OPEN;
+    }
+    else
+    {
+        return MULTIPLE;
+    }
+}
+
+
+// ************************************************************************* //