The ExodusII_IO class implements reading meshes in the ExodusII file format from Sandia National Labs. More...

#include <exodusII_io.h>

Inheritance diagram for libMesh::ExodusII_IO:

Public Member Functions
	ExodusII_IO (MeshBase &mesh, bool single_precision=false)
	Constructor. More...

	ExodusII_IO (const MeshBase &mesh, bool single_precision=false)
	Constructor. More...

	ExodusII_IO (ExodusII_IO &&)=default
	ExodusII_IO special functions: More...

	ExodusII_IO (const ExodusII_IO &)=delete

ExodusII_IO &	operator= (const ExodusII_IO &)=delete

ExodusII_IO &	operator= (ExodusII_IO &&)=delete

virtual	~ExodusII_IO ()

virtual void	read (const std::string &name) override
	This method implements reading a mesh from a specified file. More...

ExodusHeaderInfo	read_header (const std::string &name)
	Read only the header information, instead of the entire mesh. More...

virtual void	write (const std::string &fname) override
	This method implements writing a mesh to a specified file. More...

void	verbose (bool set_verbosity)
	Set the flag indicating if we should be verbose. More...

void	write_complex_magnitude (bool val)
	Set the flag indicating whether the complex modulus should be written when complex numbers are enabled. More...

void	write_added_sides (bool val)
	By default, we only write out the elements physically stored in the mesh. More...

virtual bool	get_add_sides () override

const std::vector< Real > &	get_time_steps ()

int	get_num_time_steps ()

void	copy_nodal_solution (System &system, std::string system_var_name, std::string exodus_var_name, unsigned int timestep=1)
	If we read in a nodal solution while reading in a mesh, we can attempt to copy that nodal solution into an EquationSystems object. More...

void	copy_elemental_solution (System &system, std::string system_var_name, std::string exodus_var_name, unsigned int timestep=1)
	If we read in a elemental solution while reading in a mesh, we can attempt to copy that elemental solution into an EquationSystems object. More...

void	copy_scalar_solution (System &system, std::vector< std::string > system_var_names, std::vector< std::string > exodus_var_names, unsigned int timestep=1)
	Copy global variables into scalar variables of a System object. More...

void	read_elemental_variable (std::string elemental_var_name, unsigned int timestep, std::map< unsigned int, Real > &unique_id_to_value_map)
	Given an elemental variable and a time step, returns a mapping from the elements (top parent) unique IDs to the value of the elemental variable at the corresponding time step index. More...

void	read_global_variable (std::vector< std::string > global_var_names, unsigned int timestep, std::vector< Real > &global_values)
	Given a vector of global variables and a time step, returns the values of the global variable at the corresponding time step index. More...

void	write_discontinuous_exodusII (const std::string &name, const EquationSystems &es, const std::set< std::string > *system_names=nullptr)
	Writes a exodusII file with discontinuous data. More...

void	write_timestep_discontinuous (const std::string &fname, const EquationSystems &es, const int timestep, const Real time, const std::set< std::string > *system_names=nullptr)
	Writes a discontinuous solution at a specific timestep. More...

void	write_element_data (const EquationSystems &es)
	Write out element solution. More...

void	write_element_data_from_discontinuous_nodal_data (const EquationSystems &es, const std::set< std::string > *system_names=nullptr, const std::string &var_suffix="_elem_node_")
	Similar to the function above, but instead of only handling (CONSTANT, MONOMIAL) data, writes out a general discontinuous solution field, e.g. More...

virtual void	write_nodal_data (const std::string &, const std::vector< Number > &, const std::vector< std::string > &) override
	Write out a nodal solution. More...

void	write_nodal_data_discontinuous (const std::string &, const std::vector< Number > &, const std::vector< std::string > &) override
	Write out a discontinuous nodal solution. More...

void	write_global_data (const std::vector< Number > &, const std::vector< std::string > &)
	Write out global variables. More...

void	write_information_records (const std::vector< std::string > &)
	Write out information records. More...

void	write_timestep (const std::string &fname, const EquationSystems &es, const int timestep, const Real time, const std::set< std::string > *system_names=nullptr)
	Writes out the solution at a specific timestep. More...

virtual void	write_equation_systems (const std::string &fname, const EquationSystems &es, const std::set< std::string > *system_names=nullptr) override
	Writes out the solution for no specific time or timestep. More...

void	write_elemsets ()
	Write elemsets stored on the Mesh to file. More...

void	write_sideset_data (int timestep, const std::vector< std::string > &var_names, const std::vector< std::set< boundary_id_type >> &side_ids, const std::vector< std::map< BoundaryInfo::BCTuple, Real >> &bc_vals)
	The Exodus format can also store values on sidesets. More...

void	read_sideset_data (int timestep, std::vector< std::string > &var_names, std::vector< std::set< boundary_id_type >> &side_ids, std::vector< std::map< BoundaryInfo::BCTuple, Real >> &bc_vals)
	Similar to write_sideset_data(), this function is used to read the data at a particular timestep. More...

void	get_sideset_data_indices (std::map< BoundaryInfo::BCTuple, unsigned int > &bc_array_indices)
	Similar to read_sideset_data(), but instead of creating one std::map per sideset per variable, creates a single map of (elem, side, boundary_id) tuples, and stores the exo file array indices for any/all sideset variables on that sideset (they are all the same). More...

void	write_nodeset_data (int timestep, const std::vector< std::string > &var_names, const std::vector< std::set< boundary_id_type >> &node_boundary_ids, const std::vector< std::map< BoundaryInfo::NodeBCTuple, Real >> &bc_vals)
	The Exodus format can also store values on nodesets. More...

void	read_nodeset_data (int timestep, std::vector< std::string > &var_names, std::vector< std::set< boundary_id_type >> &node_boundary_ids, std::vector< std::map< BoundaryInfo::NodeBCTuple, Real >> &bc_vals)
	Read all the nodeset data at a particular timestep. More...

void	get_nodeset_data_indices (std::map< BoundaryInfo::NodeBCTuple, unsigned int > &bc_array_indices)
	Similar to read_nodeset_data(), but instead of creating one std::map per nodeset per variable, creates a single map of (node_id, boundary_id) tuples, and stores the exo file array indices for any/all nodeset variables on that nodeset (they are all the same). More...

void	write_elemset_data (int timestep, const std::vector< std::string > &var_names, const std::vector< std::set< elemset_id_type >> &elemset_ids_in, const std::vector< std::map< std::pair< dof_id_type, elemset_id_type >, Real >> &elemset_vals)
	The Exodus format can also store values on elemsets. More...

void	read_elemset_data (int timestep, std::vector< std::string > &var_names, std::vector< std::set< elemset_id_type >> &elemset_ids_in, std::vector< std::map< std::pair< dof_id_type, elemset_id_type >, Real >> &elemset_vals)
	Read all the elemset data at a particular timestep. More...

void	get_elemset_data_indices (std::map< std::pair< dof_id_type, elemset_id_type >, unsigned int > &elemset_array_indices)
	Similar to read_elemset_data(), but instead of creating one std::map per nodeset per variable, creates a single map of (elem_id, elemset_id) pairs -> exo file array indices for any/all variables on that elemset (they are all the same). More...

void	set_extra_integer_vars (const std::vector< std::string > &extra_integer_vars)
	Set the elemental variables in the Exodus file to be read into extra element integers. More...

void	set_output_variables (const std::vector< std::string > &output_variables, bool allow_empty=true)
	Sets the list of variable names to be included in the output. More...

void	use_mesh_dimension_instead_of_spatial_dimension (bool val)
	In the general case, meshes containing 2D elements can be manifolds living in 3D space, thus by default we write all meshes with the Exodus dimension set to LIBMESH_DIM = mesh.spatial_dimension(). More...

void	write_as_dimension (unsigned dim)
	Directly control the num_dim which is written to the Exodus file. More...

void	set_coordinate_offset (Point p)
	Allows you to set a vector that is added to the coordinates of all of the nodes. More...

void	append (bool val)
	If true, this flag will cause the ExodusII_IO object to attempt to open an existing file for writing, rather than creating a new file. More...

const std::vector< std::string > &	get_elem_var_names ()
	Return list of the elemental variable names. More...

const std::vector< std::string > &	get_nodal_var_names ()
	Return list of the nodal variable names. More...

const std::vector< std::string > &	get_global_var_names ()
	Return list of the global variable names. More...

const std::vector< int > &	get_elem_num_map () const
	Returns a const reference to the elem_num_map, which is a vector that is created when a Mesh is read from file. More...

const std::vector< int > &	get_node_num_map () const
	Identical to the behavior of get_elem_num_map(), but for the node_num_map instead. More...

ExodusII_IO_Helper &	get_exio_helper ()
	Return a reference to the ExodusII_IO_Helper object. More...

void	set_hdf5_writing (bool write_hdf5)
	Set to true (the default) to write files in an HDF5-based file format (when HDF5 is available), or to false to write files in the old NetCDF3-based format. More...

void	set_discontinuous_bex (bool disc_bex)
	Set to true (false is the default) to generate independent nodes for every Bezier Extraction element in an input file containing them, causing those elements to be disconnected from each other. More...

void	write_nodal_data_common (std::string fname, const std::vector< std::string > &names, bool continuous=true)
	This function factors out a bunch of code which is common to the write_nodal_data() and write_nodal_data_discontinuous() functions. More...

bool	is_parallel_format () const
	Returns true iff this mesh file format and input class are parallelized, so that all processors can read their share of the data at once. 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 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...

unsigned int &	ascii_precision ()
	Return/set the precision to use when writing ASCII files. More...

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static int	get_exodus_version ()

Protected Member Functions
MeshBase &	mesh ()

void	set_n_partitions (unsigned int n_parts)
	Sets the number of partitions in the mesh. More...

void	skip_comment_lines (std::istream &in, const char comment_start)
	Reads input from `in`, skipping all the lines that start with the character `comment_start`. More...

const MeshBase &	mesh () const

Protected Attributes
std::vector< bool >	elems_of_dimension
	A vector of bools describing what dimension elements have been encountered when reading a mesh. More...

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...

const Parallel::Communicator &	_communicator

Private Attributes
std::unique_ptr< ExodusII_IO_Helper >	exio_helper
	Only attempt to instantiate an ExodusII helper class if the Exodus API is defined. More...

int	_timestep
	Stores the current value of the timestep when calling ExodusII_IO::write_timestep(). More...

bool	_verbose
	should we be verbose? More...

bool	_append
	Default false. More...

std::vector< std::string >	_extra_integer_vars
	An optional list of variables in the EXODUS file that are to be used to set extra integers when loading the file into a mesh. More...

std::vector< std::string >	_output_variables
	The names of the variables to be output. More...

bool	_allow_empty_variables
	Flag which controls the behavior of _output_variables: .) If true, _output_variables is allowed to remain empty. More...

bool	_write_complex_abs
	By default, when complex numbers are enabled, for each variable we write out three values: the real part, "r_u" the imaginary part, "i_u", and the complex modulus, a_u := sqrt(r_ur_u + i_ui_u), which is also the value returned by std::abs(std::complex). More...

bool	_disc_bex
	Set to true (false is the default) to generate independent nodes for every Bezier Extraction element. More...

Detailed Description

The ExodusII_IO class implements reading meshes in the ExodusII file format from Sandia National Labs.

By default, LibMesh expects ExodusII files to have a ".exd" or ".e" file extension.

Author: Benjamin Kirk; John Peterson

Date: 2004 Handles reading and writing of Exodus binary files.

Definition at line 52 of file exodusII_io.h.

Constructor & Destructor Documentation

◆ ExodusII_IO() [1/4]

libMesh::ExodusII_IO::ExodusII_IO	(	MeshBase &	mesh,
		bool	single_precision = `false`
	)

explicit

Constructor.

Takes a writable reference to a mesh object. This is the constructor required to read a mesh.

Definition at line 122 of file exodusII_io.C.

References libMesh::libmesh_ignore().

                                                  :
   MeshInput<MeshBase> (mesh),
   MeshOutput<MeshBase> (mesh,
                         /* is_parallel_format = */ false,
                         /* serial_only_needed_on_proc_0 = */ true),
   ParallelObject(mesh),
 #ifdef LIBMESH_HAVE_EXODUS_API
   exio_helper(std::make_unique<ExodusII_IO_Helper>(*this, false, true, single_precision)),
   _timestep(1),
   _verbose(false),
   _append(false),
 #endif
   _allow_empty_variables(false),
   _write_complex_abs(true),
   _disc_bex(false)
 {
   // if !LIBMESH_HAVE_EXODUS_API, we didn't use this
   libmesh_ignore(single_precision);
 }

◆ ExodusII_IO() [2/4]

libMesh::ExodusII_IO::ExodusII_IO	(	const MeshBase &	mesh,
		bool	single_precision = `false`
	)

explicit

Constructor.

Takes a const reference to a mesh object. This constructor is acceptable to write meshes.

Definition at line 145 of file exodusII_io.C.

References libMesh::libmesh_ignore().

                                                  :
   MeshInput<MeshBase> (),
   MeshOutput<MeshBase> (mesh,
                         /* is_parallel_format = */ false,
                         /* serial_only_needed_on_proc_0 = */ true),
   ParallelObject(mesh),
 #ifdef LIBMESH_HAVE_EXODUS_API
   exio_helper(std::make_unique<ExodusII_IO_Helper>(*this, false, true, single_precision)),
   _timestep(1),
   _verbose(false),
   _append(false),
 #endif
   _allow_empty_variables(false),
   _write_complex_abs(true),
   _disc_bex(false)
 {
   // if !LIBMESH_HAVE_EXODUS_API, we didn't use this
   libmesh_ignore(single_precision);
 }

◆ ExodusII_IO() [3/4]

libMesh::ExodusII_IO::ExodusII_IO ( ExodusII_IO && )

default

ExodusII_IO special functions:

Can't be (default) copy constructed or assigned since they contain a unique_ptr member.
Can't be (default) move assigned because the helper class contains references.
The destructor is responsible for closing the file

◆ ExodusII_IO() [4/4]

libMesh::ExodusII_IO::ExodusII_IO ( const ExodusII_IO & )

delete

◆ ~ExodusII_IO()

libMesh::ExodusII_IO::~ExodusII_IO ( )

virtualdefault

Definition at line 238 of file exodusII_io.C.

References exio_helper.

 {
   exio_helper->close();
 }

Member Function Documentation

◆ append()

void libMesh::ExodusII_IO::append ( bool val )

If true, this flag will cause the ExodusII_IO object to attempt to open an existing file for writing, rather than creating a new file.

Obviously this will only work if the file already exists.

Definition at line 998 of file exodusII_io.C.

References _append.

Referenced by main().

 {
   _append = val;
 }

◆ 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::STLIO::write(), libMesh::MEDITIO::write_ascii(), libMesh::TecplotIO::write_ascii(), libMesh::GMVIO::write_ascii_new_impl(), and libMesh::GMVIO::write_ascii_old_impl().

 {
   return _ascii_precision;
 }

◆ comm()

const Parallel::Communicator& libMesh::ParallelObject::comm ( ) const

inlineinherited

Returns: A reference to the Parallel::Communicator object used by this mesh.

Definition at line 97 of file parallel_object.h.

References libMesh::ParallelObject::_communicator.

98 { return _communicator; }

libMesh::ParallelObject::_communicator

const Parallel::Communicator & _communicator

Definition: parallel_object.h:120

◆ copy_elemental_solution()

void libMesh::ExodusII_IO::copy_elemental_solution	(	System &	system,
		std::string	system_var_name,
		std::string	exodus_var_name,
		unsigned int	timestep = `1`
	)

If we read in a elemental solution while reading in a mesh, we can attempt to copy that elemental solution into an EquationSystems object.

Definition at line 1147 of file exodusII_io.C.

References libMesh::ParallelObject::comm(), libMesh::CONSTANT, libMesh::DofObject::dof_number(), exio_helper, libMesh::System::get_dof_map(), libMesh::MeshBase::is_serial_on_zero(), libMesh::libmesh_isnan(), libMesh::DofMap::local_index(), libMesh::MeshInput< MeshBase >::mesh(), libMesh::MeshInput< MT >::mesh(), libMesh::MONOMIAL, libMesh::MONOMIAL_VEC, libMesh::DofObject::n_comp(), libMesh::System::number(), libMesh::ParallelObject::processor_id(), libMesh::MeshBase::query_elem_ptr(), TIMPI::Communicator::rank(), libMesh::Real, libMesh::System::solution, libMesh::System::update(), libMesh::System::variable_number(), and libMesh::System::variable_type().

Referenced by MeshInputTest::testExodusWriteElementDataFromDiscontinuousNodalData().

 {
   LOG_SCOPE("copy_elemental_solution()", "ExodusII_IO");
 
   const unsigned int var_num = system.variable_number(system_var_name);
   // Assert that variable is an elemental one.
   //
   // NOTE: Currently, this reader is capable of reading only individual components of MONOMIAL_VEC
   //       types, and each must be written out to its own CONSTANT MONOMIAL variable
   libmesh_error_msg_if((system.variable_type(var_num) != FEType(CONSTANT, MONOMIAL))
                        && (system.variable_type(var_num) != FEType(CONSTANT, MONOMIAL_VEC)),
                        "Error! Trying to copy elemental solution into a variable that is not of CONSTANT MONOMIAL nor CONSTANT MONOMIAL_VEC type.");
 
   const MeshBase & mesh = MeshInput<MeshBase>::mesh();
   const DofMap & dof_map = system.get_dof_map();
 
   // Map from element ID to elemental variable value.  We need to use
   // a map here rather than a vector (e.g. elem_var_values) since the
   // libmesh element numbering can contain "holes".  This is the case
   // if we are reading elemental var values from an adaptively refined
   // mesh that has not been sequentially renumbered.
   std::map<dof_id_type, Real> elem_var_value_map;
 
   // With Exodus files we only open them on processor 0, so that's the
   // where we have to do the data read too.
   if (system.comm().rank() == 0)
     {
       libmesh_error_msg_if(!exio_helper->opened_for_reading,
                            "ERROR, ExodusII file must be opened for reading before copying an elemental solution!");
 
       exio_helper->read_elemental_var_values(exodus_var_name, timestep, elem_var_value_map);
     }
 
   const bool serial_on_zero = mesh.is_serial_on_zero();
 
   // If our mesh isn't serial, then non-root processors need to
   // request the data for their parts of the mesh and insert it
   // themselves.
   if (!serial_on_zero)
     {
       std::unordered_map<processor_id_type, std::vector<dof_id_type>> elem_ids_to_request;
       if (this->processor_id() != 0)
         {
           std::vector<dof_id_type> elem_ids;
           for (auto & elem : mesh.active_local_element_ptr_range())
             elem_ids.push_back(elem->id());
 
           if (!elem_ids.empty())
             elem_ids_to_request[0] = std::move(elem_ids);
         }
 
       auto value_gather_functor =
         [& elem_var_value_map]
         (processor_id_type,
          const std::vector<dof_id_type> & ids,
          std::vector<Real> & values)
         {
           const std::size_t query_size = ids.size();
           values.resize(query_size);
           for (std::size_t i=0; i != query_size; ++i)
             {
               if (const auto it = elem_var_value_map.find(ids[i]);
                   it != elem_var_value_map.end())
                 {
                   values[i] = it->second;
                   elem_var_value_map.erase(it);
                 }
               else
                 values[i] = std::numeric_limits<Real>::quiet_NaN();
             }
         };
 
       auto value_action_functor =
         [& elem_var_value_map]
         (processor_id_type,
          const std::vector<dof_id_type> & ids,
          const std::vector<Real> & values)
         {
           const std::size_t query_size = ids.size();
           for (std::size_t i=0; i != query_size; ++i)
             if (!libmesh_isnan(values[i]))
               elem_var_value_map[ids[i]] = values[i];
         };
 
       Real * value_ex = nullptr;
       Parallel::pull_parallel_vector_data
         (system.comm(), elem_ids_to_request, value_gather_functor,
          value_action_functor, value_ex);
     }
 
   std::map<dof_id_type, Real>::iterator
     it = elem_var_value_map.begin(),
     end = elem_var_value_map.end();
 
   // Everybody inserts the data they've received.  If we're
   // serial_on_zero then proc 0 inserts everybody's data and other
   // procs have empty map ranges.
   for (; it!=end; ++it)
     {
       const Elem * elem = mesh.query_elem_ptr(it->first);
 
       if (elem && elem->n_comp(system.number(), var_num) > 0)
         {
           dof_id_type dof_index = elem->dof_number(system.number(), var_num, 0);
           if (serial_on_zero || dof_map.local_index(dof_index ))
             system.solution->set (dof_index, it->second);
         }
     }
 
   system.solution->close();
   system.update();
 }

◆ copy_nodal_solution()

void libMesh::ExodusII_IO::copy_nodal_solution	(	System &	system,
		std::string	system_var_name,
		std::string	exodus_var_name,
		unsigned int	timestep = `1`
	)

If we read in a nodal solution while reading in a mesh, we can attempt to copy that nodal solution into an EquationSystems object.

Definition at line 1043 of file exodusII_io.C.

References libMesh::ParallelObject::comm(), libMesh::DofObject::dof_number(), exio_helper, libMesh::MeshBase::is_serial_on_zero(), libMesh::libmesh_isnan(), libMesh::MeshInput< MeshBase >::mesh(), libMesh::MeshInput< MT >::mesh(), libMesh::DofObject::n_comp(), libMesh::System::number(), libMesh::ParallelObject::processor_id(), libMesh::DofObject::processor_id(), TIMPI::Communicator::rank(), libMesh::Real, libMesh::System::solution, libMesh::System::update(), and libMesh::System::variable_number().

Referenced by WriteVecAndScalar::testWriteExodus().

 {
   LOG_SCOPE("copy_nodal_solution()", "ExodusII_IO");
 
   const unsigned int var_num = system.variable_number(system_var_name);
 
   const MeshBase & mesh = MeshInput<MeshBase>::mesh();
 
   // With Exodus files we only open them on processor 0, so that's the
   // where we have to do the data read too.
   if (system.comm().rank() == 0)
     {
       libmesh_error_msg_if(!exio_helper->opened_for_reading,
                            "ERROR, ExodusII file must be opened for reading before copying a nodal solution!");
 
       exio_helper->read_nodal_var_values(exodus_var_name, timestep);
     }
 
   auto & node_var_value_map = exio_helper->nodal_var_values;
 
   const bool serial_on_zero = mesh.is_serial_on_zero();
 
   // If our mesh isn't serial, then non-root processors need to
   // request the data for their parts of the mesh and insert it
   // themselves.
   if (!serial_on_zero)
     {
       std::unordered_map<processor_id_type, std::vector<dof_id_type>> node_ids_to_request;
       if (this->processor_id() != 0)
         {
           std::vector<dof_id_type> node_ids;
           for (auto & node : mesh.local_node_ptr_range())
             node_ids.push_back(node->id());
           if (!node_ids.empty())
             node_ids_to_request[0] = std::move(node_ids);
         }
 
       auto value_gather_functor =
         [& node_var_value_map]
         (processor_id_type,
          const std::vector<dof_id_type> & ids,
          std::vector<Real> & values)
         {
           const std::size_t query_size = ids.size();
           values.resize(query_size);
           for (std::size_t i=0; i != query_size; ++i)
             {
               if (const auto it = node_var_value_map.find(ids[i]);
                   it != node_var_value_map.end())
                 {
                   values[i] = it->second;
                   node_var_value_map.erase(it);
                 }
               else
                 values[i] = std::numeric_limits<Real>::quiet_NaN();
             }
         };
 
       auto value_action_functor =
         [& node_var_value_map]
         (processor_id_type,
          const std::vector<dof_id_type> & ids,
          const std::vector<Real> & values)
         {
           const std::size_t query_size = ids.size();
           for (std::size_t i=0; i != query_size; ++i)
             if (!libmesh_isnan(values[i]))
               node_var_value_map[ids[i]] = values[i];
         };
 
       Real * value_ex = nullptr;
       Parallel::pull_parallel_vector_data
         (system.comm(), node_ids_to_request, value_gather_functor,
          value_action_functor, value_ex);
     }
 
   // Everybody inserts the data they've received.  If we're
   // serial_on_zero then proc 0 inserts everybody's data and other
   // procs have empty map ranges.
   for (auto p : exio_helper->nodal_var_values)
     {
       dof_id_type i = p.first;
       const Node * node = MeshInput<MeshBase>::mesh().query_node_ptr(i);
 
       if (node &&
           (serial_on_zero || node->processor_id() == system.processor_id()) &&
           node->n_comp(system.number(), var_num) > 0)
         {
           dof_id_type dof_index = node->dof_number(system.number(), var_num, 0);
 
           // If the dof_index is local to this processor, set the value
           system.solution->set (dof_index, p.second);
         }
     }
 
   system.solution->close();
   system.update();
 }

◆ copy_scalar_solution()

void libMesh::ExodusII_IO::copy_scalar_solution	(	System &	system,
		std::vector< std::string >	system_var_names,
		std::vector< std::string >	exodus_var_names,
		unsigned int	timestep = `1`
	)

Copy global variables into scalar variables of a System object.

Definition at line 1263 of file exodusII_io.C.

References libMesh::ParallelObject::comm(), exio_helper, libMesh::System::get_dof_map(), TIMPI::Communicator::get_unique_tag(), libMesh::index_range(), libMesh::ParallelObject::n_processors(), libMesh::ParallelObject::processor_id(), read_global_variable(), TIMPI::Communicator::receive(), libMesh::DofMap::SCALAR_dof_indices(), TIMPI::Communicator::send(), libMesh::System::solution, libMesh::System::update(), and libMesh::System::variable_scalar_number().

 {
   LOG_SCOPE("copy_scalar_solution()", "ExodusII_IO");
 
   libmesh_error_msg_if(!exio_helper->opened_for_reading,
                        "ERROR, ExodusII file must be opened for reading before copying a scalar solution!");
 
   libmesh_error_msg_if(system_var_names.size() != exodus_var_names.size(),
                        "ERROR, the number of system_var_names must match exodus_var_names.");
 
   std::vector<Real> values_from_exodus;
   read_global_variable(exodus_var_names, timestep, values_from_exodus);
 
 #ifdef LIBMESH_HAVE_MPI
   if (this->n_processors() > 1)
   {
     const Parallel::MessageTag tag = this->comm().get_unique_tag(1);
     if (this->processor_id() == this->n_processors()-1)
       this->comm().receive(0, values_from_exodus, tag);
     if (this->processor_id() == 0)
       this->comm().send(this->n_processors()-1, values_from_exodus, tag);
   }
 #endif
 
   if (system.processor_id() == (system.n_processors()-1))
   {
     const DofMap & dof_map = system.get_dof_map();
 
     for (auto i : index_range(system_var_names))
     {
       const unsigned int var_num = system.variable_scalar_number(system_var_names[i], 0);
 
       std::vector<dof_id_type> SCALAR_dofs;
       dof_map.SCALAR_dof_indices(SCALAR_dofs, var_num);
 
       system.solution->set (SCALAR_dofs[0], values_from_exodus[i]);
     }
   }
 
   system.solution->close();
   system.update();
 }

◆ get_add_sides()

bool libMesh::ExodusII_IO::get_add_sides ( )

overridevirtual

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 from libMesh::MeshOutput< MeshBase >.

Definition at line 1012 of file exodusII_io.C.

References exio_helper.

Referenced by write_discontinuous_exodusII(), and write_nodal_data_discontinuous().

 {
   return exio_helper->get_add_sides();
 }

◆ get_elem_num_map()

const std::vector< int > & libMesh::ExodusII_IO::get_elem_num_map ( ) const

Returns a const reference to the elem_num_map, which is a vector that is created when a Mesh is read from file.

LibMesh will number its mesh elements consistently with the elem_num_map array, except that the indices in this array are 1-based, and libmesh always uses a 0-based numbering. For example, given: elem_num_map = [4,2,3,1] libmesh will assign the first element it reads from the Exodus file elem->id() == 3, the second will get elem->id() == 1, and so on. We note that not all Exodus files contain an elem_num_map, and in that case, calling this function will return a reference to a vector containing the 1-based identity array, [1,2,3,...]

Definition at line 2381 of file exodusII_io.C.

References exio_helper.

Referenced by ExtraIntegersTest::testExtraIntegersExodusReading().

 {
   // We could make this function non-const and have it call
   // exio_helper->read_elem_num_map() before returning a reference...
   // but the intention is that this will be called some time after a
   // mesh is read in, in which case it would be doing extra work to
   // read in the elem_num_map twice.
   return exio_helper->elem_num_map;
 }

◆ get_elem_var_names()

const std::vector< std::string > & libMesh::ExodusII_IO::get_elem_var_names ( )

Return list of the elemental variable names.

Definition at line 2369 of file exodusII_io.C.

References libMesh::ExodusII_IO_Helper::ELEMENTAL, and exio_helper.

 {
   exio_helper->read_var_names(ExodusII_IO_Helper::ELEMENTAL);
   return exio_helper->elem_var_names;
 }

◆ get_elemset_data_indices()

void libMesh::ExodusII_IO::get_elemset_data_indices ( std::map< std::pair< dof_id_type, elemset_id_type >, unsigned int > & elemset_array_indices )

Similar to read_elemset_data(), but instead of creating one std::map per nodeset per variable, creates a single map of (elem_id, elemset_id) pairs -> exo file array indices for any/all variables on that elemset (they are all the same).

In cases where there are hundreds of elemset variables on a single elemset, it is more efficient to store the array indices in a quickly searchable data structure than to repeat the indexing once per variable as is done in the read_elemset_data() case.

Definition at line 2170 of file exodusII_io.C.

References exio_helper.

 {
   libmesh_error_msg_if(!exio_helper->opened_for_reading,
                        "ERROR, ExodusII file must be opened for reading "
                        "before calling ExodusII_IO::get_elemset_data_indices()!");
 
   exio_helper->get_elemset_data_indices(elemset_array_indices);
 }

◆ get_exio_helper()

ExodusII_IO_Helper & libMesh::ExodusII_IO::get_exio_helper ( )

Return a reference to the ExodusII_IO_Helper object.

Definition at line 2396 of file exodusII_io.C.

References exio_helper.

 {
   // Provide a warning when accessing the helper object
   // since it is a non-public API and is likely to see
   // future API changes
   libmesh_experimental();
 
   return *exio_helper;
 }

◆ get_exodus_version()

int libMesh::ExodusII_IO::get_exodus_version ( )

static

Returns: The ExodusII API version, in "nodot" format (e.g. 822 for 8.22), or 0 if ExodusII support is not available.

Definition at line 168 of file exodusII_io.C.

References libMesh::ExodusII_IO_Helper::get_exodus_version().

Referenced by MeshInputTest::testExodusFileMappings(), and MeshInputTest::testExodusIGASidesets().

 {
 #ifdef LIBMESH_HAVE_EXODUS_API
   return ExodusII_IO_Helper::get_exodus_version();
 #else
   return 0;
 #endif
 }

◆ get_global_var_names()

const std::vector< std::string > & libMesh::ExodusII_IO::get_global_var_names ( )

Return list of the global variable names.

Definition at line 2375 of file exodusII_io.C.

References exio_helper, and libMesh::ExodusII_IO_Helper::GLOBAL.

 {
   exio_helper->read_var_names(ExodusII_IO_Helper::GLOBAL);
   return exio_helper->global_var_names;
 }

◆ get_nodal_var_names()

const std::vector< std::string > & libMesh::ExodusII_IO::get_nodal_var_names ( )

Return list of the nodal variable names.

Definition at line 2363 of file exodusII_io.C.

References exio_helper, and libMesh::ExodusII_IO_Helper::NODAL.

Referenced by WriteVecAndScalar::testWriteExodus().

 {
   exio_helper->read_var_names(ExodusII_IO_Helper::NODAL);
   return exio_helper->nodal_var_names;
 }

◆ get_node_num_map()

const std::vector< int > & libMesh::ExodusII_IO::get_node_num_map ( ) const

Identical to the behavior of get_elem_num_map(), but for the node_num_map instead.

Definition at line 2391 of file exodusII_io.C.

References exio_helper.

Referenced by ExtraIntegersTest::testExtraIntegersExodusReading().

 {
   return exio_helper->node_num_map;
 }

◆ get_nodeset_data_indices()

void libMesh::ExodusII_IO::get_nodeset_data_indices ( std::map< BoundaryInfo::NodeBCTuple, unsigned int > & bc_array_indices )

Similar to read_nodeset_data(), but instead of creating one std::map per nodeset per variable, creates a single map of (node_id, boundary_id) tuples, and stores the exo file array indices for any/all nodeset variables on that nodeset (they are all the same).

In cases where there are hundreds of nodeset variables on a single nodeset, it is more efficient to store the array indices in a quickly searchable data structure than to repeat the indexing once per variable as is done in the read_nodeset_data() case.

Definition at line 2100 of file exodusII_io.C.

References exio_helper.

 {
   libmesh_error_msg_if(!exio_helper->opened_for_reading,
                        "ERROR, ExodusII file must be opened for reading "
                        "before calling ExodusII_IO::get_nodeset_data_indices()!");
 
   exio_helper->get_nodeset_data_indices(bc_array_indices);
 }

◆ get_num_time_steps()

int libMesh::ExodusII_IO::get_num_time_steps ( )

Returns: The number of timesteps currently stored in the Exodus file.

Knowing the number of time steps currently stored in the file is sometimes necessary when appending, so we can know where to start writing new data. Throws an error if the file is not currently open for reading or writing.

Definition at line 1032 of file exodusII_io.C.

References exio_helper.

 {
   libmesh_error_msg_if(!exio_helper->opened_for_reading && !exio_helper->opened_for_writing,
                        "ERROR, ExodusII file must be opened for reading or writing before calling ExodusII_IO::get_num_time_steps()!");
 
   exio_helper->read_num_time_steps();
   return exio_helper->num_time_steps;
 }

◆ get_sideset_data_indices()

void libMesh::ExodusII_IO::get_sideset_data_indices ( std::map< BoundaryInfo::BCTuple, unsigned int > & bc_array_indices )

Similar to read_sideset_data(), but instead of creating one std::map per sideset per variable, creates a single map of (elem, side, boundary_id) tuples, and stores the exo file array indices for any/all sideset variables on that sideset (they are all the same).

In cases where there are hundreds of sideset variables on a single sideset, it is more efficient to store the array indices in a quickly searchable data structure than to repeat the indexing once per variable as is done in the read_sideset_data() case.

Definition at line 2087 of file exodusII_io.C.

References exio_helper, libMesh::MeshInput< MeshBase >::mesh(), and libMesh::MeshOutput< MT >::mesh().

 {
   libmesh_error_msg_if(!exio_helper->opened_for_reading,
                        "ERROR, ExodusII file must be opened for reading "
                        "before calling ExodusII_IO::get_sideset_data_indices()!");
 
   const MeshBase & mesh = MeshOutput<MeshBase>::mesh();
   exio_helper->get_sideset_data_indices(mesh, bc_array_indices);
 }

◆ get_time_steps()

const std::vector< Real > & libMesh::ExodusII_IO::get_time_steps ( )

Returns: An array containing the timesteps in the file.

Definition at line 1020 of file exodusII_io.C.

References exio_helper.

 {
   libmesh_error_msg_if
     (!exio_helper->opened_for_reading,
      "ERROR, ExodusII file must be opened for reading before calling ExodusII_IO::get_time_steps()!");
 
   exio_helper->read_time_steps();
   return exio_helper->time_steps;
 }

◆ is_parallel_format()

bool libMesh::MeshInput< MeshBase >::is_parallel_format ( ) const

inlineinherited

Returns true iff this mesh file format and input class are parallelized, so that all processors can read their share of the data at once.

Definition at line 87 of file mesh_input.h.

References libMesh::MeshInput< MT >::_is_parallel_format.

87 { return this->_is_parallel_format; }

libMesh::MeshInput< MeshBase >::_is_parallel_format

const bool _is_parallel_format

Flag specifying whether this format is parallel-capable.

Definition: mesh_input.h:130

◆ mesh() [1/2]

MeshBase & libMesh::MeshInput< MeshBase >::mesh ( )

inlineprotectedinherited

Returns: The object as a writable reference.

Definition at line 178 of file mesh_input.h.

 {
   libmesh_error_msg_if(_obj == nullptr, "ERROR: _obj should not be nullptr!");
   return *_obj;
 }

◆ mesh() [2/2]

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

inlineprotectedinherited

Returns: The object as a read-only reference.

Definition at line 259 of file mesh_output.h.

References libMesh::libmesh_assert().

Referenced by libMesh::FroIO::write(), libMesh::TecplotIO::write(), libMesh::MEDITIO::write(), libMesh::PostscriptIO::write(), libMesh::EnsightIO::write(), libMesh::TecplotIO::write_ascii(), libMesh::TecplotIO::write_binary(), libMesh::TecplotIO::write_nodal_data(), libMesh::MEDITIO::write_nodal_data(), and libMesh::GnuPlotIO::write_solution().

 {
   libmesh_assert(_obj);
   return *_obj;
 }

◆ n_processors()

processor_id_type libMesh::ParallelObject::n_processors ( ) const

inlineinherited

Returns: The number of processors in the group.

Definition at line 103 of file parallel_object.h.

References libMesh::ParallelObject::_communicator, libMesh::libmesh_assert(), and TIMPI::Communicator::size().

   {
     processor_id_type returnval =
       cast_int<processor_id_type>(_communicator.size());
     libmesh_assert(returnval); // We never have an empty comm
     return returnval;
   }

◆ operator=() [1/2]

ExodusII_IO& libMesh::ExodusII_IO::operator= ( const ExodusII_IO & )

delete

◆ operator=() [2/2]

ExodusII_IO& libMesh::ExodusII_IO::operator= ( ExodusII_IO && )

delete

◆ processor_id()

processor_id_type libMesh::ParallelObject::processor_id ( ) const

inlineinherited

Returns: The rank of this processor in the group.

Definition at line 114 of file parallel_object.h.

References libMesh::ParallelObject::_communicator, and TIMPI::Communicator::rank().

115 { return cast_int<processor_id_type>(_communicator.rank()); }

TIMPI::Communicator::rank

processor_id_type rank() const

libMesh::ParallelObject::_communicator

const Parallel::Communicator & _communicator

Definition: parallel_object.h:120

◆ read()

void libMesh::ExodusII_IO::read ( const std::string & name )

overridevirtual

This method implements reading a mesh from a specified file.

Open the file named name and read the mesh in Sandia National Lab's ExodusII format. This is the method to use for reading in meshes generated by cubit. Works in 2D for TRIs, TRI6s, QUAD s, and QUAD9s. Works in 3D for TET4s, TET10s, HEX8s, and HEX27s.

Implements libMesh::MeshInput< MeshBase >.

Definition at line 244 of file exodusII_io.C.

Referenced by libMesh::NameBasedIO::read(), ExodusTest< elem_type >::test_read_gold(), ExodusTest< elem_type >::test_write(), ExtraIntegersTest::testBadExtraIntegersExodusReading(), MeshInputTest::testExodusFileMappings(), MeshInputTest::testExodusIGASidesets(), MeshInputTest::testExodusWriteElementDataFromDiscontinuousNodalData(), ExtraIntegersTest::testExtraIntegersExodusReading(), MeshInputTest::testLowOrderEdgeBlocks(), and WriteVecAndScalar::testWriteExodus().

 {
   LOG_SCOPE("read()", "ExodusII_IO");
 
   // Get a reference to the mesh we are reading
   MeshBase & mesh = MeshInput<MeshBase>::mesh();
 
   // Add extra integers into the mesh
   std::vector<unsigned int> extra_ids;
   for (auto & name : _extra_integer_vars)
     extra_ids.push_back(mesh.add_elem_integer(name));
 
   // Clear any existing mesh data
   mesh.clear();
 
   // Keep track of what kinds of elements this file contains
   elems_of_dimension.clear();
   elems_of_dimension.resize(4, false);
 
   // Open the exodus file in EX_READ mode
   exio_helper->open(fname.c_str(), /*read_only=*/true);
 
   // Get header information from exodus file
   exio_helper->read_and_store_header_info();
 
   // Read the QA records
   exio_helper->read_qa_records();
 
   // Print header information
   exio_helper->print_header();
 
   // Read nodes from the exodus file
   exio_helper->read_nodes();
 
   // Reserve space for the nodes.
   mesh.reserve_nodes(exio_helper->num_nodes);
 
   // Read the node number map from the Exodus file.  This is
   // required if we want to preserve the numbering of nodes as it
   // exists in the Exodus file.  If the Exodus file does not contain
   // a node_num_map, the identity map is returned by this call.
   exio_helper->read_node_num_map();
 
   // Read the element number map from the Exodus file.  This is
   // required if we want to preserve the numbering of elements as it
   // exists in the Exodus file.  If the Exodus file does not contain
   // an elem_num_map, the identity map is returned by this call.
   //
   // We now do this before creating nodes, so if we have any spline
   // nodes that need a NodeElem attached we can give them unused elem
   // ids.
   exio_helper->read_elem_num_map();
 
   // Read any Bezier Extraction coefficient vectors from the file,
   // such as might occur in an IsoGeometric Analysis (IGA) mesh.
   exio_helper->read_bex_cv_blocks();
 
   // If we have Rational Bezier weights, we'll need to
   // store them.
   unsigned char weight_index = 0;
   const bool weights_exist = !exio_helper->w.empty();
 
   // If we have Bezier extraction coefficients, we'll need to put
   // NODEELEM elements on spline nodes, since our Rational Bezier
   // elements will be connected to nodes derived from those; nothing
   // else will be directly connected to the spline nodes.
   const bool bex_cv_exist = !exio_helper->bex_dense_constraint_vecs.empty();
 
   // Even if weights don't exist, we still use RATIONAL_BERNSTEIN for
   // Bezier-Bernstein BEX elements, we just use 1.0 as the weight on
   // every node.
   if (bex_cv_exist)
     {
       const Real default_weight = 1.0;
       weight_index = cast_int<unsigned char>
         (mesh.add_node_datum<Real>("rational_weight", true,
                                    &default_weight));
       mesh.set_default_mapping_type(RATIONAL_BERNSTEIN_MAP);
       mesh.set_default_mapping_data(weight_index);
     }
 
   std::unordered_map<const Node *, Elem *> spline_nodeelem_ptrs;
 
   // Loop over the nodes, create Nodes with local processor_id 0.
   for (int i=0; i<exio_helper->num_nodes; i++)
     {
       // Use the node_num_map to get the correct ID for Exodus
       int exodus_id = exio_helper->node_num_map[i];
 
       // Catch the node that was added to the mesh
       Node * added_node = mesh.add_point (Point(exio_helper->x[i], exio_helper->y[i], exio_helper->z[i]), exodus_id-1);
 
       // If the Mesh assigned an ID different from what is in the
       // Exodus file, we should probably error.
       libmesh_error_msg_if(added_node->id() != static_cast<unsigned>(exodus_id-1),
                            "Error!  Mesh assigned node ID "
                            << added_node->id()
                            << " which is different from the (zero-based) Exodus ID "
                            << exodus_id-1
                            << "!");
 
       // If we have a set of spline weights, these nodes are going to
       // be used as control points for Bezier elements, and we need
       // to attach a NodeElem to each to make sure it doesn't get
       // flagged as an unused node.
       if (weights_exist)
         {
           const auto w = exio_helper->w[i];
           Point & p = *added_node;
           p /= w;  // Exodus Bezier Extraction stores spline nodes in projective space
 
           added_node->set_extra_datum<Real>(weight_index, exio_helper->w[i]);
         }
 
       if (bex_cv_exist)
         {
           std::unique_ptr<Elem> elem = Elem::build(NODEELEM);
 
           // Give the NodeElem ids at the end, so we can match any
           // existing ids in the file for other elements
           elem->set_id() = exio_helper->end_elem_id() + i;
 
           elem->set_node(0, added_node);
           Elem * added_elem = mesh.add_elem(std::move(elem));
           spline_nodeelem_ptrs[added_node] = added_elem;
         }
     }
 
   // This assert is no longer valid if the nodes are not numbered
   // sequentially starting from 1 in the Exodus file.
   // libmesh_assert_equal_to (static_cast<unsigned int>(exio_helper->num_nodes), mesh.n_nodes());
 
   // Get information about all the element and edge blocks
   exio_helper->read_block_info();
 
   // Reserve space for the elements.  Account for any NodeElem that
   // have already been attached to spline control nodes.
   mesh.reserve_elem(exio_helper->w.size() + exio_helper->num_elem);
 
   // Read variables for extra integer IDs
   std::vector<std::map<dof_id_type, Real>> elem_ids(extra_ids.size());
   // We use the last time step to load the IDs
   exio_helper->read_num_time_steps();
   unsigned int last_step = exio_helper->num_time_steps;
   for (auto i : index_range(extra_ids))
     exio_helper->read_elemental_var_values(_extra_integer_vars[i], last_step, elem_ids[i]);
 
   // Read in the element connectivity for each block.
   int nelem_last_block = 0;
 
   // If we're building Bezier elements from spline nodes, we need to
   // calculate those elements' local nodes on the fly, and we'll be
   // calculating them from constraint matrix columns, and we'll need
   // to make sure that the same node is found each time it's
   // calculated from multiple neighboring elements.
   std::map<std::vector<std::pair<dof_id_type, Real>>, Node *> local_nodes;
 
   // We'll set any spline NodeElem subdomain_id() values to exceed the
   // maximum of subdomain_id() values set via Exodus block ids.
   int max_subdomain_id = std::numeric_limits<int>::min();
 
   // We've already added all the nodes explicitly specified in the
   // file, but if we have spline nodes we may need to add assembly
   // element nodes based on them.  Add them contiguously so we're
   // compatible with any subsequent code paths (including our
   // ExodusII_IO::write()!) that don't support sparse ids.
   dof_id_type n_nodes = mesh.n_nodes();
 
   // Loop over all the element blocks
   for (int i=0; i<exio_helper->num_elem_blk; i++)
     {
       // Read the information for block i
       exio_helper->read_elem_in_block (i);
       const int subdomain_id = exio_helper->get_block_id(i);
       max_subdomain_id = std::max(max_subdomain_id, subdomain_id);
 
       // populate the map of names
       std::string subdomain_name = exio_helper->get_block_name(i);
       if (!subdomain_name.empty())
         mesh.subdomain_name(static_cast<subdomain_id_type>(subdomain_id)) = subdomain_name;
 
       // Set any relevant node/edge maps for this element
       const std::string type_str (exio_helper->get_elem_type());
       const auto & conv = exio_helper->get_conversion(type_str);
 
       // Loop over all the faces in this block
       int jmax = nelem_last_block+exio_helper->num_elem_this_blk;
       for (int j=nelem_last_block; j<jmax; j++)
         {
           auto uelem = Elem::build(conv.libmesh_elem_type());
 
           const int elem_num = j - nelem_last_block;
 
           // Make sure that Exodus's number of nodes per Elem matches
           // the number of Nodes for this type of Elem. We only check
           // this for the first Elem in each block, since these values
           // are the same for every Elem in the block.
           if (!elem_num)
             libmesh_error_msg_if(exio_helper->num_nodes_per_elem != static_cast<int>(uelem->n_nodes()),
                                  "Error: Exodus file says "
                                  << exio_helper->num_nodes_per_elem
                                  << " nodes per Elem, but Elem type "
                                  << Utility::enum_to_string(uelem->type())
                                  << " has " << uelem->n_nodes() << " nodes.");
 
           // Assign the current subdomain to this Elem
           uelem->subdomain_id() = static_cast<subdomain_id_type>(subdomain_id);
 
           // Use the elem_num_map to obtain the ID of this element in
           // the Exodus file.  Make sure we aren't reading garbage if
           // the file is corrupt.
           libmesh_error_msg_if(std::size_t(j) >= exio_helper->elem_num_map.size(),
                                "Error: Trying to read Exodus file with more elements than elem_num_map entries.\n");
           int exodus_id = exio_helper->elem_num_map[j];
 
           // Assign this element the same ID it had in the Exodus
           // file, but make it zero-based by subtracting 1.  Note:
           // some day we could use 1-based numbering in libmesh and
           // thus match the Exodus numbering exactly, but at the
           // moment libmesh is zero-based.
           uelem->set_id(exodus_id-1);
 
           // Record that we have seen an element of dimension uelem->dim()
           elems_of_dimension[uelem->dim()] = true;
 
           // Catch the Elem pointer that the Mesh throws back
           Elem * elem = mesh.add_elem(std::move(uelem));
 
           // If the Mesh assigned an ID different from what is in the
           // Exodus file, we should probably error.
           libmesh_error_msg_if(elem->id() != static_cast<unsigned>(exodus_id-1),
                                "Error!  Mesh assigned ID "
                                << elem->id()
                                << " which is different from the (zero-based) Exodus ID "
                                << exodus_id-1
                                << "!");
 
           // Assign extra integer IDs
           for (auto & id : extra_ids)
           {
             const Real v = elem_ids[id][elem->id()];
 
             if (v == Real(-1))
               {
                 elem->set_extra_integer(id, DofObject::invalid_id);
                 continue;
               }
 
             // Ignore FE_INVALID here even if we've enabled FPEs; a
             // thrown exception is preferred over an FPE signal.
             FPEDisabler disable_fpes;
             const long long iv = std::llround(v);
 
             // Check if the real number is outside of the range we can
             // convert exactly
 
             long long max_representation = 1;
             max_representation = (max_representation << std::min(std::numeric_limits<Real>::digits,
                                                                  std::numeric_limits<double>::digits));
             libmesh_error_msg_if(iv > max_representation,
                                  "Error! An element integer value higher than "
                                  << max_representation
                                  << " was found! Exodus uses real numbers for storing element "
                                  " integers, which can only represent integers from 0 to "
                                  << max_representation
                                  << ".");
 
             libmesh_error_msg_if(iv < 0,
                                  "Error! An element integer value less than -1"
                                  << " was found! Exodus uses real numbers for storing element "
                                  " integers, which can only represent integers from 0 to "
                                  << max_representation
                                  << ".");
 
 
             elem->set_extra_integer(id, cast_int<dof_id_type>(iv));
           }
 
           // Set all the nodes for this element
           //
           // If we don't have any Bezier extraction operators, this
           // is easy: we've already built all our nodes and just need
           // to link to them.
           if (exio_helper->bex_cv_conn.empty())
             {
               for (int k=0; k<exio_helper->num_nodes_per_elem; k++)
                 {
                   // global index
                   int gi = (elem_num)*exio_helper->num_nodes_per_elem + conv.get_node_map(k);
 
                   // The entries in 'connect' are actually (1-based)
                   // indices into the node_num_map, so to get the right
                   // node ID we:
                   // 1.) Subtract 1 from connect[gi]
                   // 2.) Pass it through node_num_map to get the corresponding Exodus ID
                   // 3.) Subtract 1 from that, since libmesh node numbering is "zero"-based,
                   //     even when the Exodus node numbering doesn't start with 1.
                   int libmesh_node_id = exio_helper->node_num_map[exio_helper->connect[gi] - 1] - 1;
 
                   // Set the node pointer in the Elem
                   elem->set_node(k, mesh.node_ptr(libmesh_node_id));
                 }
             }
           else // We have Bezier Extraction data
             {
               auto & constraint_rows = mesh.get_constraint_rows();
 
               const DynaIO::ElementDefinition & dyna_elem_defn =
                 DynaIO::find_elem_definition(elem->type(),
                                              elem->dim(),
                                              int(elem->default_order()));
 
               std::vector<std::vector<Real>>
                 my_constraint_mat(exio_helper->bex_num_elem_cvs);
               for (auto spline_node_index :
                    make_range(exio_helper->bex_num_elem_cvs))
                 {
                   my_constraint_mat[spline_node_index].resize(elem->n_nodes());
 
                   const auto & my_constraint_rows = exio_helper->bex_cv_conn[elem_num];
                   const unsigned long elem_coef_vec_index =
                     my_constraint_rows[spline_node_index] - 1; // Exodus isn't 0-based
                   const auto & my_vec = bex_constraint_vec(elem_coef_vec_index, *exio_helper);
                   for (auto elem_node_index :
                        make_range(elem->n_nodes()))
                     {
                       my_constraint_mat[spline_node_index][elem_node_index] =
                         my_vec[elem_node_index];
                     }
 
                 }
 
               // The tailing entries in each element's connectivity
               // vector are indices to Exodus constraint coefficient
               // rows.
 
               // Concatenating these rows gives a matrix with
               // all the constraints for the element nodes: each
               // column of that matrix is the constraint coefficients
               // for the node associated with that column (via the
               // Exodus numbering, not the libMesh numbering).
               const auto & my_constraint_rows = exio_helper->bex_cv_conn[elem_num];
 
               for (auto elem_node_index :
                    make_range(elem->n_nodes()))
                 {
                   // New finite element node data = dot product of
                   // constraint matrix columns with spline node data.
                   // Store each column's non-zero entries, along with
                   // the global spline node indices, as a key to
                   // identify shared finite element nodes.
                   std::vector<std::pair<dof_id_type, Real>> key;
 
                   for (auto spline_node_index :
                        make_range(exio_helper->bex_num_elem_cvs))
                     {
                       // Pick out a row of the element constraint matrix
                       const unsigned long elem_coef_vec_index =
                         my_constraint_rows[spline_node_index] - 1; // Exodus isn't 0-based
 
                       auto & coef_vec =
                         bex_constraint_vec(elem_coef_vec_index, *exio_helper);
 
                       // Get coef from this node's column intersect that row
                       const Real coef =
                         libmesh_vector_at(coef_vec, elem_node_index);
 
                       // Get the libMesh node corresponding to that row
                       const int gi = (elem_num)*exio_helper->bex_num_elem_cvs +
                         spline_node_index;
                       const dof_id_type libmesh_node_id =
                         exio_helper->node_num_map[exio_helper->connect[gi] - 1] - 1;
 
                       if (coef != 0) // Ignore irrelevant spline nodes
                         key.emplace_back(libmesh_node_id, coef);
                     }
 
                   // Have we already created this node?  Connect it.
                   if (const auto local_node_it = local_nodes.find(key);
                       local_node_it != local_nodes.end())
                     elem->set_node(dyna_elem_defn.nodes[elem_node_index], local_node_it->second);
                   // Have we not yet created this node?  Construct it,
                   // along with its weight and libMesh constraint row,
                   // then connect it.
                   else
                     {
                       Point p(0);
                       Real w = 0;
                       std::vector<std::pair<std::pair<const Elem *, unsigned int>, Real>> constraint_row;
 
                       for (auto [libmesh_spline_node_id, coef] : key)
                         {
                           const Node & spline_node = mesh.node_ref(libmesh_spline_node_id);
 
                           p.add_scaled(spline_node, coef);
                           const Real spline_w = weights_exist ?
                             spline_node.get_extra_datum<Real>(weight_index) : 1;
                           w += coef * spline_w;
 
                           const Elem * nodeelem =
                             libmesh_map_find(spline_nodeelem_ptrs, &spline_node);
                           constraint_row.emplace_back(std::make_pair(nodeelem, 0), coef);
                         }
 
                       Node *n = mesh.add_point(p, n_nodes++);
                       if (weights_exist)
                         n->set_extra_datum<Real>(weight_index, w);
 
                       // If we're building disconnected Bezier
                       // extraction elements then we don't want to
                       // find the new nodes to reuse later; each
                       // finite element node will connect to only one
                       // element.
                       if (!_disc_bex)
                         local_nodes[key] = n;
                       elem->set_node(dyna_elem_defn.nodes[elem_node_index], n);
 
                       constraint_rows[n] = constraint_row;
                     }
                 }
             }
         }
 
       // running sum of # of elements per block,
       // (should equal total number of elements in the end)
       nelem_last_block += exio_helper->num_elem_this_blk;
     }
 
   // Now we know enough to fix any spline NodeElem subdomains
   max_subdomain_id++;
   for (auto p : spline_nodeelem_ptrs)
     p.second->subdomain_id() = max_subdomain_id;
 
   // Read in edge blocks, storing information in the BoundaryInfo object.
   // Edge blocks are treated as BCs.
   exio_helper->read_edge_blocks(mesh);
 
   // Set the mesh dimension to the largest encountered for an element
   for (unsigned char i=0; i!=4; ++i)
     if (elems_of_dimension[i])
       mesh.set_mesh_dimension(i);
 
   // Read in sideset information -- this is useful for applying boundary conditions
   {
     // Get basic information about all sidesets
     exio_helper->read_sideset_info();
     int offset=0;
     for (int i=0; i<exio_helper->num_side_sets; i++)
       {
         // Compute new offset
         offset += (i > 0 ? exio_helper->num_sides_per_set[i-1] : 0);
         exio_helper->read_sideset (i, offset);
 
         std::string sideset_name = exio_helper->get_side_set_name(i);
         if (!sideset_name.empty())
           mesh.get_boundary_info().sideset_name
             (cast_int<boundary_id_type>(exio_helper->get_side_set_id(i)))
             = sideset_name;
       }
 
     for (auto e : index_range(exio_helper->elem_list))
       {
         // The numbers in the Exodus file sidesets should be thought
         // of as (1-based) indices into the elem_num_map array.  So,
         // to get the right element ID we have to:
         // 1.) Subtract 1 from elem_list[e] (to get a zero-based index)
         // 2.) Pass it through elem_num_map (to get the corresponding Exodus ID)
         // 3.) Subtract 1 from that, since libmesh is "zero"-based,
         //     even when the Exodus numbering doesn't start with 1.
         dof_id_type libmesh_elem_id =
           cast_int<dof_id_type>(exio_helper->elem_num_map[exio_helper->elem_list[e] - 1] - 1);
 
         // Set any relevant node/edge maps for this element
         Elem & elem = mesh.elem_ref(libmesh_elem_id);
 
         const auto & conv = exio_helper->get_conversion(elem.type());
 
         // Map the zero-based Exodus side numbering to the libmesh side numbering
         unsigned int raw_side_index = exio_helper->side_list[e]-1;
         std::size_t side_index_offset = conv.get_shellface_index_offset();
 
         if (raw_side_index < side_index_offset)
           {
             // We assume this is a "shell face"
             int mapped_shellface = raw_side_index;
 
             // Check for errors
             libmesh_error_msg_if(mapped_shellface < 0 || mapped_shellface >= 2,
                                  "Bad 0-based shellface id: "
                                  << mapped_shellface
                                  << " detected in Exodus file "
                                  << exio_helper->current_filename);
 
             // Add this (elem,shellface,id) triplet to the BoundaryInfo object.
             mesh.get_boundary_info().add_shellface (libmesh_elem_id,
                                                     cast_int<unsigned short>(mapped_shellface),
                                                     cast_int<boundary_id_type>(exio_helper->id_list[e]));
           }
         else
           {
             unsigned int side_index = static_cast<unsigned int>(raw_side_index - side_index_offset);
             int mapped_side = conv.get_side_map(side_index);
 
             // Check for errors
             libmesh_error_msg_if(mapped_side == ExodusII_IO_Helper::Conversion::invalid_id,
                                  "Invalid 1-based side id: "
                                  << side_index
                                  << " detected for "
                                  << Utility::enum_to_string(elem.type())
                                  << " in Exodus file "
                                  << exio_helper->current_filename);
 
             libmesh_error_msg_if(mapped_side < 0 ||
                                  cast_int<unsigned int>(mapped_side) >= elem.n_sides(),
                                  "Bad 0-based side id: "
                                  << mapped_side
                                  << " detected for "
                                  << Utility::enum_to_string(elem.type())
                                  << " in Exodus file "
                                  << exio_helper->current_filename);
 
             // Add this (elem,side,id) triplet to the BoundaryInfo object.
             mesh.get_boundary_info().add_side (libmesh_elem_id,
                                                cast_int<unsigned short>(mapped_side),
                                                cast_int<boundary_id_type>(exio_helper->id_list[e]));
           }
       } // end for (elem_list)
   } // end read sideset info
 
   // Read in elemset information and apply to Mesh elements if present
   {
     exio_helper->read_elemset_info();
 
     // Mimic behavior of sideset case where we store all the set
     // information in a single array with offsets.
     int offset=0;
     for (int i=0; i<exio_helper->num_elem_sets; i++)
       {
         // Compute new offset
         offset += (i > 0 ? exio_helper->num_elems_per_set[i-1] : 0);
         exio_helper->read_elemset (i, offset);
 
         // TODO: add support for elemset names
         // std::string elemset_name = exio_helper->get_elem_set_name(i);
         // if (!elemset_name.empty())
         //   mesh.get_boundary_info().elemset_name(cast_int<boundary_id_type>(exio_helper->get_elem_set_id(i))) = elemset_name;
       }
 
     // Debugging: print the concatenated list of elemset ids
     // libMesh::out << "Concatenated list of elemset Elem ids (Exodus numbering):" << std::endl;
     // for (const auto & id : exio_helper->elemset_list)
     //   libMesh::out << id << " ";
     // libMesh::out << std::endl;
 
     // Next we need to assign the elemset ids to the mesh using the
     // Elem's "extra_integers" support, if we have any.
     if (exio_helper->num_elem_all_elemsets)
       {
         // Build map from Elem -> {elemsets}. This is needed only
         // temporarily to determine a unique set of elemset codes.
         std::map<Elem *, MeshBase::elemset_type> elem_to_elemsets;
         for (auto e : index_range(exio_helper->elemset_list))
           {
             // Follow standard (see sideset case above) approach for
             // converting the ids stored in the elemset_list to
             // libmesh Elem ids.
             //
             // TODO: this should be moved to a helper function so we
             // don't duplicate the code.
             dof_id_type libmesh_elem_id =
               cast_int<dof_id_type>(exio_helper->elem_num_map[exio_helper->elemset_list[e] - 1] - 1);
 
             // Get a pointer to this Elem
             Elem * elem = mesh.elem_ptr(libmesh_elem_id);
 
             // Debugging:
             // libMesh::out << "Elem " << elem->id() << " is in elemset " << exio_helper->elemset_id_list[e] << std::endl;
 
             // Store elemset id in the map
             elem_to_elemsets[elem].insert(exio_helper->elemset_id_list[e]);
           }
 
         // Create a set of unique elemsets
         std::set<MeshBase::elemset_type> unique_elemsets;
         for (const auto & pr : elem_to_elemsets)
           unique_elemsets.insert(pr.second);
 
         // Debugging: print the unique elemsets
         // libMesh::out << "The set of unique elemsets which exist on the Mesh:" << std::endl;
         // for (const auto & s : unique_elemsets)
         //   {
         //     for (const auto & elemset_id : s)
         //       libMesh::out << elemset_id << " ";
         //     libMesh::out << std::endl;
         //   }
 
         // Enumerate the unique_elemsets and tell the mesh about them
         dof_id_type code = 0;
         for (const auto & s : unique_elemsets)
           mesh.add_elemset_code(code++, s);
 
         // Sanity check: make sure that MeshBase::n_elemsets() reports
         // the expected value after calling MeshBase::add_elemset_code()
         // one or more times.
         libmesh_assert_msg(exio_helper->num_elem_sets == cast_int<int>(mesh.n_elemsets()),
                            "Error: mesh.n_elemsets() is " << mesh.n_elemsets()
                            << ", but mesh should have " << exio_helper->num_elem_sets << " elemsets.");
 
         // Create storage for the extra integer on all Elems. Elems which
         // are not in any set will use the default value of DofObject::invalid_id
         unsigned int elemset_index =
           mesh.add_elem_integer("elemset_code",
                                 /*allocate_data=*/true);
 
         // Store the appropriate extra_integer value on all Elems that need it.
         for (const auto & [elem, s] : elem_to_elemsets)
           elem->set_extra_integer(elemset_index, mesh.get_elemset_code(s));
       }
   } // done reading elemset info
 
   // Read nodeset info
   {
     // This fills in the following fields of the helper for later use:
     // nodeset_ids
     // num_nodes_per_set
     // num_node_df_per_set
     // node_sets_node_index
     // node_sets_dist_index
     // node_sets_node_list
     // node_sets_dist_fact
     exio_helper->read_all_nodesets();
 
     for (int nodeset=0; nodeset<exio_helper->num_node_sets; nodeset++)
       {
         boundary_id_type nodeset_id =
           cast_int<boundary_id_type>(exio_helper->nodeset_ids[nodeset]);
 
         std::string nodeset_name = exio_helper->get_node_set_name(nodeset);
         if (!nodeset_name.empty())
           mesh.get_boundary_info().nodeset_name(nodeset_id) = nodeset_name;
 
         // Get starting index of node ids for current nodeset.
         unsigned int offset = exio_helper->node_sets_node_index[nodeset];
 
         for (int i=0; i<exio_helper->num_nodes_per_set[nodeset]; ++i)
           {
             int exodus_id = exio_helper->node_sets_node_list[i + offset];
 
             // It's possible for nodesets to have invalid ids in them
             // by accident.  Instead of possibly accessing past the
             // end of node_num_map, let's make sure we have that many
             // entries.
             libmesh_error_msg_if(static_cast<std::size_t>(exodus_id - 1) >= exio_helper->node_num_map.size(),
                                  "Invalid Exodus node id " << exodus_id
                                  << " found in nodeset " << nodeset_id);
 
             // As before, the entries in 'node_list' are 1-based
             // indices into the node_num_map array, so we have to map
             // them.  See comment above.
             int libmesh_node_id = exio_helper->node_num_map[exodus_id - 1] - 1;
             mesh.get_boundary_info().add_node(cast_int<dof_id_type>(libmesh_node_id),
                                               nodeset_id);
           }
       }
   }
 
 #if LIBMESH_DIM < 3
   libmesh_error_msg_if(mesh.mesh_dimension() > LIBMESH_DIM,
                        "Cannot open dimension "
                        << mesh.mesh_dimension()
                        << " mesh file when configured without "
                        << mesh.mesh_dimension()
                        << "D support.");
 #endif
 }

◆ read_elemental_variable()

void libMesh::ExodusII_IO::read_elemental_variable	(	std::string	elemental_var_name,
		unsigned int	timestep,
		std::map< unsigned int, Real > &	unique_id_to_value_map
	)

Given an elemental variable and a time step, returns a mapping from the elements (top parent) unique IDs to the value of the elemental variable at the corresponding time step index.

Note that this function MUST only be called before renumbering! This function is essentially a wrapper for read_elemental_var_values from the exodus helper (which is not accessible outside this class).

Parameters

elemental_var_name	Name of an elemental variable
timestep	The corresponding time step index
unique_id_to_value_map	The map to be filled

Definition at line 1309 of file exodusII_io.C.

References exio_helper, libMesh::MeshInput< MT >::mesh(), libMesh::Elem::top_parent(), and libMesh::DofObject::unique_id().

 {
   LOG_SCOPE("read_elemental_variable()", "ExodusII_IO");
 
   // Note that this function MUST be called before renumbering
   std::map<dof_id_type, Real> elem_var_value_map;
 
   exio_helper->read_elemental_var_values(elemental_var_name, timestep, elem_var_value_map);
   for (auto & pr : elem_var_value_map)
     {
       const Elem * elem = MeshInput<MeshBase>::mesh().query_elem_ptr(pr.first);
       unique_id_to_value_map.emplace(elem->top_parent()->unique_id(), pr.second);
     }
 }

◆ read_elemset_data()

void libMesh::ExodusII_IO::read_elemset_data	(	int	timestep,
		std::vector< std::string > &	var_names,
		std::vector< std::set< elemset_id_type >> &	elemset_ids_in,
		std::vector< std::map< std::pair< dof_id_type, elemset_id_type >, Real >> &	elemset_vals
	)

Read all the elemset data at a particular timestep.

Definition at line 2157 of file exodusII_io.C.

References exio_helper.

 {
   libmesh_error_msg_if(!exio_helper->opened_for_reading,
                        "ERROR, ExodusII file must be opened for reading "
                        "before calling ExodusII_IO::read_elemset_data()!");
 
   exio_helper->read_elemset_data(timestep, var_names, elemset_ids_in, elemset_vals);
 }

◆ read_global_variable()

void libMesh::ExodusII_IO::read_global_variable	(	std::vector< std::string >	global_var_names,
		unsigned int	timestep,
		std::vector< Real > &	global_values
	)

Given a vector of global variables and a time step, returns the values of the global variable at the corresponding time step index.

Parameters

global_var_names	Vector of names of global variables
timestep	The corresponding time step index
global_values	The vector to be filled

Definition at line 1326 of file exodusII_io.C.

References exio_helper, and libMesh::ExodusII_IO_Helper::GLOBAL.

Referenced by copy_scalar_solution().

 {
   LOG_SCOPE("read_global_variable()", "ExodusII_IO");
 
   std::size_t size = global_var_names.size();
   libmesh_error_msg_if(size == 0, "ERROR, empty list of global variables to read from the Exodus file.");
 
   // read the values for all global variables
   std::vector<Real> values_from_exodus;
   exio_helper->read_var_names(ExodusII_IO_Helper::GLOBAL);
   exio_helper->read_global_values(values_from_exodus, timestep);
   std::vector<std::string> global_var_names_exodus = exio_helper->global_var_names;
 
   if (values_from_exodus.size() == 0)
     return;   // This will happen in parallel on procs that are not 0
 
   global_values.clear();
   for (std::size_t i = 0; i != size; ++i)
     {
       // for each global variable in global_var_names, look the corresponding one in global_var_names_from_exodus
       // and fill global_values accordingly
       auto it = find(global_var_names_exodus.begin(), global_var_names_exodus.end(), global_var_names[i]);
       if (it != global_var_names_exodus.end())
         global_values.push_back(values_from_exodus[it - global_var_names_exodus.begin()]);
       else
         libmesh_error_msg("ERROR, Global variable " << global_var_names[i] << \
                           " not found in Exodus file.");
     }
 
 }

◆ read_header()

ExodusHeaderInfo libMesh::ExodusII_IO::read_header ( const std::string & name )

Read only the header information, instead of the entire mesh.

After the header is read, the file is closed and the MeshBase object remains unchanged. This capability is useful if you only need to know the mesh "metadata" and should be faster than reading in all the nodes and elems.

The header information currently includes: .) Title string .) Mesh spatial dimension .) Number of nodes .) Number of elements .) Number of element blocks .) Number of node sets .) Number of side sets .) Number of edge blocks/edges

Definition at line 923 of file exodusII_io.C.

References libMesh::ExodusHeaderInfo::broadcast(), libMesh::ParallelObject::comm(), exio_helper, libMesh::MeshInput< MT >::mesh(), libMesh::MeshInput< MeshBase >::mesh(), and libMesh::ParallelObject::processor_id().

Referenced by MeshInputTest::testExodusReadHeader().

 {
   // We will need the Communicator of the Mesh we were created with.
   MeshBase & mesh = MeshInput<MeshBase>::mesh();
 
   // Eventual return value
   ExodusHeaderInfo header_info;
 
   // File I/O is done on processor 0, then broadcast to other procs
   if (mesh.processor_id() == 0)
     {
       // Open the exodus file in EX_READ mode
       exio_helper->open(fname.c_str(), /*read_only=*/true);
 
       // Get header information from exodus file without updating the
       // Helper object's internal data structures.
       header_info = exio_helper->read_header();
 
       // Close the file, we are now done with it. The goal is to keep the
       // exio_helper object unchanged while calling this function,
       // although it can't quite be marked "const" because we do have to
       // actually open/close the file. This way, it should be possible to
       // use the same ExodusII_IO object to read the headers of multiple
       // different mesh files.
       exio_helper->close();
     }
 
   // Broadcast header_info to other procs before returning
   header_info.broadcast(mesh.comm());
 
   // Return the information we read back to the user.
   return header_info;
 }

◆ read_nodeset_data()

void libMesh::ExodusII_IO::read_nodeset_data	(	int	timestep,
		std::vector< std::string > &	var_names,
		std::vector< std::set< boundary_id_type >> &	node_boundary_ids,
		std::vector< std::map< BoundaryInfo::NodeBCTuple, Real >> &	bc_vals
	)

Read all the nodeset data at a particular timestep.

TODO: currently all the nodeset variables are read, but we might want to change this to only read the requested ones.

Definition at line 2127 of file exodusII_io.C.

References exio_helper.

 {
   libmesh_error_msg_if(!exio_helper->opened_for_reading,
                        "ERROR, ExodusII file must be opened for reading "
                        "before calling ExodusII_IO::read_nodeset_data()!");
 
   exio_helper->read_nodeset_data(timestep, var_names, node_boundary_ids, bc_vals);
 }

◆ read_sideset_data()

void libMesh::ExodusII_IO::read_sideset_data	(	int	timestep,
		std::vector< std::string > &	var_names,
		std::vector< std::set< boundary_id_type >> &	side_ids,
		std::vector< std::map< BoundaryInfo::BCTuple, Real >> &	bc_vals
	)

Similar to write_sideset_data(), this function is used to read the data at a particular timestep.

TODO: currently all the sideset variables are read, but we might want to change this to only read the requested ones.

Definition at line 2070 of file exodusII_io.C.

References exio_helper, libMesh::MeshInput< MeshBase >::mesh(), and libMesh::MeshOutput< MT >::mesh().

 {
   libmesh_error_msg_if(!exio_helper->opened_for_reading,
                        "ERROR, ExodusII file must be opened for reading "
                        "before calling ExodusII_IO::read_sideset_data()!");
 
   const MeshBase & mesh = MeshOutput<MeshBase>::mesh();
   exio_helper->read_sideset_data(mesh, timestep, var_names, side_ids, bc_vals);
 }

◆ set_coordinate_offset()

void libMesh::ExodusII_IO::set_coordinate_offset ( Point p )

Allows you to set a vector that is added to the coordinates of all of the nodes.

Effectively, this "moves" the mesh to a particular position.

Deprecated:: As requested by Roy in libmesh PR #90, this function was "deprecated on arrival". There is not really a suitable replacement for it in the works, however. The same effect could be achieved by calling MeshTools::Modification::translate() twice, but that approach seems inefficient in the case of very large problems with millions of nodes. That said, this should probably become a base class API so that it works for all the different IO subclasses.

Definition at line 990 of file exodusII_io.C.

References exio_helper.

 {
   libmesh_warning("This method may be deprecated in the future");
   exio_helper->set_coordinate_offset(p);
 }

◆ set_discontinuous_bex()

void libMesh::ExodusII_IO::set_discontinuous_bex ( bool disc_bex )

Set to true (false is the default) to generate independent nodes for every Bezier Extraction element in an input file containing them, causing those elements to be disconnected from each other.

(Their vertices will be multiple distinct nodes at overlapping points.) If this is false, Bezier Extraction elements will be connected where possible, thereby using fewer redundancies and less memory, but the input mesh will need to be conforming.

Definition at line 2413 of file exodusII_io.C.

References _disc_bex.

Referenced by MeshInputTest::testExodusFileMappings().

 {
   _disc_bex = disc_bex;
 }

◆ set_extra_integer_vars()

void libMesh::ExodusII_IO::set_extra_integer_vars ( const std::vector< std::string > & extra_integer_vars )

Set the elemental variables in the Exodus file to be read into extra element integers.

The names of these elemental variables will be used to name the extra element integers.

Definition at line 179 of file exodusII_io.C.

References _extra_integer_vars.

Referenced by ExtraIntegersTest::testBadExtraIntegersExodusReading(), and ExtraIntegersTest::testExtraIntegersExodusReading().

 {
   _extra_integer_vars = extra_integer_vars;
 }

◆ set_hdf5_writing()

void libMesh::ExodusII_IO::set_hdf5_writing ( bool write_hdf5 )

Set to true (the default) to write files in an HDF5-based file format (when HDF5 is available), or to false to write files in the old NetCDF3-based format.

If HDF5 is unavailable, this setting does nothing.

Definition at line 2407 of file exodusII_io.C.

References exio_helper.

 {
   exio_helper->set_hdf5_writing(write_hdf5);
 }

◆ set_n_partitions()

void libMesh::MeshInput< MeshBase >::set_n_partitions ( unsigned int n_parts )

inlineprotectedinherited

Sets the number of partitions in the mesh.

Typically this gets done by the partitioner, but some parallel file formats begin "pre-partitioned".

Definition at line 101 of file mesh_input.h.

References libMesh::MeshInput< MT >::mesh().

Referenced by libMesh::Nemesis_IO::read(), and libMesh::XdrIO::read_header().

101 { this->mesh().set_n_partitions() = n_parts; }

libMesh::MeshInput< MeshBase >::mesh

MeshBase & mesh()

Definition: mesh_input.h:178

libMesh::MeshBase::set_n_partitions

unsigned int & set_n_partitions()

Definition: mesh_base.h:1865

◆ set_output_variables()

void libMesh::ExodusII_IO::set_output_variables	(	const std::vector< std::string > &	output_variables,
		bool	allow_empty = `true`
	)

Sets the list of variable names to be included in the output.

This is optional. If this is never called then all variables will be present. If this is called and an empty vector is supplied no variables will be output. Setting the allow_empty = false will result in empty vectors supplied here to also be populated with all variables.

Definition at line 184 of file exodusII_io.C.

References _allow_empty_variables, and _output_variables.

Referenced by main().

 {
   _output_variables = output_variables;
   _allow_empty_variables = allow_empty;
 }

◆ skip_comment_lines()

void libMesh::MeshInput< MeshBase >::skip_comment_lines	(	std::istream &	in,
		const char	comment_start
	)

protectedinherited

Reads input from in, skipping all the lines that start with the character comment_start.

Definition at line 187 of file mesh_input.h.

Referenced by libMesh::TetGenIO::read(), and libMesh::UCDIO::read_implementation().

 {
   char c, line[256];
 
   while (in.get(c), c==comment_start)
     in.getline (line, 255);
 
   // put back first character of
   // first non-comment line
   in.putback (c);
 }

◆ use_mesh_dimension_instead_of_spatial_dimension()

void libMesh::ExodusII_IO::use_mesh_dimension_instead_of_spatial_dimension ( bool val )

In the general case, meshes containing 2D elements can be manifolds living in 3D space, thus by default we write all meshes with the Exodus dimension set to LIBMESH_DIM = mesh.spatial_dimension().

In certain cases, however, the user may know his 2D mesh actually lives in the z=0 plane, and therefore wants to write a truly 2D Exodus mesh. In such a case, he should call this function with val=true.

Definition at line 976 of file exodusII_io.C.

References exio_helper.

 {
   exio_helper->use_mesh_dimension_instead_of_spatial_dimension(val);
 }

◆ verbose()

void libMesh::ExodusII_IO::verbose ( bool set_verbosity )

Set the flag indicating if we should be verbose.

Definition at line 959 of file exodusII_io.C.

References _verbose, and exio_helper.

 {
   _verbose = set_verbosity;
 
   // Set the verbose flag in the helper object as well.
   exio_helper->verbose = _verbose;
 }

◆ write()

void libMesh::ExodusII_IO::write ( const std::string & fname )

overridevirtual

This method implements writing a mesh to a specified file.

Note that writes may be buffered for efficiency, and so may not reach disk until after the file has been closed, which happens when the ExodusII_IO object is destructed.

Implements libMesh::MeshOutput< MeshBase >.

Definition at line 2180 of file exodusII_io.C.

References _append, _verbose, exio_helper, libMesh::MeshBase::get_boundary_info(), libMesh::libmesh_assert(), libMesh::MeshInput< MeshBase >::mesh(), libMesh::MeshOutput< MT >::mesh(), and libMesh::BoundaryInfo::n_edge_conds().

Referenced by main(), libMesh::ErrorVector::plot_error(), ExodusTest< elem_type >::test_write(), MeshInputTest::testExodusFileMappings(), MeshInputTest::testExodusIGASidesets(), and libMesh::NameBasedIO::write().

 {
   LOG_SCOPE("write()", "ExodusII_IO");
 
   const MeshBase & mesh = MeshOutput<MeshBase>::mesh();
 
   // We may need to gather a DistributedMesh to output it, making that
   // const qualifier in our constructor a dirty lie
   // The "true" specifies that we only need the mesh serialized to processor 0
   MeshSerializer serialize
     (const_cast<MeshBase &>(mesh),
      !MeshOutput<MeshBase>::_is_parallel_format, true);
 
   libmesh_assert( !exio_helper->opened_for_writing );
 
   // If the user has set the append flag here, it doesn't really make
   // sense: the intent of this function is to write a Mesh with no
   // data, while "appending" is really intended to add data to an
   // existing file.  If we're verbose, print a message to this effect.
   if (_append && _verbose)
     libmesh_warning("Warning: Appending in ExodusII_IO::write() does not make sense.\n"
                     "Creating a new file instead!");
 
   exio_helper->create(fname);
   exio_helper->initialize(fname,mesh);
   exio_helper->write_nodal_coordinates(mesh);
   exio_helper->write_elements(mesh);
   exio_helper->write_sidesets(mesh);
   exio_helper->write_nodesets(mesh);
   exio_helper->write_elemsets(mesh);
 
   if ((mesh.get_boundary_info().n_edge_conds() > 0) && _verbose)
     libmesh_warning("Warning: Mesh contains edge boundary IDs, but these "
                     "are not supported by the ExodusII format.");
 }

◆ write_added_sides()

void libMesh::ExodusII_IO::write_added_sides ( bool val )

By default, we only write out the elements physically stored in the mesh.

If we have any SIDE_DISCONTINUOUS variables, however, we cannot easily output them on elements with sides that share vertices (and in 3D, edges). We can set this flag to instead create extra "side elements" on which to visualize such variables.

By default this flag is set to false.

Definition at line 1005 of file exodusII_io.C.

References exio_helper.

 {
   exio_helper->set_add_sides(val);
 }

◆ write_as_dimension()

void libMesh::ExodusII_IO::write_as_dimension ( unsigned dim )

Directly control the num_dim which is written to the Exodus file.

If non-zero, this value supersedes all other dimensions, including: 1.) MeshBase::spatial_dimension() 2.) MeshBase::mesh_dimension() 3.) Any value passed to use_mesh_dimension_instead_of_spatial_dimension() This is useful/necessary for working around a bug in Paraview which prevents the "Plot Over Line" filter from working on 1D meshes.

Definition at line 983 of file exodusII_io.C.

References dim, and exio_helper.

 {
   exio_helper->write_as_dimension(dim);
 }

◆ write_complex_magnitude()

void libMesh::ExodusII_IO::write_complex_magnitude ( bool val )

Set the flag indicating whether the complex modulus should be written when complex numbers are enabled.

By default this flag is set to true.

Definition at line 969 of file exodusII_io.C.

References _write_complex_abs.

 {
   _write_complex_abs = val;
 }

◆ 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 write_timestep_discontinuous().

 {
   LOG_SCOPE("write_discontinuous_equation_systems()", "MeshOutput");
 
   // We may need to gather and/or renumber a DistributedMesh to output
   // it, making that const qualifier in our constructor a dirty lie
   MT & my_mesh = const_cast<MT &>(*_obj);
 
   // If we're asked to write data that's associated with a different
   // mesh, output files full of garbage are the result.
   libmesh_assert_equal_to(&es.get_mesh(), _obj);
 
   // A non-renumbered mesh may not have a contiguous numbering, and
   // that needs to be fixed before we can build a solution vector.
   if (my_mesh.max_elem_id() != my_mesh.n_elem() ||
       my_mesh.max_node_id() != my_mesh.n_nodes())
     {
       // If we were allowed to renumber then we should have already
       // been properly renumbered...
       libmesh_assert(!my_mesh.allow_renumbering());
 
       libmesh_do_once(libMesh::out <<
                       "Warning:  This MeshOutput subclass only supports meshes which are contiguously renumbered!"
                       << std::endl;);
 
       my_mesh.allow_renumbering(true);
 
       my_mesh.renumber_nodes_and_elements();
 
       // Not sure what good going back to false will do here, the
       // renumbering horses have already left the barn...
       my_mesh.allow_renumbering(false);
     }
 
   MeshSerializer serialize(const_cast<MT &>(*_obj), !_is_parallel_format, _serial_only_needed_on_proc_0);
 
   // Build the list of variable names that will be written.
   std::vector<std::string> names;
   es.build_variable_names  (names, nullptr, system_names);
 
   if (!_is_parallel_format)
     {
       // Build the nodal solution values & get the variable
       // names from the EquationSystems object
       std::vector<Number> soln;
       es.build_discontinuous_solution_vector (soln, system_names,
                                               nullptr, false, /* defaults */
                                               this->get_add_sides());
 
       this->write_nodal_data_discontinuous (fname, soln, names);
     }
   else // _is_parallel_format
     {
       libmesh_not_implemented();
     }
 }

◆ write_discontinuous_exodusII()

void libMesh::ExodusII_IO::write_discontinuous_exodusII	(	const std::string &	name,
		const EquationSystems &	es,
		const std::set< std::string > *	system_names = `nullptr`
	)

Writes a exodusII file with discontinuous data.

Definition at line 193 of file exodusII_io.C.

References libMesh::EquationSystems::build_discontinuous_solution_vector(), libMesh::EquationSystems::build_variable_names(), get_add_sides(), and write_nodal_data_discontinuous().

Referenced by main().

 {
   std::vector<std::string> solution_names;
   std::vector<Number>      v;
 
   es.build_variable_names  (solution_names, nullptr, system_names);
   es.build_discontinuous_solution_vector (v, system_names,
                                           nullptr, false, /* defaults */
                                           this->get_add_sides());
   this->write_nodal_data_discontinuous(name, v, solution_names);
 }

◆ write_element_data()

void libMesh::ExodusII_IO::write_element_data ( const EquationSystems & es )

Write out element solution.

Definition at line 1359 of file exodusII_io.C.

References _output_variables, _timestep, _write_complex_abs, libMesh::EquationSystems::build_elemental_solution_vector(), libMesh::EquationSystems::build_variable_names(), libMesh::CONSTANT, exio_helper, libMesh::EquationSystems::get_vars_active_subdomains(), libMesh::MeshInput< MeshBase >::mesh(), libMesh::MeshOutput< MT >::mesh(), libMesh::MONOMIAL, libMesh::MONOMIAL_VEC, and libMesh::ParallelObject::processor_id().

Referenced by main(), and libMesh::ErrorVector::plot_error().

 {
   LOG_SCOPE("write_element_data()", "ExodusII_IO");
 
   // Be sure the file has been opened for writing!
   libmesh_error_msg_if(MeshOutput<MeshBase>::mesh().processor_id() == 0 && !exio_helper->opened_for_writing,
                        "ERROR, ExodusII file must be initialized before outputting element variables.");
 
   // This function currently only works on serialized meshes. We rely
   // on having a reference to a non-const MeshBase object from our
   // MeshInput parent class to construct a MeshSerializer object,
   // similar to what is done in ExodusII_IO::write().  Note that
   // calling ExodusII_IO::write_timestep() followed by
   // ExodusII_IO::write_element_data() when the underlying Mesh is a
   // DistributedMesh will result in an unnecessary additional
   // serialization/re-parallelization step.
   // The "true" specifies that we only need the mesh serialized to processor 0
   MeshSerializer serialize(MeshInput<MeshBase>::mesh(), !MeshOutput<MeshBase>::_is_parallel_format, true);
 
   // To be (possibly) filled with a filtered list of variable names to output.
   std::vector<std::string> names;
 
   // If _output_variables is populated, only output the monomials which are
   // also in the _output_variables vector.
   if (_output_variables.size() > 0)
     {
       // Create a list of CONSTANT MONOMIAL variable names
       std::vector<std::string> monomials;
       FEType type(CONSTANT, MONOMIAL);
       es.build_variable_names(monomials, &type);
 
       // Now concatenate a list of CONSTANT MONOMIAL_VEC variable names
       type = FEType(CONSTANT, MONOMIAL_VEC);
       es.build_variable_names(monomials, &type);
 
       // Filter that list against the _output_variables list.  Note: if names is still empty after
       // all this filtering, all the monomial variables will be gathered
       for (const auto & var : monomials)
         if (std::find(_output_variables.begin(), _output_variables.end(), var) != _output_variables.end())
           names.push_back(var);
     }
 
   // If we pass in a list of names to "build_elemental_solution_vector()"
   // it'll filter the variables coming back.
   std::vector<Number> soln;
   es.build_elemental_solution_vector(soln, names);
 
   // Also, store the list of subdomains on which each variable is active
   std::vector<std::set<subdomain_id_type>> vars_active_subdomains;
   es.get_vars_active_subdomains(names, vars_active_subdomains);
 
   if (soln.empty()) // If there is nothing to write just return
     return;
 
   // The data must ultimately be written block by block.  This means that this data
   // must be sorted appropriately.
   if (MeshOutput<MeshBase>::mesh().processor_id())
     return;
 
   const MeshBase & mesh = MeshOutput<MeshBase>::mesh();
 
 #ifdef LIBMESH_USE_COMPLEX_NUMBERS
 
   std::vector<std::string> complex_names =
     exio_helper->get_complex_names(names, _write_complex_abs);
 
   std::vector<std::set<subdomain_id_type>>
     complex_vars_active_subdomains =
     exio_helper->get_complex_vars_active_subdomains(vars_active_subdomains,
                                                     _write_complex_abs);
   exio_helper->initialize_element_variables(complex_names, complex_vars_active_subdomains);
 
   const std::vector<Real> complex_soln =
     complex_soln_components(soln, names.size(), _write_complex_abs);
 
   exio_helper->write_element_values(mesh, complex_soln, _timestep, complex_vars_active_subdomains);
 
 #else
   exio_helper->initialize_element_variables(names, vars_active_subdomains);
   exio_helper->write_element_values(mesh, soln, _timestep, vars_active_subdomains);
 #endif
 }

◆ write_element_data_from_discontinuous_nodal_data()

void libMesh::ExodusII_IO::write_element_data_from_discontinuous_nodal_data	(	const EquationSystems &	es,
		const std::set< std::string > *	system_names = `nullptr`,
		const std::string &	var_suffix = `"_elem_node_"`
	)

Similar to the function above, but instead of only handling (CONSTANT, MONOMIAL) data, writes out a general discontinuous solution field, e.g.

(FIRST, L2_LAGRANGE) or (SECOND, MONOMIAL) as a number of elemental fields equal to the number of vertices in each element. For example, if you have a (FIRST, L2_LAGRANGE) variable "u" defined on HEX8 elements, calling this function would by default write 8 elemental fields named u_elem_node_0, u_elem_node_1, u_elem_node_2, etc.

This may be useful if you have a viz tool which is capable of interpreting this element data as a discontinuous solution field. Note that (CONSTANT, MONOMIAL) data is still written as a single value per element, as it makes no sense to write n_vertices copies of the same value.

The 'var_suffix' parameter, which defaults to "_elem_node_", is used to generate the elemental variable names, and is inserted between the base variable name and the node id which the variable applies to, e.g. "u_elem_node_0", "u_elem_node_1", etc.

Definition at line 1446 of file exodusII_io.C.

References libMesh::EquationSystems::build_discontinuous_solution_vector(), libMesh::EquationSystems::build_variable_names(), libMesh::CONSTANT, distance(), libMesh::EquationSystems::get_vars_active_subdomains(), libMesh::index_range(), libMesh::MeshOutput< MT >::mesh(), mesh, libMesh::MONOMIAL, and libMesh::Quality::name().

Referenced by MeshInputTest::testExodusWriteElementDataFromDiscontinuousNodalData().

 {
   LOG_SCOPE("write_element_data_from_discontinuous_nodal_data()", "ExodusII_IO");
 
   // Be sure that some other function has already opened the file and prepared it
   // for writing. This is the same behavior as the write_element_data() function
   // which we are trying to mimic.
   libmesh_error_msg_if(MeshOutput<MeshBase>::mesh().processor_id() == 0 && !exio_helper->opened_for_writing,
                        "ERROR, ExodusII file must be initialized before outputting element variables.");
 
   // This function currently only works on serialized meshes.  The
   // "true" flag specifies that we only need the mesh serialized to
   // processor 0
   MeshSerializer serialize(MeshInput<MeshBase>::mesh(),
                            !MeshOutput<MeshBase>::_is_parallel_format,
                            true);
 
   // Note: in general we want to respect the contents of
   // _output_variables, only building a solution vector with values
   // from the requested variables. First build a list of all variable
   // names, then throw out ones that aren't in _output_variables, if
   // any.
   std::vector<std::string> var_names;
   es.build_variable_names (var_names, /*fetype=*/nullptr, system_names);
 
   // Get a subset of all variable names that are CONSTANT,
   // MONOMIALs. We treat those slightly differently since they can
   // truly only have a single value per Elem.
   //
   // Should the same apply here for CONSTANT MONOMIAL_VECs? [CW]
   // That is, get rid of 'const' on 'fe_type' and rerun:
   //    fe_type = FEType(CONSTANT, MONOMIAL_VEC);
   //    es.build_variable_names(monomial_var_names, &fe_type);
   // Then, es.find_variable_numbers() can be used without a type
   // (since we know for sure they're monomials) like:
   //    var_nums = es.find_variable_numbers(monomial_var_names)
   // for which the DOF indices for 'var_nums' have to be resolved
   // manually like in build_elemental_solution_vector()
   std::vector<std::string> monomial_var_names;
   const FEType fe_type(CONSTANT, MONOMIAL);
   es.build_variable_names(monomial_var_names, &fe_type);
 
   // Remove all names from var_names that are not in _output_variables.
   // Note: This approach avoids errors when the user provides invalid
   // variable names in _output_variables, as the code will not try to
   // write a variable that doesn't exist.
   if (!_output_variables.empty())
     {
       var_names.erase
         (std::remove_if
          (var_names.begin(),
           var_names.end(),
           [this](const std::string & name)
           {return !std::count(_output_variables.begin(),
                               _output_variables.end(),
                               name);}),
          var_names.end());
 
       // Also filter the monomial variable names.
       monomial_var_names.erase
         (std::remove_if
          (monomial_var_names.begin(),
           monomial_var_names.end(),
           [this](const std::string & name)
           {return !std::count(_output_variables.begin(),
                               _output_variables.end(),
                               name);}),
          monomial_var_names.end());
     }
 
   // Build a solution vector, limiting the results to the variables in
   // var_names and the Systems in system_names, and only computing values
   // at the vertices.
   std::vector<Number> v;
   es.build_discontinuous_solution_vector
     (v, system_names, &var_names, /*vertices_only=*/true);
 
   // Get active subdomains for each variable in var_names.
   std::vector<std::set<subdomain_id_type>> vars_active_subdomains;
   es.get_vars_active_subdomains(var_names, vars_active_subdomains);
 
   // Determine names of variables to write based on the number of
   // nodes/vertices the elements in different subdomains have.
   const MeshBase & mesh = MeshOutput<MeshBase>::mesh();
   std::map<subdomain_id_type, unsigned int> subdomain_id_to_vertices_per_elem;
   for (const auto & elem : mesh.active_element_ptr_range())
     {
       // Try to insert key/value pair into the map. If this returns
       // false, check the returned iterator's value to make sure it
       // matches. It shouldn't actually be possible for this to fail
       // (since if the Mesh was like this it would have already
       // failed) but it doesn't hurt to be on the safe side.
       auto pr2 = subdomain_id_to_vertices_per_elem.emplace
         (elem->subdomain_id(), elem->n_vertices());
       libmesh_error_msg_if(!pr2.second && pr2.first->second != elem->n_vertices(),
                            "Elem with different number of vertices found.");
     }
 
   // Determine "derived" variable names. These names are created by
   // starting with the base variable name and appending the user's
   // variable_suffix (default: "_elem_node_") followed by a node id.
   //
   // Not every derived variable will be active on every subdomain,
   // even if the original variable _is_ active.  Subdomains can have
   // different geometric element types (with differing numbers of
   // nodes), so some of the derived variable names will be inactive on
   // those subdomains.
   //
   // Since we would otherwise generate the same name once per
   // subdomain, we keep the list of names unique as we are creating
   // it. We can't use a std::set for this because we don't want the
   // variables names to be in a different order from the order
   // they were written in the call to: build_discontinuous_solution_vector()
   //
   // The list of derived variable names includes one for each vertex,
   // for higher-order elements we currently only write out vertex
   // values, but this could be changed in the future without too much
   // trouble.
   std::vector<std::string> derived_var_names;
 
   // Keep track of mapping from derived_name to (orig_name, node_id)
   // pair. We will use this later to determine whether a given
   // variable is active on a given subdomain.
   std::map<std::string, std::pair<std::string, unsigned int>>
     derived_name_to_orig_name_and_node_id;
 
   for (const auto & pr : subdomain_id_to_vertices_per_elem)
     {
       const subdomain_id_type sbd_id = pr.first;
       const unsigned int vertices_per_elem =
         subdomain_id_to_vertices_per_elem[sbd_id];
 
       std::ostringstream oss;
       for (unsigned int n=0; n<vertices_per_elem; ++n)
         for (const auto & orig_var_name : var_names)
           {
             oss.str("");
             oss.clear();
             oss << orig_var_name << var_suffix << n;
             std::string derived_name = oss.str();
 
             // Only add this var name if it's not already in the list.
             if (!std::count(derived_var_names.begin(), derived_var_names.end(), derived_name))
               {
                 derived_var_names.push_back(derived_name);
                 // Add entry for derived_name -> (orig_name, node_id) mapping.
                 derived_name_to_orig_name_and_node_id[derived_name] =
                   std::make_pair(orig_var_name, n);
               }
           }
     }
 
   // For each derived variable name, determine whether it is active
   // based on how many nodes/vertices the elements in a given subdomain have,
   // and whether they were active on the subdomain to begin with.
   std::vector<std::set<subdomain_id_type>>
     derived_vars_active_subdomains(derived_var_names.size());
 
   // A new data structure for keeping track of a list of variable names
   // that are in the discontinuous solution vector on each subdomain. Used
   // for indexing. Note: if a variable was inactive at the System level,
   // an entry for it will still be in the discontinuous solution vector,
   // but it will just have a value of zero. On the other hand, when we
   // create the derived variable names some of them are "inactive" on
   // different subdomains in the sense that they don't exist at all, i.e.
   // there is no zero padding for them. We need to be able to distinguish
   // between these two types in order to do the indexing into this vector
   // correctly.
   std::map<subdomain_id_type, std::vector<std::string>>
     subdomain_to_var_names;
 
   for (auto derived_var_id : index_range(derived_var_names))
     {
       const auto & derived_name = derived_var_names[derived_var_id];
       const auto & [orig_name, node_id] =
         libmesh_map_find (derived_name_to_orig_name_and_node_id,
                   derived_name);
 
       // For each subdomain, determine whether the current variable
       // should be active on that subdomain.
       for (const auto & pr : subdomain_id_to_vertices_per_elem)
         {
           // Convenience variables for the current subdomain and the
           // number of nodes elements in this subdomain have.
           subdomain_id_type sbd_id = pr.first;
           unsigned int vertices_per_elem_this_sbd =
             subdomain_id_to_vertices_per_elem[sbd_id];
 
           // Check whether variable orig_name was active on this
           // subdomain to begin with by looking in the
           // vars_active_subdomains container. We assume that the
           // location of orig_name in the var_names vector matches its
           // index in the vars_active_subdomains container.
           auto var_loc = std::find(var_names.begin(), var_names.end(), orig_name);
           libmesh_error_msg_if(var_loc == var_names.end(),
                                "Variable " << orig_name << " somehow not found in var_names array.");
           auto var_id = std::distance(var_names.begin(), var_loc);
 
           // The derived_var will only be active if this subdomain has
           // enough vertices for that to be the case.
           if (node_id < vertices_per_elem_this_sbd)
             {
               // Regardless of whether the original variable was not active on this subdomain,
               // the discontinuous solution vector will have zero padding for it, and
               // we will need to account for it. Therefore it should still be added to
               // the subdomain_to_var_names data structure!
               subdomain_to_var_names[sbd_id].push_back(derived_name);
 
               // If the original variable was not active on the
               // current subdomain, it should not be added to the
               // derived_vars_active_subdomains data structure, since
               // it will not be written to the Exodus file.
 
               // Determine if the original variable was active on the
               // current subdomain.
               bool orig_var_active =
                 (vars_active_subdomains[var_id].empty() ||
                  vars_active_subdomains[var_id].count(sbd_id));
 
               // And only if it was, add it to the
               // derived_vars_active_subdomains data structure.
               if (orig_var_active)
                 derived_vars_active_subdomains[derived_var_id].insert(sbd_id);
             }
         } // end loop over subdomain_id_to_vertices_per_elem
     } // end loop over derived_var_names
 
   // At this point we've built the "true" list of derived names, but
   // if there are any CONSTANT MONOMIALS in this list, we now want to
   // remove all but one copy of them from the derived_var_names list,
   // and rename them in (but not remove them from) the
   // subdomain_to_var_names list, and then update the
   // derived_vars_active_subdomains containers before finally calling
   // the Exodus helper functions.
   for (auto & derived_var_name : derived_var_names)
     {
       // Get the original name associated with this derived name.
       const auto & name_and_id =
         libmesh_map_find (derived_name_to_orig_name_and_node_id,
                   derived_var_name);
 
       // Convenience variables for the map entry's contents.
       const std::string & orig_name = name_and_id.first;
 
       // Was the original name a constant monomial?
       if (std::count(monomial_var_names.begin(),
                      monomial_var_names.end(),
                      orig_name))
         {
           // Rename this variable in the subdomain_to_var_names vectors.
           for (auto & pr : subdomain_to_var_names)
             {
               // Reference to ordered list of variable names on this subdomain.
               auto & name_vec = pr.second;
 
               auto name_vec_it =
                 std::find(name_vec.begin(),
                           name_vec.end(),
                           derived_var_name);
 
               if (name_vec_it != name_vec.end())
                 {
                   // Actually rename it back to the orig_name, dropping
                   // the "_elem_corner_" stuff.
                   *name_vec_it = orig_name;
                 }
             }
 
           // Finally, rename the variable in the derived_var_names vector itself.
           derived_var_name = orig_name;
         } // if (monomial)
     } // end loop over derived names
 
   // Now remove duplicate entries from derived_var_names after the first.
   // Also update the derived_vars_active_subdomains container in a consistent way.
   {
     std::vector<std::string> derived_var_names_edited;
     std::vector<std::set<subdomain_id_type>> derived_vars_active_subdomains_edited;
     std::vector<unsigned int> found_first(monomial_var_names.size());
 
     for (auto i : index_range(derived_var_names))
       {
         const auto & derived_var_name = derived_var_names[i];
         const auto & active_set = derived_vars_active_subdomains[i];
 
         // Determine whether we will keep this derived variable name in
         // the final container.
         bool keep = true;
         for (auto j : index_range(monomial_var_names))
           if (derived_var_name == monomial_var_names[j])
             {
               if (!found_first[j])
                 found_first[j] = 1;
 
               else
                 keep = false;
             }
 
         // We also don't keep variables that are not active on any subdomains.
         // Contrary to other uses of the var_active_subdomains container where
         // the empty set means "all" subdomains, here it really means "none".
         if (active_set.empty())
           keep = false;
 
         if (keep)
           {
             derived_var_names_edited.push_back(derived_var_name);
             derived_vars_active_subdomains_edited.push_back(active_set);
           }
       }
 
     // We built the filtered ranges, now swap them with the originals.
     derived_var_names.swap(derived_var_names_edited);
     derived_vars_active_subdomains.swap(derived_vars_active_subdomains_edited);
   }
 
 #ifdef LIBMESH_USE_COMPLEX_NUMBERS
   // Build complex variable names "r_foo", "i_foo", "a_foo" and the lists of
   // subdomains on which they are active.
   auto complex_var_names =
     exio_helper->get_complex_names(derived_var_names,
                                    _write_complex_abs);
   auto complex_vars_active_subdomains =
     exio_helper->get_complex_vars_active_subdomains(derived_vars_active_subdomains,
                                                     _write_complex_abs);
   auto complex_subdomain_to_var_names =
     exio_helper->get_complex_subdomain_to_var_names(subdomain_to_var_names,
                                                     _write_complex_abs);
 
   // Make expanded version of vector "v" in which each entry in the
   // original expands to an ("r_", "i_", "a_") triple.
   // "nco" is the number of complex outputs, which depends on whether
   // or not we are writing out the complex magnitudes.
   std::vector<Real> complex_v;
   int nco = _write_complex_abs ? 3 : 2;
   complex_v.reserve(nco * v.size());
   for (const auto & val : v)
     {
       complex_v.push_back(val.real());
       complex_v.push_back(val.imag());
       if (_write_complex_abs)
         complex_v.push_back(std::abs(val));
     }
 
   // Finally, initialize storage for the variables and write them to file.
   exio_helper->initialize_element_variables
     (complex_var_names, complex_vars_active_subdomains);
   exio_helper->write_element_values_element_major
     (mesh, complex_v, _timestep,
      complex_vars_active_subdomains,
      complex_var_names,
      complex_subdomain_to_var_names);
 #else
 
   // Call function which writes the derived variable names to the
   // Exodus file.
   exio_helper->initialize_element_variables(derived_var_names, derived_vars_active_subdomains);
 
   // ES::build_discontinuous_solution_vector() creates a vector with
   // an element-major ordering, so call Helper::write_element_values()
   // passing false for the last argument.
   exio_helper->write_element_values_element_major
     (mesh, v, _timestep,
      derived_vars_active_subdomains,
      derived_var_names,
      subdomain_to_var_names);
 #endif
 }

◆ write_elemset_data()

void libMesh::ExodusII_IO::write_elemset_data	(	int	timestep,
		const std::vector< std::string > &	var_names,
		const std::vector< std::set< elemset_id_type >> &	elemset_ids_in,
		const std::vector< std::map< std::pair< dof_id_type, elemset_id_type >, Real >> &	elemset_vals
	)

The Exodus format can also store values on elemsets.

The inputs to the function are: .) var_names[i] is the name of the ith elemset variable to be written to file. .) elemset_ids_in[i] is a set of elemset ids where var_names[i] is active. .) elemset_vals[i] is a map from (elem-id, elemset-id) pairs to the corresponding real-valued data.

Note: You must have already written the mesh by calling e.g. write() (and write_elemsets()) before calling this function, because it uses the ordering of the elemsets that have already been written to the Exodus file.

Definition at line 2141 of file exodusII_io.C.

References exio_helper.

 {
   libmesh_error_msg_if(!exio_helper->opened_for_writing,
                        "ERROR, ExodusII file must be opened for writing "
                        "before calling ExodusII_IO::write_elemset_data()!");
 
   exio_helper->write_elemset_data(timestep, var_names, elemset_ids_in, elemset_vals);
 }

◆ write_elemsets()

void libMesh::ExodusII_IO::write_elemsets ( )

Write elemsets stored on the Mesh to file.

Note: An elemset is a concept which is related to (but distinct from) both elem blocks and sidesets/nodesets. Elements in an elemset can be from multiple different blocks. In addition, one can define an elemset variable, which is like an elemental variable but exists only on the elements defined in the set.

Definition at line 2041 of file exodusII_io.C.

References exio_helper, libMesh::MeshInput< MeshBase >::mesh(), and libMesh::MeshOutput< MT >::mesh().

 {
   libmesh_error_msg_if(!exio_helper->opened_for_writing,
                        "ERROR, ExodusII file must be opened for writing "
                        "before calling ExodusII_IO::write_elemsets()!");
 
   const MeshBase & mesh = MeshOutput<MeshBase>::mesh();
   exio_helper->write_elemsets(mesh);
 }

◆ write_equation_systems()

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

overridevirtual

Writes out the solution for no specific time or timestep.

Parameters

fname	Name of the file to write to
es	EquationSystems object which contains the solution vector.
system_names	Optional list of systems to write solutions for.

Reimplemented from libMesh::MeshOutput< MeshBase >.

Definition at line 2033 of file exodusII_io.C.

References write_timestep().

Referenced by libMesh::RBConstruction::enrich_basis_from_rhs_terms(), main(), scale_mesh_and_plot(), MeshInputTest::testExodusWriteElementDataFromDiscontinuousNodalData(), PeriodicBCTest::testPeriodicBC(), WriteVecAndScalar::testWriteExodus(), transform_mesh_and_plot(), libMesh::TransientRBConstruction::truth_solve(), and libMesh::RBConstruction::truth_solve().

 {
   write_timestep(fname, es, 1, 0, system_names);
 }

◆ write_global_data()

void libMesh::ExodusII_IO::write_global_data	(	const std::vector< Number > &	soln,
		const std::vector< std::string > &	names
	)

Write out global variables.

Definition at line 1985 of file exodusII_io.C.

References _timestep, _write_complex_abs, exio_helper, and libMesh::ParallelObject::processor_id().

 {
   LOG_SCOPE("write_global_data()", "ExodusII_IO");
 
   if (MeshOutput<MeshBase>::mesh().processor_id())
     return;
 
   libmesh_error_msg_if(!exio_helper->opened_for_writing,
                        "ERROR, ExodusII file must be initialized before outputting global variables.");
 
 #ifdef LIBMESH_USE_COMPLEX_NUMBERS
 
   std::vector<std::string> complex_names =
     exio_helper->get_complex_names(names,
                                    _write_complex_abs);
 
   exio_helper->initialize_global_variables(complex_names);
 
   const std::vector<Real> complex_soln =
     complex_soln_components(soln, names.size(), _write_complex_abs);
 
   exio_helper->write_global_values(complex_soln, _timestep);
 
 #else
   exio_helper->initialize_global_variables(names);
   exio_helper->write_global_values(soln, _timestep);
 #endif
 }

◆ write_information_records()

void libMesh::ExodusII_IO::write_information_records ( const std::vector< std::string > & records )

Write out information records.

Definition at line 1972 of file exodusII_io.C.

References exio_helper, and libMesh::ParallelObject::processor_id().

 {
   if (MeshOutput<MeshBase>::mesh().processor_id())
     return;
 
   libmesh_error_msg_if(!exio_helper->opened_for_writing,
                        "ERROR, ExodusII file must be initialized before outputting information records.");
 
   exio_helper->write_information_records(records);
 }

◆ write_nodal_data() [1/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().

 {
   // This is the fallback implementation for parallel I/O formats that
   // do not yet implement proper writing in parallel, and instead rely
   // on the full solution vector being available on all processors.
   std::vector<Number> soln;
   parallel_soln.localize(soln);
   this->write_nodal_data(fname, soln, names);
 }

◆ write_nodal_data() [2/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().

 {
   std::vector<std::string> names;
   es.build_variable_names  (names, nullptr, system_names);
 
   std::unique_ptr<NumericVector<Number>> parallel_soln =
     es.build_parallel_solution_vector(system_names);
 
   this->write_nodal_data (fname, *parallel_soln, names);
 }

◆ write_nodal_data() [3/3]

void libMesh::ExodusII_IO::write_nodal_data	(	const std::string &	fname,
		const std::vector< Number > &	soln,
		const std::vector< std::string > &	names
	)

overridevirtual

Write out a nodal solution.

Reimplemented from libMesh::MeshOutput< MeshBase >.

Definition at line 1819 of file exodusII_io.C.

References _allow_empty_variables, _output_variables, _timestep, _write_complex_abs, exio_helper, libMesh::MeshTools::Generation::Private::idx(), std::imag(), libMesh::index_range(), libMesh::libmesh_ignore(), libMesh::MeshInput< MeshBase >::mesh(), libMesh::MeshOutput< MT >::mesh(), libMesh::MeshBase::n_nodes(), libMesh::ParallelObject::n_processors(), libMesh::ParallelObject::processor_id(), std::real(), libMesh::EquationSystems::redundant_added_side(), and write_nodal_data_common().

Referenced by libMesh::NameBasedIO::write_nodal_data().

 {
   LOG_SCOPE("write_nodal_data()", "ExodusII_IO");
 
   const MeshBase & mesh = MeshOutput<MeshBase>::mesh();
 
   int num_vars = cast_int<int>(names.size());
   dof_id_type num_nodes = mesh.n_nodes();
 
   // The names of the variables to be output
   std::vector<std::string> output_names;
 
   if (_allow_empty_variables || !_output_variables.empty())
     output_names = _output_variables;
   else
     output_names = names;
 
 #ifdef LIBMESH_USE_COMPLEX_NUMBERS
   std::vector<std::string> complex_names =
     exio_helper->get_complex_names(output_names,
                                    _write_complex_abs);
 
   // Call helper function for opening/initializing data, giving it the
   // complex variable names
   this->write_nodal_data_common(fname, complex_names, /*continuous=*/true);
 #else
   // Call helper function for opening/initializing data
   this->write_nodal_data_common(fname, output_names, /*continuous=*/true);
 #endif
 
   if (mesh.processor_id())
     return;
 
   // This will count the number of variables actually output
   for (int c=0; c<num_vars; c++)
     {
       std::stringstream name_to_find;
 
       std::vector<std::string>::iterator pos =
         std::find(output_names.begin(), output_names.end(), names[c]);
       if (pos == output_names.end())
         continue;
 
       unsigned int variable_name_position =
         cast_int<unsigned int>(pos - output_names.begin());
 
       // Set up temporary vectors to be passed to Exodus to write the
       // nodal values for a single variable at a time.
 #ifdef LIBMESH_USE_REAL_NUMBERS
       std::vector<Number> cur_soln;
 
       // num_nodes is either exactly how much space we will need for
       // each vector, or a safe upper bound for the amount of memory
       // we will require when there are gaps in the numbering.
       cur_soln.reserve(num_nodes);
 #else
       std::vector<Real> real_parts;
       std::vector<Real> imag_parts;
       std::vector<Real> magnitudes;
       real_parts.reserve(num_nodes);
       imag_parts.reserve(num_nodes);
       if (_write_complex_abs)
         magnitudes.reserve(num_nodes);
 #endif
 
       // There could be gaps in soln based on node numbering, but in
       // serial the empty numbers are left empty.
       // There could also be offsets in soln based on "fake" nodes
       // inserted on each processor (because NumericVector indices
       // have to be contiguous); the helper keeps track of those.
       // We now copy the proper solution values contiguously into
       // "cur_soln", removing the gaps.
       for (const auto & node : mesh.node_ptr_range())
         {
           const dof_id_type idx =
             (exio_helper->node_id_to_vec_id(node->id()))
             * num_vars + c;
 #ifdef LIBMESH_USE_REAL_NUMBERS
           cur_soln.push_back(soln[idx]);
 #else
           real_parts.push_back(soln[idx].real());
           imag_parts.push_back(soln[idx].imag());
           if (_write_complex_abs)
             magnitudes.push_back(std::abs(soln[idx]));
 #endif
         }
 
       // If we're adding extra sides, we need to add their data too.
       //
       // Because soln was created from a parallel NumericVector, its
       // numbering was contiguous on each processor; we need to use
       // the same offsets here, and we need to loop through elements
       // from earlier ranks first.
       if (exio_helper->get_add_sides())
         {
           std::vector<std::vector<const Elem *>>
             elems_by_pid(mesh.n_processors());
 
           for (const auto & elem : mesh.active_element_ptr_range())
             elems_by_pid[elem->processor_id()].push_back(elem);
 
           for (auto p : index_range(elems_by_pid))
             {
               dof_id_type global_idx =
                 exio_helper->added_node_offset_on(p) * num_vars + c;
               for (const Elem * elem : elems_by_pid[p])
                 {
                   for (auto s : elem->side_index_range())
                     {
                       if (EquationSystems::redundant_added_side(*elem,s))
                         continue;
 
                       const std::vector<unsigned int> side_nodes =
                         elem->nodes_on_side(s);
 
                       for (auto n : index_range(side_nodes))
                         {
                           libmesh_ignore(n);
                           libmesh_assert_less(global_idx, soln.size());
 #ifdef LIBMESH_USE_REAL_NUMBERS
                           cur_soln.push_back(soln[global_idx]);
 #else
                           real_parts.push_back(soln[global_idx].real());
                           imag_parts.push_back(soln[global_idx].imag());
                           if (_write_complex_abs)
                             magnitudes.push_back(std::abs(soln[global_idx]));
 #endif
                           global_idx += num_vars;
                         }
                     }
                 }
             }
         }
 
       // Finally, actually call the Exodus API to write to file.
 #ifdef LIBMESH_USE_REAL_NUMBERS
       exio_helper->write_nodal_values(variable_name_position+1, cur_soln, _timestep);
 #else
       int nco = _write_complex_abs ? 3 : 2;
       exio_helper->write_nodal_values(nco*variable_name_position+1, real_parts, _timestep);
       exio_helper->write_nodal_values(nco*variable_name_position+2, imag_parts, _timestep);
       if (_write_complex_abs)
         exio_helper->write_nodal_values(3*variable_name_position+3, magnitudes, _timestep);
 #endif
 
     }
 }

◆ write_nodal_data_common()

void libMesh::ExodusII_IO::write_nodal_data_common	(	std::string	fname,
		const std::vector< std::string > &	names,
		bool	continuous = `true`
	)

This function factors out a bunch of code which is common to the write_nodal_data() and write_nodal_data_discontinuous() functions.

Definition at line 2296 of file exodusII_io.C.

References _append, exio_helper, libMesh::MeshInput< MeshBase >::mesh(), libMesh::MeshOutput< MT >::mesh(), and libMesh::ParallelObject::processor_id().

Referenced by write_nodal_data(), and write_nodal_data_discontinuous().

 {
   const MeshBase & mesh = MeshOutput<MeshBase>::mesh();
 
   // This function can be called multiple times, we only want to open
   // the ExodusII file the first time it's called.
   if (!exio_helper->opened_for_writing)
     {
       // If we're appending, open() the file with read_only=false,
       // otherwise create() it and write the contents of the mesh to
       // it.
       if (_append)
         {
           // We do our writing only from proc 0, to avoid race
           // conditions with Exodus 8
           if (!MeshOutput<MeshBase>::mesh().processor_id())
             {
               exio_helper->open(fname.c_str(), /*read_only=*/false);
               // If we're appending, it's not valid to call exio_helper->initialize()
               // or exio_helper->initialize_nodal_variables(), but we do need to set up
               // certain aspects of the Helper object itself, such as the number of nodes
               // and elements.  We do that by reading the header...
               exio_helper->read_and_store_header_info();
 
               // ...and reading the block info
               exio_helper->read_block_info();
             }
             // Keep other processors aware of what we've done on root
           else
             {
               exio_helper->opened_for_writing = true;
               exio_helper->current_filename = fname;
             }
         }
       else
         {
           exio_helper->create(fname);
 
           // But some of our write calls are parallel-only, due to
           // calls to parallel-only getter functions.
           exio_helper->initialize(fname, mesh, !continuous);
 
           exio_helper->write_nodal_coordinates(mesh, !continuous);
           exio_helper->write_elements(mesh, !continuous);
 
           exio_helper->write_sidesets(mesh);
           exio_helper->write_nodesets(mesh);
           exio_helper->write_elemsets(mesh);
 
           exio_helper->initialize_nodal_variables(names);
         }
     }
   else
     {
       // We are already open for writing, so check that the filename
       // passed to this function matches the filename currently in use
       // by the helper.
       libmesh_error_msg_if(fname != exio_helper->current_filename,
                            "Error! This ExodusII_IO object is already associated with file: "
                            << exio_helper->current_filename
                            << ", cannot use it with requested file: "
                            << fname);
     }
 }

◆ write_nodal_data_discontinuous()

void libMesh::ExodusII_IO::write_nodal_data_discontinuous	(	const std::string &	fname,
		const std::vector< Number > &	soln,
		const std::vector< std::string > &	names
	)

overridevirtual

Write out a discontinuous nodal solution.

Reimplemented from libMesh::MeshOutput< MeshBase >.

Definition at line 2218 of file exodusII_io.C.

References _timestep, _write_complex_abs, exio_helper, get_add_sides(), libMesh::MeshInput< MeshBase >::mesh(), libMesh::MeshOutput< MT >::mesh(), libMesh::ParallelObject::processor_id(), and write_nodal_data_common().

Referenced by write_discontinuous_exodusII().

 {
   LOG_SCOPE("write_nodal_data_discontinuous()", "ExodusII_IO");
 
   const MeshBase & mesh = MeshOutput<MeshBase>::mesh();
 
 #ifdef LIBMESH_USE_COMPLEX_NUMBERS
 
   std::vector<std::string> complex_names =
     exio_helper->get_complex_names(names,
                                    _write_complex_abs);
 
   // Call helper function for opening/initializing data, giving it the
   // complex variable names
   this->write_nodal_data_common(fname, complex_names, /*continuous=*/false);
 #else
   // Call helper function for opening/initializing data
   this->write_nodal_data_common(fname, names, /*continuous=*/false);
 #endif
 
   if (mesh.processor_id())
     return;
 
   int num_vars = cast_int<int>(names.size());
   libmesh_assert_equal_to(soln.size() % num_vars, 0);
   int num_nodes = soln.size() / num_vars;
   libmesh_assert_equal_to(exio_helper->num_nodes, num_nodes);
 
 #ifndef NDEBUG
   if (!this->get_add_sides())
     {
       int num_real_nodes = 0;
       for (const auto & elem : mesh.active_element_ptr_range())
         num_real_nodes += elem->n_nodes();
       libmesh_assert_equal_to(num_real_nodes, num_nodes);
     }
 #endif
 
   for (int c=0; c<num_vars; c++)
     {
 #ifdef LIBMESH_USE_COMPLEX_NUMBERS
       std::vector<Real> real_parts(num_nodes);
       std::vector<Real> imag_parts(num_nodes);
       std::vector<Real> magnitudes;
       if (_write_complex_abs)
         magnitudes.resize(num_nodes);
 
       // The number of complex outputs depends on whether or not we are
       // writing out the absolute values.
       int nco = _write_complex_abs ? 3 : 2;
 
       for (int i=0; i<num_nodes; ++i)
         {
           real_parts[i] = soln[i*num_vars + c].real();
           imag_parts[i] = soln[i*num_vars + c].imag();
           if (_write_complex_abs)
             magnitudes[i] = std::abs(soln[i*num_vars + c]);
         }
       exio_helper->write_nodal_values(nco*c+1, real_parts, _timestep);
       exio_helper->write_nodal_values(nco*c+2, imag_parts, _timestep);
       if (_write_complex_abs)
         exio_helper->write_nodal_values(3*c+3, magnitudes, _timestep);
 #else
       // Copy out this variable's solution
       std::vector<Number> cur_soln(num_nodes);
 
       for (int i=0; i<num_nodes; i++)
         cur_soln[i] = soln[i*num_vars + c];
 
       exio_helper->write_nodal_values(c+1,cur_soln,_timestep);
 #endif
     }
 }

◆ write_nodeset_data()

void libMesh::ExodusII_IO::write_nodeset_data	(	int	timestep,
		const std::vector< std::string > &	var_names,
		const std::vector< std::set< boundary_id_type >> &	node_boundary_ids,
		const std::vector< std::map< BoundaryInfo::NodeBCTuple, Real >> &	bc_vals
	)

The Exodus format can also store values on nodesets.

This can be thought of as an alternative to defining a nodal variable field on lower-dimensional elements making up a part of the boundary. The inputs to the function are: .) var_names[i] is the name of the ith sideset variable to be written to file. .) node_boundary_ids[i] is a set of node_ids where var_names[i] is active. .) bc_vals[i] is a map from (node-id, boundary-id) NodeBCTuple objects to the corresponding real-valued data.

Note: You must have already written the mesh by calling e.g. write() before calling this function, because it uses the existing ordering of the Exodus nodesets.

Definition at line 2111 of file exodusII_io.C.

References exio_helper.

 {
   libmesh_error_msg_if(!exio_helper->opened_for_writing,
                        "ERROR, ExodusII file must be opened for writing "
                        "before calling ExodusII_IO::write_nodeset_data()!");
 
   exio_helper->write_nodeset_data(timestep, var_names, node_boundary_ids, bc_vals);
 }

◆ write_sideset_data()

void libMesh::ExodusII_IO::write_sideset_data	(	int	timestep,
		const std::vector< std::string > &	var_names,
		const std::vector< std::set< boundary_id_type >> &	side_ids,
		const std::vector< std::map< BoundaryInfo::BCTuple, Real >> &	bc_vals
	)

The Exodus format can also store values on sidesets.

This can be thought of as an alternative to defining an elemental variable field on lower-dimensional elements making up a part of the boundary. The inputs to the function are: .) var_names[i] is the name of the ith sideset variable to be written to file. .) side_ids[i] is a set of side_ids where var_names[i] is active. .) bc_vals[i] is a map from (elem,side,id) BCTuple objects to the corresponding real-valued data.

Note: You must have already written the mesh by calling e.g. write() before calling this function, because it uses the existing ordering of the Exodus sidesets.

Definition at line 2053 of file exodusII_io.C.

References exio_helper, libMesh::MeshInput< MeshBase >::mesh(), and libMesh::MeshOutput< MT >::mesh().

 {
   libmesh_error_msg_if(!exio_helper->opened_for_writing,
                        "ERROR, ExodusII file must be opened for writing "
                        "before calling ExodusII_IO::write_sideset_data()!");
 
   const MeshBase & mesh = MeshOutput<MeshBase>::mesh();
   exio_helper->write_sideset_data(mesh, timestep, var_names, side_ids, bc_vals);
 }

◆ write_timestep()

void libMesh::ExodusII_IO::write_timestep	(	const std::string &	fname,
		const EquationSystems &	es,
		const int	timestep,
		const Real	time,
		const std::set< std::string > *	system_names = `nullptr`
	)

Writes out the solution at a specific timestep.

Parameters

fname	Name of the file to write to
es	EquationSystems object which contains the solution vector.
timestep	The timestep to write out, should be 1 indexed.
time	The current simulation time.
system_names	Optional list of systems to write solutions for.

Definition at line 2017 of file exodusII_io.C.

References _timestep, exio_helper, libMesh::ParallelObject::processor_id(), and libMesh::MeshOutput< MT >::write_equation_systems().

Referenced by main(), libMesh::ClawSystem::solve_conservation_law(), write_equation_systems(), and write_output().

 {
   _timestep = timestep;
   MeshOutput<MeshBase>::write_equation_systems(fname,es,system_names);
 
   if (MeshOutput<MeshBase>::mesh().processor_id())
     return;
 
   exio_helper->write_timestep(timestep, time);
 }

◆ write_timestep_discontinuous()

void libMesh::ExodusII_IO::write_timestep_discontinuous	(	const std::string &	fname,
		const EquationSystems &	es,
		const int	timestep,
		const Real	time,
		const std::set< std::string > *	system_names = `nullptr`
	)

Writes a discontinuous solution at a specific timestep.

Parameters

fname	Name of the file to be written
es	EquationSystems object which contains the solution vector
timestep	The timestep to write out. (should be 1 indexed)
time	The current simulation time
system_names	Optional list of systems to write solutions for.

Definition at line 209 of file exodusII_io.C.

References _timestep, exio_helper, libMesh::ParallelObject::processor_id(), and libMesh::MeshOutput< MeshBase >::write_discontinuous_equation_systems().

 {
   _timestep = timestep;
   write_discontinuous_equation_systems (fname,es,system_names);
 
   if (MeshOutput<MeshBase>::mesh().processor_id())
     return;
 
   exio_helper->write_timestep(timestep, time);
 }

Member Data Documentation

◆ _allow_empty_variables

bool libMesh::ExodusII_IO::_allow_empty_variables

private

Flag which controls the behavior of _output_variables: .) If true, _output_variables is allowed to remain empty.

.) If false, if _output_variables is empty it will be populated with a complete list of all variables. .) By default, calling set_output_variables() sets this flag to true, but it provides an override.

Definition at line 660 of file exodusII_io.h.

Referenced by set_output_variables(), and write_nodal_data().

◆ _append

bool libMesh::ExodusII_IO::_append

private

Default false.

If true, files will be opened with EX_WRITE rather than created from scratch when writing.

Definition at line 636 of file exodusII_io.h.

Referenced by append(), write(), and write_nodal_data_common().

◆ _communicator

const Parallel::Communicator& libMesh::ParallelObject::_communicator

protectedinherited

Definition at line 120 of file parallel_object.h.

Referenced by libMesh::EquationSystems::build_parallel_elemental_solution_vector(), libMesh::EquationSystems::build_parallel_solution_vector(), libMesh::StaticCondensation::close(), libMesh::ParallelObject::comm(), libMesh::CondensedEigenSystem::initialize_condensed_matrices(), libMesh::ParallelObject::n_processors(), libMesh::ParallelObject::operator=(), libMesh::ParallelObject::processor_id(), libMesh::BoundaryInfo::regenerate_id_sets(), and libMesh::StaticCondensationDofMap::reinit().

◆ _disc_bex

bool libMesh::ExodusII_IO::_disc_bex

private

Set to true (false is the default) to generate independent nodes for every Bezier Extraction element.

Definition at line 677 of file exodusII_io.h.

Referenced by read(), and set_discontinuous_bex().

◆ _extra_integer_vars

std::vector<std::string> libMesh::ExodusII_IO::_extra_integer_vars

private

An optional list of variables in the EXODUS file that are to be used to set extra integers when loading the file into a mesh.

The variable names will be used to name the extra integers.

Definition at line 644 of file exodusII_io.h.

Referenced by read(), and set_extra_integer_vars().

◆ _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(), libMesh::PostscriptIO::write(), and libMesh::EnsightIO::write().

◆ _output_variables

std::vector<std::string> libMesh::ExodusII_IO::_output_variables

private

The names of the variables to be output.

If this is empty then all variables are output.

Definition at line 650 of file exodusII_io.h.

Referenced by set_output_variables(), write_element_data(), and write_nodal_data().

◆ _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.

◆ _timestep

int libMesh::ExodusII_IO::_timestep

private

Stores the current value of the timestep when calling ExodusII_IO::write_timestep().

Definition at line 625 of file exodusII_io.h.

Referenced by write_element_data(), write_global_data(), write_nodal_data(), write_nodal_data_discontinuous(), write_timestep(), and write_timestep_discontinuous().

◆ _verbose

bool libMesh::ExodusII_IO::_verbose

private

should we be verbose?

Definition at line 630 of file exodusII_io.h.

Referenced by verbose(), and write().

◆ _write_complex_abs

bool libMesh::ExodusII_IO::_write_complex_abs

private

By default, when complex numbers are enabled, for each variable we write out three values: the real part, "r_u" the imaginary part, "i_u", and the complex modulus, a_u := sqrt(r_u*r_u + i_u*i_u), which is also the value returned by std::abs(std::complex).

Since the modulus is not an independent quantity, we can set this flag to false and save some file space by not writing out.

Definition at line 671 of file exodusII_io.h.

Referenced by write_complex_magnitude(), write_element_data(), write_global_data(), write_nodal_data(), and write_nodal_data_discontinuous().

◆ elems_of_dimension

std::vector<bool> libMesh::MeshInput< MeshBase >::elems_of_dimension

protectedinherited

A vector of bools describing what dimension elements have been encountered when reading a mesh.

Definition at line 107 of file mesh_input.h.

Referenced by libMesh::GMVIO::_read_one_cell(), libMesh::UNVIO::elements_in(), libMesh::UNVIO::max_elem_dimension_seen(), libMesh::AbaqusIO::max_elem_dimension_seen(), libMesh::AbaqusIO::read(), libMesh::GMVIO::read(), libMesh::Nemesis_IO::read(), read(), libMesh::VTKIO::read(), libMesh::AbaqusIO::read_elements(), libMesh::UCDIO::read_implementation(), libMesh::UNVIO::read_implementation(), and libMesh::XdrIO::read_serialized_connectivity().

◆ exio_helper

std::unique_ptr<ExodusII_IO_Helper> libMesh::ExodusII_IO::exio_helper

private

Only attempt to instantiate an ExodusII helper class if the Exodus API is defined.

This class will have no functionality when LIBMESH_HAVE_EXODUS_API is not defined.

Definition at line 619 of file exodusII_io.h.

Referenced by copy_elemental_solution(), copy_nodal_solution(), copy_scalar_solution(), get_add_sides(), get_elem_num_map(), get_elem_var_names(), get_elemset_data_indices(), get_exio_helper(), get_global_var_names(), get_nodal_var_names(), get_node_num_map(), get_nodeset_data_indices(), get_num_time_steps(), get_sideset_data_indices(), get_time_steps(), read(), read_elemental_variable(), read_elemset_data(), read_global_variable(), read_header(), read_nodeset_data(), read_sideset_data(), set_coordinate_offset(), set_hdf5_writing(), use_mesh_dimension_instead_of_spatial_dimension(), verbose(), write(), write_added_sides(), write_as_dimension(), write_element_data(), write_elemset_data(), write_elemsets(), write_global_data(), write_information_records(), write_nodal_data(), write_nodal_data_common(), write_nodal_data_discontinuous(), write_nodeset_data(), write_sideset_data(), write_timestep(), write_timestep_discontinuous(), and ~ExodusII_IO().

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

include/mesh/exodusII_io.h
src/mesh/exodusII_io.C

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ExodusII_IO() [1/4]

◆ ExodusII_IO() [2/4]

◆ ExodusII_IO() [3/4]

◆ ExodusII_IO() [4/4]

◆ ~ExodusII_IO()

Member Function Documentation

◆ append()

◆ ascii_precision()

◆ comm()

◆ copy_elemental_solution()

◆ copy_nodal_solution()

◆ copy_scalar_solution()

◆ get_add_sides()

◆ get_elem_num_map()

◆ get_elem_var_names()

◆ get_elemset_data_indices()

◆ get_exio_helper()

◆ get_exodus_version()

◆ get_global_var_names()

◆ get_nodal_var_names()

◆ get_node_num_map()

◆ get_nodeset_data_indices()

◆ get_num_time_steps()

◆ get_sideset_data_indices()

◆ get_time_steps()

◆ is_parallel_format()

◆ mesh() [1/2]

◆ mesh() [2/2]

◆ n_processors()

◆ operator=() [1/2]

◆ operator=() [2/2]

◆ processor_id()

◆ read()

◆ read_elemental_variable()

◆ read_elemset_data()

◆ read_global_variable()

◆ read_header()

◆ read_nodeset_data()

◆ read_sideset_data()

◆ set_coordinate_offset()

◆ set_discontinuous_bex()

◆ set_extra_integer_vars()

◆ set_hdf5_writing()

◆ set_n_partitions()

◆ set_output_variables()

◆ skip_comment_lines()

◆ use_mesh_dimension_instead_of_spatial_dimension()

◆ verbose()

◆ write()

◆ write_added_sides()

◆ write_as_dimension()

◆ write_complex_magnitude()

◆ write_discontinuous_equation_systems()

◆ write_discontinuous_exodusII()

◆ write_element_data()

◆ write_element_data_from_discontinuous_nodal_data()

◆ write_elemset_data()

◆ write_elemsets()

◆ write_equation_systems()

◆ write_global_data()

◆ write_information_records()

◆ write_nodal_data() [1/3]

◆ write_nodal_data() [2/3]

◆ write_nodal_data() [3/3]

◆ write_nodal_data_common()

◆ write_nodal_data_discontinuous()

◆ write_nodeset_data()

◆ write_sideset_data()

◆ write_timestep()

◆ write_timestep_discontinuous()

Member Data Documentation

◆ _allow_empty_variables

◆ _append