src/dynamicMesh/dynamicTopoFvMesh/fieldMapping/topoCellMapperTemplates.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 \*---------------------------------------------------------------------------*/
  25
  26 namespace Foam
  27 {
  28
  29 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  30
  31 //- Conservatively map the internal field
  32 template <class Type, class gradType>
  33 void topoCellMapper::mapInternalField
  34 (
  35     const word& fieldName,
  36     const Field<gradType>& gF,
  37     Field<Type>& iF
  38 ) const
  39 {
  40     if (iF.size() != sizeBeforeMapping() || gF.size() != sizeBeforeMapping())
  41     {
  42         FatalErrorIn
  43         (
  44             "\n\n"
  45             "void topoCellMapper::mapInternalField<Type>\n"
  46             "(\n"
  47             "    const word& fieldName,\n"
  48             "    const Field<gradType>& gF,\n"
  49             "    Field<Type>& iF\n"
  50             ") const\n"
  51         )  << "Incompatible size before mapping." << nl
  52            << " Field: " << fieldName << nl
  53            << " Field size: " << iF.size() << nl
  54            << " Gradient Field size: " << gF.size() << nl
  55            << " map size: " << sizeBeforeMapping() << nl
  56            << abort(FatalError);
  57     }
  58
  59     // If we have direct addressing, map and bail out
  60     if (direct())
  61     {
  62         iF.autoMap(*this);
  63         return;
  64     }
  65
  66     // Fetch addressing
  67     const labelListList& cAddressing = addressing();
  68     const List<scalarField>& wC = intersectionWeights();
  69     const List<vectorField>& xC = intersectionCentres();
  70
  71     // Fetch geometry
  72     const vectorField& centres = tMapper_.internalCentres();
  73
  74     // Copy the original field
  75     Field<Type> fieldCpy(iF);
  76
  77     // Resize to current dimensions
  78     iF.setSize(size());
  79
  80     // Map the internal field
  81     forAll(iF, cellI)
  82     {
  83         const labelList& addr = cAddressing[cellI];
  84
  85         iF[cellI] = pTraits<Type>::zero;
  86
  87         forAll(addr, cellJ)
  88         {
  89             // Accumulate volume-weighted Taylor-series interpolate
  90             iF[cellI] +=
  91             (
  92                 wC[cellI][cellJ] *
  93                 (
  94                     fieldCpy[addr[cellJ]]
  95                   + (
  96                         gF[addr[cellJ]] &
  97                         (
  98                             xC[cellI][cellJ]
  99                           - centres[addr[cellJ]]
 100                         )
 101                     )
 102                 )
 103             );
 104         }
 105     }
 106 }
 107
 108
 109 } // End namespace Foam
 110
 111 // ************************************************************************* //