src/foam/interpolations/patchToPatchInterpolation/PatchToPatchInterpolationTemplate.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Description
  25     Interpolation class dealing with transfer of data between two
  26     primitivePatches
  27
  28 \*---------------------------------------------------------------------------*/
  29
  30 #include "PatchToPatchInterpolationTemplate.H"
  31 #include "demandDrivenData.H"
  32
  33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  34
  35 namespace Foam
  36 {
  37
  38 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  39
  40 template<class FromPatch, class ToPatch>
  41 const Foam::debug::tolerancesSwitch
  42 PatchToPatchInterpolation<FromPatch, ToPatch>::directHitTol_
  43 (
  44     "patchToPatchDirectHit",
  45     1e-5
  46 );
  47
  48 template<class FromPatch, class ToPatch>
  49 const Foam::debug::tolerancesSwitch
  50 PatchToPatchInterpolation<FromPatch, ToPatch>::projectionTol_
  51 (
  52     "patchToPatchProjectionTol",
  53     0.05
  54 );
  55
  56 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  57
  58 template<class FromPatch, class ToPatch>
  59 void PatchToPatchInterpolation<FromPatch, ToPatch>::clearOut()
  60 {
  61     deleteDemandDrivenData(pointAddressingPtr_);
  62     deleteDemandDrivenData(pointWeightsPtr_);
  63     deleteDemandDrivenData(pointDistancePtr_);
  64     deleteDemandDrivenData(faceAddressingPtr_);
  65     deleteDemandDrivenData(faceWeightsPtr_);
  66     deleteDemandDrivenData(faceDistancePtr_);
  67 }
  68
  69
  70 template<class FromPatch, class ToPatch>
  71 void PatchToPatchInterpolation<FromPatch, ToPatch>::setWeights
  72 (
  73     labelList* paPtr,
  74     FieldField<Field, scalar>* pwPtr,
  75     scalarField* pdPtr,
  76     labelList* faPtr,
  77     FieldField<Field, scalar>* fwPtr,
  78     scalarField* fdPtr
  79 )
  80 {
  81     clearOut();
  82
  83     pointAddressingPtr_ = paPtr;
  84     pointWeightsPtr_ = pwPtr;
  85     pointDistancePtr_ = pdPtr;
  86     faceAddressingPtr_ = faPtr;;
  87     faceWeightsPtr_ = fwPtr;
  88     faceDistancePtr_ = fdPtr;;
  89 }
  90
  91
  92 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  93
  94 // Construct from components
  95 template<class FromPatch, class ToPatch>
  96 PatchToPatchInterpolation<FromPatch, ToPatch>::PatchToPatchInterpolation
  97 (
  98     const FromPatch& fromPatch,
  99     const ToPatch& toPatch,
 100     intersection::algorithm alg,
 101     const intersection::direction dir
 102 )
 103 :
 104     fromPatch_(fromPatch),
 105     toPatch_(toPatch),
 106     alg_(alg),
 107     dir_(dir),
 108     pointAddressingPtr_(NULL),
 109     pointWeightsPtr_(NULL),
 110     pointDistancePtr_(NULL),
 111     faceAddressingPtr_(NULL),
 112     faceWeightsPtr_(NULL),
 113     faceDistancePtr_(NULL)
 114 {
 115 }
 116
 117
 118 // Construct as copy
 119 template<class FromPatch, class ToPatch>
 120 PatchToPatchInterpolation<FromPatch, ToPatch>::PatchToPatchInterpolation
 121 (
 122     const PatchToPatchInterpolation<FromPatch, ToPatch>& ppi
 123 )
 124 :
 125     PatchToPatchInterpolationName(),
 126     fromPatch_(ppi.fromPatch_),
 127     toPatch_(ppi.toPatch_),
 128     alg_(ppi.alg_),
 129     dir_(ppi.dir_),
 130     pointAddressingPtr_(NULL),
 131     pointWeightsPtr_(NULL),
 132     pointDistancePtr_(NULL),
 133     faceAddressingPtr_(NULL),
 134     faceWeightsPtr_(NULL),
 135     faceDistancePtr_(NULL)
 136 {}
 137
 138
 139 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
 140
 141 template<class FromPatch, class ToPatch>
 142 PatchToPatchInterpolation<FromPatch, ToPatch>::~PatchToPatchInterpolation()
 143 {
 144     clearOut();
 145 }
 146
 147
 148 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 149
 150 template<class FromPatch, class ToPatch>
 151 const labelList&
 152 PatchToPatchInterpolation<FromPatch, ToPatch>::pointAddr() const
 153 {
 154     if (!pointAddressingPtr_)
 155     {
 156         calcPointAddressing();
 157     }
 158
 159     return *pointAddressingPtr_;
 160 }
 161
 162
 163 template<class FromPatch, class ToPatch>
 164 const FieldField<Field, scalar>&
 165 PatchToPatchInterpolation<FromPatch, ToPatch>::pointWeights() const
 166 {
 167     if (!pointWeightsPtr_)
 168     {
 169         calcPointAddressing();
 170     }
 171
 172     return *pointWeightsPtr_;
 173 }
 174
 175
 176 template<class FromPatch, class ToPatch>
 177 const labelList&
 178 PatchToPatchInterpolation<FromPatch, ToPatch>::faceAddr() const
 179 {
 180     if (!faceAddressingPtr_)
 181     {
 182         calcFaceAddressing();
 183     }
 184
 185     return *faceAddressingPtr_;
 186 }
 187
 188
 189 template<class FromPatch, class ToPatch>
 190 const FieldField<Field, scalar>&
 191 PatchToPatchInterpolation<FromPatch, ToPatch>::faceWeights() const
 192 {
 193     if (!faceWeightsPtr_)
 194     {
 195         calcFaceAddressing();
 196     }
 197
 198     return *faceWeightsPtr_;
 199 }
 200
 201
 202 template<class FromPatch, class ToPatch>
 203 const scalarField&
 204 PatchToPatchInterpolation<FromPatch, ToPatch>
 205 ::pointDistanceToIntersection() const
 206 {
 207     if (!pointDistancePtr_)
 208     {
 209         calcPointAddressing();
 210     }
 211
 212     return *pointDistancePtr_;
 213 }
 214
 215
 216 template<class FromPatch, class ToPatch>
 217 const scalarField&
 218 PatchToPatchInterpolation<FromPatch, ToPatch>
 219 ::faceDistanceToIntersection() const
 220 {
 221     if (!faceDistancePtr_)
 222     {
 223         calcFaceAddressing();
 224     }
 225
 226     return *faceDistancePtr_;
 227 }
 228
 229
 230 template<class FromPatch, class ToPatch>
 231 bool PatchToPatchInterpolation<FromPatch, ToPatch>::movePoints()
 232 {
 233     clearOut();
 234
 235     return true;
 236 }
 237
 238
 239 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 240
 241 } // End namespace Foam
 242
 243 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 244
 245 #   include "CalcPatchToPatchWeights.C"
 246 #   include "PatchToPatchInterpolate.C"
 247
 248 // ************************************************************************* //