writeVTK.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2011-2016 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 "triSurface.H"
27 
28 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
29 
30 namespace Foam
31 {
32 
33 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
34 
35 void triSurface::writeVTK(const bool writeSorted, Ostream& os) const
36 {
37  // Write header
38  os << "# vtk DataFile Version 2.0" << nl
39  << "triSurface" << nl
40  << "ASCII" << nl
41  << "DATASET POLYDATA"
42  << nl;
43 
44  const pointField& ps = points();
45 
46  os << "POINTS " << ps.size() << " float" << nl;
47 
48  // Write vertex coords
49  forAll(ps, pointi)
50  {
51  if (pointi > 0 && (pointi % 10) == 0)
52  {
53  os << nl;
54  }
55  else
56  {
57  os << ' ';
58  }
59  os << ps[pointi].x() << ' '
60  << ps[pointi].y() << ' '
61  << ps[pointi].z();
62  }
63  os << nl;
64 
65  os << "POLYGONS " << size() << ' ' << 4*size() << nl;
66 
68  surfacePatchList patches(calcPatches(faceMap));
69 
70  if (writeSorted)
71  {
72  label faceIndex = 0;
73 
75  {
76  // Print all faces belonging to this patch
77 
78  for
79  (
80  label patchFacei = 0;
81  patchFacei < patches[patchi].size();
82  patchFacei++
83  )
84  {
85  if (faceIndex > 0 && (faceIndex % 10) == 0)
86  {
87  os << nl;
88  }
89  else
90  {
91  os << ' ';
92  }
93 
94  const label facei = faceMap[faceIndex++];
95 
96  os << "3 "
97  << operator[](facei)[0] << ' '
98  << operator[](facei)[1] << ' '
99  << operator[](facei)[2];
100  }
101  }
102  os << nl;
103 
104 
105  // Print region numbers
106 
107  os << "CELL_DATA " << size() << nl;
108  os << "FIELD attributes 1" << nl;
109  os << "region 1 " << size() << " float" << nl;
110 
111  faceIndex = 0;
112 
114  {
115  for
116  (
117  label patchFacei = 0;
118  patchFacei < patches[patchi].size();
119  patchFacei++
120  )
121  {
122  if (faceIndex > 0 && (faceIndex % 10) == 0)
123  {
124  os << nl;
125  }
126  else
127  {
128  os << ' ';
129  }
130 
131  const label facei = faceMap[faceIndex++];
132 
133  os << operator[](facei).region();
134  }
135  }
136  os << nl;
137  }
138  else
139  {
140  forAll(*this, facei)
141  {
142  if (facei > 0 && (facei % 10) == 0)
143  {
144  os << nl;
145  }
146  else
147  {
148  os << ' ';
149  }
150  os << "3 "
151  << operator[](facei)[0] << ' '
152  << operator[](facei)[1] << ' '
153  << operator[](facei)[2];
154  }
155  os << nl;
156 
157  os << "CELL_DATA " << size() << nl;
158  os << "FIELD attributes 1" << nl;
159  os << "region 1 " << size() << " float" << nl;
160 
161  forAll(*this, facei)
162  {
163  if (facei > 0 && (facei % 10) == 0)
164  {
165  os << nl;
166  }
167  else
168  {
169  os << ' ';
170  }
171  os << operator[](facei).region();
172  }
173  os << nl;
174  }
175 }
176 
177 
178 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
179 
180 } // End namespace Foam
181 
182 // ************************************************************************* //
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:428
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
T & operator[](const label)
Return element of UList.
Definition: UListI.H:167
void size(const label)
Override size to be inconsistent with allocated storage.
Definition: ListI.H:76
const Field< point > & points() const
Return reference to global points.
List< surfacePatch > surfacePatchList
Pair< int > faceMap(const label facePi, const face &faceP, const label faceNi, const face &faceN)
vectorField pointField
pointField is a vectorField.
Definition: pointFieldFwd.H:42
const geometricSurfacePatchList & patches() const
Definition: triSurface.H:314
List< label > labelList
A List of labels.
Definition: labelList.H:56
static const char nl
Definition: Ostream.H:262
label size() const
Return the number of elements in the UList.
label patchi
Namespace for OpenFOAM.