libMesh
Public Member Functions | Public Attributes | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes | Static Private Attributes | List of all members
libMesh::PostscriptIO Class Reference

This class implements writing 2D meshes in Postscript. More...

#include <postscript_io.h>

Inheritance diagram for libMesh::PostscriptIO:
[legend]

Public Member Functions

 PostscriptIO (const MeshBase &mesh)
 Constructor. More...
 
virtual ~PostscriptIO ()
 Destructor. More...
 
virtual void write (const std::string &) override
 This method implements writing a mesh to a specified file. More...
 
void plot_quadratic_elem (const Elem *elem)
 Draws an element with Bezier curves. More...
 
void plot_linear_elem (const Elem *elem)
 Draws an element with straight lines. More...
 
virtual void write_equation_systems (const std::string &, const EquationSystems &, const std::set< std::string > *system_names=nullptr)
 This method implements writing a mesh with data to a specified file where the data is taken from the EquationSystems object. More...
 
virtual void write_discontinuous_equation_systems (const std::string &, const EquationSystems &, const std::set< std::string > *system_names=nullptr)
 This method implements writing a mesh with discontinuous data to a specified file where the data is taken from the EquationSystems object. More...
 
virtual void write_nodal_data (const std::string &, const std::vector< Number > &, const std::vector< std::string > &)
 This method implements writing a mesh with nodal data to a specified file where the nodal data and variable names are provided. More...
 
virtual void write_nodal_data (const std::string &, const NumericVector< Number > &, const std::vector< std::string > &)
 This method may be overridden by "parallel" output formats for writing nodal data. More...
 
virtual void write_nodal_data (const std::string &, const EquationSystems &, const std::set< std::string > *)
 This method should be overridden by "parallel" output formats for writing nodal data. More...
 
virtual void write_nodal_data_discontinuous (const std::string &, const std::vector< Number > &, const std::vector< std::string > &)
 This method implements writing a mesh with discontinuous data to a specified file where the nodal data and variables names are provided. More...
 
unsigned intascii_precision ()
 Return/set the precision to use when writing ASCII files. More...
 

Public Attributes

Real shade_value
 Controls greyscale shading of cells. More...
 
Real line_width
 Control the thickness of the lines used. More...
 

Protected Member Functions

const MeshBasemesh () const
 
virtual bool get_add_sides ()
 

Protected Attributes

const bool _is_parallel_format
 Flag specifying whether this format is parallel-capable. More...
 
const bool _serial_only_needed_on_proc_0
 Flag specifying whether this format can be written by only serializing the mesh to processor zero. More...
 

Private Member Functions

void _compute_edge_bezier_coeffs (const Elem *elem)
 Given a quadratic edge Elem which lies in the x-y plane, computes the Bezier coefficients. More...
 

Private Attributes

std::vector< Point_bezier_coeffs
 Vector containing 3 points corresponding to Bezier coefficients, as computed by _compute_edge_bezier_coeffs. More...
 
Point _offset
 Amount to add to every (x,y) point to place it in Postscript coordinates. More...
 
Real _scale
 Amount by which to stretch each point to place it in Postscript coordinates. More...
 
Point _current_point
 A point object used for temporary calculations. More...
 
std::ostringstream _cell_string
 Drawing style-independent data for a single cell. More...
 
std::ofstream _out
 Output file stream which will be opened when the file name is known. More...
 

Static Private Attributes

static const float _bezier_transform [3][3]
 Coefficients of the transformation from physical-space edge coordinates to Bezier basis coefficients. More...
 

Detailed Description

This class implements writing 2D meshes in Postscript.

It borrows several ideas from, and is a more simple-minded version of, the DataOutBase::write_eps() function from Deal II. Only output is supported here, and only the Mesh (none of the data) is written. The main use I imagined for this class is creating nice Mesh images for publications, since I didn't find/don't know of a free visualization program which would do this.

Author
John W. Peterson
Date
2008

Definition at line 53 of file postscript_io.h.

Constructor & Destructor Documentation

◆ PostscriptIO()

libMesh::PostscriptIO::PostscriptIO ( const MeshBase mesh)
explicit

Constructor.

Definition at line 43 of file postscript_io.C.

References _bezier_coeffs.

43  :
44  MeshOutput<MeshBase> (mesh_in),
45  shade_value(0.0),
46  line_width(0.5),
47  //_M(3,3),
48  _offset(0., 0.),
49  _scale(1.0),
50  _current_point(0., 0.)
51 {
52  // This code is still undergoing some development.
53  libmesh_experimental();
54 
55  // Entries of transformation matrix from physical to Bezier coords.
56  // _M(0,0) = -1./6.; _M(0,1) = 1./6.; _M(0,2) = 1.;
57  // _M(1,0) = -1./6.; _M(1,1) = 0.5 ; _M(1,2) = 1./6.;
58  // _M(2,0) = 0. ; _M(2,1) = 1. ; _M(2,2) = 0.;
59 
60  // Make sure there is enough room to store Bezier coefficients.
61  _bezier_coeffs.resize(3);
62 }
Point _offset
Amount to add to every (x,y) point to place it in Postscript coordinates.
Point _current_point
A point object used for temporary calculations.
template class LIBMESH_EXPORT MeshOutput< MeshBase >
Definition: mesh_output.C:180
Real shade_value
Controls greyscale shading of cells.
Definition: postscript_io.h:80
Real line_width
Control the thickness of the lines used.
Definition: postscript_io.h:88
Real _scale
Amount by which to stretch each point to place it in Postscript coordinates.
std::vector< Point > _bezier_coeffs
Vector containing 3 points corresponding to Bezier coefficients, as computed by _compute_edge_bezier_...

◆ ~PostscriptIO()

libMesh::PostscriptIO::~PostscriptIO ( )
virtualdefault

Destructor.

Member Function Documentation

◆ _compute_edge_bezier_coeffs()

void libMesh::PostscriptIO::_compute_edge_bezier_coeffs ( const Elem elem)
private

Given a quadratic edge Elem which lies in the x-y plane, computes the Bezier coefficients.

These may be passed to the Postscript routine "curveto".

Definition at line 263 of file postscript_io.C.

References _bezier_coeffs, _bezier_transform, _offset, _scale, libMesh::EDGE3, libMesh::Elem::point(), and libMesh::Elem::type().

Referenced by plot_quadratic_elem().

264 {
265  // I only know how to do this for an Edge3!
266  libmesh_assert_equal_to (elem->type(), EDGE3);
267 
268  // Get x-coordinates into an array, transform them,
269  // and repeat for y.
270  float phys_coords[3] = {0., 0., 0.};
271  float bez_coords[3] = {0., 0., 0.};
272 
273  for (unsigned int i=0; i<2; ++i)
274  {
275  // Initialize vectors. Physical coordinates are initialized
276  // by their postscript-scaled values.
277  for (unsigned int j=0; j<3; ++j)
278  {
279  phys_coords[j] = static_cast<float>
280  ((elem->point(j)(i) - _offset(i)) * _scale);
281  bez_coords[j] = 0.; // zero out result vector
282  }
283 
284  // Multiply matrix times vector
285  for (unsigned int j=0; j<3; ++j)
286  for (unsigned int k=0; k<3; ++k)
287  bez_coords[j] += _bezier_transform[j][k]*phys_coords[k];
288 
289  // Store result in _bezier_coeffs
290  for (unsigned int j=0; j<3; ++j)
291  _bezier_coeffs[j](i) = phys_coords[j];
292  }
293 }
Point _offset
Amount to add to every (x,y) point to place it in Postscript coordinates.
static const float _bezier_transform[3][3]
Coefficients of the transformation from physical-space edge coordinates to Bezier basis coefficients...
Real _scale
Amount by which to stretch each point to place it in Postscript coordinates.
std::vector< Point > _bezier_coeffs
Vector containing 3 points corresponding to Bezier coefficients, as computed by _compute_edge_bezier_...

◆ ascii_precision()

unsigned int & libMesh::MeshOutput< MeshBase >::ascii_precision ( )
inlineinherited

Return/set the precision to use when writing ASCII files.

By default we use numeric_limits<Real>::max_digits10, which should be enough to write out to ASCII and get the exact same Real back when reading in.

Definition at line 269 of file mesh_output.h.

Referenced by libMesh::UNVIO::nodes_out(), libMesh::FroIO::write(), libMesh::MEDITIO::write_ascii(), libMesh::TecplotIO::write_ascii(), libMesh::GMVIO::write_ascii_new_impl(), and libMesh::GMVIO::write_ascii_old_impl().

270 {
271  return _ascii_precision;
272 }
unsigned int _ascii_precision
Precision to use when writing ASCII files.
Definition: mesh_output.h:207

◆ get_add_sides()

virtual bool libMesh::MeshOutput< MeshBase >::get_add_sides ( )
inlineprotectedvirtualinherited
Returns
Whether or not added sides are expected to be output, to plot SIDE_DISCONTINUOUS data. Subclasses should override this if they are capable of plotting such data.

Reimplemented in libMesh::ExodusII_IO.

Definition at line 176 of file mesh_output.h.

176 { return false; }

◆ mesh()

const MeshBase & libMesh::MeshOutput< MeshBase >::mesh ( ) const
inlineprotectedinherited

◆ plot_linear_elem()

void libMesh::PostscriptIO::plot_linear_elem ( const Elem elem)

Draws an element with straight lines.

Definition at line 189 of file postscript_io.C.

References _cell_string, _current_point, _offset, _out, _scale, libMesh::Elem::n_vertices(), libMesh::Elem::point(), and shade_value.

Referenced by write().

190 {
191  // Clear the string contents. Yes, this really is how you do that...
192  _cell_string.str("");
193 
194  // The general strategy is:
195  // 1.) Use m := {moveto} to go to vertex 0.
196  // 2.) Use l := {lineto} commands to draw lines to vertex 1, 2, ... N-1.
197  // 3a.) Use lx := {lineto closepath stroke} command at vertex N to draw the last line.
198  // 3b.) lf := {lineto closepath fill} command to shade the cell just drawn
199  // All of our 2D elements' vertices are numbered in counterclockwise order,
200  // so we can just draw them in the same order.
201 
202  // 1.)
203  _current_point = (elem->point(0) - _offset) * _scale;
204  _cell_string << _current_point(0) << " " << _current_point(1) << " "; // write x y
205  _cell_string << "m ";
206 
207  // 2.)
208  const unsigned int nv=elem->n_vertices();
209  for (unsigned int v=1; v<nv-1; ++v)
210  {
211  _current_point = (elem->point(v) - _offset) * _scale;
212  _cell_string << _current_point(0) << " " << _current_point(1) << " "; // write x y
213  _cell_string << "l ";
214  }
215 
216  // 3.)
217  _current_point = (elem->point(nv-1) - _offset) * _scale;
218  _cell_string << _current_point(0) << " " << _current_point(1) << " "; // write x y
219 
220  // We draw the shaded (interior) parts first, if applicable.
221  if (shade_value > 0.0)
222  _out << shade_value << " sg " << _cell_string.str() << "lf\n";
223 
224  // Draw the black lines (I guess we will always do this)
225  _out << "0 sg " << _cell_string.str() << "lx\n";
226 }
std::ofstream _out
Output file stream which will be opened when the file name is known.
Point _offset
Amount to add to every (x,y) point to place it in Postscript coordinates.
Point _current_point
A point object used for temporary calculations.
Real shade_value
Controls greyscale shading of cells.
Definition: postscript_io.h:80
Real _scale
Amount by which to stretch each point to place it in Postscript coordinates.
std::ostringstream _cell_string
Drawing style-independent data for a single cell.

◆ plot_quadratic_elem()

void libMesh::PostscriptIO::plot_quadratic_elem ( const Elem elem)

Draws an element with Bezier curves.

Definition at line 231 of file postscript_io.C.

References _bezier_coeffs, _compute_edge_bezier_coeffs(), _current_point, _offset, _out, _scale, libMesh::Elem::build_side_ptr(), libMesh::EDGE3, libMesh::index_range(), and libMesh::Elem::side_index_range().

232 {
233  for (auto ns : elem->side_index_range())
234  {
235  // Build the quadratic side
236  std::unique_ptr<const Elem> side = elem->build_side_ptr(ns);
237 
238  // Be sure it's quadratic (Edge2). Eventually we could
239  // handle cubic elements as well...
240  libmesh_assert_equal_to ( side->type(), EDGE3 );
241 
242  _out << "0 sg ";
243 
244  // Move to the first point on this side.
245  _current_point = (side->point(0) - _offset) * _scale;
246  _out << _current_point(0) << " " << _current_point(1) << " "; // write x y
247  _out << "m ";
248 
249  // Compute _bezier_coeffs for this edge. This fills up
250  // the _bezier_coeffs vector.
251  this->_compute_edge_bezier_coeffs(side.get());
252 
253  // Print curveto path to file
254  for (auto i : index_range(_bezier_coeffs))
255  _out << _bezier_coeffs[i](0) << " " << _bezier_coeffs[i](1) << " ";
256  _out << " cs\n";
257  }
258 }
std::ofstream _out
Output file stream which will be opened when the file name is known.
Point _offset
Amount to add to every (x,y) point to place it in Postscript coordinates.
void _compute_edge_bezier_coeffs(const Elem *elem)
Given a quadratic edge Elem which lies in the x-y plane, computes the Bezier coefficients.
Point _current_point
A point object used for temporary calculations.
Real _scale
Amount by which to stretch each point to place it in Postscript coordinates.
std::vector< Point > _bezier_coeffs
Vector containing 3 points corresponding to Bezier coefficients, as computed by _compute_edge_bezier_...
auto index_range(const T &sizable)
Helper function that returns an IntRange<std::size_t> representing all the indices of the passed-in v...
Definition: int_range.h:111

◆ write()

void libMesh::PostscriptIO::write ( const std::string &  fname)
overridevirtual

This method implements writing a mesh to a specified file.

Implements libMesh::MeshOutput< MeshBase >.

Definition at line 70 of file postscript_io.C.

References libMesh::MeshOutput< MeshBase >::_is_parallel_format, _offset, _out, _scale, libMesh::MeshTools::create_bounding_box(), line_width, libMesh::MeshOutput< MeshBase >::mesh(), libMesh::MeshOutput< MT >::mesh(), libMesh::MeshBase::mesh_dimension(), plot_linear_elem(), and libMesh::Real.

71 {
72  // We may need to gather a DistributedMesh to output it, making that
73  // const qualifier in our constructor a dirty lie
74  MeshSerializer serialize(const_cast<MeshBase &>(this->mesh()), !_is_parallel_format);
75 
76  if (this->mesh().processor_id() == 0)
77  {
78  // Get a constant reference to the mesh.
79  const MeshBase & the_mesh = MeshOutput<MeshBase>::mesh();
80 
81  // Only works in 2D
82  libmesh_assert_equal_to (the_mesh.mesh_dimension(), 2);
83 
84  // Create output file stream.
85  // _out is now a private member of the class.
86  _out.open(fname.c_str());
87 
88  // Make sure it opened correctly
89  if (!_out.good())
90  libmesh_file_error(fname.c_str());
91 
92  // The mesh bounding box gives us info about what the
93  // Postscript bounding box should be.
94  BoundingBox bbox = MeshTools::create_bounding_box(the_mesh);
95 
96  // Add a little extra padding to the "true" bounding box so
97  // that we can still see the boundary
98  const Real percent_padding = 0.01;
99  const Real dx=bbox.second(0)-bbox.first(0); libmesh_assert_greater (dx, 0.0);
100  const Real dy=bbox.second(1)-bbox.first(1); libmesh_assert_greater (dy, 0.0);
101 
102  const Real x_min = bbox.first(0) - percent_padding*dx;
103  const Real y_min = bbox.first(1) - percent_padding*dy;
104  const Real x_max = bbox.second(0) + percent_padding*dx;
105  const Real y_max = bbox.second(1) + percent_padding*dy;
106 
107  // Width of the output as given in postscript units.
108  // This usually is given by the strange unit 1/72 inch.
109  // A width of 300 represents a size of roughly 10 cm.
110  const Real width = 300;
111  _scale = width / (x_max-x_min);
112  _offset(0) = x_min;
113  _offset(1) = y_min;
114 
115  // Header writing stuff stolen from Deal.II
116  std::time_t time1= std::time (0);
117  std::tm * time = std::localtime(&time1);
118  _out << "%!PS-Adobe-2.0 EPSF-1.2" << '\n'
119  //<< "%!PS-Adobe-1.0" << '\n' // Lars' PS version
120  << "%%Filename: " << fname << '\n'
121  << "%%Title: LibMesh Output" << '\n'
122  << "%%Creator: LibMesh: A C++ finite element library" << '\n'
123  << "%%Creation Date: "
124  << time->tm_year+1900 << "/"
125  << time->tm_mon+1 << "/"
126  << time->tm_mday << " - "
127  << time->tm_hour << ":"
128  << std::setw(2) << time->tm_min << ":"
129  << std::setw(2) << time->tm_sec << '\n'
130  << "%%BoundingBox: "
131  // lower left corner
132  << "0 0 "
133  // upper right corner
134  << static_cast<unsigned int>( rint(double((x_max-x_min) * _scale )))
135  << ' '
136  << static_cast<unsigned int>( rint(double((y_max-y_min) * _scale )))
137  << '\n';
138 
139  // define some abbreviations to keep
140  // the output small:
141  // m=move turtle to
142  // l=define a line
143  // s=set rgb color
144  // sg=set gray value
145  // lx=close the line and plot the line
146  // lf=close the line and fill the interior
147  _out << "/m {moveto} bind def" << '\n'
148  << "/l {lineto} bind def" << '\n'
149  << "/s {setrgbcolor} bind def" << '\n'
150  << "/sg {setgray} bind def" << '\n'
151  << "/cs {curveto stroke} bind def" << '\n'
152  << "/lx {lineto closepath stroke} bind def" << '\n'
153  << "/lf {lineto closepath fill} bind def" << '\n';
154 
155  _out << "%%EndProlog" << '\n';
156  // << '\n';
157 
158  // Set line width in the postscript file.
159  _out << line_width << " setlinewidth" << '\n';
160 
161  // Set line cap and join options
162  _out << "1 setlinecap" << '\n';
163  _out << "1 setlinejoin" << '\n';
164 
165  // allow only five digits for output (instead of the default
166  // six); this should suffice even for fine grids, but reduces
167  // the file size significantly
168  _out << std::setprecision (5);
169 
170  // Loop over the active elements, draw lines for the edges. We
171  // draw even quadratic elements with straight sides, i.e. a straight
172  // line sits between each pair of vertices. Also we draw every edge
173  // for an element regardless of the fact that it may overlap with
174  // another. This would probably be a useful optimization...
175  for (const auto & elem : the_mesh.active_element_ptr_range())
176  this->plot_linear_elem(elem);
177 
178  // Issue the showpage command, and we're done.
179  _out << "showpage" << std::endl;
180 
181  } // end if (this->mesh().processor_id() == 0)
182 }
const MeshBase & mesh() const
Definition: mesh_output.h:259
void plot_linear_elem(const Elem *elem)
Draws an element with straight lines.
std::ofstream _out
Output file stream which will be opened when the file name is known.
Point _offset
Amount to add to every (x,y) point to place it in Postscript coordinates.
libMesh::BoundingBox create_bounding_box(const MeshBase &mesh)
Definition: mesh_tools.C:558
const bool _is_parallel_format
Flag specifying whether this format is parallel-capable.
Definition: mesh_output.h:184
Real line_width
Control the thickness of the lines used.
Definition: postscript_io.h:88
Real _scale
Amount by which to stretch each point to place it in Postscript coordinates.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ write_discontinuous_equation_systems()

void libMesh::MeshOutput< MeshBase >::write_discontinuous_equation_systems ( const std::string &  fname,
const EquationSystems es,
const std::set< std::string > *  system_names = nullptr 
)
virtualinherited

This method implements writing a mesh with discontinuous data to a specified file where the data is taken from the EquationSystems object.

Definition at line 89 of file mesh_output.C.

References libMesh::EquationSystems::build_discontinuous_solution_vector(), libMesh::EquationSystems::build_variable_names(), libMesh::EquationSystems::get_mesh(), libMesh::libmesh_assert(), and libMesh::out.

Referenced by libMesh::ExodusII_IO::write_timestep_discontinuous().

92 {
93  LOG_SCOPE("write_discontinuous_equation_systems()", "MeshOutput");
94 
95  // We may need to gather and/or renumber a DistributedMesh to output
96  // it, making that const qualifier in our constructor a dirty lie
97  MT & my_mesh = const_cast<MT &>(*_obj);
98 
99  // If we're asked to write data that's associated with a different
100  // mesh, output files full of garbage are the result.
101  libmesh_assert_equal_to(&es.get_mesh(), _obj);
102 
103  // A non-renumbered mesh may not have a contiguous numbering, and
104  // that needs to be fixed before we can build a solution vector.
105  if (my_mesh.max_elem_id() != my_mesh.n_elem() ||
106  my_mesh.max_node_id() != my_mesh.n_nodes())
107  {
108  // If we were allowed to renumber then we should have already
109  // been properly renumbered...
110  libmesh_assert(!my_mesh.allow_renumbering());
111 
112  libmesh_do_once(libMesh::out <<
113  "Warning: This MeshOutput subclass only supports meshes which are contiguously renumbered!"
114  << std::endl;);
115 
116  my_mesh.allow_renumbering(true);
117 
118  my_mesh.renumber_nodes_and_elements();
119 
120  // Not sure what good going back to false will do here, the
121  // renumbering horses have already left the barn...
122  my_mesh.allow_renumbering(false);
123  }
124 
125  MeshSerializer serialize(const_cast<MT &>(*_obj), !_is_parallel_format, _serial_only_needed_on_proc_0);
126 
127  // Build the list of variable names that will be written.
128  std::vector<std::string> names;
129  es.build_variable_names (names, nullptr, system_names);
130 
131  if (!_is_parallel_format)
132  {
133  // Build the nodal solution values & get the variable
134  // names from the EquationSystems object
135  std::vector<Number> soln;
136  es.build_discontinuous_solution_vector (soln, system_names,
137  nullptr, false, /* defaults */
138  this->get_add_sides());
139 
140  this->write_nodal_data_discontinuous (fname, soln, names);
141  }
142  else // _is_parallel_format
143  {
144  libmesh_not_implemented();
145  }
146 }
virtual void write_nodal_data_discontinuous(const std::string &, const std::vector< Number > &, const std::vector< std::string > &)
This method implements writing a mesh with discontinuous data to a specified file where the nodal dat...
Definition: mesh_output.h:118
const MeshBase *const _obj
A pointer to a constant object.
Definition: mesh_output.h:202
const bool _is_parallel_format
Flag specifying whether this format is parallel-capable.
Definition: mesh_output.h:184
libmesh_assert(ctx)
OStreamProxy out
const bool _serial_only_needed_on_proc_0
Flag specifying whether this format can be written by only serializing the mesh to processor zero...
Definition: mesh_output.h:193

◆ write_equation_systems()

void libMesh::MeshOutput< MeshBase >::write_equation_systems ( const std::string &  fname,
const EquationSystems es,
const std::set< std::string > *  system_names = nullptr 
)
virtualinherited

This method implements writing a mesh with data to a specified file where the data is taken from the EquationSystems object.

Reimplemented in libMesh::NameBasedIO.

Definition at line 31 of file mesh_output.C.

References libMesh::EquationSystems::build_solution_vector(), libMesh::EquationSystems::build_variable_names(), libMesh::EquationSystems::get_mesh(), libMesh::libmesh_assert(), and libMesh::out.

Referenced by libMesh::Nemesis_IO::write_timestep(), and libMesh::ExodusII_IO::write_timestep().

34 {
35  LOG_SCOPE("write_equation_systems()", "MeshOutput");
36 
37  // We may need to gather and/or renumber a DistributedMesh to output
38  // it, making that const qualifier in our constructor a dirty lie
39  MT & my_mesh = const_cast<MT &>(*_obj);
40 
41  // If we're asked to write data that's associated with a different
42  // mesh, output files full of garbage are the result.
43  libmesh_assert_equal_to(&es.get_mesh(), _obj);
44 
45  // A non-parallel format, non-renumbered mesh may not have a contiguous
46  // numbering, and that needs to be fixed before we can build a solution vector.
47  if (!_is_parallel_format &&
48  (my_mesh.max_elem_id() != my_mesh.n_elem() ||
49  my_mesh.max_node_id() != my_mesh.n_nodes()))
50  {
51  // If we were allowed to renumber then we should have already
52  // been properly renumbered...
53  libmesh_assert(!my_mesh.allow_renumbering());
54 
55  libmesh_do_once(libMesh::out <<
56  "Warning: This MeshOutput subclass only supports meshes which are contiguously renumbered!"
57  << std::endl;);
58 
59  my_mesh.allow_renumbering(true);
60 
61  my_mesh.renumber_nodes_and_elements();
62 
63  // Not sure what good going back to false will do here, the
64  // renumbering horses have already left the barn...
65  my_mesh.allow_renumbering(false);
66  }
67 
69  {
70  MeshSerializer serialize(const_cast<MT &>(*_obj), !_is_parallel_format, _serial_only_needed_on_proc_0);
71 
72  // Build the list of variable names that will be written.
73  std::vector<std::string> names;
74  es.build_variable_names (names, nullptr, system_names);
75 
76  // Build the nodal solution values & get the variable
77  // names from the EquationSystems object
78  std::vector<Number> soln;
79  es.build_solution_vector (soln, system_names,
80  this->get_add_sides());
81 
82  this->write_nodal_data (fname, soln, names);
83  }
84  else // _is_parallel_format
85  this->write_nodal_data (fname, es, system_names);
86 }
virtual void write_nodal_data(const std::string &, const std::vector< Number > &, const std::vector< std::string > &)
This method implements writing a mesh with nodal data to a specified file where the nodal data and va...
Definition: mesh_output.h:109
const MeshBase *const _obj
A pointer to a constant object.
Definition: mesh_output.h:202
const bool _is_parallel_format
Flag specifying whether this format is parallel-capable.
Definition: mesh_output.h:184
libmesh_assert(ctx)
OStreamProxy out
const bool _serial_only_needed_on_proc_0
Flag specifying whether this format can be written by only serializing the mesh to processor zero...
Definition: mesh_output.h:193

◆ write_nodal_data() [1/3]

virtual void libMesh::MeshOutput< MeshBase >::write_nodal_data ( const std::string &  ,
const std::vector< Number > &  ,
const std::vector< std::string > &   
)
inlinevirtualinherited

This method implements writing a mesh with nodal data to a specified file where the nodal data and variable names are provided.

Reimplemented in libMesh::ExodusII_IO, libMesh::Nemesis_IO, libMesh::GmshIO, libMesh::NameBasedIO, libMesh::VTKIO, libMesh::UCDIO, libMesh::GMVIO, libMesh::MEDITIO, libMesh::GnuPlotIO, and libMesh::TecplotIO.

Definition at line 109 of file mesh_output.h.

112  { libmesh_not_implemented(); }

◆ write_nodal_data() [2/3]

void libMesh::MeshOutput< MeshBase >::write_nodal_data ( const std::string &  fname,
const NumericVector< Number > &  parallel_soln,
const std::vector< std::string > &  names 
)
virtualinherited

This method may be overridden by "parallel" output formats for writing nodal data.

Instead of getting a localized copy of the nodal solution vector, it is passed a NumericVector of type=PARALLEL which is in node-major order i.e. (u0,v0,w0, u1,v1,w1, u2,v2,w2, u3,v3,w3, ...) and contains n_nodes*n_vars total entries. Then, it is up to the individual I/O class to extract the required solution values from this vector and write them in parallel.

If not implemented, localizes the parallel vector into a std::vector and calls the other version of this function.

Reimplemented in libMesh::Nemesis_IO.

Definition at line 149 of file mesh_output.C.

References libMesh::NumericVector< T >::localize().

152 {
153  // This is the fallback implementation for parallel I/O formats that
154  // do not yet implement proper writing in parallel, and instead rely
155  // on the full solution vector being available on all processors.
156  std::vector<Number> soln;
157  parallel_soln.localize(soln);
158  this->write_nodal_data(fname, soln, names);
159 }
virtual void write_nodal_data(const std::string &, const std::vector< Number > &, const std::vector< std::string > &)
This method implements writing a mesh with nodal data to a specified file where the nodal data and va...
Definition: mesh_output.h:109
virtual void localize(std::vector< T > &v_local) const =0
Creates a copy of the global vector in the local vector v_local.

◆ write_nodal_data() [3/3]

void libMesh::MeshOutput< MeshBase >::write_nodal_data ( const std::string &  fname,
const EquationSystems es,
const std::set< std::string > *  system_names 
)
virtualinherited

This method should be overridden by "parallel" output formats for writing nodal data.

Instead of getting a localized copy of the nodal solution vector, it directly uses EquationSystems current_local_solution vectors to look up nodal values.

If not implemented, reorders the solutions into a nodal-only NumericVector and calls the above version of this function.

Reimplemented in libMesh::Nemesis_IO.

Definition at line 162 of file mesh_output.C.

References libMesh::EquationSystems::build_parallel_solution_vector(), and libMesh::EquationSystems::build_variable_names().

165 {
166  std::vector<std::string> names;
167  es.build_variable_names (names, nullptr, system_names);
168 
169  std::unique_ptr<NumericVector<Number>> parallel_soln =
170  es.build_parallel_solution_vector(system_names);
171 
172  this->write_nodal_data (fname, *parallel_soln, names);
173 }
virtual void write_nodal_data(const std::string &, const std::vector< Number > &, const std::vector< std::string > &)
This method implements writing a mesh with nodal data to a specified file where the nodal data and va...
Definition: mesh_output.h:109

◆ write_nodal_data_discontinuous()

virtual void libMesh::MeshOutput< MeshBase >::write_nodal_data_discontinuous ( const std::string &  ,
const std::vector< Number > &  ,
const std::vector< std::string > &   
)
inlinevirtualinherited

This method implements writing a mesh with discontinuous data to a specified file where the nodal data and variables names are provided.

Reimplemented in libMesh::ExodusII_IO.

Definition at line 118 of file mesh_output.h.

121  { libmesh_not_implemented(); }

Member Data Documentation

◆ _bezier_coeffs

std::vector<Point> libMesh::PostscriptIO::_bezier_coeffs
private

Vector containing 3 points corresponding to Bezier coefficients, as computed by _compute_edge_bezier_coeffs.

Definition at line 120 of file postscript_io.h.

Referenced by _compute_edge_bezier_coeffs(), plot_quadratic_elem(), and PostscriptIO().

◆ _bezier_transform

const float libMesh::PostscriptIO::_bezier_transform
staticprivate
Initial value:
=
{
{-1.f/6.f, 1.f/6.f, 1.},
{-1.f/6.f, 0.5, 1.f/6.f},
{0., 1., 0.}
}

Coefficients of the transformation from physical-space edge coordinates to Bezier basis coefficients.

Transforms x and y separately.

Definition at line 114 of file postscript_io.h.

Referenced by _compute_edge_bezier_coeffs().

◆ _cell_string

std::ostringstream libMesh::PostscriptIO::_cell_string
private

Drawing style-independent data for a single cell.

This can be used as a temporary buffer for storing data which may be sent to the output stream multiple times.

Definition at line 142 of file postscript_io.h.

Referenced by plot_linear_elem().

◆ _current_point

Point libMesh::PostscriptIO::_current_point
private

A point object used for temporary calculations.

Definition at line 135 of file postscript_io.h.

Referenced by plot_linear_elem(), and plot_quadratic_elem().

◆ _is_parallel_format

const bool libMesh::MeshOutput< MeshBase >::_is_parallel_format
protectedinherited

Flag specifying whether this format is parallel-capable.

If this is false (default) I/O is only permitted when the mesh has been serialized.

Definition at line 184 of file mesh_output.h.

Referenced by libMesh::FroIO::write(), write(), and libMesh::EnsightIO::write().

◆ _offset

Point libMesh::PostscriptIO::_offset
private

Amount to add to every (x,y) point to place it in Postscript coordinates.

Definition at line 125 of file postscript_io.h.

Referenced by _compute_edge_bezier_coeffs(), plot_linear_elem(), plot_quadratic_elem(), and write().

◆ _out

std::ofstream libMesh::PostscriptIO::_out
private

Output file stream which will be opened when the file name is known.

Definition at line 147 of file postscript_io.h.

Referenced by plot_linear_elem(), plot_quadratic_elem(), and write().

◆ _scale

Real libMesh::PostscriptIO::_scale
private

Amount by which to stretch each point to place it in Postscript coordinates.

Definition at line 130 of file postscript_io.h.

Referenced by _compute_edge_bezier_coeffs(), plot_linear_elem(), plot_quadratic_elem(), and write().

◆ _serial_only_needed_on_proc_0

const bool libMesh::MeshOutput< MeshBase >::_serial_only_needed_on_proc_0
protectedinherited

Flag specifying whether this format can be written by only serializing the mesh to processor zero.

If this is false (default) the mesh will be serialized to all processors

Definition at line 193 of file mesh_output.h.

◆ line_width

Real libMesh::PostscriptIO::line_width

Control the thickness of the lines used.

0.5 is a reasonable default for printed images, but you may need to decrease this value (or choose it adaptively) when there are very slim cells present in the mesh.

Definition at line 88 of file postscript_io.h.

Referenced by write().

◆ shade_value

Real libMesh::PostscriptIO::shade_value

Controls greyscale shading of cells.

By default this value is 0.0 (which actually corresponds to black) and this indicates "no shading" i.e. only mesh lines will be drawn. Any other value in (0,1] will cause the cells to be grey-shaded to some degree, with higher values being lighter. A value of 0.75 gives decent results.

Definition at line 80 of file postscript_io.h.

Referenced by plot_linear_elem().


The documentation for this class was generated from the following files: