src/meshTools/cellDist/wallPoint/wallPointDataI.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 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
  27 \*---------------------------------------------------------------------------*/
  28
  29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33
  34 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  35
  36
  37 // Update this with w2 if w2 nearer to pt.
  38 template <class Type>
  39 inline bool wallPointData<Type>::update
  40 (
  41     const point& pt,
  42     const wallPointData<Type>& w2,
  43     const scalar tol
  44 )
  45 {
  46     scalar dist2 = magSqr(pt - w2.origin());
  47
  48     if (!valid())
  49     {
  50         // current not yet set so use any value
  51         distSqr() = dist2;
  52         origin() = w2.origin();
  53         data_ = w2.data();
  54
  55         return true;
  56     }
  57
  58     scalar diff = distSqr() - dist2;
  59
  60     if (diff < 0)
  61     {
  62         // already nearer to pt
  63         return false;
  64     }
  65
  66     if ((diff < SMALL) || ((distSqr() > SMALL) && (diff/distSqr() < tol)))
  67     {
  68         // don't propagate small changes
  69         return false;
  70     }
  71     else
  72     {
  73         // update with new values
  74         distSqr() = dist2;
  75         origin() = w2.origin();
  76         data_ = w2.data();
  77
  78         return true;
  79     }
  80 }
  81
  82
  83 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  84
  85 // Null constructor
  86 template <class Type>
  87 inline wallPointData<Type>::wallPointData()
  88 :
  89     wallPoint(),
  90     data_()
  91 {}
  92
  93
  94 // Construct from components
  95 template <class Type>
  96 inline wallPointData<Type>::wallPointData
  97 (
  98     const point& origin,
  99     const Type& data,
 100     const scalar distSqr
 101 )
 102 :
 103     wallPoint(origin, distSqr),
 104     data_(data)
 105 {}
 106
 107
 108 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 109
 110 template <class Type>
 111 inline const Type& wallPointData<Type>::data() const
 112 {
 113     return data_;
 114 }
 115
 116
 117 template <class Type>
 118 inline Type& wallPointData<Type>::data()
 119 {
 120     return data_;
 121 }
 122
 123
 124 // Update this with w2 if w2 nearer to pt.
 125 template <class Type>
 126 inline bool wallPointData<Type>::updateCell
 127 (
 128     const polyMesh& mesh,
 129     const label thisCellI,
 130     const label,
 131     const wallPointData<Type>& neighbourWallInfo,
 132     const scalar tol
 133 )
 134 {
 135     const vectorField& cellCentres = mesh.primitiveMesh::cellCentres();
 136
 137     return update
 138     (
 139         cellCentres[thisCellI],
 140         neighbourWallInfo,
 141         tol
 142     );
 143 }
 144
 145
 146 // Update this with w2 if w2 nearer to pt.
 147 template <class Type>
 148 inline bool wallPointData<Type>::updateFace
 149 (
 150     const polyMesh& mesh,
 151     const label thisFaceI,
 152     const label,
 153     const wallPointData<Type>& neighbourWallInfo,
 154     const scalar tol
 155 )
 156 {
 157     const vectorField& faceCentres = mesh.faceCentres();
 158
 159     return update
 160     (
 161         faceCentres[thisFaceI],
 162         neighbourWallInfo,
 163         tol
 164     );
 165 }
 166
 167
 168 // Update this with w2 if w2 nearer to pt.
 169 template <class Type>
 170 inline bool wallPointData<Type>::updateFace
 171 (
 172     const polyMesh& mesh,
 173     const label thisFaceI,
 174     const wallPointData<Type>& neighbourWallInfo,
 175     const scalar tol
 176 )
 177 {
 178     const vectorField& faceCentres = mesh.faceCentres();
 179
 180     return update
 181     (
 182         faceCentres[thisFaceI],
 183         neighbourWallInfo,
 184         tol
 185     );
 186 }
 187
 188
 189 // ************************************************************************* //
 190
 191 } // End namespace Foam
 192
 193 // ************************************************************************* //