src/dynamicMesh/dynamicFvMesh/dynamicTopoFvMesh/fieldMapping/topoCellMapperTemplates.C

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