FreeFOAM The Cross-Platform CFD Toolkit
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) 1991-2010 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
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/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 
67  labelList faceMap;
68  surfacePatchList myPatches(calcPatches(faceMap));
69 
70  if (writeSorted)
71  {
72  label faceIndex = 0;
73 
74  forAll(myPatches, patchI)
75  {
76  // Print all faces belonging to this patch
77 
78  for
79  (
80  label patchFaceI = 0;
81  patchFaceI < myPatches[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 
113  forAll(myPatches, patchI)
114  {
115  for
116  (
117  label patchFaceI = 0;
118  patchFaceI < myPatches[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 // ************************ vim: set sw=4 sts=4 et: ************************ //