42 const point& pt = points[pointI];
44 os << pt.
x() <<
' ' << pt.y() <<
' ' << pt.z() <<
' ';
59 os << ct.
x() <<
' ' << ct.y() <<
' ' << ct.z() <<
' ';
67 const fileName& fieldName,
75 os <<
"# x y z " << fieldName <<
endl;
78 if (values.size() == points.size())
82 writeGeometry(points, elemI, os);
83 os << values[elemI] <<
nl;
90 writeGeometry(points, faces, elemI, os);
91 os << values[elemI] <<
nl;
103 const fileName& fieldName,
112 << fieldName <<
"_x "
113 << fieldName <<
"_y "
114 << fieldName <<
"_z "
118 if (values.size() == points.size())
122 writeGeometry(points, elemI, os);
124 const vector& v = values[elemI];
125 os << v[0] <<
' ' << v[1] <<
' ' << v[2] <<
nl;
132 writeGeometry(points, faces, elemI, os);
134 const vector& v = values[elemI];
135 os << v[0] <<
' ' << v[1] <<
' ' << v[2] <<
nl;
146 const fileName& fieldName,
155 os << fieldName <<
"_ii" <<
endl;
158 if (values.size() == points.size())
162 writeGeometry(points, elemI, os);
172 writeGeometry(points, faces, elemI, os);
185 const fileName& fieldName,
193 os <<
"# xx xy xz yy yz ";
194 for(
int i=0; i<6; i++)
196 os << fieldName <<
"_" << i <<
" ";
201 if (values.size() == points.size())
205 writeGeometry(points, elemI, os);
209 os << v[0] <<
' ' << v[1] <<
' ' << v[2] <<
' '
210 << v[3] <<
' ' << v[4] <<
' ' << v[5] <<
' '
218 writeGeometry(points, faces, elemI, os);
222 os << v[0] <<
' ' << v[1] <<
' ' << v[2] <<
' '
223 << v[3] <<
' ' << v[4] <<
' ' << v[5] <<
' '
234 const fileName& fieldName,
242 os <<
"# xx xy xz yx yy yz zx zy zz";
243 for (
int i=0; i<9; ++i)
245 os << fieldName <<
"_" << i <<
" ";
250 if (values.size() == points.size())
254 writeGeometry(points, elemI, os);
256 const tensor& v = values[elemI];
257 os << v[0] <<
' ' << v[1] <<
' ' << v[2] <<
' '
258 << v[3] <<
' ' << v[4] <<
' ' << v[5] <<
' '
259 << v[6] <<
' ' << v[7] <<
' ' << v[8] <<
nl;
266 writeGeometry(points, faces, elemI, os);
268 const tensor& v = values[elemI];
269 os << v[0] <<
' ' << v[1] <<
' ' << v[2] <<
' '
270 << v[3] <<
' ' << v[4] <<
' ' << v[5] <<
' '
271 << v[6] <<
' ' << v[7] <<
' ' << v[8] <<
nl;
305 if (!
isDir(outputDir))
312 outputDir/surfaceName +
".raw"
322 os <<
"# geometry NO_DATA " << faces.
size() <<
nl
323 <<
"# x y z" <<
endl;
328 writeGeometry(points, faces, elemI, os);
366 if (!
isDir(outputDir))
373 outputDir/fieldName +
'_' + surfaceName +
".raw"
378 Info<<
"Writing field " << fieldName <<
" to " << os.
name() <<
endl;
383 os <<
"# " << fieldName;
386 os <<
" POINT_DATA ";
393 os << values.
size() <<
nl;
395 writeData(fieldName, points, faces, values, os);