applications/utilities/postProcessing/dataConversion/foamToVTK/surfaceMeshWriterTemplates.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 "surfaceMeshWriter.H"
  27 #include "writeFuns.H"
  28
  29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  30
  31 template<class Type>
  32 Foam::tmp<Field<Type> > Foam::surfaceMeshWriter::getFaceField
  33 (
  34     const GeometricField<Type, fvsPatchField, surfaceMesh>& sfld
  35 ) const
  36 {
  37     const polyBoundaryMesh& patches = sfld.mesh().boundaryMesh();
  38
  39     tmp<Field<Type> > tfld(new Field<Type>(pp_.size()));
  40     Field<Type>& fld = tfld();
  41
  42     forAll(pp_.addressing(), i)
  43     {
  44         label faceI = pp_.addressing()[i];
  45
  46         label patchI = patches.whichPatch(faceI);
  47
  48         if (patchI == -1)
  49         {
  50             fld[i] = sfld[faceI];
  51         }
  52         else
  53         {
  54             label localFaceI = faceI - patches[patchI].start();
  55             fld[i] = sfld.boundaryField()[patchI][localFaceI];
  56         }
  57     }
  58
  59     return tfld;
  60 }
  61
  62
  63 template<class Type>
  64 void Foam::surfaceMeshWriter::write
  65 (
  66     const PtrList<GeometricField<Type, fvsPatchField, surfaceMesh> >& sflds
  67 )
  68 {
  69     forAll(sflds, fieldI)
  70     {
  71         const GeometricField<Type, fvsPatchField, surfaceMesh>& fld =
  72             sflds[fieldI];
  73
  74         os_ << fld.name() << ' ' << pTraits<Type>::nComponents << ' '
  75             << pp_.size() << " float" << std::endl;
  76
  77         DynamicList<floatScalar> fField(pTraits<Type>::nComponents*pp_.size());
  78         writeFuns::insert(getFaceField(fld)(), fField);
  79         writeFuns::write(os_, binary_, fField);
  80     }
  81 }
  82
  83
  84 // ************************************************************************* //