src/dynamicMesh/dynamicTopoFvMesh/fieldMapping/topoSurfaceMapperTemplates.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 //- Map the internal field
  32 template <class Type>
  33 void topoSurfaceMapper::mapInternalField
  34 (
  35     const word& fieldName,
  36     Field<Type>& iF
  37 ) const
  38 {
  39     if (iF.size() != sizeBeforeMapping())
  40     {
  41         FatalErrorIn
  42         (
  43             "\n\n"
  44             "void topoSurfaceMapper::mapInternalField<Type>\n"
  45             "(\n"
  46             "    Field<Type>& iF\n"
  47             ") const\n"
  48         )  << "Incompatible size before mapping." << nl
  49            << " Field: " << fieldName << nl
  50            << " Field size: " << iF.size() << nl
  51            << " map size: " << sizeBeforeMapping() << nl
  52            << abort(FatalError);
  53     }
  54
  55     // Map the internal field
  56     iF.autoMap(*this);
  57
  58     // Flip the flux
  59     const labelList flipFaces = flipFaceFlux().toc();
  60
  61     forAll (flipFaces, i)
  62     {
  63         if (flipFaces[i] < iF.size())
  64         {
  65             iF[flipFaces[i]] *= -1.0;
  66         }
  67         else
  68         {
  69             FatalErrorIn
  70             (
  71                 "\n\n"
  72                 "void topoSurfaceMapper::mapInternalField<Type>\n"
  73                 "(\n"
  74                 "    Field<Type>& iF\n"
  75                 ") const\n"
  76             )  << "Cannot flip boundary face fluxes." << nl
  77                << " Field: " << fieldName << nl
  78                << " Field size: " << iF.size() << nl
  79                << " Face flip index: " << flipFaces[i] << nl
  80                << abort(FatalError);
  81         }
  82     }
  83 }
  84
  85
  86 } // End namespace Foam
  87
  88 // ************************************************************************* //