src/postProcessing/functionObjects/field/surfaceInterpolateFields/surfaceInterpolateFieldsTemplates.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "surfaceInterpolateFields.H"
  27 #include "volFields.H"
  28 #include "linear.H"
  29
  30 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  31
  32 template<class Type>
  33 void Foam::surfaceInterpolateFields::interpolateFields
  34 (
  35     PtrList<GeometricField<Type, fvsPatchField, surfaceMesh> >& sflds
  36 ) const
  37 {
  38     typedef GeometricField<Type, fvPatchField, volMesh> vfType;
  39     typedef GeometricField<Type, fvsPatchField, surfaceMesh> sfType;
  40
  41     // Convert field to map
  42     HashTable<word> fieldMap(2*fieldSet_.size());
  43     forAll(fieldSet_, i)
  44     {
  45         fieldMap.insert(fieldSet_[i].first(), fieldSet_[i].second());
  46     }
  47
  48
  49     HashTable<const vfType*> flds(obr_.lookupClass<vfType>());
  50
  51     forAllConstIter(typename HashTable<const vfType*>, flds, iter)
  52     {
  53         const vfType& fld = *iter();
  54
  55         if (fieldMap.found(fld.name()))
  56         {
  57             //const word sName = "interpolate(" + fld.name() + ')';
  58             const word& sName = fieldMap[fld.name()];
  59
  60             if (obr_.found(sName))
  61             {
  62                 Info<< "    a surfaceField " << sName << " already exists"
  63                     << endl;
  64             }
  65             else
  66             {
  67                 label sz = sflds.size();
  68                 sflds.setSize(sz+1);
  69                 sflds.set(sz, new sfType(sName, linearInterpolate(fld)));
  70
  71                 Info<< "    interpolated " << fld.name() << " to create "
  72                     << sflds[sz].name() << endl;
  73             }
  74         }
  75     }
  76 }
  77
  78
  79 // ************************************************************************* //