43 <<
"# vtk DataFile Version 2.0" <<
nl
44 <<
"sampleSurface" <<
nl
46 <<
"DATASET POLYDATA" <<
nl;
49 os <<
"POINTS " << points.size() <<
" float" <<
nl;
52 const point& pt = points[pointI];
53 os << float(pt.x()) <<
' '
54 <<
float(pt.y()) <<
' '
55 <<
float(pt.z()) << nl;
64 nNodes += faces[faceI].size();
67 os <<
"POLYGONS " << faces.size() <<
' '
68 << faces.size() + nNodes <<
nl;
72 const face&
f = faces[faceI];
92 const Field<Foam::scalar>& values
95 os <<
"1 " << values.size() <<
" float" <<
nl;
111 const scalar& v = values[elemI];
122 const Field<Foam::vector>& values
125 os <<
"3 " << values.size() <<
" float" <<
nl;
129 const vector& v = values[elemI];
130 os << float(v[0]) <<
' ' << float(v[1]) <<
' ' << float(v[2])
141 const Field<sphericalTensor>& values
144 os <<
"1 " << values.size() <<
" float" <<
nl;
149 os << float(v[0]) <<
nl;
159 const Field<symmTensor>& values
162 os <<
"6 " << values.size() <<
" float" <<
nl;
167 os << float(v[0]) <<
' ' << float(v[1]) <<
' ' << float(v[2])
168 << float(v[3]) <<
' ' << float(v[4]) <<
' ' << float(v[5])
180 const Field<tensor>& values
183 os <<
"9 " << values.size() <<
" float" <<
nl;
187 const tensor& v = values[elemI];
188 os << float(v[0]) <<
' ' << float(v[1]) <<
' ' << float(v[2])
189 << float(v[3]) <<
' ' << float(v[4]) <<
' ' << float(v[5])
190 << float(v[6]) <<
' ' << float(v[7]) <<
' ' << float(v[8])
203 const Field<Type>& values
206 os <<
"1 " << values.size() <<
" float" <<
nl;
210 os << float(0) <<
nl;
244 if (!
isDir(outputDir))
249 fileName fName(outputDir/surfaceName +
".vtk");
253 Info<<
"Writing geometry to " << fName <<
endl;
257 writeGeometry(os, points, faces);
273 if (!
isDir(outputDir))
280 outputDir/fieldName +
'_' + surfaceName +
".vtk"
285 Info<<
"Writing field " << fieldName <<
" to " << os.
name() <<
endl;
288 writeGeometry(os, points, faces);
300 os << values.
size() << nl
301 <<
"FIELD attributes 1" << nl
302 << fieldName.c_str() <<
" ";
305 writeData(os, values);