src/mesh/cfMesh/utilities/triSurfaceTools/triSurfaceCurvatureEstimator/triSurfaceCurvatureEstimator.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | cfMesh: A library for mesh generation
   4    \\    /   O peration     |
   5     \\  /    A nd           | Author: Franjo Juretic (franjo.juretic@c-fields.com)
   6      \\/     M anipulation  | Copyright (C) Creative Fields, Ltd.
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of cfMesh.
  10
  11     cfMesh 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 3 of the License, or (at your
  14     option) any later version.
  15
  16     cfMesh 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 cfMesh.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Description
  25
  26 \*---------------------------------------------------------------------------*/
  27
  28 #include "triSurfaceCurvatureEstimator.H"
  29 #include "demandDrivenData.H"
  30
  31 //#define DEBUGMorph
  32
  33 # ifdef DEBUGMorph
  34 #include <sstream>
  35 # endif
  36
  37 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  38
  39 namespace Foam
  40 {
  41
  42 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  43
  44 triSurfaceCurvatureEstimator::triSurfaceCurvatureEstimator
  45 (
  46     const triSurf& surface
  47 )
  48 :
  49     surface_(surface),
  50     edgePointCurvature_(),
  51     patchPositions_(),
  52     gaussianCurvature_(),
  53     meanCurvature_(),
  54     maxCurvature_(),
  55     minCurvature_(),
  56     maxCurvatureVector_(),
  57     minCurvatureVector_()
  58 {
  59     calculateEdgeCurvature();
  60     calculateSurfaceCurvatures();
  61     //calculateGaussianCurvature();
  62     //calculateMeanCurvature();
  63     //calculateMinAndMaxCurvature();
  64 }
  65
  66 triSurfaceCurvatureEstimator::~triSurfaceCurvatureEstimator()
  67 {}
  68
  69 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  70
  71 scalar triSurfaceCurvatureEstimator::edgePointCurvature(const label pI) const
  72 {
  73     return edgePointCurvature_[pI];
  74 }
  75
  76 scalar triSurfaceCurvatureEstimator::curvatureAtEdge(const label edgeI) const
  77 {
  78     const edge& edg = surface_.edges()[edgeI];
  79
  80     const scalar k1 = edgePointCurvature_[edg.start()];
  81     const scalar k2 = edgePointCurvature_[edg.end()];
  82
  83     return (k1 + k2) / 2.0;
  84 }
  85
  86 scalar triSurfaceCurvatureEstimator::gaussianCurvatureAtTriangle
  87 (
  88     const label triI
  89 ) const
  90 {
  91     const labelledTri& tri = surface_[triI];
  92
  93     scalar curv(0.0);
  94     for(label i=0;i<3;++i)
  95         curv += gaussianCurvature_[tri[i]][patchPositions_(triI, i)];
  96     curv /= 3.0;
  97
  98     return curv;
  99 }
 100
 101 scalar triSurfaceCurvatureEstimator::meanCurvatureAtTriangle
 102 (
 103     const label triI
 104 ) const
 105 {
 106     const labelledTri& tri = surface_[triI];
 107
 108     scalar curv(0.0);
 109     for(label i=0;i<3;++i)
 110         curv += meanCurvature_[tri[i]][patchPositions_(triI, i)];
 111     curv /= 3.0;
 112
 113     return curv;
 114 }
 115
 116 scalar triSurfaceCurvatureEstimator::maxCurvatureAtTriangle
 117 (
 118     const label triI
 119 ) const
 120 {
 121     const labelledTri& tri = surface_[triI];
 122
 123     scalar curv(0.0);
 124     for(label i=0;i<3;++i)
 125         curv += maxCurvature_[tri[i]][patchPositions_(triI, i)];
 126     curv /= 3.0;
 127
 128     return curv;
 129 }
 130
 131 scalar triSurfaceCurvatureEstimator::minCurvatureAtTriangle
 132 (
 133     const label triI
 134 ) const
 135 {
 136     const labelledTri& tri = surface_[triI];
 137
 138     scalar curv(0.0);
 139     for(label i=0;i<3;++i)
 140         curv += minCurvature_[tri[i]][patchPositions_(triI, i)];
 141     curv /= 3.0;
 142
 143     return curv;
 144 }
 145
 146 vector triSurfaceCurvatureEstimator::maxCurvatureVectorAtTriangle
 147 (
 148     const label triI
 149 ) const
 150 {
 151     const labelledTri& tri = surface_[triI];
 152
 153     vector curv(vector::zero);
 154     for(label i=0;i<3;++i)
 155         curv += maxCurvatureVector_[tri[i]][patchPositions_(triI, i)];
 156     curv /= 3.0;
 157
 158     return curv;
 159 }
 160
 161 vector triSurfaceCurvatureEstimator::minCurvatureVectorAtTriangle
 162 (
 163     const label triI
 164 ) const
 165 {
 166     const labelledTri& tri = surface_[triI];
 167
 168     vector curv(vector::zero);
 169     for(label i=0;i<3;++i)
 170         curv += minCurvatureVector_[tri[i]][patchPositions_(triI, i)];
 171     curv /= 3.0;
 172
 173     return curv;
 174 }
 175
 176 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 177
 178 } // End namespace Foam
 179
 180 // ************************************************************************* //