src/surfMesh/surfaceFormats/vtk/VTKsurfaceFormatCore.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   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 #include "VTKsurfaceFormatCore.H"
  28 #include "clock.H"
  29
  30 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  31
  32 void Foam::fileFormats::VTKsurfaceFormatCore::writeHeader
  33 (
  34     Ostream& os,
  35     const pointField& pointLst
  36 )
  37 {
  38     // Write header
  39     os  << "# vtk DataFile Version 2.0" << nl
  40         << "surface written " << clock::dateTime().c_str() << nl
  41         << "ASCII" << nl
  42         << nl
  43         << "DATASET POLYDATA" << nl;
  44
  45     // Write vertex coords
  46     os  << "POINTS " << pointLst.size() << " float" << nl;
  47     forAll(pointLst, ptI)
  48     {
  49         const point& pt = pointLst[ptI];
  50
  51         os  << pt.x() << ' ' << pt.y() << ' ' << pt.z() << nl;
  52     }
  53 }
  54
  55
  56 void Foam::fileFormats::VTKsurfaceFormatCore::writeTail
  57 (
  58     Ostream& os,
  59     const UList<surfZone>& zoneLst
  60 )
  61 {
  62     label nFaces = 0;
  63     forAll(zoneLst, zoneI)
  64     {
  65         nFaces += zoneLst[zoneI].size();
  66     }
  67
  68     // Print zone numbers
  69     os  << nl
  70         << "CELL_DATA " << nFaces << nl
  71         << "FIELD attributes 1" << nl
  72         << "zone 1 " << nFaces << " float" << nl;
  73
  74
  75     forAll(zoneLst, zoneI)
  76     {
  77         forAll(zoneLst[zoneI], localFaceI)
  78         {
  79             if (localFaceI)
  80             {
  81                 if (localFaceI % 20)
  82                 {
  83                     os << ' ';
  84                 }
  85                 else
  86                 {
  87                     os << nl;
  88                 }
  89             }
  90             os  << zoneI + 1;
  91         }
  92         os  << nl;
  93     }
  94 }
  95
  96
  97 void Foam::fileFormats::VTKsurfaceFormatCore::writeTail
  98 (
  99     Ostream& os,
 100     const UList<label>& zoneIds
 101 )
 102 {
 103     // Print zone numbers
 104     os  << nl
 105         << "CELL_DATA " << zoneIds.size() << nl
 106         << "FIELD attributes 1" << nl
 107         << "zone 1 " << zoneIds.size() << " float" << nl;
 108
 109     forAll(zoneIds, faceI)
 110     {
 111         if (faceI)
 112         {
 113             if (faceI % 20)
 114             {
 115                 os << ' ';
 116             }
 117             else
 118             {
 119                 os << nl;
 120             }
 121         }
 122         os  << zoneIds[faceI] + 1;
 123     }
 124     os  << nl;
 125 }
 126
 127
 128 // ************************************************************************* //