1 #include "rheolef/geo.h"
22 #include "rheolef/geo_domain.h"
23 #include "rheolef/rheostream.h"
24 #include "rheolef/iorheo.h"
29 template <
class T> odiststream&
visu_gnuplot (odiststream&,
const geo_basic<T,sequential>&);
30 template <
class T> odiststream&
visu_vtk_paraview (odiststream&,
const geo_basic<T,sequential>&);
31 template <
class T> odiststream&
geo_put_gmsh (odiststream&,
const geo_basic<T,sequential>&);
32 template <
class T> odiststream&
geo_put_bamg (odiststream&,
const geo_basic<T,sequential>&);
37 iorheo::flag_type format = iorheo::flags(ops.
os()) & iorheo::format_field;
42 return base::data().put (ops);
48 size_type n_vert = base::_node.dis_size ();
49 size_type n_elt = base::_gs.ownership_by_dimension[base::_gs._map_dimension].dis_size ();
53 ops << setprecision(numeric_limits<T>::digits10)
59 if (base::_version == 4) {
61 h.dimension = base::_dimension;
62 h.sys_coord = base::_sys_coord;
64 h.dis_size_by_variant [0] = base::_node.dis_size();
65 if (base::_have_connectivity) {
68 h.dis_size_by_variant [
variant] = base::_geo_element [
variant].dis_size();
74 h.dis_size_by_variant [
variant] = base::_geo_element [
variant].dis_size();
77 ops << endl <<
h << endl;
82 T rounding_prec = iorheo::getrounding_precision(ops.
os());
83 if (rounding_prec == 0) {
84 base::_node.put_values (ops,
_point_put<T>(base::_dimension));
92 if (base::_gs._map_dimension > 0) {
95 base::_geo_element [
variant].put_values (ops);
102 if (base::_gs._map_dimension > 1 && base::_have_connectivity) {
103 for (
size_type side_dim = base::_gs._map_dimension - 1; side_dim >= 1; side_dim--) {
106 base::_geo_element [
variant].put_values (ops);
114 iter = base::_domains.begin(), last = base::_domains.end();
115 iter != last; ++iter) {