This is the Nemesis_IO_Helper class. More...

#include <nemesis_io_helper.h>

Inheritance diagram for libMesh::Nemesis_IO_Helper:

Public Types
enum	ExodusVarType { NODAL =0, ELEMENTAL =1, GLOBAL =2, SIDESET =3 }
	Wraps calls to exII::ex_get_var_names() and exII::ex_get_var_param(). More...

Public Member Functions
	Nemesis_IO_Helper (const ParallelObject &parent, bool verbose=false, bool single_precision=false)
	Constructor. More...

virtual	~Nemesis_IO_Helper ()
	Destructor. More...

void	get_init_global ()
	Reading functions. More...

void	get_ss_param_global ()
	Fills: global_sideset_ids, num_global_side_counts, num_global_side_df_counts Call after: get_init_global() More...

void	get_ns_param_global ()

void	get_eb_info_global ()

void	get_init_info ()

void	get_loadbal_param ()

void	get_elem_map ()

void	get_node_map ()

void	get_cmap_params ()

void	get_node_cmap ()

void	get_elem_cmap ()

void	put_init_info (unsigned num_proc, unsigned num_proc_in_file, const char *ftype)
	Writing functions. More...

void	put_init_global (dof_id_type num_nodes_global, dof_id_type num_elems_global, unsigned num_elem_blks_global, unsigned num_node_sets_global, unsigned num_side_sets_global)
	Writes global information including: .) global number of nodes .) global number of elems .) global number of element blocks .) global number of node sets .) global number of side sets. More...

void	put_eb_info_global (std::vector< int > &global_elem_blk_ids, std::vector< int > &global_elem_blk_cnts)
	Writes global block information to the file .) global_elem_blk_ids - list of block IDs for all blocks present in the mesh .) global_elem_blk_cnts - number of elements in each block for the global mesh. More...

void	put_ns_param_global (std::vector< int > &global_nodeset_ids, std::vector< int > &num_global_node_counts, std::vector< int > &num_global_node_df_counts)
	This function writes information about global node sets. More...

void	put_ss_param_global (std::vector< int > &global_sideset_ids, std::vector< int > &num_global_side_counts, std::vector< int > &num_global_side_df_counts)
	This function writes information about global side sets. More...

void	put_loadbal_param (unsigned num_internal_nodes, unsigned num_border_nodes, unsigned num_external_nodes, unsigned num_internal_elems, unsigned num_border_elems, unsigned num_node_cmaps, unsigned num_elem_cmaps)
	Writes load balance parameters, some of which are described below: .) num_internal_nodes - nodes "wholly" owned by the current processor .) num_border_nodes - nodes local to a processor but residing in an element which also has nodes on other processors .) num_external_nodes - nodes that reside on other processors but whose element "partially" resides on the current processor – we assert this should be zero on reading! .) num_border_elems - elements local to this processor but whose nodes reside on other processors as well. More...

void	put_cmap_params (std::vector< int > &node_cmap_ids, std::vector< int > &node_cmap_node_cnts, std::vector< int > &elem_cmap_ids, std::vector< int > &elem_cmap_elem_cnts)
	Outputs initial information for communication maps. More...

void	put_node_cmap (std::vector< std::vector< int >> &node_cmap_node_ids, std::vector< std::vector< int >> &node_cmap_proc_ids)
	Outputs all of the nodal communication maps for this processor. More...

void	put_node_map (std::vector< int > &node_mapi, std::vector< int > &node_mapb, std::vector< int > &node_mape)
	Outputs IDs of internal, border, and external nodes. More...

void	put_elem_cmap (std::vector< std::vector< int >> &elem_cmap_elem_ids, std::vector< std::vector< int >> &elem_cmap_side_ids, std::vector< std::vector< int >> &elem_cmap_proc_ids)
	Writes information about elemental communication map. More...

void	put_elem_map (std::vector< int > &elem_mapi, std::vector< int > &elem_mapb)
	Outputs IDs of internal and border elements. More...

void	put_n_coord (unsigned start_node_num, unsigned num_nodes, std::vector< Real > &x_coor, std::vector< Real > &y_coor, std::vector< Real > &z_coor)
	Writes the specified number of coordinate values starting at the specified index. More...

virtual void	write_nodal_coordinates (const MeshBase &mesh, bool use_discontinuous=false) override
	This function is specialized from ExodusII_IO_Helper to write only the nodal coordinates stored on the local piece of the Mesh. More...

virtual void	write_elements (const MeshBase &mesh, bool use_discontinuous=false) override
	This function is specialized to write the connectivity. More...

virtual void	write_sidesets (const MeshBase &mesh) override
	Writes the sidesets for this processor. More...

virtual void	write_nodesets (const MeshBase &mesh) override
	Writes the nodesets for this processor. More...

virtual void	create (std::string filename) override
	This function is specialized from ExodusII_IO_Helper to create the nodal coordinates stored on the local piece of the Mesh. More...

virtual void	initialize (std::string title, const MeshBase &mesh, bool use_discontinuous=false) override
	Specialization of the initialize function from ExodusII_IO_Helper that also writes global initial data to file. More...

void	write_nodal_solution (const NumericVector< Number > &parallel_soln, const std::vector< std::string > &names, int timestep, const std::vector< std::string > &output_names)
	Takes a parallel solution vector containing the node-major solution vector for all variables and outputs it to the files. More...

void	write_nodal_solution (const EquationSystems &es, const std::vector< std::pair< unsigned int, unsigned int >> &var_nums, int timestep, const std::vector< std::string > &output_names)
	Outputs EquationSystems current_local_solution nodal values. More...

void	write_nodal_solution (const std::vector< Number > &values, const std::vector< std::string > &names, int timestep)
	Takes a solution vector containing the solution for all variables and outputs it to the files. More...

virtual void	initialize_element_variables (std::vector< std::string > names, const std::vector< std::set< subdomain_id_type >> &vars_active_subdomains) override
	Override the Exodus Helper's implementation of this function so that it works correctly in parallel. More...

void	write_element_values (const MeshBase &mesh, const EquationSystems &es, const std::vector< std::pair< unsigned int, unsigned int >> &var_nums, int timestep, const std::vector< std::set< subdomain_id_type >> &vars_active_subdomains)
	Writes the vector of elemental variable values, one variable and one subdomain at a time. More...

std::string	construct_nemesis_filename (const std::string &base_filename)
	Given base_filename, foo.e, constructs the Nemesis filename foo.e.X.Y, where X=n. More...

const char *	get_elem_type () const

void	open (const char *filename, bool read_only)
	Opens an `ExodusII` mesh file named `filename`. More...

void	read_header ()
	Reads an `ExodusII` mesh file header. More...

void	read_qa_records ()
	Reads the QA records from an ExodusII file. More...

void	print_header ()
	Prints the `ExodusII` mesh file header, which includes the mesh title, the number of nodes, number of elements, mesh dimension, number of sidesets, and number of nodesets. More...

void	read_nodes ()
	Reads the nodal data (x,y,z coordinates) from the `ExodusII` mesh file. More...

void	read_node_num_map ()
	Reads the optional `node_num_map` from the `ExodusII` mesh file. More...

void	print_nodes (std::ostream &out=libMesh::out)
	Prints the nodal information, by default to `libMesh::out`. More...

void	read_block_info ()
	Reads information for all of the blocks in the `ExodusII` mesh file. More...

int	get_block_id (int index)
	Get the block number for the given block index. More...

std::string	get_block_name (int index)
	Get the block name for the given block index if supplied in the mesh file. More...

int	get_side_set_id (int index)
	Get the side set id for the given side set index. More...

std::string	get_side_set_name (int index)
	Get the side set name for the given side set index if supplied in the mesh file. More...

int	get_node_set_id (int index)
	Get the node set id for the given node set index. More...

std::string	get_node_set_name (int index)
	Get the node set name for the given node set index if supplied in the mesh file. More...

void	read_elem_in_block (int block)
	Reads all of the element connectivity for block `block` in the `ExodusII` mesh file. More...

void	read_edge_blocks (MeshBase &mesh)
	Read in edge blocks, storing information in the BoundaryInfo object. More...

void	read_elem_num_map ()
	Reads the optional `node_num_map` from the `ExodusII` mesh file. More...

void	read_sideset_info ()
	Reads information about all of the sidesets in the `ExodusII` mesh file. More...

void	read_nodeset_info ()
	Reads information about all of the nodesets in the `ExodusII` mesh file. More...

void	read_sideset (int id, int offset)
	Reads information about sideset `id` and inserts it into the global sideset array at the position `offset`. More...

void	read_nodeset (int id)
	Reads information about nodeset `id` and inserts it into the global nodeset array at the position `offset`. More...

void	read_all_nodesets ()
	New API that reads all nodesets simultaneously. More...

void	close ()
	Closes the `ExodusII` mesh file. More...

int	inquire (int req_info, std::string error_msg="")

void	read_time_steps ()
	Reads and stores the timesteps in the 'time_steps' array. More...

void	read_num_time_steps ()
	Reads the number of timesteps currently stored in the Exodus file and stores it in the num_time_steps variable. More...

void	read_nodal_var_values (std::string nodal_var_name, int time_step)
	Reads the nodal values for the variable 'nodal_var_name' at the specified time into the 'nodal_var_values' array. More...

void	read_elemental_var_values (std::string elemental_var_name, int time_step, std::map< dof_id_type, Real > &elem_var_value_map)
	Reads elemental values for the variable 'elemental_var_name' at the specified timestep into the 'elem_var_value_map' which is passed in. More...

void	initialize_nodal_variables (std::vector< std::string > names)
	Sets up the nodal variables. More...

void	initialize_global_variables (std::vector< std::string > names)
	Sets up the global variables. More...

void	write_timestep (int timestep, Real time)
	Writes the time for the timestep. More...

void	write_sideset_data (const MeshBase &mesh, 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)
	Write sideset data for the requested timestep. More...

void	read_sideset_data (const MeshBase &mesh, 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)
	Read sideset variables, if any, into the provided data structures. More...

void	write_element_values (const MeshBase &mesh, const std::vector< Real > &values, int timestep, const std::vector< std::set< subdomain_id_type >> &vars_active_subdomains)
	Writes the vector of values to the element variables. More...

void	write_element_values_element_major (const MeshBase &mesh, const std::vector< Real > &values, int timestep, const std::vector< std::set< subdomain_id_type >> &vars_active_subdomains, const std::vector< std::string > &derived_var_names, const std::map< subdomain_id_type, std::vector< std::string >> &subdomain_to_var_names)
	Same as the function above, but assume the input 'values' vector is in element-major order, i.e. More...

void	write_nodal_values (int var_id, const std::vector< Real > &values, int timestep)
	Writes the vector of values to a nodal variable. More...

void	write_information_records (const std::vector< std::string > &records)
	Writes the vector of information records. More...

void	write_global_values (const std::vector< Real > &values, int timestep)
	Writes the vector of global variables. More...

void	read_global_values (std::vector< Real > &values, int timestep)
	Reads the vector of global variables. More...

void	use_mesh_dimension_instead_of_spatial_dimension (bool val)
	Sets the underlying value of the boolean flag _use_mesh_dimension_instead_of_spatial_dimension. More...

void	write_as_dimension (unsigned dim)
	Sets the value of _write_as_dimension. 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...

std::vector< std::string >	get_complex_names (const std::vector< std::string > &names) const

std::vector< std::set< subdomain_id_type > >	get_complex_vars_active_subdomains (const std::vector< std::set< subdomain_id_type >> &vars_active_subdomains) const
	returns a "tripled" copy of `vars_active_subdomains`, which is necessary in the complex-valued case. More...

std::map< subdomain_id_type, std::vector< std::string > >	get_complex_subdomain_to_var_names (const std::map< subdomain_id_type, std::vector< std::string >> &subdomain_to_var_names) const
	Takes a map from subdomain id -> vector of active variable names as input and returns a corresponding map where the original variable names have been replaced by their complex counterparts. More...

void	message (const std::string &msg)
	Prints the message defined in `msg`. More...

void	message (const std::string &msg, int i)
	Prints the message defined in `msg`, and appends the number `i` to the end of the message. More...

void	read_var_names (ExodusVarType type)

const ExodusII_IO_Helper::Conversion &	get_conversion (const ElemType type) const

const ExodusII_IO_Helper::Conversion &	get_conversion (std::string type_str) const

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Public Attributes
int	nemesis_err_flag
	Member data. More...

int	num_nodes_global
	Global initial information. More...

int	num_elems_global

int	num_elem_blks_global

int	num_node_sets_global

int	num_side_sets_global

int	num_proc
	The number of processors for which the NEMESIS I file was created. More...

int	num_proc_in_file
	The number of processors for which the NEMESIS I file stores information. More...

char	ftype
	The type of file to be written. More...

std::vector< int >	global_sideset_ids
	Containers for reading global sideset (boundary conditions) information. More...

std::vector< int >	num_global_side_counts

std::vector< int >	num_global_side_df_counts

std::vector< int >	global_nodeset_ids
	Containers for reading global nodeset information. More...

std::vector< int >	num_global_node_counts

std::vector< int >	num_global_node_df_counts

std::vector< int >	global_elem_blk_ids
	Read the global element block IDs and counts. More...

std::vector< int >	global_elem_blk_cnts

std::set< int >	nodes_attached_to_local_elems
	libMesh numbered node ids attached to local elems. More...

std::map< subdomain_id_type, std::vector< dof_id_type > >	subdomain_map
	Map of subdomains to element numbers. More...

std::map< int, std::vector< int > >	block_id_to_elem_connectivity
	This is the block connectivity, i.e. More...

int	num_internal_nodes
	To be used with the Nemesis::ne_get_loadbal_param() routine. More...

int	num_border_nodes
	The number of FEM nodes local to a processor but residing in an element which also has FEM nodes on other processors. More...

int	num_external_nodes
	The number of FEM nodes that reside on another processor but whose element partially resides on the current processor. More...

int	num_internal_elems
	The number of internal FEM elements. More...

int	num_border_elems
	The number of border FEM elements. More...

int	num_node_cmaps
	The number of nodal communication maps for this processor. More...

int	num_elem_cmaps
	The number of elemental communication maps for this processor. More...

std::vector< int >	elem_mapi
	Vector which stores internal element IDs. More...

std::vector< int >	elem_mapb
	Vector which stores border element IDs. More...

std::vector< int >	node_mapi
	Vector which stores internal node IDs. More...

std::vector< int >	node_mapb
	Vector which stores border node IDs. More...

std::vector< int >	node_mape
	Vector which stores external node IDs. More...

std::vector< int >	node_cmap_ids
	Vectors for storing the communication map parameters. More...

std::vector< int >	node_cmap_node_cnts

std::vector< int >	elem_cmap_ids

std::vector< int >	elem_cmap_elem_cnts

std::vector< std::vector< int > >	node_cmap_node_ids
	2 vectors of vectors for storing the node communication IDs for this processor. More...

std::vector< std::vector< int > >	node_cmap_proc_ids

std::vector< std::vector< int > >	elem_cmap_elem_ids
	3 vectors of vectors for storing element communication IDs for this processor. More...

std::vector< std::vector< int > >	elem_cmap_side_ids

std::vector< std::vector< int > >	elem_cmap_proc_ids

int	ex_id

int	ex_err

int	num_dim

int	num_global_vars

int	num_sideset_vars

int	num_nodes

int	num_elem

int	num_elem_blk

int	num_edge

int	num_edge_blk

int	num_node_sets

int	num_side_sets

int	num_elem_this_blk

int	num_nodes_per_elem

int	num_attr

int	num_elem_all_sidesets

std::vector< int >	block_ids

std::vector< int >	edge_block_ids

std::vector< int >	connect

std::vector< int >	ss_ids

std::vector< int >	nodeset_ids

std::vector< int >	num_sides_per_set

std::vector< int >	num_nodes_per_set

std::vector< int >	num_df_per_set

std::vector< int >	num_node_df_per_set

std::vector< int >	node_sets_node_index

std::vector< int >	node_sets_dist_index

std::vector< int >	node_sets_node_list

std::vector< Real >	node_sets_dist_fact

std::vector< int >	elem_list

std::vector< int >	side_list

std::vector< int >	node_list

std::vector< int >	id_list

std::vector< int >	node_num_map

std::vector< int >	elem_num_map

std::vector< Real >	x

std::vector< Real >	y

std::vector< Real >	z

std::vector< char >	title

std::vector< char >	elem_type

std::map< int, int >	libmesh_elem_num_to_exodus

std::vector< int >	exodus_elem_num_to_libmesh

std::map< int, int >	libmesh_node_num_to_exodus

std::vector< int >	exodus_node_num_to_libmesh

int	num_time_steps

std::vector< Real >	time_steps

int	num_nodal_vars

std::vector< std::string >	nodal_var_names

std::vector< Real >	nodal_var_values

int	num_elem_vars

std::vector< std::string >	elem_var_names

std::vector< Real >	elem_var_values

std::vector< std::string >	global_var_names

std::vector< std::string >	sideset_var_names

std::map< int, std::string >	id_to_block_names

std::map< int, std::string >	id_to_edge_block_names

std::map< int, std::string >	id_to_ss_names

std::map< int, std::string >	id_to_ns_names

bool	verbose

bool	opened_for_writing

bool	opened_for_reading

std::string	current_filename

Protected Member Functions
void	check_existing_vars (ExodusVarType type, std::vector< std::string > &names, std::vector< std::string > &names_from_file)
	When appending: during initialization, check that variable names in the file match those you attempt to initialize with. More...

void	write_var_names (ExodusVarType type, const std::vector< std::string > &names)
	Wraps calls to exII::ex_put_var_names() and exII::ex_put_var_param(). More...

Protected Attributes
bool	_run_only_on_proc0

bool	_elem_vars_initialized

bool	_global_vars_initialized

bool	_nodal_vars_initialized

bool	_use_mesh_dimension_instead_of_spatial_dimension

unsigned	_write_as_dimension

Point	_coordinate_offset

bool	_single_precision

const Parallel::Communicator &	_communicator

Private Types
typedef std::map< unsigned, std::set< unsigned > >::iterator	proc_nodes_touched_iterator
	Typedef for an iterator into the data structure above. More...

typedef std::map< unsigned, std::set< std::pair< unsigned, unsigned > > >::iterator	proc_border_elem_sets_iterator
	Typedef for an iterator into the data structure above. More...

Private Member Functions
void	compute_num_global_elem_blocks (const MeshBase &pmesh)
	This function uses global communication routines to determine the number of element blocks across the entire mesh. More...

void	compute_num_global_nodesets (const MeshBase &pmesh)
	This function uses global communication routines to determine the number of nodesets across the entire mesh. More...

void	compute_num_global_sidesets (const MeshBase &pmesh)
	This function uses global communication routines to determine the number of sidesets across the entire mesh. More...

void	build_element_and_node_maps (const MeshBase &pmesh)
	This function builds the libmesh -> exodus and exodus -> libmesh node and element maps. More...

void	compute_border_node_ids (const MeshBase &pmesh)
	This function constructs the set of border node IDs present on the current mesh. More...

void	compute_internal_and_border_elems_and_internal_nodes (const MeshBase &pmesh)
	This function constructs the set of border and internal element IDs and internal node IDs present on the current mesh. More...

void	compute_communication_map_parameters ()
	This function determines the communication map parameters which will eventually be written to file. More...

void	compute_node_communication_maps ()
	Compute the node communication maps (really just pack vectors) in preparation for writing them to file. More...

void	compute_node_maps ()
	Compute the node maps (really just pack vectors) which map the nodes to internal, border, and external nodes in the file. More...

void	compute_elem_communication_maps ()
	This function computes element communication maps (really just packs vectors) in preparation for writing them to file. More...

void	compute_element_maps ()
	This function computes element maps (really just packs vectors) which map the elements to internal and border elements. More...

void	write_exodus_initialization_info (const MeshBase &pmesh, const std::string &title)
	This function writes exodus-specific initialization information. More...

void	read_var_names_impl (const char *var_type, int &count, std::vector< std::string > &result)
	read_var_names() dispatches to this function. More...

void	write_var_names_impl (const char *var_type, int &count, const std::vector< std::string > &names)
	write_var_names() dispatches to this function. More...

void	init_element_equivalence_map ()

void	init_conversion_map ()

Private Attributes
std::map< subdomain_id_type, unsigned >	local_subdomain_counts
	This map keeps track of the number of elements in each subdomain (block) for this processor. More...

std::set< unsigned >	border_node_ids
	The set which will eventually contain the IDs of "border nodes". More...

std::map< unsigned, std::set< unsigned > >	proc_nodes_touched_intersections
	Another map to store sets of intersections with each other processor (other than ourself, of course). More...

std::map< unsigned, std::set< std::pair< unsigned, unsigned > > >	proc_border_elem_sets
	Map between processor ID and (element,side) pairs bordering that processor ID. More...

std::set< unsigned >	internal_node_ids
	A set of internal node IDs for this processor. More...

std::set< unsigned >	internal_elem_ids
	A set of internal elem IDs for this processor. More...

std::set< unsigned >	border_elem_ids
	A set of border elem IDs for this processor. More...

std::map< std::string, ElemType >	element_equivalence_map
	Defines equivalence classes of Exodus element types that map to libmesh ElemTypes. More...

std::map< ElemType, ExodusII_IO_Helper::Conversion >	conversion_map
	Associates libMesh ElemTypes with node/face/edge/etc. More...

Detailed Description

This is the Nemesis_IO_Helper class.

Think of it as a big struct with storage for all the stuff one might want to pull from a Nemesis file. Derived from ExodusII_IO_Helper object, since Nemesis is based on the same file format.

Author: John W. Peterson

Date: 2008

Definition at line 65 of file nemesis_io_helper.h.

Member Typedef Documentation

◆ proc_border_elem_sets_iterator

typedef std::map<unsigned, std::set<std::pair<unsigned,unsigned> > >::iterator libMesh::Nemesis_IO_Helper::proc_border_elem_sets_iterator

private

Typedef for an iterator into the data structure above.

Definition at line 603 of file nemesis_io_helper.h.

◆ proc_nodes_touched_iterator

typedef std::map<unsigned, std::set<unsigned> >::iterator libMesh::Nemesis_IO_Helper::proc_nodes_touched_iterator

private

Typedef for an iterator into the data structure above.

Definition at line 593 of file nemesis_io_helper.h.

Member Enumeration Documentation

◆ ExodusVarType

enum libMesh::ExodusII_IO_Helper::ExodusVarType

inherited

Wraps calls to exII::ex_get_var_names() and exII::ex_get_var_param().

The enumeration controls whether nodal, elemental, or global variable names are read and which class members are filled in. NODAL: num_nodal_vars nodal_var_names ELEMENTAL: num_elem_vars elem_var_names GLOBAL: num_global_vars global_var_names SIDESET: num_sideset_vars sideset_var_names

Enumerator
NODAL
ELEMENTAL
GLOBAL
SIDESET

Definition at line 692 of file exodusII_io_helper.h.

692 {NODAL=0, ELEMENTAL=1, GLOBAL=2, SIDESET=3};

Constructor & Destructor Documentation

◆ Nemesis_IO_Helper()

libMesh::Nemesis_IO_Helper::Nemesis_IO_Helper	(	const ParallelObject &	parent,
		bool	verbose = `false`,
		bool	single_precision = `false`
	)

explicit

Constructor.

Definition at line 48 of file nemesis_io_helper.C.

                                                                              :
   ExodusII_IO_Helper(parent, verbose_in, /*run_only_on_proc0=*/false, /*single_precision=*/single_precision),
   nemesis_err_flag(0),
   num_nodes_global(0),
   num_elems_global(0),
   num_elem_blks_global(0),
   num_node_sets_global(0),
   num_side_sets_global(0),
   num_proc(0),
   num_proc_in_file(0),
   ftype('\0'),
   num_internal_nodes(0),
   num_border_nodes(0),
   num_external_nodes(0),
   num_internal_elems(0),
   num_border_elems(0),
   num_node_cmaps(0),
   num_elem_cmaps(0)
 {
   // Warn about using untested code!
   libmesh_experimental();
 }

◆ ~Nemesis_IO_Helper()

libMesh::Nemesis_IO_Helper::~Nemesis_IO_Helper ( )

virtual

Destructor.

Definition at line 73 of file nemesis_io_helper.C.

 {
   // Our destructor is called from Nemesis_IO.  We close the Exodus file here since we have
   // responsibility for managing the file's lifetime.  Only call ex_update() if the file was
   // opened for writing!
   if (this->opened_for_writing)
     {
       this->ex_err = exII::ex_update(this->ex_id);
       EX_EXCEPTIONLESS_CHECK_ERR(ex_err, "Error flushing buffers to file.");
     }
   this->close();
 }

References libMesh::ExodusII_IO_Helper::close(), libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, and libMesh::ExodusII_IO_Helper::opened_for_writing.

Member Function Documentation

◆ build_element_and_node_maps()

void libMesh::Nemesis_IO_Helper::build_element_and_node_maps ( const MeshBase & pmesh )

private

This function builds the libmesh -> exodus and exodus -> libmesh node and element maps.

These maps allow us to have a consistent numbering scheme within an Exodus file, given an existing globally consistent numbering scheme from LibMesh.

Definition at line 1667 of file nemesis_io_helper.C.

 {
   // If we don't have any local subdomains, it had better be because
   // we don't have any local elements
 #ifdef DEBUG
   if (local_subdomain_counts.empty())
     {
       libmesh_assert(pmesh.active_local_elements_begin() ==
                      pmesh.active_local_elements_end());
       libmesh_assert(this->nodes_attached_to_local_elems.empty());
     }
 #endif
  
   // Elements have to be numbered contiguously based on what block
   // number they are in.  Therefore we have to do a bit of work to get
   // the block (ie subdomain) numbers first and store them off as
   // block_ids.
  
   // Make sure there is no leftover information in the subdomain_map, and reserve
   // enough space to store the elements we need.
   this->subdomain_map.clear();
   for (const auto & pr : local_subdomain_counts)
     {
       if (verbose)
         {
           libMesh::out << "[" << this->processor_id() << "] "
                        << "local_subdomain_counts [" << static_cast<unsigned>(pr.first) << "]= "
                        << pr.second
                        << std::endl;
         }
  
       // pr.first is the subdomain ID, pr.second is the number of elements it contains
       this->subdomain_map[pr.first].reserve(pr.second);
     }
  
  
   // First loop over the elements to figure out which elements are in which subdomain
   for (const auto & elem : pmesh.active_local_element_ptr_range())
     {
       // Grab the nodes while we're here.
       for (auto n : elem->node_index_range())
         this->nodes_attached_to_local_elems.insert( elem->node_id(n) );
  
       subdomain_id_type cur_subdomain = elem->subdomain_id();
  
       this->subdomain_map[cur_subdomain].push_back(elem->id());
     }
  
   // Set num_nodes which is used by exodusII_io_helper
   this->num_nodes =
     cast_int<int>(this->nodes_attached_to_local_elems.size());
  
   // Now come up with a 1-based numbering for these nodes
   this->exodus_node_num_to_libmesh.clear(); // Make sure it's empty
   this->exodus_node_num_to_libmesh.reserve(this->nodes_attached_to_local_elems.size());
  
   // Also make sure there's no leftover information in the map which goes the
   // other direction.
   this->libmesh_node_num_to_exodus.clear();
  
   // Set the map for nodes
   for (const auto & id : nodes_attached_to_local_elems)
     {
       // I.e. given exodus_node_id,
       // exodus_node_num_to_libmesh[ exodus_node_id ] returns the libmesh ID for that node.
       // Note that even though most of Exodus is 1-based, this code will map an Exodus ID of
       // zero to some libmesh node ID.  Is that a problem?
       this->exodus_node_num_to_libmesh.push_back(id);
  
       // Likewise, given libmesh_node_id,
       // libmesh_node_num_to_exodus[ libmesh_node_id ] returns the *Exodus* ID for that node.
       // Unlike the exodus_node_num_to_libmesh vector above, this one is a std::map
       this->libmesh_node_num_to_exodus[id] =
         cast_int<int>(this->exodus_node_num_to_libmesh.size()); // should never be zero...
     }
  
   // Now we're going to loop over the subdomain map and build a few things right
   // now that we'll use later.
  
   // First make sure our data structures don't have any leftover data...
   this->exodus_elem_num_to_libmesh.clear();
   this->block_ids.clear();
   this->libmesh_elem_num_to_exodus.clear();
  
   // Now loop over each subdomain and get a unique numbering for the elements
   for (auto & pr : subdomain_map)
     {
       block_ids.push_back(pr.first);
  
       // Vector of element IDs for this subdomain
       std::vector<dof_id_type> & elem_ids_this_subdomain = pr.second;
  
       // The code below assumes this subdomain block is not empty, make sure that's the case!
       if (elem_ids_this_subdomain.size() == 0)
         libmesh_error_msg("Error, no element IDs found in subdomain " << pr.first);
  
       // Use the first element in this block to get representative information.
       // Note that Exodus assumes all elements in a block are of the same type!
       // We are using that same assumption here!
       const auto & conv = get_conversion
         (pmesh.elem_ref(elem_ids_this_subdomain[0]).type());
       this->num_nodes_per_elem =
         pmesh.elem_ref(elem_ids_this_subdomain[0]).n_nodes();
  
       // Get a reference to the connectivity vector for this subdomain.  This vector
       // is most likely empty, we are going to fill it up now.
       std::vector<int> & current_block_connectivity = this->block_id_to_elem_connectivity[pr.first];
  
       // Just in case it's not already empty...
       current_block_connectivity.clear();
       current_block_connectivity.resize(elem_ids_this_subdomain.size() * this->num_nodes_per_elem);
  
       for (auto i : index_range(elem_ids_this_subdomain))
         {
           auto elem_id = elem_ids_this_subdomain[i];
  
           // Set the number map for elements
           this->exodus_elem_num_to_libmesh.push_back(elem_id);
           this->libmesh_elem_num_to_exodus[elem_id] =
             cast_int<int>(this->exodus_elem_num_to_libmesh.size());
  
           const Elem & elem = pmesh.elem_ref(elem_id);
  
           // Exodus/Nemesis want every block to have the same element type
           // libmesh_assert_equal_to (elem->type(), conv.libmesh_elem_type());
  
           // But we can get away with writing e.g. HEX8 and INFHEX8 in
           // the same block...
           libmesh_assert_equal_to (elem.n_nodes(), Elem::build(conv.libmesh_elem_type(), nullptr)->n_nodes());
  
           for (auto j : IntRange<int>(0, this->num_nodes_per_elem))
             {
               const unsigned int connect_index   = (i*this->num_nodes_per_elem)+j;
               const unsigned int elem_node_index = conv.get_node_map(j);
  
               current_block_connectivity[connect_index] =
                 libmesh_map_find(libmesh_node_num_to_exodus,
                                  elem.node_id(elem_node_index));
             }
         } // End loop over elems in this subdomain
     } // end loop over subdomain_map
 }

Referenced by initialize().

◆ check_existing_vars()

void libMesh::ExodusII_IO_Helper::check_existing_vars	(	ExodusVarType	type,
		std::vector< std::string > &	names,
		std::vector< std::string > &	names_from_file
	)

protectedinherited

When appending: during initialization, check that variable names in the file match those you attempt to initialize with.

Definition at line 2522 of file exodusII_io_helper.C.

 {
   // There may already be global variables in the file (for example,
   // if we're appending) and in that case, we
   // 1.) Cannot initialize them again.
   // 2.) Should check to be sure that the global variable names are the same.
  
   // Fills up names_from_file for us
   this->read_var_names(type);
  
   // Both the number of variables and their names (up to the first
   // MAX_STR_LENGTH characters) must match for the names we are
   // planning to write and the names already in the file.
   bool match =
     std::equal(names.begin(), names.end(),
                names_from_file.begin(),
                [](const std::string & a,
                   const std::string & b) -> bool
                {
                  return a.compare(/*pos=*/0, /*len=*/MAX_STR_LENGTH, b) == 0;
                });
  
   if (!match)
     {
       libMesh::err << "Error! The Exodus file already contains the variables:" << std::endl;
       for (const auto & name : names_from_file)
         libMesh::err << name << std::endl;
  
       libMesh::err << "And you asked to write:" << std::endl;
       for (const auto & name : names)
         libMesh::err << name << std::endl;
  
       libmesh_error_msg("Cannot overwrite existing variables in Exodus II file.");
     }
 }

References equal(), libMesh::err, libMesh::Quality::name(), and libMesh::ExodusII_IO_Helper::read_var_names().

Referenced by libMesh::ExodusII_IO_Helper::initialize_element_variables(), initialize_element_variables(), libMesh::ExodusII_IO_Helper::initialize_global_variables(), and libMesh::ExodusII_IO_Helper::initialize_nodal_variables().

◆ close()

void libMesh::ExodusII_IO_Helper::close ( )

inherited

Closes the ExodusII mesh file.

Definition at line 1228 of file exodusII_io_helper.C.

 {
   // Always call close on processor 0.
   // If we're running on multiple processors, i.e. as one of several Nemesis files,
   // we call close on all processors...
   if ((this->processor_id() == 0) || (!_run_only_on_proc0))
     {
       // Don't close the file if it was never opened, this raises an Exodus error
       if (opened_for_writing || opened_for_reading)
         {
           ex_err = exII::ex_close(ex_id);
           EX_CHECK_ERR(ex_err, "Error closing Exodus file.");
           message("Exodus file closed successfully.");
         }
     }
 }

References libMesh::ExodusII_IO_Helper::_run_only_on_proc0, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::message(), libMesh::ExodusII_IO_Helper::opened_for_reading, libMesh::ExodusII_IO_Helper::opened_for_writing, and libMesh::ParallelObject::processor_id().

Referenced by ~Nemesis_IO_Helper().

◆ comm()

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

inlineinherited

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

Definition at line 94 of file parallel_object.h.

95 { return _communicator; }

References libMesh::ParallelObject::_communicator.

◆ compute_border_node_ids()

void libMesh::Nemesis_IO_Helper::compute_border_node_ids ( const MeshBase & pmesh )

private

This function constructs the set of border node IDs present on the current mesh.

These are nodes which live on the "border" between elements which live on different processors.

Definition at line 1814 of file nemesis_io_helper.C.

 {
   // The set which will eventually contain the IDs of "border nodes".  These are nodes
   // that lie on the boundary between one or more processors.
   //std::set<unsigned> border_node_ids;
  
   // map from processor ID to set of nodes which elements from this processor "touch",
   // that is,
   // proc_nodes_touched[p] = (set all node IDs found in elements owned by processor p)
   std::map<unsigned, std::set<unsigned>> proc_nodes_touched;
  
  
   // We are going to create a lot of intermediate data structures here, so make sure
   // as many as possible all cleaned up by creating scope!
   {
     // Loop over active (not just active local) elements, make sets of node IDs for each
     // processor which has an element that "touches" a node.
     for (const auto & elem : pmesh.active_element_ptr_range())
       {
         // Get reference to the set for this processor.  If it does not exist
         // it will be created.
         std::set<unsigned> & set_p = proc_nodes_touched[ elem->processor_id() ];
  
         // Insert all nodes touched by this element into the set
         for (auto node : elem->node_index_range())
           set_p.insert(elem->node_id(node));
       }
  
     // The number of node communication maps is the number of other processors
     // with which we share nodes. (I think.) This is just the size of the map we just
     // created, minus 1 unless this processor has no nodes of its own.
     this->num_node_cmaps =
       cast_int<int>(proc_nodes_touched.size() -
                     proc_nodes_touched.count(this->processor_id()));
  
     // We can't be connecting to more processors than exist outside
     // ourselves
     libmesh_assert_less (this->num_node_cmaps, this->n_processors());
  
     if (verbose)
       {
         libMesh::out << "[" << this->processor_id()
                      << "] proc_nodes_touched contains "
                      << proc_nodes_touched.size()
                      << " sets of nodes."
                      << std::endl;
  
         for (const auto & pr : proc_nodes_touched)
           libMesh::out << "[" << this->processor_id()
                        << "] proc_nodes_touched[" << pr.first << "] has "
                        << pr.second.size()
                        << " entries."
                        << std::endl;
       }
  
  
     // Loop over all the sets we just created and compute intersections with the
     // this processor's set.  Obviously, don't intersect with ourself.
     for (auto & pr : proc_nodes_touched)
       {
         // Don't compute intersections with ourself
         if (pr.first == this->processor_id())
           continue;
  
         // Otherwise, compute intersection with other processor and ourself
         std::set<unsigned> & my_set = proc_nodes_touched[this->processor_id()];
         std::set<unsigned> & other_set = pr.second;
         std::set<unsigned> & result_set = this->proc_nodes_touched_intersections[pr.first]; // created if does not exist
  
         std::set_intersection(my_set.begin(), my_set.end(),
                               other_set.begin(), other_set.end(),
                               std::inserter(result_set, result_set.end()));
       }
  
     if (verbose)
       {
         for (const auto & pr : proc_nodes_touched_intersections)
           libMesh::out << "[" << this->processor_id()
                        << "] this->proc_nodes_touched_intersections[" << pr.first << "] has "
                        << pr.second.size()
                        << " entries."
                        << std::endl;
       }
  
     // Compute the set_union of all the preceding intersections.  This will be the set of
     // border node IDs for this processor.
     for (auto & pr : proc_nodes_touched_intersections)
       {
         std::set<unsigned> & other_set = pr.second;
         std::set<unsigned> intermediate_result; // Don't think we can insert into one of the sets we're unioning...
  
         std::set_union(this->border_node_ids.begin(), this->border_node_ids.end(),
                        other_set.begin(), other_set.end(),
                        std::inserter(intermediate_result, intermediate_result.end()));
  
         // Swap our intermediate result into the final set
         this->border_node_ids.swap(intermediate_result);
       }
  
     libmesh_assert_less_equal
       (this->proc_nodes_touched_intersections.size(),
        std::size_t(this->num_node_cmaps));
  
     if (verbose)
       {
         libMesh::out << "[" << this->processor_id()
                      << "] border_node_ids.size()=" << this->border_node_ids.size()
                      << std::endl;
       }
   } // end scope for border node ID creation
  
   // Store the number of border node IDs to be written to Nemesis file
   this->num_border_nodes = cast_int<int>(this->border_node_ids.size());
 }

References libMesh::MeshBase::active_element_ptr_range(), border_node_ids, libMesh::ParallelObject::n_processors(), num_border_nodes, num_node_cmaps, libMesh::out, proc_nodes_touched_intersections, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

Referenced by initialize().

◆ compute_communication_map_parameters()

void libMesh::Nemesis_IO_Helper::compute_communication_map_parameters ( )

private

This function determines the communication map parameters which will eventually be written to file.

Definition at line 1147 of file nemesis_io_helper.C.

 {
   // For the nodes, these are the number of entries in the sets in proc_nodes_touched_intersections
   // map computed above.  Note: this map does not contain self-intersections so we can loop over it
   // directly.
   this->node_cmap_node_cnts.clear(); // Make sure we don't have any leftover information...
   this->node_cmap_ids.clear();       // Make sure we don't have any leftover information...
   this->node_cmap_node_cnts.resize(this->num_node_cmaps);
   this->node_cmap_ids.resize(this->num_node_cmaps);
  
   {
     unsigned cnt=0; // Index into the vector
     proc_nodes_touched_iterator
       it = this->proc_nodes_touched_intersections.begin(),
       end = this->proc_nodes_touched_intersections.end();
  
     for (; it != end; ++it)
       {
         this->node_cmap_ids[cnt] = it->first; // The ID of the proc we communicate with
         this->node_cmap_node_cnts[cnt] = cast_int<int>(it->second.size()); // The number of nodes we communicate
         cnt++; // increment vector index!
       }
   }
  
   // Print the packed vectors we just filled
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] node_cmap_node_cnts = ";
       for (const auto & node_cnt : node_cmap_node_cnts)
         libMesh::out << node_cnt << ", ";
       libMesh::out << std::endl;
  
       libMesh::out << "[" << this->processor_id() << "] node_cmap_ids = ";
       for (const auto & node_id : node_cmap_ids)
         libMesh::out << node_id << ", ";
       libMesh::out << std::endl;
     }
  
   // For the elements, we have not yet computed all this information..
   this->elem_cmap_elem_cnts.clear(); // Make sure we don't have any leftover information...
   this->elem_cmap_ids.clear();       // Make sure we don't have any leftover information...
   this->elem_cmap_elem_cnts.resize(this->num_elem_cmaps);
   this->elem_cmap_ids.resize(this->num_elem_cmaps);
  
   // Pack the elem_cmap_ids and elem_cmap_elem_cnts vectors
   {
     unsigned cnt=0; // Index into the vectors we're filling
     proc_border_elem_sets_iterator
       it = this->proc_border_elem_sets.begin(),
       end = this->proc_border_elem_sets.end();
  
     for (; it != end; ++it)
       {
         this->elem_cmap_ids[cnt] = it->first; // The ID of the proc we communicate with
         this->elem_cmap_elem_cnts[cnt] = cast_int<int>(it->second.size()); // The number of elems we communicate to/from that proc
         cnt++; // increment vector index!
       }
   }
  
   // Print the packed vectors we just filled
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] elem_cmap_elem_cnts = ";
       for (const auto & elem_cnt : elem_cmap_elem_cnts)
         libMesh::out << elem_cnt << ", ";
       libMesh::out << std::endl;
  
       libMesh::out << "[" << this->processor_id() << "] elem_cmap_ids = ";
       for (const auto & elem_id : elem_cmap_ids)
         libMesh::out << elem_id << ", ";
       libMesh::out << std::endl;
     }
 }

References elem_cmap_elem_cnts, elem_cmap_ids, end, node_cmap_ids, node_cmap_node_cnts, num_elem_cmaps, num_node_cmaps, libMesh::out, proc_border_elem_sets, proc_nodes_touched_intersections, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

Referenced by initialize().

◆ compute_elem_communication_maps()

void libMesh::Nemesis_IO_Helper::compute_elem_communication_maps ( )

private

This function computes element communication maps (really just packs vectors) in preparation for writing them to file.

Definition at line 995 of file nemesis_io_helper.C.

 {
   // Make sure there is no leftover information
   this->elem_cmap_elem_ids.clear();
   this->elem_cmap_side_ids.clear();
   this->elem_cmap_proc_ids.clear();
  
   // Allocate enough space for all our element maps
   this->elem_cmap_elem_ids.resize(this->num_elem_cmaps);
   this->elem_cmap_side_ids.resize(this->num_elem_cmaps);
   this->elem_cmap_proc_ids.resize(this->num_elem_cmaps);
   {
     unsigned cnt=0; // Index into vectors
     proc_border_elem_sets_iterator
       it = this->proc_border_elem_sets.begin(),
       end = this->proc_border_elem_sets.end();
  
     for (; it != end; ++it)
       {
         // Make sure the current elem_cmap_id matches the index in our map of node intersections
         libmesh_assert_equal_to (static_cast<unsigned>(this->elem_cmap_ids[cnt]), it->first);
  
         // Get reference to the set of IDs to be packed into the vector
         std::set<std::pair<unsigned,unsigned>> & elem_set = it->second;
  
         // Resize the vectors to receive their payload
         this->elem_cmap_elem_ids[cnt].resize(elem_set.size());
         this->elem_cmap_side_ids[cnt].resize(elem_set.size());
         this->elem_cmap_proc_ids[cnt].resize(elem_set.size());
  
         std::set<std::pair<unsigned,unsigned>>::iterator elem_set_iter = elem_set.begin();
  
         // Pack the vectors with elem IDs, side IDs, and processor IDs.
         for (std::size_t j=0, eceis=this->elem_cmap_elem_ids[cnt].size(); j<eceis; ++j, ++elem_set_iter)
           {
             this->elem_cmap_elem_ids[cnt][j] =
               libmesh_map_find(libmesh_elem_num_to_exodus, elem_set_iter->first);
             this->elem_cmap_side_ids[cnt][j] = elem_set_iter->second;     // Side ID, this has already been converted above
             this->elem_cmap_proc_ids[cnt][j] = it->first; // All have the same processor ID
           }
  
         // increment vector index to go to next processor
         cnt++;
       }
   } // end scope for packing
 }

References elem_cmap_elem_ids, elem_cmap_ids, elem_cmap_proc_ids, elem_cmap_side_ids, end, libMesh::ExodusII_IO_Helper::libmesh_elem_num_to_exodus, num_elem_cmaps, and proc_border_elem_sets.

Referenced by initialize().

◆ compute_element_maps()

void libMesh::Nemesis_IO_Helper::compute_element_maps ( )

private

This function computes element maps (really just packs vectors) which map the elements to internal and border elements.

Definition at line 970 of file nemesis_io_helper.C.

 {
   // Make sure we don't have any leftover info
   this->elem_mapi.clear();
   this->elem_mapb.clear();
  
   // Copy set contents into vectors
   this->elem_mapi.resize(this->internal_elem_ids.size());
   this->elem_mapb.resize(this->border_elem_ids.size());
  
   {
     unsigned cnt = 0;
     for (const auto & id : this->internal_elem_ids)
       this->elem_mapi[cnt++] = libmesh_map_find(libmesh_elem_num_to_exodus, id);
   }
  
   {
     unsigned cnt = 0;
     for (const auto & id : this->border_elem_ids)
       this->elem_mapb[cnt++] = libmesh_map_find(libmesh_elem_num_to_exodus, id);
   }
 }

References border_elem_ids, elem_mapb, elem_mapi, internal_elem_ids, and libMesh::ExodusII_IO_Helper::libmesh_elem_num_to_exodus.

Referenced by initialize().

◆ compute_internal_and_border_elems_and_internal_nodes()

void libMesh::Nemesis_IO_Helper::compute_internal_and_border_elems_and_internal_nodes ( const MeshBase & pmesh )

private

This function constructs the set of border and internal element IDs and internal node IDs present on the current mesh.

Definition at line 1225 of file nemesis_io_helper.C.

 {
   // Set of all local, active element IDs.  After we have identified border element
   // IDs, the set_difference between this set and the border_elem_ids set will give us
   // the set of internal_elem_ids.
   std::set<unsigned> all_elem_ids;
  
   // A set of processor IDs which elements on this processor have as
   // neighbors.  The size of this set will determine the number of
   // element communication maps in Exodus.
   std::set<unsigned> neighboring_processor_ids;
  
   for (const auto & elem : pmesh.active_local_element_ptr_range())
     {
       // Add this Elem's ID to all_elem_ids, later we will take the difference
       // between this set and the set of border_elem_ids, to get the set of
       // internal_elem_ids.
       all_elem_ids.insert(elem->id());
  
       // Will be set to true if element is determined to be a border element
       bool is_border_elem = false;
  
       // Construct a conversion object for this Element.  This will help us map
       // Libmesh numberings into Nemesis numberings for sides.
       const auto & conv = get_conversion(elem->type());
  
       // Add all this element's node IDs to the set of all node IDs.
       // The set of internal_node_ids will be the set difference between
       // the set of all nodes and the set of border nodes.
       //
       // In addition, if any node of a local node is listed in the
       // border nodes list, then this element goes into the proc_border_elem_sets.
       // Note that there is not a 1:1 correspondence between
       // border_elem_ids and the entries which go into proc_border_elem_sets.
       // The latter is for communication purposes, ie determining which elements
       // should be shared between processors.
       for (auto node : elem->node_index_range())
         this->nodes_attached_to_local_elems.insert(elem->node_id(node));
  
       // Loop over element's neighbors, see if it has a neighbor which is off-processor
       for (auto n : elem->side_index_range())
         {
           if (elem->neighbor_ptr(n) != nullptr)
             {
               unsigned neighbor_proc_id = elem->neighbor_ptr(n)->processor_id();
  
               // If my neighbor has a different processor ID, I must be a border element.
               // Also track the neighboring processor ID if it is are different from our processor ID
               if (neighbor_proc_id != this->processor_id())
                 {
                   is_border_elem = true;
                   neighboring_processor_ids.insert(neighbor_proc_id);
  
                   // Convert libmesh side(n) of this element into a side ID for Nemesis
                   unsigned nemesis_side_id = conv.get_inverse_side_map(n);
  
                   if (verbose)
                     libMesh::out << "[" << this->processor_id() << "] LibMesh side "
                                  << n
                                  << " mapped to (1-based) Exodus side "
                                  << nemesis_side_id
                                  << std::endl;
  
                   // Add this element's ID and the ID of the side which is on the boundary
                   // to the set of border elements for this processor.
                   // Note: if the set does not already exist, this creates it.
                   this->proc_border_elem_sets[ neighbor_proc_id ].insert( std::make_pair(elem->id(), nemesis_side_id) );
                 }
             }
         } // end for loop over neighbors
  
       // If we're on a border element, add it to the set
       if (is_border_elem)
         this->border_elem_ids.insert( elem->id() );
  
     } // end for loop over active local elements
  
   // Take the set_difference between all elements and border elements to get internal
   // element IDs
   std::set_difference(all_elem_ids.begin(), all_elem_ids.end(),
                       this->border_elem_ids.begin(), this->border_elem_ids.end(),
                       std::inserter(this->internal_elem_ids, this->internal_elem_ids.end()));
  
   // Take the set_difference between all nodes and border nodes to get internal nodes
   std::set_difference(this->nodes_attached_to_local_elems.begin(), this->nodes_attached_to_local_elems.end(),
                       this->border_node_ids.begin(), this->border_node_ids.end(),
                       std::inserter(this->internal_node_ids, this->internal_node_ids.end()));
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] neighboring_processor_ids = ";
       for (const auto & id : neighboring_processor_ids)
         libMesh::out << id << " ";
       libMesh::out << std::endl;
     }
  
   // The size of the neighboring_processor_ids set should be the number of element communication maps
   this->num_elem_cmaps =
     cast_int<int>(neighboring_processor_ids.size());
  
   if (verbose)
     libMesh::out << "[" << this->processor_id() << "] "
                  << "Number of neighboring processor IDs="
                  << this->num_elem_cmaps
                  << std::endl;
  
   if (verbose)
     {
       // Print out counts of border elements for each processor
       for (const auto & pr : proc_border_elem_sets)
         {
           libMesh::out << "[" << this->processor_id() << "] "
                        << "Proc "
                        << pr.first << " communicates "
                        << pr.second.size() << " elements." << std::endl;
         }
     }
  
   // Store the number of internal and border elements, and the number of internal nodes,
   // to be written to the Nemesis file.
   this->num_internal_elems =
     cast_int<int>(this->internal_elem_ids.size());
   this->num_border_elems   =
     cast_int<int>(this->border_elem_ids.size());
   this->num_internal_nodes =
     cast_int<int>(this->internal_node_ids.size());
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] num_internal_nodes=" << this->num_internal_nodes << std::endl;
       libMesh::out << "[" << this->processor_id() << "] num_border_nodes=" << this->num_border_nodes << std::endl;
       libMesh::out << "[" << this->processor_id() << "] num_border_elems=" << this->num_border_elems << std::endl;
       libMesh::out << "[" << this->processor_id() << "] num_internal_elems=" << this->num_internal_elems << std::endl;
     }
 }

References libMesh::MeshBase::active_local_element_ptr_range(), border_elem_ids, border_node_ids, libMesh::ExodusII_IO_Helper::get_conversion(), internal_elem_ids, internal_node_ids, nodes_attached_to_local_elems, num_border_elems, num_border_nodes, num_elem_cmaps, num_internal_elems, num_internal_nodes, libMesh::out, proc_border_elem_sets, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

Referenced by initialize().

◆ compute_node_communication_maps()

void libMesh::Nemesis_IO_Helper::compute_node_communication_maps ( )

private

Compute the node communication maps (really just pack vectors) in preparation for writing them to file.

Definition at line 1075 of file nemesis_io_helper.C.

 {
   // Make sure there's no left-over information
   this->node_cmap_node_ids.clear();
   this->node_cmap_proc_ids.clear();
  
   libmesh_assert_less_equal
     (this->proc_nodes_touched_intersections.size(),
      std::size_t(this->num_node_cmaps));
  
   // Allocate enough space for all our node maps
   this->node_cmap_node_ids.resize(this->num_node_cmaps);
   this->node_cmap_proc_ids.resize(this->num_node_cmaps);
   {
     unsigned cnt=0; // Index into vectors
     proc_nodes_touched_iterator
       it = this->proc_nodes_touched_intersections.begin(),
       end = this->proc_nodes_touched_intersections.end();
  
     for (; it != end; ++it)
       {
         // Make sure the current node_cmap_id matches the index in our map of node intersections
         libmesh_assert_equal_to (static_cast<unsigned>(this->node_cmap_ids[cnt]), it->first);
  
         // Get reference to the set of IDs to be packed into the vector.
         std::set<unsigned> & node_set = it->second;
  
         // Resize the vectors to receive their payload
         this->node_cmap_node_ids[cnt].resize(node_set.size());
         this->node_cmap_proc_ids[cnt].resize(node_set.size());
  
         std::set<unsigned>::iterator node_set_iter = node_set.begin();
  
         // Pack the vectors with node IDs and processor IDs.
         for (std::size_t j=0, nceis=this->node_cmap_node_ids[cnt].size(); j<nceis; ++j, ++node_set_iter)
           {
             this->node_cmap_node_ids[cnt][j] =
               libmesh_map_find(libmesh_node_num_to_exodus, *node_set_iter);
             this->node_cmap_proc_ids[cnt][j] = it->first;
           }
  
         // increment vector index to go to next processor
         cnt++;
       }
   } // end scope for packing
  
   // Print out the vectors we just packed
   if (verbose)
     {
       for (auto i : index_range(this->node_cmap_node_ids))
         {
           libMesh::out << "[" << this->processor_id() << "] nodes communicated to proc "
                        << this->node_cmap_ids[i]
                        << " = ";
           for (const auto & node_id : this->node_cmap_node_ids[i])
             libMesh::out << node_id << " ";
           libMesh::out << std::endl;
         }
  
       for (const auto & id_vec : this->node_cmap_node_ids)
         {
           libMesh::out << "[" << this->processor_id() << "] processor ID node communicated to = ";
           for (const auto & proc_id : id_vec)
             libMesh::out << proc_id << " ";
           libMesh::out << std::endl;
         }
     }
 }

References end, libMesh::index_range(), libMesh::ExodusII_IO_Helper::libmesh_node_num_to_exodus, node_cmap_ids, node_cmap_node_ids, node_cmap_proc_ids, node_set, num_node_cmaps, libMesh::out, proc_nodes_touched_intersections, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

Referenced by initialize().

◆ compute_node_maps()

void libMesh::Nemesis_IO_Helper::compute_node_maps ( )

private

Compute the node maps (really just pack vectors) which map the nodes to internal, border, and external nodes in the file.

Definition at line 1046 of file nemesis_io_helper.C.

 {
   // Make sure we don't have any leftover information
   this->node_mapi.clear();
   this->node_mapb.clear();
   this->node_mape.clear();
  
   // Make sure there's enough space to hold all our node IDs
   this->node_mapi.resize(this->internal_node_ids.size());
   this->node_mapb.resize(this->border_node_ids.size());
  
   // Copy set contents into vectors
   {
     unsigned cnt = 0;
     for (const auto & id : this->internal_node_ids)
       this->node_mapi[cnt++] = libmesh_map_find(libmesh_node_num_to_exodus, id);
   }
  
   {
     unsigned cnt=0;
     for (const auto & id : this->border_node_ids)
       this->node_mapb[cnt++] = libmesh_map_find(libmesh_node_num_to_exodus, id);
   }
 }

References border_node_ids, internal_node_ids, libMesh::ExodusII_IO_Helper::libmesh_node_num_to_exodus, node_mapb, node_mape, and node_mapi.

Referenced by initialize().

◆ compute_num_global_elem_blocks()

void libMesh::Nemesis_IO_Helper::compute_num_global_elem_blocks ( const MeshBase & pmesh )

private

This function uses global communication routines to determine the number of element blocks across the entire mesh.

Definition at line 1567 of file nemesis_io_helper.C.

 {
   // 1.) Loop over active local elements, build up set of subdomain IDs.
   std::set<subdomain_id_type> global_subdomain_ids;
  
   // This map keeps track of the number of elements in each subdomain over all processors
   std::map<subdomain_id_type, unsigned> global_subdomain_counts;
  
   for (const auto & elem : pmesh.active_local_element_ptr_range())
     {
       subdomain_id_type cur_subdomain = elem->subdomain_id();
  
       /*
       // We can't have a zero subdomain ID in Exodus (for some reason?)
       // so map zero subdomains to a max value...
       if (cur_subdomain == 0)
       cur_subdomain = std::numeric_limits<subdomain_id_type>::max();
       */
  
       global_subdomain_ids.insert(cur_subdomain);
  
       // Increment the count of elements in this subdomain
       global_subdomain_counts[cur_subdomain]++;
     }
  
   // We're next going to this->comm().sum the subdomain counts, so save the local counts
   this->local_subdomain_counts = global_subdomain_counts;
  
   {
     // 2.) Copy local subdomain IDs into a vector for communication
     std::vector<subdomain_id_type> global_subdomain_ids_vector(global_subdomain_ids.begin(),
                                                                global_subdomain_ids.end());
  
     // 3.) Gather them into an enlarged vector
     this->comm().allgather(global_subdomain_ids_vector);
  
     // 4.) Insert any new IDs into the set (any duplicates will be dropped)
     global_subdomain_ids.insert(global_subdomain_ids_vector.begin(),
                                 global_subdomain_ids_vector.end());
   }
  
   // 5.) Now global_subdomain_ids actually contains a global list of all subdomain IDs
   this->num_elem_blks_global =
     cast_int<int>(global_subdomain_ids.size());
  
   // Print the number of elements found locally in each subdomain
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] ";
       for (const auto & pr : global_subdomain_counts)
         {
           libMesh::out << "ID: "
                        << static_cast<unsigned>(pr.first)
                        << ", Count: " << pr.second << ", ";
         }
       libMesh::out << std::endl;
     }
  
   // 6.) this->comm().sum up the number of elements in each block.  We know the global
   // subdomain IDs, so pack them into a vector one by one.  Use a vector of int since
   // that is what Nemesis wants
   this->global_elem_blk_cnts.resize(global_subdomain_ids.size());
  
   unsigned cnt=0;
   // Find the entry in the local map, note: if not found, will be created with 0 default value, which is OK...
   for (const auto & id : global_subdomain_ids)
     this->global_elem_blk_cnts[cnt++] = global_subdomain_counts[id];
  
   // Sum up subdomain counts from all processors
   this->comm().sum(this->global_elem_blk_cnts);
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] global_elem_blk_cnts = ";
       for (const auto & bc : this->global_elem_blk_cnts)
         libMesh::out << bc << ", ";
       libMesh::out << std::endl;
     }
  
   // 7.) Create a vector<int> from the global_subdomain_ids set, for passing to Nemesis
   this->global_elem_blk_ids.clear();
   this->global_elem_blk_ids.insert(this->global_elem_blk_ids.end(), // pos
                                    global_subdomain_ids.begin(),
                                    global_subdomain_ids.end());
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] global_elem_blk_ids = ";
       for (const auto & id : this->global_elem_blk_ids)
         libMesh::out << id << ", ";
       libMesh::out << std::endl;
     }
  
  
   // 8.) We will call put_eb_info_global later, it must be called after this->put_init_global().
 }

References libMesh::MeshBase::active_local_element_ptr_range(), libMesh::ParallelObject::comm(), global_elem_blk_cnts, global_elem_blk_ids, local_subdomain_counts, num_elem_blks_global, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

Referenced by initialize().

◆ compute_num_global_nodesets()

void libMesh::Nemesis_IO_Helper::compute_num_global_nodesets ( const MeshBase & pmesh )

private

This function uses global communication routines to determine the number of nodesets across the entire mesh.

Definition at line 1456 of file nemesis_io_helper.C.

 {
   std::set<boundary_id_type> local_node_boundary_ids;
  
   // 1.) Get reference to the set of node boundary IDs *for this processor*
   std::set<boundary_id_type> global_node_boundary_ids
     (pmesh.get_boundary_info().get_node_boundary_ids().begin(),
      pmesh.get_boundary_info().get_node_boundary_ids().end());
  
   // Save a copy of the local_node_boundary_ids...
   local_node_boundary_ids = global_node_boundary_ids;
  
   // 2.) Gather boundary node IDs from other processors
   this->comm().set_union(global_node_boundary_ids);
  
   // 3.) Now global_node_boundary_ids actually contains a global list of all node boundary IDs
   this->num_node_sets_global =
     cast_int<int>(global_node_boundary_ids.size());
  
   // 4.) Create a vector<int> from the global_node_boundary_ids set
   this->global_nodeset_ids.clear();
   this->global_nodeset_ids.insert(this->global_nodeset_ids.end(),
                                   global_node_boundary_ids.begin(),
                                   global_node_boundary_ids.end());
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] global_nodeset_ids = ";
       for (const auto & id : global_nodeset_ids)
         libMesh::out << id << ", ";
       libMesh::out << std::endl;
  
       libMesh::out << "[" << this->processor_id() << "] local_node_boundary_ids = ";
       for (const auto & id : local_node_boundary_ids)
         libMesh::out << id << ", ";
       libMesh::out << std::endl;
     }
  
   // 7.) We also need to know the number of nodes which is in each of the nodesets, globally.
  
   // Build list of (node-id, bc-id) tuples.
   typedef std::tuple<dof_id_type, boundary_id_type> Tuple;
   std::vector<Tuple> bc_tuples = pmesh.get_boundary_info().build_node_list();
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] boundary_node_list.size()="
                    << bc_tuples.size() << std::endl;
       libMesh::out << "[" << this->processor_id() << "] (boundary_node_id, boundary_id) = ";
       for (const auto & t : bc_tuples)
         libMesh::out << "(" << std::get<0>(t) << ", " << std::get<1>(t) << ") ";
       libMesh::out << std::endl;
     }
  
   // Now get the global information.  In this case, we only want to count boundary
   // information for nodes *owned* by this processor, so there are no duplicates.
  
   // Make sure we don't have any left over information
   this->num_global_node_counts.clear();
   this->num_global_node_counts.resize(this->global_nodeset_ids.size());
  
   // Unfortunately, we can't just count up all occurrences of a given id,
   // that would give us duplicate entries when we do the parallel summation.
   // So instead, only count entries for nodes owned by this processor.
   // Start by getting rid of all non-local node entries from the vectors.
   std::vector<Tuple>::iterator
     it = bc_tuples.begin(),
     new_end = bc_tuples.end();
  
   while (it != new_end)
     {
       if (pmesh.node_ptr(std::get<0>(*it))->processor_id() != this->processor_id())
         {
           // Back up the new end iterators to prepare for swap
           --new_end;
  
           // Swap places, the non-local node will now be "past-the-end"
           std::swap(*it, *new_end);
         }
       else // node is local, go to next
         ++it;
     }
  
   // Erase from "new" end to old end.
   bc_tuples.erase(new_end, bc_tuples.end());
  
   // Now we can do the local count for each ID...
   for (auto i : index_range(global_nodeset_ids))
     {
       int id = this->global_nodeset_ids[i];
       this->num_global_node_counts[i] =
         cast_int<int>(std::count_if(bc_tuples.begin(),
                                     bc_tuples.end(),
                                     [id](const Tuple & t)->bool { return std::get<1>(t) == id; }));
     }
  
   // And finally we can sum them up
   this->comm().sum(this->num_global_node_counts);
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] num_global_node_counts = ";
       for (const auto & cnt : num_global_node_counts)
         libMesh::out << cnt << ", ";
       libMesh::out << std::endl;
     }
 }

References libMesh::BoundaryInfo::build_node_list(), libMesh::ParallelObject::comm(), libMesh::MeshBase::get_boundary_info(), libMesh::BoundaryInfo::get_node_boundary_ids(), global_nodeset_ids, libMesh::index_range(), libMesh::MeshBase::node_ptr(), num_global_node_counts, num_node_sets_global, libMesh::out, libMesh::ParallelObject::processor_id(), libMesh::DofObject::processor_id(), swap(), and libMesh::ExodusII_IO_Helper::verbose.

Referenced by initialize().

◆ compute_num_global_sidesets()

void libMesh::Nemesis_IO_Helper::compute_num_global_sidesets ( const MeshBase & pmesh )

private

This function uses global communication routines to determine the number of sidesets across the entire mesh.

Definition at line 1363 of file nemesis_io_helper.C.

 {
   // 1.) Get reference to the set of side boundary IDs
   std::set<boundary_id_type> global_side_boundary_ids
     (pmesh.get_boundary_info().get_side_boundary_ids().begin(),
      pmesh.get_boundary_info().get_side_boundary_ids().end());
  
   // 2.) Gather boundary side IDs from other processors
   this->comm().set_union(global_side_boundary_ids);
  
   // 3.) Now global_side_boundary_ids actually contains a global list of all side boundary IDs
   this->num_side_sets_global =
     cast_int<int>(global_side_boundary_ids.size());
  
   // 4.) Pack these sidesets into a vector so they can be written by Nemesis
   this->global_sideset_ids.clear(); // Make sure there is no leftover information
   this->global_sideset_ids.insert(this->global_sideset_ids.end(),
                                   global_side_boundary_ids.begin(),
                                   global_side_boundary_ids.end());
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] global_sideset_ids = ";
       for (const auto & id : this->global_sideset_ids)
         libMesh::out << id << ", ";
       libMesh::out << std::endl;
     }
  
   // We also need global counts of sides in each of the sidesets.
   // Build a list of (elem, side, bc) tuples.
   typedef std::tuple<dof_id_type, unsigned short int, boundary_id_type> Tuple;
   std::vector<Tuple> bc_triples = pmesh.get_boundary_info().build_side_list();
  
   // Iterators to the beginning and end of the current range.
   std::vector<Tuple>::iterator
     it = bc_triples.begin(),
     new_end = bc_triples.end();
  
   while (it != new_end)
     {
       if (pmesh.elem_ref(std::get<0>(*it)).processor_id() != this->processor_id())
         {
           // Back up the new end iterators to prepare for swap
           --new_end;
  
           // Swap places, the non-local elem will now be "past-the-end"
           std::swap (*it, *new_end);
         }
       else // elem is local, go to next
         ++it;
     }
  
   // Erase from "new" end to old.
   bc_triples.erase(new_end, bc_triples.end());
  
   this->num_global_side_counts.clear(); // Make sure we don't have any leftover information
   this->num_global_side_counts.resize(this->global_sideset_ids.size());
  
   // Get the count for each global sideset ID
   for (auto i : index_range(global_sideset_ids))
     {
       int id = global_sideset_ids[i];
       this->num_global_side_counts[i] =
         cast_int<int>(std::count_if(bc_triples.begin(),
                                     bc_triples.end(),
                                     [id](const Tuple & t)->bool { return std::get<2>(t) == id; }));
     }
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] num_global_side_counts = ";
       for (const auto & cnt : this->num_global_side_counts)
         libMesh::out << cnt << ", ";
       libMesh::out << std::endl;
     }
  
   // Finally sum up the result
   this->comm().sum(this->num_global_side_counts);
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] num_global_side_counts = ";
       for (const auto & cnt : this->num_global_side_counts)
         libMesh::out << cnt << ", ";
       libMesh::out << std::endl;
     }
 }

References libMesh::BoundaryInfo::build_side_list(), libMesh::ParallelObject::comm(), libMesh::MeshBase::elem_ref(), libMesh::MeshBase::get_boundary_info(), libMesh::BoundaryInfo::get_side_boundary_ids(), global_sideset_ids, libMesh::index_range(), num_global_side_counts, num_side_sets_global, libMesh::out, libMesh::ParallelObject::processor_id(), libMesh::DofObject::processor_id(), swap(), and libMesh::ExodusII_IO_Helper::verbose.

Referenced by initialize().

◆ construct_nemesis_filename()

std::string libMesh::Nemesis_IO_Helper::construct_nemesis_filename ( const std::string & base_filename )

Given base_filename, foo.e, constructs the Nemesis filename foo.e.X.Y, where X=n.

CPUs and Y=processor ID

Definition at line 2663 of file nemesis_io_helper.C.

 {
   // Build a filename for this processor.  This code is cut-n-pasted from the read function
   // and should probably be put into a separate function...
   std::ostringstream file_oss;
  
   // We have to be a little careful here: Nemesis left pads its file
   // numbers based on the number of processors, so for example on 10
   // processors, we'd have:
   // mesh.e.10.00
   // mesh.e.10.01
   // mesh.e.10.02
   // ...
   // mesh.e.10.09
  
   // And on 100 you'd have:
   // mesh.e.100.000
   // mesh.e.100.001
   // ...
   // mesh.e.128.099
  
   // Find the length of the highest processor ID
   file_oss << (this->n_processors());
   unsigned int field_width = cast_int<unsigned int>(file_oss.str().size());
  
   if (verbose)
     libMesh::out << "field_width=" << field_width << std::endl;
  
   file_oss.str(""); // reset the string stream
   file_oss << base_filename
            << '.' << this->n_processors()
            << '.' << std::setfill('0') << std::setw(field_width) << this->processor_id();
  
   // Return the resulting string
   return file_oss.str();
 }

References libMesh::ParallelObject::n_processors(), libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ create()

void libMesh::Nemesis_IO_Helper::create ( std::string filename )

overridevirtual

This function is specialized from ExodusII_IO_Helper to create the nodal coordinates stored on the local piece of the Mesh.

Reimplemented from libMesh::ExodusII_IO_Helper.

Definition at line 730 of file nemesis_io_helper.C.

 {
   // Fall back on double precision when necessary since ExodusII
   // doesn't seem to support long double
   int
     comp_ws = 0,
     io_ws = 0;
  
   if (_single_precision)
     {
       comp_ws = sizeof(float);
       io_ws = sizeof(float);
     }
   else
     {
       comp_ws = cast_int<int>(std::min(sizeof(Real), sizeof(double)));
       io_ws = cast_int<int>(std::min(sizeof(Real), sizeof(double)));
     }
  
   this->ex_id = exII::ex_create(filename.c_str(), EX_CLOBBER, &comp_ws, &io_ws);
  
   EX_CHECK_ERR(ex_id, "Error creating Nemesis mesh file.");
  
   if (verbose)
     libMesh::out << "File created successfully." << std::endl;
  
   this->opened_for_writing = true;
 }

References libMesh::ExodusII_IO_Helper::_single_precision, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::opened_for_writing, libMesh::out, libMesh::Real, and libMesh::ExodusII_IO_Helper::verbose.

◆ get_block_id()

int libMesh::ExodusII_IO_Helper::get_block_id ( int index )

inherited

Get the block number for the given block index.

Definition at line 742 of file exodusII_io_helper.C.

 {
   libmesh_assert_less (index, block_ids.size());
  
   return block_ids[index];
 }

References libMesh::ExodusII_IO_Helper::block_ids.

Referenced by libMesh::ExodusII_IO_Helper::write_element_values(), and libMesh::ExodusII_IO_Helper::write_element_values_element_major().

◆ get_block_name()

std::string libMesh::ExodusII_IO_Helper::get_block_name ( int index )

inherited

Get the block name for the given block index if supplied in the mesh file.

Otherwise an empty string is returned.

Definition at line 751 of file exodusII_io_helper.C.

 {
   libmesh_assert_less (index, block_ids.size());
  
   return id_to_block_names[block_ids[index]];
 }

References libMesh::ExodusII_IO_Helper::block_ids, and libMesh::ExodusII_IO_Helper::id_to_block_names.

◆ get_cmap_params()

void libMesh::Nemesis_IO_Helper::get_cmap_params ( )

Definition at line 325 of file nemesis_io_helper.C.

 {
   node_cmap_ids.resize(num_node_cmaps);
   node_cmap_node_cnts.resize(num_node_cmaps);
   elem_cmap_ids.resize(num_elem_cmaps);
   elem_cmap_elem_cnts.resize(num_elem_cmaps);
  
   nemesis_err_flag =
     Nemesis::ne_get_cmap_params(ex_id,
                                 node_cmap_ids.empty()       ? nullptr : node_cmap_ids.data(),
                                 node_cmap_node_cnts.empty() ? nullptr : node_cmap_node_cnts.data(),
                                 elem_cmap_ids.empty()       ? nullptr : elem_cmap_ids.data(),
                                 elem_cmap_elem_cnts.empty() ? nullptr : elem_cmap_elem_cnts.data(),
                                 this->processor_id());
   EX_CHECK_ERR(nemesis_err_flag, "Error reading cmap parameters!");
  
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] ";
       for (auto i : index_range(node_cmap_ids))
         libMesh::out << "node_cmap_ids[" << i << "]=" << node_cmap_ids[i] << " ";
       libMesh::out << std::endl;
  
       libMesh::out << "[" << this->processor_id() << "] ";
       for (auto i : index_range(node_cmap_node_cnts))
         libMesh::out << "node_cmap_node_cnts[" << i << "]=" << node_cmap_node_cnts[i] << " ";
       libMesh::out << std::endl;
  
       libMesh::out << "[" << this->processor_id() << "] ";
       for (auto i : index_range(elem_cmap_ids))
         libMesh::out << "elem_cmap_ids[" << i << "]=" << elem_cmap_ids[i] << " ";
       libMesh::out << std::endl;
  
       libMesh::out << "[" << this->processor_id() << "] ";
       for (auto i : index_range(elem_cmap_elem_cnts))
         libMesh::out << "elem_cmap_elem_cnts[" << i << "]=" << elem_cmap_elem_cnts[i] << " ";
       libMesh::out << std::endl;
     }
 }

References elem_cmap_elem_cnts, elem_cmap_ids, libMesh::ExodusII_IO_Helper::ex_id, libMesh::index_range(), nemesis_err_flag, node_cmap_ids, node_cmap_node_cnts, num_elem_cmaps, num_node_cmaps, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ get_complex_names()

std::vector< std::string > libMesh::ExodusII_IO_Helper::get_complex_names ( const std::vector< std::string > & names ) const

inherited

Returns: A vector with three copies of each element in the provided name vector, starting with r_, i_ and a_ respectively.

Definition at line 3188 of file exodusII_io_helper.C.

 {
   std::vector<std::string> complex_names;
  
   // This will loop over all names and create new "complex" names
   // (i.e. names that start with r_, i_ or a_
   for (const auto & name : names)
     {
       complex_names.push_back("r_" + name);
       complex_names.push_back("i_" + name);
       complex_names.push_back("a_" + name);
     }
  
   return complex_names;
 }

References libMesh::Quality::name().

◆ get_complex_subdomain_to_var_names()

std::map< subdomain_id_type, std::vector< std::string > > libMesh::ExodusII_IO_Helper::get_complex_subdomain_to_var_names ( const std::map< subdomain_id_type, std::vector< std::string >> & subdomain_to_var_names ) const

inherited

Takes a map from subdomain id -> vector of active variable names as input and returns a corresponding map where the original variable names have been replaced by their complex counterparts.

Used by the ExodusII_IO::write_element_data_from_discontinuous_nodal_data() function.

Definition at line 3227 of file exodusII_io_helper.C.

 {
   // Eventual return value
   std::map<subdomain_id_type, std::vector<std::string>> ret;
  
   for (const auto & pr : subdomain_to_var_names)
     {
       // Initialize entry for current subdomain
       auto & vec = ret[pr.first];
  
       // Get list of non-complex variable names active on this subdomain.
       const auto & varnames = pr.second;
  
       // Allocate space for 3x the number of entries
       vec.reserve(3 * varnames.size());
  
       // For each varname in the input map, write three variable names
       // to the output formed by prepending "r_", "i_", and "a_",
       // respectively.
       for (const auto & varname : varnames)
         {
           vec.push_back("r_" + varname);
           vec.push_back("i_" + varname);
           vec.push_back("a_" + varname);
         }
     }
   return ret;
 }

◆ get_complex_vars_active_subdomains()

std::vector< std::set< subdomain_id_type > > libMesh::ExodusII_IO_Helper::get_complex_vars_active_subdomains ( const std::vector< std::set< subdomain_id_type >> & vars_active_subdomains ) const

inherited

returns a "tripled" copy of vars_active_subdomains, which is necessary in the complex-valued case.

Definition at line 3206 of file exodusII_io_helper.C.

 {
   std::vector<std::set<subdomain_id_type>> complex_vars_active_subdomains;
  
   for (auto & s : vars_active_subdomains)
     {
       // Push back the same data three times to match the tripling of the variables
       // for the real, imag, and modulus for the complex-valued solution.
       complex_vars_active_subdomains.push_back(s);
       complex_vars_active_subdomains.push_back(s);
       complex_vars_active_subdomains.push_back(s);
     }
  
   return complex_vars_active_subdomains;
 }

◆ get_conversion() [1/2]

const ExodusII_IO_Helper::Conversion & libMesh::ExodusII_IO_Helper::get_conversion ( const ElemType type ) const

inherited

Definition at line 383 of file exodusII_io_helper.C.

 {
   return libmesh_map_find(conversion_map, type);
 }

References libMesh::ExodusII_IO_Helper::conversion_map.

Referenced by build_element_and_node_maps(), compute_internal_and_border_elems_and_internal_nodes(), libMesh::ExodusII_IO_Helper::get_conversion(), libMesh::ExodusII_IO_Helper::read_edge_blocks(), libMesh::ExodusII_IO_Helper::read_sideset_data(), write_elements(), libMesh::ExodusII_IO_Helper::write_elements(), libMesh::ExodusII_IO_Helper::write_sideset_data(), write_sidesets(), and libMesh::ExodusII_IO_Helper::write_sidesets().

◆ get_conversion() [2/2]

const ExodusII_IO_Helper::Conversion & libMesh::ExodusII_IO_Helper::get_conversion ( std::string type_str ) const

inherited

Definition at line 389 of file exodusII_io_helper.C.

 {
   // Do only upper-case comparisons
   std::transform(type_str.begin(), type_str.end(), type_str.begin(), ::toupper);
   return get_conversion (libmesh_map_find(element_equivalence_map, type_str));
 }

References libMesh::ExodusII_IO_Helper::element_equivalence_map, and libMesh::ExodusII_IO_Helper::get_conversion().

◆ get_eb_info_global()

void libMesh::Nemesis_IO_Helper::get_eb_info_global ( )

Definition at line 179 of file nemesis_io_helper.C.

 {
   global_elem_blk_ids.resize(num_elem_blks_global);
   global_elem_blk_cnts.resize(num_elem_blks_global);
  
   if (num_elem_blks_global > 0)
     {
       nemesis_err_flag =
         Nemesis::ne_get_eb_info_global(ex_id,
                                        global_elem_blk_ids.data(),
                                        global_elem_blk_cnts.data());
       EX_CHECK_ERR(nemesis_err_flag, "Error reading global element block info!");
     }
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] " << "Global Element Block IDs and Counts:" << std::endl;
       for (auto bn : index_range(global_elem_blk_ids))
         {
           libMesh::out << "  [" << this->processor_id() << "] "
                        << "global_elem_blk_ids[" << bn << "]=" << global_elem_blk_ids[bn]
                        << ", global_elem_blk_cnts[" << bn << "]=" << global_elem_blk_cnts[bn]
                        << std::endl;
         }
     }
 }

References libMesh::ExodusII_IO_Helper::ex_id, global_elem_blk_cnts, global_elem_blk_ids, libMesh::index_range(), nemesis_err_flag, num_elem_blks_global, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ get_elem_cmap()

void libMesh::Nemesis_IO_Helper::get_elem_cmap ( )

Definition at line 410 of file nemesis_io_helper.C.

 {
   elem_cmap_elem_ids.resize(num_elem_cmaps);
   elem_cmap_side_ids.resize(num_elem_cmaps);
   elem_cmap_proc_ids.resize(num_elem_cmaps);
  
   for (auto i : index_range(elem_cmap_elem_ids))
     {
       elem_cmap_elem_ids[i].resize(elem_cmap_elem_cnts[i]);
       elem_cmap_side_ids[i].resize(elem_cmap_elem_cnts[i]);
       elem_cmap_proc_ids[i].resize(elem_cmap_elem_cnts[i]);
  
       if (elem_cmap_elem_cnts[i] > 0)
         {
           nemesis_err_flag =
             Nemesis::ne_get_elem_cmap(ex_id,
                                       elem_cmap_ids[i],
                                       elem_cmap_elem_ids[i].data(),
                                       elem_cmap_side_ids[i].data(),
                                       elem_cmap_proc_ids[i].data(),
                                       this->processor_id());
           EX_CHECK_ERR(nemesis_err_flag, "Error reading elem cmap elem, side, and processor ids!");
         }
  
       if (verbose)
         {
           libMesh::out << "[" << this->processor_id() << "] elem_cmap_elem_ids[" << i << "]=";
           for (const auto & dof : elem_cmap_elem_ids[i])
             libMesh::out << dof << " ";
           libMesh::out << std::endl;
  
           // These must be the (local) side IDs (in the ExodusII face numbering scheme)
           // of the sides shared across processors.
           libMesh::out << "[" << this->processor_id() << "] elem_cmap_side_ids[" << i << "]=";
           for (const auto & dof : elem_cmap_side_ids[i])
             libMesh::out << dof << " ";
           libMesh::out << std::endl;
  
           // This is basically a vector, all entries of which are = elem_cmap_ids[i]
           // Not sure if it's always guaranteed to be that or what...
           libMesh::out << "[" << this->processor_id() << "] elem_cmap_proc_ids[" << i << "]=";
           for (const auto & dof : elem_cmap_proc_ids[i])
             libMesh::out << dof << " ";
           libMesh::out << std::endl;
         }
     }
 }

References data, elem_cmap_elem_cnts, elem_cmap_elem_ids, elem_cmap_ids, elem_cmap_proc_ids, elem_cmap_side_ids, libMesh::ExodusII_IO_Helper::ex_id, libMesh::index_range(), nemesis_err_flag, num_elem_cmaps, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ get_elem_map()

void libMesh::Nemesis_IO_Helper::get_elem_map ( )

Definition at line 257 of file nemesis_io_helper.C.

 {
   elem_mapi.resize(num_internal_elems);
   elem_mapb.resize(num_border_elems);
  
   nemesis_err_flag =
     Nemesis::ne_get_elem_map(ex_id,
                              elem_mapi.empty() ? nullptr : elem_mapi.data(),
                              elem_mapb.empty() ? nullptr : elem_mapb.data(),
                              this->processor_id()
                              );
   EX_CHECK_ERR(nemesis_err_flag, "Error reading element maps!");
  
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] elem_mapi[i] = ";
       for (auto i : IntRange<int>(0, num_internal_elems-1))
         libMesh::out << elem_mapi[i] << ", ";
       libMesh::out << "... " << elem_mapi.back() << std::endl;
  
       libMesh::out << "[" << this->processor_id() << "] elem_mapb[i] = ";
       for (auto i : IntRange<int>(0, std::min(10, num_border_elems-1)))
         libMesh::out << elem_mapb[i] << ", ";
       libMesh::out << "... " << elem_mapb.back() << std::endl;
     }
 }

References elem_mapb, elem_mapi, libMesh::ExodusII_IO_Helper::ex_id, nemesis_err_flag, num_border_elems, num_internal_elems, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ get_elem_type()

const char * libMesh::ExodusII_IO_Helper::get_elem_type ( ) const

inherited

Returns: The current element type.

Note: The default behavior is for this value to be in all capital letters, e.g. HEX27.

Definition at line 396 of file exodusII_io_helper.C.

 {
   return elem_type.data();
 }

References libMesh::ExodusII_IO_Helper::elem_type.

◆ get_init_global()

void libMesh::Nemesis_IO_Helper::get_init_global ( )

Reading functions.

These just allocate memory for you and call the Nemesis routines of the same name. They also handle error checking for the Nemesis return value. Be careful calling these at random, some depend on others being called first... Fills: num_nodes_global, num_elems_global, num_elem_blks_global, num_node_sets_global, num_side_sets_global Call after: read_header() Call before: Any other get_* function from this class

Definition at line 88 of file nemesis_io_helper.C.

 {
   nemesis_err_flag =
     Nemesis::ne_get_init_global(ex_id,
                                 &num_nodes_global,
                                 &num_elems_global,
                                 &num_elem_blks_global,
                                 &num_node_sets_global,
                                 &num_side_sets_global);
   EX_CHECK_ERR(nemesis_err_flag, "Error reading initial global data!");
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] " << "num_nodes_global=" << num_nodes_global << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "num_elems_global=" << num_elems_global << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "num_elem_blks_global=" << num_elem_blks_global << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "num_node_sets_global=" << num_node_sets_global << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "num_side_sets_global=" << num_side_sets_global << std::endl;
     }
 }

References libMesh::ExodusII_IO_Helper::ex_id, nemesis_err_flag, num_elem_blks_global, num_elems_global, num_node_sets_global, num_nodes_global, num_side_sets_global, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ get_init_info()

void libMesh::Nemesis_IO_Helper::get_init_info ( )

Definition at line 208 of file nemesis_io_helper.C.

 {
   nemesis_err_flag =
     Nemesis::ne_get_init_info(ex_id,
                               &num_proc,
                               &num_proc_in_file,
                               &ftype);
   EX_CHECK_ERR(nemesis_err_flag, "Error reading initial info!");
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] " << "num_proc=" << num_proc << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "num_proc_in_file=" << num_proc_in_file << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "ftype=" << ftype << std::endl;
     }
 }

References libMesh::ExodusII_IO_Helper::ex_id, ftype, nemesis_err_flag, num_proc, num_proc_in_file, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ get_loadbal_param()

void libMesh::Nemesis_IO_Helper::get_loadbal_param ( )

Definition at line 227 of file nemesis_io_helper.C.

 {
   nemesis_err_flag =
     Nemesis::ne_get_loadbal_param(ex_id,
                                   &num_internal_nodes,
                                   &num_border_nodes,
                                   &num_external_nodes,
                                   &num_internal_elems,
                                   &num_border_elems,
                                   &num_node_cmaps,
                                   &num_elem_cmaps,
                                   this->processor_id() // The ID of the processor for which info is to be read
                                   );
   EX_CHECK_ERR(nemesis_err_flag, "Error reading load balance parameters!");
  
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] " << "num_internal_nodes=" << num_internal_nodes << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "num_border_nodes=" << num_border_nodes << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "num_external_nodes=" << num_external_nodes << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "num_internal_elems=" << num_internal_elems << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "num_border_elems=" << num_border_elems << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "num_node_cmaps=" << num_node_cmaps << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "num_elem_cmaps=" << num_elem_cmaps << std::endl;
     }
 }

References libMesh::ExodusII_IO_Helper::ex_id, nemesis_err_flag, num_border_elems, num_border_nodes, num_elem_cmaps, num_external_nodes, num_internal_elems, num_internal_nodes, num_node_cmaps, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ get_node_cmap()

void libMesh::Nemesis_IO_Helper::get_node_cmap ( )

Definition at line 368 of file nemesis_io_helper.C.

 {
   node_cmap_node_ids.resize(num_node_cmaps);
   node_cmap_proc_ids.resize(num_node_cmaps);
  
   for (auto i : index_range(node_cmap_node_ids))
     {
       node_cmap_node_ids[i].resize(node_cmap_node_cnts[i]);
       node_cmap_proc_ids[i].resize(node_cmap_node_cnts[i]);
  
       // Don't call ne_get_node_cmap() if there is nothing there to
       // get, Nemesis throws an error in this case.
       if (node_cmap_node_cnts[i] > 0)
         {
           nemesis_err_flag =
             Nemesis::ne_get_node_cmap(ex_id,
                                       node_cmap_ids[i],
                                       node_cmap_node_ids[i].data(),
                                       node_cmap_proc_ids[i].data(),
                                       this->processor_id());
           EX_CHECK_ERR(nemesis_err_flag, "Error reading node cmap node and processor ids!");
         }
  
       if (verbose)
         {
           libMesh::out << "[" << this->processor_id() << "] node_cmap_node_ids[" << i << "]=";
           for (const auto & dof : node_cmap_node_ids[i])
             libMesh::out << dof << " ";
           libMesh::out << std::endl;
  
           // This is basically a vector, all entries of which are = node_cmap_ids[i]
           // Not sure if it's always guaranteed to be that or what...
           libMesh::out << "[" << this->processor_id() << "] node_cmap_proc_ids[" << i << "]=";
           for (const auto & dof : node_cmap_proc_ids[i])
             libMesh::out << dof << " ";
           libMesh::out << std::endl;
         }
     }
 }

References data, libMesh::ExodusII_IO_Helper::ex_id, libMesh::index_range(), nemesis_err_flag, node_cmap_ids, node_cmap_node_cnts, node_cmap_node_ids, node_cmap_proc_ids, num_node_cmaps, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ get_node_map()

void libMesh::Nemesis_IO_Helper::get_node_map ( )

Definition at line 288 of file nemesis_io_helper.C.

 {
   node_mapi.resize(num_internal_nodes);
   node_mapb.resize(num_border_nodes);
   node_mape.resize(num_external_nodes);
  
   nemesis_err_flag =
     Nemesis::ne_get_node_map(ex_id,
                              node_mapi.empty() ? nullptr : node_mapi.data(),
                              node_mapb.empty() ? nullptr : node_mapb.data(),
                              node_mape.empty() ? nullptr : node_mape.data(),
                              this->processor_id()
                              );
   EX_CHECK_ERR(nemesis_err_flag, "Error reading node maps!");
  
   if (verbose)
     {
       // Remark: The Exodus/Nemesis node numbering is always (?) 1-based!  This means the first interior node id will
       // always be == 1.
       libMesh::out << "[" << this->processor_id() << "] " << "first interior node id=" << node_mapi[0] << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "last interior node id=" << node_mapi.back() << std::endl;
  
       libMesh::out << "[" << this->processor_id() << "] " << "first boundary node id=" << node_mapb[0] << std::endl;
       libMesh::out << "[" << this->processor_id() << "] " << "last boundary node id=" << node_mapb.back() << std::endl;
  
       // The number of external nodes is sometimes zero, don't try to access
       // node_mape.back() in this case!
       if (num_external_nodes > 0)
         {
           libMesh::out << "[" << this->processor_id() << "] " << "first external node id=" << node_mape[0] << std::endl;
           libMesh::out << "[" << this->processor_id() << "] " << "last external node id=" << node_mape.back() << std::endl;
         }
     }
 }

References libMesh::ExodusII_IO_Helper::ex_id, nemesis_err_flag, node_mapb, node_mape, node_mapi, num_border_nodes, num_external_nodes, num_internal_nodes, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ get_node_set_id()

int libMesh::ExodusII_IO_Helper::get_node_set_id ( int index )

inherited

Get the node set id for the given node set index.

Definition at line 778 of file exodusII_io_helper.C.

 {
   libmesh_assert_less (index, nodeset_ids.size());
  
   return nodeset_ids[index];
 }

References libMesh::ExodusII_IO_Helper::nodeset_ids.

◆ get_node_set_name()

std::string libMesh::ExodusII_IO_Helper::get_node_set_name ( int index )

inherited

Get the node set name for the given node set index if supplied in the mesh file.

Otherwise an empty string is returned.

Definition at line 787 of file exodusII_io_helper.C.

 {
   libmesh_assert_less (index, nodeset_ids.size());
  
   return id_to_ns_names[nodeset_ids[index]];
 }

References libMesh::ExodusII_IO_Helper::id_to_ns_names, and libMesh::ExodusII_IO_Helper::nodeset_ids.

◆ get_ns_param_global()

void libMesh::Nemesis_IO_Helper::get_ns_param_global ( )

Definition at line 147 of file nemesis_io_helper.C.

 {
   if (num_node_sets_global > 0)
     {
       global_nodeset_ids.resize(num_node_sets_global);
       num_global_node_counts.resize(num_node_sets_global);
       num_global_node_df_counts.resize(num_node_sets_global);
  
       nemesis_err_flag =
         Nemesis::ne_get_ns_param_global(ex_id,
                                         global_nodeset_ids.data(),
                                         num_global_node_counts.data(),
                                         num_global_node_df_counts.data());
       EX_CHECK_ERR(nemesis_err_flag, "Error reading global nodeset parameters!");
  
       if (verbose)
         {
           libMesh::out << "[" << this->processor_id() << "] " << "Global Nodeset IDs, Node Counts, and DF counts:" << std::endl;
           for (auto bn : index_range(global_nodeset_ids))
             {
               libMesh::out << "  [" << this->processor_id() << "] "
                            << "global_nodeset_ids["<<bn<<"]=" << global_nodeset_ids[bn]
                            << ", num_global_node_counts["<<bn<<"]=" << num_global_node_counts[bn]
                            << ", num_global_node_df_counts["<<bn<<"]=" << num_global_node_df_counts[bn]
                            << std::endl;
             }
         }
     }
 }

References libMesh::ExodusII_IO_Helper::ex_id, global_nodeset_ids, libMesh::index_range(), nemesis_err_flag, num_global_node_counts, num_global_node_df_counts, num_node_sets_global, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ get_side_set_id()

int libMesh::ExodusII_IO_Helper::get_side_set_id ( int index )

inherited

Get the side set id for the given side set index.

Definition at line 760 of file exodusII_io_helper.C.

 {
   libmesh_assert_less (index, ss_ids.size());
  
   return ss_ids[index];
 }

References libMesh::ExodusII_IO_Helper::ss_ids.

◆ get_side_set_name()

std::string libMesh::ExodusII_IO_Helper::get_side_set_name ( int index )

inherited

Get the side set name for the given side set index if supplied in the mesh file.

Otherwise an empty string is returned.

Definition at line 769 of file exodusII_io_helper.C.

 {
   libmesh_assert_less (index, ss_ids.size());
  
   return id_to_ss_names[ss_ids[index]];
 }

References libMesh::ExodusII_IO_Helper::id_to_ss_names, and libMesh::ExodusII_IO_Helper::ss_ids.

◆ get_ss_param_global()

void libMesh::Nemesis_IO_Helper::get_ss_param_global ( )

Fills: global_sideset_ids, num_global_side_counts, num_global_side_df_counts Call after: get_init_global()

Definition at line 111 of file nemesis_io_helper.C.

 {
   if (num_side_sets_global > 0)
     {
       global_sideset_ids.resize(num_side_sets_global);
       num_global_side_counts.resize(num_side_sets_global);
  
       // df = "distribution factor", not really sure what that is.  I don't yet have a file
       // which has distribution factors so I guess we'll worry about it later...
       num_global_side_df_counts.resize(num_side_sets_global);
  
       nemesis_err_flag =
         Nemesis::ne_get_ss_param_global(ex_id,
                                         global_sideset_ids.data(),
                                         num_global_side_counts.data(),
                                         num_global_side_df_counts.data());
       EX_CHECK_ERR(nemesis_err_flag, "Error reading global sideset parameters!");
  
       if (verbose)
         {
           libMesh::out << "[" << this->processor_id() << "] " << "Global Sideset IDs, Side Counts, and DF counts:" << std::endl;
           for (auto bn : index_range(global_sideset_ids))
             {
               libMesh::out << "  [" << this->processor_id() << "] "
                            << "global_sideset_ids["<<bn<<"]=" << global_sideset_ids[bn]
                            << ", num_global_side_counts["<<bn<<"]=" << num_global_side_counts[bn]
                            << ", num_global_side_df_counts["<<bn<<"]=" << num_global_side_df_counts[bn]
                            << std::endl;
             }
         }
     }
 }

References libMesh::ExodusII_IO_Helper::ex_id, global_sideset_ids, libMesh::index_range(), nemesis_err_flag, num_global_side_counts, num_global_side_df_counts, num_side_sets_global, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ init_conversion_map()

void libMesh::ExodusII_IO_Helper::init_conversion_map ( )

privateinherited

Definition at line 146 of file exodusII_io_helper.C.

 {
   {
     auto & conv = conversion_map[NODEELEM];
     conv.libmesh_type = NODEELEM;
     conv.exodus_type = "SPHERE";
   }
   {
     auto & conv = conversion_map[EDGE2];
     conv.libmesh_type = EDGE2;
     conv.exodus_type = "EDGE2";
   }
   {
     auto & conv = conversion_map[EDGE3];
     conv.libmesh_type = EDGE3;
     conv.exodus_type = "EDGE3";
   }
   {
     auto & conv = conversion_map[QUAD4];
     conv.libmesh_type = QUAD4;
     conv.exodus_type = "QUAD4";
   }
   {
     auto & conv = conversion_map[QUADSHELL4];
     conv.inverse_side_map = &quadshell4_inverse_edge_map;
     conv.shellface_index_offset = 2;
     conv.libmesh_type = QUADSHELL4;
     conv.exodus_type = "SHELL4";
   }
   {
     auto & conv = conversion_map[QUAD8];
     conv.libmesh_type = QUAD8;
     conv.exodus_type = "QUAD8";
   }
   {
     auto & conv = conversion_map[QUADSHELL8];
     conv.inverse_side_map = &quadshell4_inverse_edge_map;
     conv.shellface_index_offset = 2;
     conv.libmesh_type = QUADSHELL8;
     conv.exodus_type = "SHELL8";
   }
   {
     auto & conv = conversion_map[QUAD9];
     conv.libmesh_type = QUAD9;
     conv.exodus_type = "QUAD9";
   }
   {
     auto & conv = conversion_map[TRI3];
     conv.libmesh_type = TRI3;
     conv.exodus_type = "TRI3";
   }
   {
     auto & conv = conversion_map[TRISHELL3];
     conv.inverse_side_map = &trishell3_inverse_edge_map;
     conv.shellface_index_offset = 2;
     conv.libmesh_type = TRISHELL3;
     conv.exodus_type = "TRISHELL3";
   }
   {
     auto & conv = conversion_map[TRI3SUBDIVISION];
     conv.libmesh_type = TRI3SUBDIVISION;
     conv.exodus_type = "TRI3";
   }
   {
     auto & conv = conversion_map[TRI6];
     conv.libmesh_type = TRI6;
     conv.exodus_type = "TRI6";
   }
   {
     auto & conv = conversion_map[HEX8];
     conv.side_map = &hex_face_map;
     conv.inverse_side_map = &hex_inverse_face_map;
     conv.libmesh_type = HEX8;
     conv.exodus_type = "HEX8";
   }
   {
     auto & conv = conversion_map[HEX20];
     conv.side_map = &hex_face_map;
     conv.inverse_side_map = &hex_inverse_face_map;
     conv.libmesh_type = HEX20;
     conv.exodus_type = "HEX20";
   }
   {
     auto & conv = conversion_map[HEX27];
     conv.node_map = &hex27_node_map;
     conv.inverse_node_map = &hex27_inverse_node_map;
     conv.side_map = &hex_face_map;
     conv.inverse_side_map = &hex_inverse_face_map;
     conv.libmesh_type = HEX27;
     conv.exodus_type = "HEX27";
   }
   {
     auto & conv = conversion_map[TET4];
     conv.side_map = &tet_face_map;
     conv.inverse_side_map = &tet_inverse_face_map;
     conv.libmesh_type = TET4;
     conv.exodus_type = "TETRA4";
   }
   {
     auto & conv = conversion_map[TET10];
     conv.side_map = &tet_face_map;
     conv.inverse_side_map = &tet_inverse_face_map;
     conv.libmesh_type = TET10;
     conv.exodus_type = "TETRA10";
   }
   {
     auto & conv = conversion_map[PRISM6];
     conv.side_map = &prism_face_map;
     conv.inverse_side_map = &prism_inverse_face_map;
     conv.libmesh_type = PRISM6;
     conv.exodus_type = "WEDGE";
   }
   {
     auto & conv = conversion_map[PRISM15];
     conv.side_map = &prism_face_map;
     conv.inverse_side_map = &prism_inverse_face_map;
     conv.libmesh_type = PRISM15;
     conv.exodus_type = "WEDGE15";
   }
   {
     auto & conv = conversion_map[PRISM18];
     conv.side_map = &prism_face_map;
     conv.inverse_side_map = &prism_inverse_face_map;
     conv.libmesh_type = PRISM18;
     conv.exodus_type = "WEDGE18";
   }
   {
     auto & conv = conversion_map[PYRAMID5];
     conv.libmesh_type = PYRAMID5;
     conv.exodus_type = "PYRAMID5";
   }
   {
     auto & conv = conversion_map[PYRAMID13];
     conv.libmesh_type = PYRAMID13;
     conv.exodus_type = "PYRAMID13";
   }
   {
     auto & conv = conversion_map[PYRAMID14];
     conv.libmesh_type = PYRAMID14;
     conv.exodus_type = "PYRAMID14";
   }
 }

References libMesh::ExodusII_IO_Helper::conversion_map, libMesh::EDGE2, libMesh::EDGE3, libMesh::HEX20, libMesh::HEX27, libMesh::HEX8, libMesh::NODEELEM, libMesh::PRISM15, libMesh::PRISM18, libMesh::PRISM6, libMesh::PYRAMID13, libMesh::PYRAMID14, libMesh::PYRAMID5, libMesh::QUAD4, libMesh::QUAD8, libMesh::QUAD9, libMesh::QUADSHELL4, libMesh::QUADSHELL8, libMesh::TET10, libMesh::TET4, libMesh::TRI3, libMesh::TRI3SUBDIVISION, libMesh::TRI6, and libMesh::TRISHELL3.

Referenced by libMesh::ExodusII_IO_Helper::ExodusII_IO_Helper().

◆ init_element_equivalence_map()

void libMesh::ExodusII_IO_Helper::init_element_equivalence_map ( )

privateinherited

Definition at line 293 of file exodusII_io_helper.C.

 {
   // We use an ExodusII SPHERE element to represent a NodeElem
   element_equivalence_map["SPHERE"] = NODEELEM;
  
   // EDGE2 equivalences
   element_equivalence_map["EDGE"]   = EDGE2;
   element_equivalence_map["EDGE2"]  = EDGE2;
   element_equivalence_map["TRUSS"]  = EDGE2;
   element_equivalence_map["BEAM"]   = EDGE2;
   element_equivalence_map["BAR"]    = EDGE2;
   element_equivalence_map["TRUSS2"] = EDGE2;
   element_equivalence_map["BEAM2"]  = EDGE2;
   element_equivalence_map["BAR2"]   = EDGE2;
  
   // EDGE3 equivalences
   element_equivalence_map["EDGE3"]  = EDGE3;
   element_equivalence_map["TRUSS3"] = EDGE3;
   element_equivalence_map["BEAM3"]  = EDGE3;
   element_equivalence_map["BAR3"]   = EDGE3;
  
   // QUAD4 equivalences
   element_equivalence_map["QUAD"]   = QUAD4;
   element_equivalence_map["QUAD4"]  = QUAD4;
  
   // QUADSHELL4 equivalences
   element_equivalence_map["SHELL"]  = QUADSHELL4;
   element_equivalence_map["SHELL4"] = QUADSHELL4;
  
   // QUAD8 equivalences
   element_equivalence_map["QUAD8"]  = QUAD8;
  
   // QUADSHELL8 equivalences
   element_equivalence_map["SHELL8"] = QUADSHELL8;
  
   // QUAD9 equivalences
   element_equivalence_map["QUAD9"]  = QUAD9;
   // element_equivalence_map["SHELL9"] = QUAD9;
  
   // TRI3 equivalences
   element_equivalence_map["TRI"]       = TRI3;
   element_equivalence_map["TRI3"]      = TRI3;
   element_equivalence_map["TRIANGLE"]  = TRI3;
  
   // TRISHELL3 equivalences
   element_equivalence_map["TRISHELL"]  = TRISHELL3;
   element_equivalence_map["TRISHELL3"] = TRISHELL3;
  
   // TRI6 equivalences
   element_equivalence_map["TRI6"]      = TRI6;
   // element_equivalence_map["TRISHELL6"] = TRI6;
  
   // HEX8 equivalences
   element_equivalence_map["HEX"]  = HEX8;
   element_equivalence_map["HEX8"] = HEX8;
  
   // HEX20 equivalences
   element_equivalence_map["HEX20"] = HEX20;
  
   // HEX27 equivalences
   element_equivalence_map["HEX27"] = HEX27;
  
   // TET4 equivalences
   element_equivalence_map["TETRA"]  = TET4;
   element_equivalence_map["TETRA4"] = TET4;
  
   // TET10 equivalences
   element_equivalence_map["TETRA10"] = TET10;
  
   // PRISM6 equivalences
   element_equivalence_map["WEDGE"] = PRISM6;
  
   // PRISM15 equivalences
   element_equivalence_map["WEDGE15"] = PRISM15;
  
   // PRISM18 equivalences
   element_equivalence_map["WEDGE18"] = PRISM18;
  
   // PYRAMID5 equivalences
   element_equivalence_map["PYRAMID"]  = PYRAMID5;
   element_equivalence_map["PYRAMID5"] = PYRAMID5;
  
   // PYRAMID13 equivalences
   element_equivalence_map["PYRAMID13"] = PYRAMID13;
  
   // PYRAMID14 equivalences
   element_equivalence_map["PYRAMID14"] = PYRAMID14;
 }

References libMesh::EDGE2, libMesh::EDGE3, libMesh::ExodusII_IO_Helper::element_equivalence_map, libMesh::HEX20, libMesh::HEX27, libMesh::HEX8, libMesh::NODEELEM, libMesh::PRISM15, libMesh::PRISM18, libMesh::PRISM6, libMesh::PYRAMID13, libMesh::PYRAMID14, libMesh::PYRAMID5, libMesh::QUAD4, libMesh::QUAD8, libMesh::QUAD9, libMesh::QUADSHELL4, libMesh::QUADSHELL8, libMesh::TET10, libMesh::TET4, libMesh::TRI3, libMesh::TRI6, and libMesh::TRISHELL3.

Referenced by libMesh::ExodusII_IO_Helper::ExodusII_IO_Helper().

◆ initialize()

void libMesh::Nemesis_IO_Helper::initialize	(	std::string	title,
		const MeshBase &	mesh,
		bool	use_discontinuous = `false`
	)

overridevirtual

Specialization of the initialize function from ExodusII_IO_Helper that also writes global initial data to file.

Reimplemented from libMesh::ExodusII_IO_Helper.

Definition at line 766 of file nemesis_io_helper.C.

 {
   // Make sure that the reference passed in is really a DistributedMesh
   // const DistributedMesh & pmesh = cast_ref<const DistributedMesh &>(mesh);
   const MeshBase & pmesh = mesh;
  
   // According to Nemesis documentation, first call when writing should be to
   // ne_put_init_info().  Our reader doesn't actually call this, but we should
   // strive to be as close to a normal nemesis file as possible...
   this->put_init_info(this->n_processors(), 1, "p");
  
  
   // Gather global "initial" information for Nemesis.  This consists of
   // three parts labeled I, II, and III below...
  
   //
   // I.) Need to compute the number of global element blocks.  To be consistent with
   // Exodus, we also incorrectly associate the number of element blocks with the
   // number of libmesh subdomains...
   //
   this->compute_num_global_elem_blocks(pmesh);
  
   //
   // II.) Determine the global number of nodesets by communication.
   // This code relies on BoundaryInfo storing side and node
   // boundary IDs separately at the time they are added to the
   // BoundaryInfo object.
   //
   this->compute_num_global_nodesets(pmesh);
  
   //
   // III.) Need to compute the global number of sidesets by communication:
   // This code relies on BoundaryInfo storing side and node
   // boundary IDs separately at the time they are added to the
   // BoundaryInfo object.
   //
   this->compute_num_global_sidesets(pmesh);
  
   // Now write the global data obtained in steps I, II, and III to the Nemesis file
   this->put_init_global(pmesh.parallel_n_nodes(),
                         pmesh.parallel_n_elem(),
                         this->num_elem_blks_global, /* I.   */
                         this->num_node_sets_global, /* II.  */
                         this->num_side_sets_global  /* III. */
                         );
  
   // Next, we'll write global element block information to the file.  This was already
   // gathered in step I. above
   this->put_eb_info_global(this->global_elem_blk_ids,
                            this->global_elem_blk_cnts);
  
  
   // Next, write global nodeset information to the file.  This was already gathered in
   // step II. above.
   this->num_global_node_df_counts.clear();
   this->num_global_node_df_counts.resize(this->global_nodeset_ids.size()); // distribution factors all zero...
   this->put_ns_param_global(this->global_nodeset_ids,
                             this->num_global_node_counts,
                             this->num_global_node_df_counts);
  
  
   // Next, write global sideset information to the file.  This was already gathered in
   // step III. above.
   this->num_global_side_df_counts.clear();
   this->num_global_side_df_counts.resize(this->global_sideset_ids.size()); // distribution factors all zero...
   this->put_ss_param_global(this->global_sideset_ids,
                             this->num_global_side_counts,
                             this->num_global_side_df_counts);
  
  
   // Before we go any further we need to derive consistent node and
   // element numbering schemes for all local elems and nodes connected
   // to local elements.
   //
   // Must be called *after* the local_subdomain_counts map has been constructed
   // by the compute_num_global_elem_blocks() function!
   this->build_element_and_node_maps(pmesh);
  
   // Next step is to write "load balance" parameters.  Several things need to
   // be computed first though...
  
   // First we'll collect IDs of border nodes.
   this->compute_border_node_ids(pmesh);
  
   // Next we'll collect numbers of internal and border elements, and internal nodes.
   // Note: "A border node does not a border element make...", that is, just because one
   // of an element's nodes has been identified as a border node, the element is not
   // necessarily a border element.  It must have a side on the boundary between processors,
   // i.e. have a face neighbor with a different processor id...
   this->compute_internal_and_border_elems_and_internal_nodes(pmesh);
  
   // Finally we are ready to write the loadbal information to the file
   this->put_loadbal_param(this->num_internal_nodes,
                           this->num_border_nodes,
                           this->num_external_nodes,
                           this->num_internal_elems,
                           this->num_border_elems,
                           this->num_node_cmaps,
                           this->num_elem_cmaps);
  
  
   // Now we need to compute the "communication map" parameters.  These are basically
   // lists of nodes and elements which need to be communicated between different processors
   // when the mesh file is read back in.
   this->compute_communication_map_parameters();
  
   // Do we have communication maps to write?  Note that
   // ne_put_cmap_params expects us to have either *both* node and elem
   // cmaps or *neither*
   if (!this->node_cmap_ids.empty() &&
       !this->node_cmap_node_cnts.empty() &&
       !this->elem_cmap_ids.empty() &&
       !this->elem_cmap_elem_cnts.empty())
     {
       // Write communication map parameters to file.
       this->put_cmap_params(this->node_cmap_ids,
                             this->node_cmap_node_cnts,
                             this->elem_cmap_ids,
                             this->elem_cmap_elem_cnts);
  
       // Ready the node communication maps.  The node IDs which
       // are communicated are the ones currently stored in
       // proc_nodes_touched_intersections.
       this->compute_node_communication_maps();
  
       // Write the packed node communication vectors to file.
       this->put_node_cmap(this->node_cmap_node_ids,
                           this->node_cmap_proc_ids);
  
       // Ready the node maps.  These have nothing to do with communication, they map
       // the nodes to internal, border, and external nodes in the file.
       this->compute_node_maps();
  
       // Call the Nemesis API to write the node maps to file.
       this->put_node_map(this->node_mapi,
                          this->node_mapb,
                          this->node_mape);
  
       // Ready the element communication maps.  This includes border
       // element IDs, sides which are on the border, and the processors to which
       // they are to be communicated...
       this->compute_elem_communication_maps();
  
       // Call the Nemesis API to write the packed element communication maps vectors to file
       this->put_elem_cmap(this->elem_cmap_elem_ids,
                           this->elem_cmap_side_ids,
                           this->elem_cmap_proc_ids);
     }
  
  
   // Ready the Nemesis element maps (internal and border) for writing to file.
   this->compute_element_maps();
  
   // Call the Nemesis API to write the internal and border element IDs.
   this->put_elem_map(this->elem_mapi,
                      this->elem_mapb);
  
   // Now write Exodus-specific initialization information, some of which is
   // different when you are using Nemesis.
   this->write_exodus_initialization_info(pmesh, title_in);
 } // end initialize()

References build_element_and_node_maps(), compute_border_node_ids(), compute_communication_map_parameters(), compute_elem_communication_maps(), compute_element_maps(), compute_internal_and_border_elems_and_internal_nodes(), compute_node_communication_maps(), compute_node_maps(), compute_num_global_elem_blocks(), compute_num_global_nodesets(), compute_num_global_sidesets(), elem_cmap_elem_cnts, elem_cmap_elem_ids, elem_cmap_ids, elem_cmap_proc_ids, elem_cmap_side_ids, elem_mapb, elem_mapi, global_elem_blk_cnts, global_elem_blk_ids, global_nodeset_ids, global_sideset_ids, mesh, libMesh::ParallelObject::n_processors(), node_cmap_ids, node_cmap_node_cnts, node_cmap_node_ids, node_cmap_proc_ids, node_mapb, node_mape, node_mapi, num_border_elems, num_border_nodes, num_elem_blks_global, num_elem_cmaps, num_external_nodes, num_global_node_counts, num_global_node_df_counts, num_global_side_counts, num_global_side_df_counts, num_internal_elems, num_internal_nodes, num_node_cmaps, num_node_sets_global, num_side_sets_global, libMesh::MeshBase::parallel_n_elem(), libMesh::MeshBase::parallel_n_nodes(), put_cmap_params(), put_eb_info_global(), put_elem_cmap(), put_elem_map(), put_init_global(), put_init_info(), put_loadbal_param(), put_node_cmap(), put_node_map(), put_ns_param_global(), put_ss_param_global(), and write_exodus_initialization_info().

◆ initialize_element_variables()

void libMesh::Nemesis_IO_Helper::initialize_element_variables	(	std::vector< std::string >	names,
		const std::vector< std::set< subdomain_id_type >> &	vars_active_subdomains
	)

overridevirtual

Override the Exodus Helper's implementation of this function so that it works correctly in parallel.

Reimplemented from libMesh::ExodusII_IO_Helper.

Definition at line 2548 of file nemesis_io_helper.C.

 {
   // Quick return if there are no element variables to write
   if (names.size() == 0)
     return;
  
   // Quick return if we have already called this function
   if (_elem_vars_initialized)
     return;
  
   // Be sure that variables in the file match what we are asking for
   if (num_elem_vars > 0)
     {
       this->check_existing_vars(ELEMENTAL, names, this->elem_var_names);
       return;
     }
  
   // Set the flag so we can skip this stuff on subsequent calls to
   // initialize_element_variables()
   _elem_vars_initialized = true;
  
   this->write_var_names(ELEMENTAL, names);
  
   // Create a truth table from global_elem_blk_ids and the information
   // in vars_active_subdomains. Create a truth table of
   // size global_elem_blk_ids.size() * names.size().
   std::vector<int> truth_tab(global_elem_blk_ids.size() * names.size());
   for (auto blk : index_range(global_elem_blk_ids))
     for (auto var : index_range(names))
       if (vars_active_subdomains[var].empty() ||
           vars_active_subdomains[var].count(cast_int<subdomain_id_type>(global_elem_blk_ids[blk])))
         truth_tab[names.size() * blk + var] = 1;
  
   // Write truth table to file.
   if (truth_tab.size())
     {
       ex_err = exII::ex_put_elem_var_tab(ex_id,
                                          cast_int<int>(global_elem_blk_ids.size()),
                                          cast_int<int>(names.size()),
                                          truth_tab.data());
       EX_CHECK_ERR(ex_err, "Error writing element truth table.");
     }
 }

References libMesh::ExodusII_IO_Helper::_elem_vars_initialized, libMesh::ExodusII_IO_Helper::check_existing_vars(), libMesh::ExodusII_IO_Helper::elem_var_names, libMesh::ExodusII_IO_Helper::ELEMENTAL, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, global_elem_blk_ids, libMesh::index_range(), libMesh::ExodusII_IO_Helper::num_elem_vars, and libMesh::ExodusII_IO_Helper::write_var_names().

◆ initialize_global_variables()

void libMesh::ExodusII_IO_Helper::initialize_global_variables ( std::vector< std::string > names )

inherited

Sets up the global variables.

Definition at line 2496 of file exodusII_io_helper.C.

 {
   if ((_run_only_on_proc0) && (this->processor_id() != 0))
     return;
  
   // Quick return if there are no global variables to write
   if (names.size() == 0)
     return;
  
   if (_global_vars_initialized)
     return;
  
   // Be sure that variables in the file match what we are asking for
   if (num_global_vars > 0)
     {
       this->check_existing_vars(GLOBAL, names, this->global_var_names);
       return;
     }
  
   _global_vars_initialized = true;
  
   this->write_var_names(GLOBAL, names);
 }

References libMesh::ExodusII_IO_Helper::_global_vars_initialized, libMesh::ExodusII_IO_Helper::_run_only_on_proc0, libMesh::ExodusII_IO_Helper::check_existing_vars(), libMesh::ExodusII_IO_Helper::GLOBAL, libMesh::ExodusII_IO_Helper::global_var_names, libMesh::ExodusII_IO_Helper::num_global_vars, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::write_var_names().

◆ initialize_nodal_variables()

void libMesh::ExodusII_IO_Helper::initialize_nodal_variables ( std::vector< std::string > names )

inherited

Sets up the nodal variables.

Definition at line 2468 of file exodusII_io_helper.C.

 {
   if ((_run_only_on_proc0) && (this->processor_id() != 0))
     return;
  
   // Quick return if there are no nodal variables to write
   if (names.size() == 0)
     return;
  
   // Quick return if we have already called this function
   if (_nodal_vars_initialized)
     return;
  
   // Be sure that variables in the file match what we are asking for
   if (num_nodal_vars > 0)
     {
       this->check_existing_vars(NODAL, names, this->nodal_var_names);
       return;
     }
  
   // Set the flag so we can skip the rest of this function on subsequent calls.
   _nodal_vars_initialized = true;
  
   this->write_var_names(NODAL, names);
 }

References libMesh::ExodusII_IO_Helper::_nodal_vars_initialized, libMesh::ExodusII_IO_Helper::_run_only_on_proc0, libMesh::ExodusII_IO_Helper::check_existing_vars(), libMesh::ExodusII_IO_Helper::NODAL, libMesh::ExodusII_IO_Helper::nodal_var_names, libMesh::ExodusII_IO_Helper::num_nodal_vars, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::write_var_names().

◆ inquire()

int libMesh::ExodusII_IO_Helper::inquire	(	int	req_info,
		std::string	error_msg = `""`
	)

inherited

Returns: The value obtained from a generic exII::ex_inquire() call.

Definition at line 1247 of file exodusII_io_helper.C.

 {
   int ret_int = 0;
   char ret_char = 0;
   float ret_float = 0.;
  
   ex_err = exII::ex_inquire(ex_id,
                             req_info_in,
                             &ret_int,
                             &ret_float,
                             &ret_char);
  
   EX_CHECK_ERR(ex_err, error_msg);
  
   return ret_int;
 }

References libMesh::ExodusII_IO_Helper::ex_err, and libMesh::ExodusII_IO_Helper::ex_id.

Referenced by libMesh::ExodusII_IO_Helper::read_all_nodesets(), libMesh::ExodusII_IO_Helper::read_num_time_steps(), libMesh::ExodusII_IO_Helper::read_qa_records(), libMesh::ExodusII_IO_Helper::read_sideset_info(), and libMesh::ExodusII_IO_Helper::write_information_records().

◆ message() [1/2]

void libMesh::ExodusII_IO_Helper::message ( const std::string & msg )

inherited

Prints the message defined in msg.

Can be turned off if verbosity is set to 0.

Definition at line 403 of file exodusII_io_helper.C.

 {
   if (verbose) libMesh::out << msg << std::endl;
 }

References libMesh::out, and libMesh::ExodusII_IO_Helper::verbose.

◆ message() [2/2]

void libMesh::ExodusII_IO_Helper::message	(	const std::string &	msg,
		int	i
	)

inherited

Prints the message defined in msg, and appends the number i to the end of the message.

Useful for printing messages in loops. Can be turned off if verbosity is set to 0.

Definition at line 410 of file exodusII_io_helper.C.

 {
   if (verbose) libMesh::out << msg << i << "." << std::endl;
 }

References libMesh::out, and libMesh::ExodusII_IO_Helper::verbose.

◆ n_processors()

processor_id_type libMesh::ParallelObject::n_processors ( ) const

inlineinherited

Returns: The number of processors in the group.

Definition at line 100 of file parallel_object.h.

101 { return cast_int<processor_id_type>(_communicator.size()); }

References libMesh::ParallelObject::_communicator.

◆ open()

void libMesh::ExodusII_IO_Helper::open	(	const char *	filename,
		bool	read_only
	)

inherited

Opens an ExodusII mesh file named filename.

If read_only==true, the file will be opened with the EX_READ flag, otherwise it will be opened with the EX_WRITE flag.

Definition at line 492 of file exodusII_io_helper.C.

 {
   // Version of Exodus you are using
   float ex_version = 0.;
  
   int comp_ws = 0;
  
   if (_single_precision)
     comp_ws = cast_int<int>(sizeof(float));
  
   // Fall back on double precision when necessary since ExodusII
   // doesn't seem to support long double
   else
     comp_ws = cast_int<int>(std::min(sizeof(Real), sizeof(double)));
  
   // Word size in bytes of the floating point data as they are stored
   // in the ExodusII file.  "If this argument is 0, the word size of the
   // floating point data already stored in the file is returned"
   int io_ws = 0;
  
   ex_id = exII::ex_open(filename,
                         read_only ? EX_READ : EX_WRITE,
                         &comp_ws,
                         &io_ws,
                         &ex_version);
  
   std::string err_msg = std::string("Error opening ExodusII mesh file: ") + std::string(filename);
   EX_CHECK_ERR(ex_id, err_msg);
   if (verbose) libMesh::out << "File opened successfully." << std::endl;
  
   if (read_only)
     opened_for_reading = true;
   else
     opened_for_writing = true;
  
   current_filename = std::string(filename);
 }

References libMesh::ExodusII_IO_Helper::_single_precision, libMesh::ExodusII_IO_Helper::current_filename, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::opened_for_reading, libMesh::ExodusII_IO_Helper::opened_for_writing, libMesh::out, libMesh::Real, and libMesh::ExodusII_IO_Helper::verbose.

◆ print_header()

void libMesh::ExodusII_IO_Helper::print_header ( )

inherited

Prints the ExodusII mesh file header, which includes the mesh title, the number of nodes, number of elements, mesh dimension, number of sidesets, and number of nodesets.

Definition at line 625 of file exodusII_io_helper.C.

 {
   if (verbose)
     libMesh::out << "Title: \t" << title.data() << std::endl
                  << "Mesh Dimension: \t"   << num_dim << std::endl
                  << "Number of Nodes: \t" << num_nodes << std::endl
                  << "Number of elements: \t" << num_elem << std::endl
                  << "Number of elt blocks: \t" << num_elem_blk << std::endl
                  << "Number of node sets: \t" << num_node_sets << std::endl
                  << "Number of side sets: \t" << num_side_sets << std::endl;
 }

References libMesh::ExodusII_IO_Helper::num_dim, libMesh::ExodusII_IO_Helper::num_elem, libMesh::ExodusII_IO_Helper::num_elem_blk, libMesh::ExodusII_IO_Helper::num_node_sets, libMesh::ExodusII_IO_Helper::num_nodes, libMesh::ExodusII_IO_Helper::num_side_sets, libMesh::out, libMesh::ExodusII_IO_Helper::title, and libMesh::ExodusII_IO_Helper::verbose.

◆ print_nodes()

void libMesh::ExodusII_IO_Helper::print_nodes ( std::ostream & out = libMesh::out )

inherited

Prints the nodal information, by default to libMesh::out.

Definition at line 684 of file exodusII_io_helper.C.

 {
   for (int i=0; i<num_nodes; i++)
     out_stream << "(" << x[i] << ", " << y[i] << ", " << z[i] << ")" << std::endl;
 }

References libMesh::ExodusII_IO_Helper::num_nodes, libMesh::ExodusII_IO_Helper::x, libMesh::ExodusII_IO_Helper::y, and libMesh::ExodusII_IO_Helper::z.

◆ processor_id()

processor_id_type libMesh::ParallelObject::processor_id ( ) const

inlineinherited

Returns: The rank of this processor in the group.

Definition at line 106 of file parallel_object.h.

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

References libMesh::ParallelObject::_communicator.

◆ put_cmap_params()

void libMesh::Nemesis_IO_Helper::put_cmap_params	(	std::vector< int > &	node_cmap_ids,
		std::vector< int > &	node_cmap_node_cnts,
		std::vector< int > &	elem_cmap_ids,
		std::vector< int > &	elem_cmap_elem_cnts
	)

Outputs initial information for communication maps.

Note: The order of the arguments specified in the Nemesis User's Manual is wrong. The correct order is (ids, counts, ids, counts). Must be called after put_loadbal_param().

Definition at line 576 of file nemesis_io_helper.C.

 {
   libmesh_assert(!node_cmap_ids_in.empty());
   libmesh_assert(!node_cmap_node_cnts_in.empty());
   libmesh_assert(!elem_cmap_ids_in.empty());
   libmesh_assert(!elem_cmap_elem_cnts_in.empty());
  
   nemesis_err_flag =
     Nemesis::ne_put_cmap_params(ex_id,
                                 node_cmap_ids_in.data(),
                                 node_cmap_node_cnts_in.data(),
                                 elem_cmap_ids_in.data(),
                                 elem_cmap_elem_cnts_in.data(),
                                 this->processor_id());
  
   EX_CHECK_ERR(nemesis_err_flag, "Error writing cmap parameters!");
 }

References libMesh::ExodusII_IO_Helper::ex_id, libMesh::libmesh_assert(), nemesis_err_flag, and libMesh::ParallelObject::processor_id().

Referenced by initialize().

◆ put_eb_info_global()

void libMesh::Nemesis_IO_Helper::put_eb_info_global	(	std::vector< int > &	global_elem_blk_ids,
		std::vector< int > &	global_elem_blk_cnts
	)

Writes global block information to the file .) global_elem_blk_ids - list of block IDs for all blocks present in the mesh .) global_elem_blk_cnts - number of elements in each block for the global mesh.

Must be called after put_init_global().

Definition at line 496 of file nemesis_io_helper.C.

 {
   nemesis_err_flag =
     Nemesis::ne_put_eb_info_global(ex_id,
                                    global_elem_blk_ids_in.data(),
                                    global_elem_blk_cnts_in.data());
  
   EX_CHECK_ERR(nemesis_err_flag, "Error writing global element block information!");
 }

References libMesh::ExodusII_IO_Helper::ex_id, and nemesis_err_flag.

Referenced by initialize().

◆ put_elem_cmap()

void libMesh::Nemesis_IO_Helper::put_elem_cmap	(	std::vector< std::vector< int >> &	elem_cmap_elem_ids,
		std::vector< std::vector< int >> &	elem_cmap_side_ids,
		std::vector< std::vector< int >> &	elem_cmap_proc_ids
	)

Writes information about elemental communication map.

Note: This class contains elem_cmap_elem_ids, elem_cmap_side_ids, abd elem_cmap_proc_ids which can be used when calling this function.

Must be called after put_cmap_params().

Definition at line 662 of file nemesis_io_helper.C.

 {
   for (auto i : index_range(elem_cmap_ids))
     {
       nemesis_err_flag =
         Nemesis::ne_put_elem_cmap(ex_id,
                                   this->elem_cmap_ids[i],
                                   elem_cmap_elem_ids_in[i].data(),
                                   elem_cmap_side_ids_in[i].data(),
                                   elem_cmap_proc_ids_in[i].data(),
                                   this->processor_id());
  
       EX_CHECK_ERR(nemesis_err_flag, "Error writing elem communication map to file!");
     }
 }

References data, elem_cmap_ids, libMesh::ExodusII_IO_Helper::ex_id, libMesh::index_range(), nemesis_err_flag, and libMesh::ParallelObject::processor_id().

Referenced by initialize().

◆ put_elem_map()

void libMesh::Nemesis_IO_Helper::put_elem_map	(	std::vector< int > &	elem_mapi,
		std::vector< int > &	elem_mapb
	)

Outputs IDs of internal and border elements.

Must be called after ne_put_loadbal_param().

Definition at line 683 of file nemesis_io_helper.C.

 {
   nemesis_err_flag =
     Nemesis::ne_put_elem_map(ex_id,
                              elem_mapi_in.empty() ? nullptr : elem_mapi_in.data(),
                              elem_mapb_in.empty() ? nullptr : elem_mapb_in.data(),
                              this->processor_id());
  
   EX_CHECK_ERR(nemesis_err_flag, "Error writing Nemesis internal and border element maps to file!");
 }

References libMesh::ExodusII_IO_Helper::ex_id, nemesis_err_flag, and libMesh::ParallelObject::processor_id().

Referenced by initialize().

◆ put_init_global()

void libMesh::Nemesis_IO_Helper::put_init_global	(	dof_id_type	num_nodes_global,
		dof_id_type	num_elems_global,
		unsigned	num_elem_blks_global,
		unsigned	num_node_sets_global,
		unsigned	num_side_sets_global
	)

Writes global information including: .) global number of nodes .) global number of elems .) global number of element blocks .) global number of node sets .) global number of side sets.

Definition at line 477 of file nemesis_io_helper.C.

 {
   nemesis_err_flag =
     Nemesis::ne_put_init_global(ex_id,
                                 num_nodes_global_in,
                                 num_elems_global_in,
                                 num_elem_blks_global_in,
                                 num_node_sets_global_in,
                                 num_side_sets_global_in);
  
   EX_CHECK_ERR(nemesis_err_flag, "Error writing initial global data!");
 }

References libMesh::ExodusII_IO_Helper::ex_id, and nemesis_err_flag.

Referenced by initialize().

◆ put_init_info()

void libMesh::Nemesis_IO_Helper::put_init_info	(	unsigned	num_proc,
		unsigned	num_proc_in_file,
		const char *	ftype
	)

Writing functions.

Writes basic info about the partitioning to file .) num_proc - number of processors .) num_proc_in_file - number of processors in the current file - generally equal to 1 .) ftype = "s" for scalar load-balance file, "p" for parallel file

Definition at line 461 of file nemesis_io_helper.C.

 {
   nemesis_err_flag =
     Nemesis::ne_put_init_info(ex_id,
                               num_proc_in,
                               num_proc_in_file_in,
                               const_cast<char *>(ftype_in));
  
   EX_CHECK_ERR(nemesis_err_flag, "Error writing initial information!");
 }

References libMesh::ExodusII_IO_Helper::ex_id, and nemesis_err_flag.

Referenced by initialize().

◆ put_loadbal_param()

void libMesh::Nemesis_IO_Helper::put_loadbal_param	(	unsigned	num_internal_nodes,
		unsigned	num_border_nodes,
		unsigned	num_external_nodes,
		unsigned	num_internal_elems,
		unsigned	num_border_elems,
		unsigned	num_node_cmaps,
		unsigned	num_elem_cmaps
	)

Writes load balance parameters, some of which are described below: .) num_internal_nodes - nodes "wholly" owned by the current processor .) num_border_nodes - nodes local to a processor but residing in an element which also has nodes on other processors .) num_external_nodes - nodes that reside on other processors but whose element "partially" resides on the current processor – we assert this should be zero on reading! .) num_border_elems - elements local to this processor but whose nodes reside on other processors as well.

.) processor - ID of the processor for which information is to be written

Definition at line 550 of file nemesis_io_helper.C.

 {
   nemesis_err_flag =
     Nemesis::ne_put_loadbal_param(ex_id,
                                   num_internal_nodes_in,
                                   num_border_nodes_in,
                                   num_external_nodes_in,
                                   num_internal_elems_in,
                                   num_border_elems_in,
                                   num_node_cmaps_in,
                                   num_elem_cmaps_in,
                                   this->processor_id());
  
   EX_CHECK_ERR(nemesis_err_flag, "Error writing loadbal parameters!");
 }

References libMesh::ExodusII_IO_Helper::ex_id, nemesis_err_flag, and libMesh::ParallelObject::processor_id().

Referenced by initialize().

◆ put_n_coord()

void libMesh::Nemesis_IO_Helper::put_n_coord	(	unsigned	start_node_num,
		unsigned	num_nodes,
		std::vector< Real > &	x_coor,
		std::vector< Real > &	y_coor,
		std::vector< Real > &	z_coor
	)

Writes the specified number of coordinate values starting at the specified index.

Definition at line 700 of file nemesis_io_helper.C.

 {
   if (num_nodes_in)
     {
       nemesis_err_flag =
         Nemesis::ne_put_n_coord(ex_id,
                                 start_node_num,
                                 num_nodes_in,
                                 x_coor.empty() ? nullptr : x_coor.data(),
                                 y_coor.empty() ? nullptr : y_coor.data(),
                                 z_coor.empty() ? nullptr : z_coor.data());
  
       EX_CHECK_ERR(nemesis_err_flag, "Error writing coords to file!");
     }
 }

References libMesh::ExodusII_IO_Helper::ex_id, and nemesis_err_flag.

◆ put_node_cmap()

void libMesh::Nemesis_IO_Helper::put_node_cmap	(	std::vector< std::vector< int >> &	node_cmap_node_ids,
		std::vector< std::vector< int >> &	node_cmap_proc_ids
	)

Outputs all of the nodal communication maps for this processor.

Internally, this function loops over all communication maps and calls Nemesis::ne_put_node_cmap() for each one.

.) node_cmap_node_ids = Nodal IDs of the FEM nodes in this communication map .) node_cmap_proc_ids = processor IDs associated with each of the nodes in node_ids

In the Nemesis file, these all appear to be written to the same chunks of data: n_comm_nids and n_comm_proc, but don't rely on these names...

Note: This class contains node_cmap_node_ids and node_cmap_proc_ids which can be used when calling this function.

Must be called after put_cmap_params().

Definition at line 600 of file nemesis_io_helper.C.

 {
  
   // Print to screen what we are about to print to Nemesis file
   if (verbose)
     {
       for (auto i : index_range(node_cmap_node_ids_in))
         {
           libMesh::out << "[" << this->processor_id() << "] put_node_cmap() : nodes communicated to proc "
                        << this->node_cmap_ids[i]
                        << " = ";
           for (const auto & node_id : node_cmap_node_ids_in[i])
             libMesh::out << node_id << " ";
           libMesh::out << std::endl;
         }
  
       for (auto i : index_range(node_cmap_node_ids_in))
         {
           libMesh::out << "[" << this->processor_id() << "] put_node_cmap() : processor IDs = ";
           for (const auto & proc_id : node_cmap_proc_ids_in[i])
             libMesh::out << proc_id << " ";
           libMesh::out << std::endl;
         }
     }
  
   for (auto i : index_range(node_cmap_node_ids_in))
     {
       int * node_ids_ptr = node_cmap_node_ids_in[i].empty() ?
         nullptr : node_cmap_node_ids_in[i].data();
       int * proc_ids_ptr = node_cmap_proc_ids_in[i].empty() ?
         nullptr : node_cmap_proc_ids_in[i].data();
  
       nemesis_err_flag =
         Nemesis::ne_put_node_cmap(ex_id, this->node_cmap_ids[i],
                                   node_ids_ptr, proc_ids_ptr,
                                   this->processor_id());
  
       EX_CHECK_ERR(nemesis_err_flag, "Error writing node communication map to file!");
     }
 }

References libMesh::ExodusII_IO_Helper::ex_id, libMesh::index_range(), nemesis_err_flag, node_cmap_ids, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

Referenced by initialize().

◆ put_node_map()

void libMesh::Nemesis_IO_Helper::put_node_map	(	std::vector< int > &	node_mapi,
		std::vector< int > &	node_mapb,
		std::vector< int > &	node_mape
	)

Outputs IDs of internal, border, and external nodes.

LibMesh asserts that the number of external nodes is zero in the Nemesis files it reads

Definition at line 645 of file nemesis_io_helper.C.

 {
   nemesis_err_flag =
     Nemesis::ne_put_node_map(ex_id,
                              node_mapi_in.empty() ? nullptr : node_mapi_in.data(),
                              node_mapb_in.empty() ? nullptr : node_mapb_in.data(),
                              node_mape_in.empty() ? nullptr : node_mape_in.data(),
                              this->processor_id());
  
   EX_CHECK_ERR(nemesis_err_flag, "Error writing Nemesis internal and border node maps to file!");
 }

References libMesh::ExodusII_IO_Helper::ex_id, nemesis_err_flag, and libMesh::ParallelObject::processor_id().

Referenced by initialize().

◆ put_ns_param_global()

void libMesh::Nemesis_IO_Helper::put_ns_param_global	(	std::vector< int > &	global_nodeset_ids,
		std::vector< int > &	num_global_node_counts,
		std::vector< int > &	num_global_node_df_counts
	)

This function writes information about global node sets.

.) global_nodeset_ids - vector of global node set IDs .) num_global_node_counts - vector of global node counts contained in each global node set .) num_global_df_count - vector of global distribution factors in each global node set

Must be called after put_init_global()

Definition at line 510 of file nemesis_io_helper.C.

 {
   // Only add nodesets if there are some
   if (global_nodeset_ids.size())
     {
       nemesis_err_flag =
         Nemesis::ne_put_ns_param_global(ex_id,
                                         global_nodeset_ids_in.data(),
                                         num_global_node_counts_in.data(),
                                         num_global_node_df_counts_in.data());
     }
  
   EX_CHECK_ERR(nemesis_err_flag, "Error writing global nodeset parameters!");
 }

References libMesh::ExodusII_IO_Helper::ex_id, global_nodeset_ids, and nemesis_err_flag.

Referenced by initialize().

◆ put_ss_param_global()

void libMesh::Nemesis_IO_Helper::put_ss_param_global	(	std::vector< int > &	global_sideset_ids,
		std::vector< int > &	num_global_side_counts,
		std::vector< int > &	num_global_side_df_counts
	)

This function writes information about global side sets.

.) global_sideset_ids - vector of global side set IDs .) num_global_side_counts - vector of global side counts contained in each global side set .) num_global_df_count - vector of global distribution factors in each global side set

Must be called after put_init_global()

Definition at line 530 of file nemesis_io_helper.C.

 {
   // Only add sidesets if there are some
   if (global_sideset_ids.size())
     {
       nemesis_err_flag =
         Nemesis::ne_put_ss_param_global(ex_id,
                                         global_sideset_ids_in.data(),
                                         num_global_side_counts_in.data(),
                                         num_global_side_df_counts_in.data());
     }
  
   EX_CHECK_ERR(nemesis_err_flag, "Error writing global sideset parameters!");
 }

References libMesh::ExodusII_IO_Helper::ex_id, global_sideset_ids, and nemesis_err_flag.

Referenced by initialize().

◆ read_all_nodesets()

void libMesh::ExodusII_IO_Helper::read_all_nodesets ( )

inherited

New API that reads all nodesets simultaneously.

This may be slightly faster than reading them one at a time. Calls ex_get_concat_node_sets() under the hood.

Definition at line 1164 of file exodusII_io_helper.C.

 {
   // Figure out how many nodesets there are in the file so we can
   // properly resize storage as necessary.
   num_node_sets =
     this->inquire
     (exII::EX_INQ_NODE_SETS,
      "Error retrieving number of node sets");
  
   // Figure out how many nodes there are in all the nodesets.
   int total_nodes_in_all_sets =
     this->inquire
     (exII::EX_INQ_NS_NODE_LEN,
      "Error retrieving number of nodes in all node sets.");
  
   // Figure out how many distribution factors there are in all the nodesets.
   int total_df_in_all_sets =
     this->inquire
     (exII::EX_INQ_NS_DF_LEN,
      "Error retrieving number of distribution factors in all node sets.");
  
   // If there are no nodesets, there's nothing to read in.
   if (num_node_sets == 0)
     return;
  
   // Allocate space to read all the nodeset data.
   // Use existing class members where possible to avoid shadowing
   nodeset_ids.clear();          nodeset_ids.resize(num_node_sets);
   num_nodes_per_set.clear();    num_nodes_per_set.resize(num_node_sets);
   num_node_df_per_set.clear();  num_node_df_per_set.resize(num_node_sets);
   node_sets_node_index.clear(); node_sets_node_index.resize(num_node_sets);
   node_sets_dist_index.clear(); node_sets_dist_index.resize(num_node_sets);
   node_sets_node_list.clear();  node_sets_node_list.resize(total_nodes_in_all_sets);
   node_sets_dist_fact.clear();  node_sets_dist_fact.resize(total_df_in_all_sets);
  
   ex_err = exII::ex_get_concat_node_sets
     (ex_id,
      nodeset_ids.data(),
      num_nodes_per_set.data(),
      num_node_df_per_set.data(),
      node_sets_node_index.data(),
      node_sets_dist_index.data(),
      node_sets_node_list.data(),
      total_df_in_all_sets ?
      MappedInputVector(node_sets_dist_fact, _single_precision).data() : nullptr);
  
   EX_CHECK_ERR(ex_err, "Error reading concatenated nodesets");
  
   // Read the nodeset names from file!
   char name_buffer[MAX_STR_LENGTH+1];
   for (int i=0; i<num_node_sets; ++i)
     {
       ex_err = exII::ex_get_name
         (ex_id,
          exII::EX_NODE_SET,
          nodeset_ids[i],
          name_buffer);
       EX_CHECK_ERR(ex_err, "Error getting node set name.");
       id_to_ns_names[nodeset_ids[i]] = name_buffer;
     }
 }

◆ read_block_info()

void libMesh::ExodusII_IO_Helper::read_block_info ( )

inherited

Reads information for all of the blocks in the ExodusII mesh file.

Definition at line 692 of file exodusII_io_helper.C.

 {
   if (num_elem_blk)
     {
       // Read all element block IDs.
       block_ids.resize(num_elem_blk);
       ex_err = exII::ex_get_ids(ex_id,
                                 exII::EX_ELEM_BLOCK,
                                 block_ids.data());
  
       EX_CHECK_ERR(ex_err, "Error getting block IDs.");
       message("All block IDs retrieved successfully.");
  
       char name_buffer[MAX_STR_LENGTH+1];
       for (int i=0; i<num_elem_blk; ++i)
         {
           ex_err = exII::ex_get_name(ex_id, exII::EX_ELEM_BLOCK,
                                      block_ids[i], name_buffer);
           EX_CHECK_ERR(ex_err, "Error getting block name.");
           id_to_block_names[block_ids[i]] = name_buffer;
         }
       message("All block names retrieved successfully.");
     }
  
   if (num_edge_blk)
     {
       // Read all edge block IDs.
       edge_block_ids.resize(num_edge_blk);
       ex_err = exII::ex_get_ids(ex_id,
                                 exII::EX_EDGE_BLOCK,
                                 edge_block_ids.data());
  
       EX_CHECK_ERR(ex_err, "Error getting edge block IDs.");
       message("All edge block IDs retrieved successfully.");
  
       // Read in edge block names
       char name_buffer[MAX_STR_LENGTH+1];
       for (int i=0; i<num_edge_blk; ++i)
         {
           ex_err = exII::ex_get_name(ex_id, exII::EX_EDGE_BLOCK,
                                      edge_block_ids[i], name_buffer);
           EX_CHECK_ERR(ex_err, "Error getting block name.");
           id_to_edge_block_names[edge_block_ids[i]] = name_buffer;
         }
       message("All edge block names retrieved successfully.");
     }
 }

References libMesh::ExodusII_IO_Helper::block_ids, libMesh::ExodusII_IO_Helper::edge_block_ids, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::id_to_block_names, libMesh::ExodusII_IO_Helper::id_to_edge_block_names, libMesh::ExodusII_IO_Helper::message(), libMesh::ExodusII_IO_Helper::num_edge_blk, and libMesh::ExodusII_IO_Helper::num_elem_blk.

◆ read_edge_blocks()

void libMesh::ExodusII_IO_Helper::read_edge_blocks ( MeshBase & mesh )

inherited

Read in edge blocks, storing information in the BoundaryInfo object.

Definition at line 852 of file exodusII_io_helper.C.

 {
   // Check for quick return if there are no edge blocks.
   if (num_edge_blk == 0)
     return;
  
   // Build data structure that we can quickly search for edges
   // and then add required BoundaryInfo information. This is a
   // map from edge->key() to a list of (elem_id, edge_id) pairs
   // for the Edge in question. Since edge->key() is edge orientation
   // invariant, this map does not distinguish different orientations
   // of the same Edge.
   typedef std::pair<dof_id_type, unsigned int> ElemEdgePair;
   std::unordered_map<dof_id_type, std::vector<ElemEdgePair>> edge_map;
   for (const auto & elem : mesh.element_ptr_range())
     for (auto e : elem->edge_index_range())
       {
         // TODO: Make various Elem::compute_key() functions
         // unprotected, this would allow us to avoid calling
         // build_edge_ptr() repeatedly in case that turns out to be
         // a bottleneck.
         std::unique_ptr<Elem> edge = elem->build_edge_ptr(e);
         dof_id_type edge_key = edge->key();
  
         // Creates vector if not already there
         auto & vec = edge_map[edge_key];
         vec.push_back(std::make_pair(elem->id(), e));
       }
  
   // Get reference to the mesh's BoundaryInfo object, as we will be
   // adding edges to this below.
   BoundaryInfo & bi = mesh.get_boundary_info();
  
   for (const auto & edge_block_id : edge_block_ids)
     {
       // exII::ex_get_block() output parameters.  Unlike the other
       // "extended" APIs, exII::ex_get_block() does not use a
       // parameter struct.
       int num_edge_this_blk = 0;
       int num_nodes_per_edge = 0;
       int num_edges_per_edge = 0;
       int num_faces_per_edge = 0;
       int num_attr_per_edge = 0;
       ex_err = exII::ex_get_block(ex_id,
                                   exII::EX_EDGE_BLOCK,
                                   edge_block_id,
                                   elem_type.data(),
                                   &num_edge_this_blk,
                                   &num_nodes_per_edge,
                                   &num_edges_per_edge, // 0 or -1 for edge blocks
                                   &num_faces_per_edge, // 0 or -1 for edge blocks
                                   &num_attr_per_edge);
  
       EX_CHECK_ERR(ex_err, "Error getting edge block info.");
       message("Info retrieved successfully for block: ", edge_block_id);
  
       // Read in the connectivity of the edges of this block,
       // watching out for the case where we actually have no
       // elements in this block (possible with parallel files)
       connect.resize(num_nodes_per_edge * num_edge_this_blk);
  
       if (!connect.empty())
         {
           ex_err = exII::ex_get_conn(ex_id,
                                      exII::EX_EDGE_BLOCK,
                                      edge_block_id,
                                      connect.data(), // node_conn
                                      nullptr,        // elem_edge_conn (unused)
                                      nullptr);       // elem_face_conn (unused)
  
           EX_CHECK_ERR(ex_err, "Error reading block connectivity.");
           message("Connectivity retrieved successfully for block: ", edge_block_id);
  
           // All edge types have an identity mapping from the corresponding
           // Exodus type, so we don't need to bother with mapping ids, but
           // we do need to know what kind of elements to build.
           const auto & conv = get_conversion(std::string(elem_type.data()));
  
           // Loop over indices in connectivity array, build edge elements,
           // look them up in the edge_map.
           for (unsigned int i=0, sz=connect.size(); i<sz; i+=num_nodes_per_edge)
             {
               auto edge = Elem::build(conv.libmesh_elem_type());
               for (int n=0; n<num_nodes_per_edge; ++n)
                 {
                   int exodus_node_id = connect[i+n];
                   int exodus_node_id_zero_based = exodus_node_id - 1;
                   int libmesh_node_id = node_num_map[exodus_node_id_zero_based] - 1;
  
                   edge->set_node(n) = mesh.node_ptr(libmesh_node_id);
                 }
  
               // Compute key for the edge Elem we just built.
               dof_id_type edge_key = edge->key();
  
               // If this key is not found in the edge_map, which is
               // supposed to include every edge in the Mesh, then we
               // need to throw an error.
               auto & elem_edge_pair_vec =
                 libmesh_map_find(edge_map, edge_key);
  
               for (const auto & elem_edge_pair : elem_edge_pair_vec)
                 {
                   // We only want to match edges which have the same
                   // orientation (node ordering) to the one in the
                   // Exodus file, otherwise we ignore this elem_edge_pair.
                   //
                   // Note: this also handles the situation where two
                   // edges have the same key (hash collision) as then
                   // this check avoids a false positive.
  
                   // Build edge indicated by elem_edge_pair
                   auto candidate_edge =
                     mesh.elem_ptr(elem_edge_pair.first)->
                     build_edge_ptr(elem_edge_pair.second);
  
                   // Determine whether this candidate edge is a "real" match,
                   // i.e. also has the same orientation.
                   bool is_match = true;
                   for (int n=0; n<num_nodes_per_edge; ++n)
                     if (candidate_edge->node_id(n) != edge->node_id(n))
                       {
                         is_match = false;
                         break;
                       }
  
                   if (is_match)
                     {
                       // Add this (elem, edge, id) combo to the BoundaryInfo object.
                       bi.add_edge(elem_edge_pair.first,
                                   elem_edge_pair.second,
                                   edge_block_id);
                     }
                 } // end loop over elem_edge_pairs
             } // end loop over connectivity array
  
           // Set edgeset name in the BoundaryInfo object.
           bi.edgeset_name(edge_block_id) = id_to_edge_block_names[edge_block_id];
         } // end if !connect.empty()
     } // end for edge_block_id : edge_block_ids
 }

◆ read_elem_in_block()

void libMesh::ExodusII_IO_Helper::read_elem_in_block ( int block )

inherited

Reads all of the element connectivity for block block in the ExodusII mesh file.

Definition at line 797 of file exodusII_io_helper.C.

 {
   libmesh_assert_less (block, block_ids.size());
  
   // Unlike the other "extended" APIs, this one does not use a parameter struct.
   int num_edges_per_elem = 0;
   int num_faces_per_elem = 0;
   ex_err = exII::ex_get_block(ex_id,
                               exII::EX_ELEM_BLOCK,
                               block_ids[block],
                               elem_type.data(),
                               &num_elem_this_blk,
                               &num_nodes_per_elem,
                               &num_edges_per_elem, // 0 or -1 if no "extended" block info
                               &num_faces_per_elem, // 0 or -1 if no "extended" block info
                               &num_attr);
  
   EX_CHECK_ERR(ex_err, "Error getting block info.");
   message("Info retrieved successfully for block: ", block);
  
   // Warn that we don't currently support reading blocks with extended info.
   // Note: the docs say -1 will be returned for this but I found that it was
   // actually 0, so not sure which it will be in general.
   if (!(num_edges_per_elem == 0) && !(num_edges_per_elem == -1))
     libmesh_warning("Exodus files with extended edge connectivity not currently supported.");
   if (!(num_faces_per_elem == 0) && !(num_faces_per_elem == -1))
     libmesh_warning("Exodus files with extended face connectivity not currently supported.");
  
   if (verbose)
     libMesh::out << "Read a block of " << num_elem_this_blk
                  << " " << elem_type.data() << "(s)"
                  << " having " << num_nodes_per_elem
                  << " nodes per element." << std::endl;
  
   // Read in the connectivity of the elements of this block,
   // watching out for the case where we actually have no
   // elements in this block (possible with parallel files)
   connect.resize(num_nodes_per_elem*num_elem_this_blk);
  
   if (!connect.empty())
     {
       ex_err = exII::ex_get_conn(ex_id,
                                  exII::EX_ELEM_BLOCK,
                                  block_ids[block],
                                  connect.data(), // node_conn
                                  nullptr,        // elem_edge_conn (unused)
                                  nullptr);       // elem_face_conn (unused)
  
       EX_CHECK_ERR(ex_err, "Error reading block connectivity.");
       message("Connectivity retrieved successfully for block: ", block);
     }
 }

References libMesh::ExodusII_IO_Helper::block_ids, libMesh::ExodusII_IO_Helper::connect, libMesh::ExodusII_IO_Helper::elem_type, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::message(), libMesh::ExodusII_IO_Helper::num_attr, libMesh::ExodusII_IO_Helper::num_elem_this_blk, libMesh::ExodusII_IO_Helper::num_nodes_per_elem, libMesh::out, and libMesh::ExodusII_IO_Helper::verbose.

◆ read_elem_num_map()

void libMesh::ExodusII_IO_Helper::read_elem_num_map ( )

inherited

Reads the optional node_num_map from the ExodusII mesh file.

Definition at line 996 of file exodusII_io_helper.C.

 {
   elem_num_map.resize(num_elem);
  
   // Note: we cannot use the exII::ex_get_num_map() here because it
   // (apparently) does not behave like ex_get_elem_num_map() when
   // there is no elem number map in the file: it throws an error
   // instead of returning a default identity array (1,2,3,...).
   ex_err = exII::ex_get_elem_num_map
     (ex_id, elem_num_map.empty() ? nullptr : elem_num_map.data());
  
   EX_CHECK_ERR(ex_err, "Error retrieving element number map.");
   message("Element numbering map retrieved successfully.");
  
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] elem_num_map[i] = ";
       for (unsigned int i=0; i<static_cast<unsigned int>(std::min(10, num_elem-1)); ++i)
         libMesh::out << elem_num_map[i] << ", ";
       libMesh::out << "... " << elem_num_map.back() << std::endl;
     }
 }

References libMesh::ExodusII_IO_Helper::elem_num_map, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::message(), libMesh::ExodusII_IO_Helper::num_elem, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ read_elemental_var_values()

void libMesh::ExodusII_IO_Helper::read_elemental_var_values	(	std::string	elemental_var_name,
		int	time_step,
		std::map< dof_id_type, Real > &	elem_var_value_map
	)

inherited

Reads elemental values for the variable 'elemental_var_name' at the specified timestep into the 'elem_var_value_map' which is passed in.

Definition at line 1476 of file exodusII_io_helper.C.

 {
   this->read_var_names(ELEMENTAL);
  
   // See if we can find the variable we are looking for
   unsigned int var_index = 0;
   bool found = false;
  
   // Do a linear search for elem_var_name in elemental_var_names
   for (; var_index != elem_var_names.size(); ++var_index)
     if (elem_var_names[var_index] == elemental_var_name)
       {
         found = true;
         break;
       }
  
   if (!found)
     {
       libMesh::err << "Available variables: " << std::endl;
       for (const auto & var_name : elem_var_names)
         libMesh::err << var_name << std::endl;
  
       libmesh_error_msg("Unable to locate variable named: " << elemental_var_name);
     }
  
   // Sequential index which we can use to look up the element ID in the elem_num_map.
   unsigned ex_el_num = 0;
  
   // Element variable truth table
   std::vector<int> var_table(block_ids.size() * elem_var_names.size());
   exII::ex_get_var_tab(ex_id, "e", block_ids.size(), elem_var_names.size(), var_table.data());
  
   for (unsigned i=0; i<static_cast<unsigned>(num_elem_blk); i++)
     {
       ex_err = exII::ex_get_elem_block(ex_id,
                                        block_ids[i],
                                        nullptr,
                                        &num_elem_this_blk,
                                        nullptr,
                                        nullptr);
       EX_CHECK_ERR(ex_err, "Error getting number of elements in block.");
  
       // If the current variable isn't active on this subdomain, advance
       // the index by the number of elements on this block and go to the
       // next loop iteration.
       if (!var_table[elem_var_names.size()*i + var_index])
         {
           ex_el_num += num_elem_this_blk;
           continue;
         }
  
       std::vector<Real> block_elem_var_values(num_elem_this_blk);
  
       ex_err = exII::ex_get_elem_var
         (ex_id,
          time_step,
          var_index+1,
          block_ids[i],
          num_elem_this_blk,
          MappedInputVector(block_elem_var_values, _single_precision).data());
       EX_CHECK_ERR(ex_err, "Error getting elemental values.");
  
       for (unsigned j=0; j<static_cast<unsigned>(num_elem_this_blk); j++)
         {
           // Use the elem_num_map to obtain the ID of this element in the Exodus file,
           // and remember to subtract 1 since libmesh is zero-based and Exodus is 1-based.
           unsigned mapped_elem_id = this->elem_num_map[ex_el_num] - 1;
  
           // Store the elemental value in the map.
           elem_var_value_map[mapped_elem_id] = block_elem_var_values[j];
  
           // Go to the next sequential element ID.
           ex_el_num++;
         }
     }
 }

References libMesh::ExodusII_IO_Helper::_single_precision, libMesh::ExodusII_IO_Helper::block_ids, data, libMesh::ExodusII_IO_Helper::elem_num_map, libMesh::ExodusII_IO_Helper::elem_var_names, libMesh::ExodusII_IO_Helper::ELEMENTAL, libMesh::err, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::num_elem_blk, libMesh::ExodusII_IO_Helper::num_elem_this_blk, and libMesh::ExodusII_IO_Helper::read_var_names().

◆ read_global_values()

void libMesh::ExodusII_IO_Helper::read_global_values	(	std::vector< Real > &	values,
		int	timestep
	)

inherited

Reads the vector of global variables.

Definition at line 3151 of file exodusII_io_helper.C.

 {
   if ((_run_only_on_proc0) && (this->processor_id() != 0))
     return;
  
   values.clear();
   values.resize(num_global_vars);
   ex_err = exII::ex_get_glob_vars
     (ex_id, timestep, num_global_vars,
      MappedInputVector(values, _single_precision).data());
  
   EX_CHECK_ERR(ex_err, "Error reading global values.");
 }

References libMesh::ExodusII_IO_Helper::_run_only_on_proc0, libMesh::ExodusII_IO_Helper::_single_precision, data, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::num_global_vars, and libMesh::ParallelObject::processor_id().

◆ read_header()

void libMesh::ExodusII_IO_Helper::read_header ( )

inherited

Reads an ExodusII mesh file header.

Definition at line 532 of file exodusII_io_helper.C.

 {
   // Read init params using newer API that reads into a struct.  For
   // backwards compatibility, assign local member values from struct
   // afterwards. Note: using the new API allows us to automatically
   // read edge and face block/set information if it's present in the
   // file.
   exII::ex_init_params params = {};
   ex_err = exII::ex_get_init_ext(ex_id, &params);
   EX_CHECK_ERR(ex_err, "Error retrieving header info.");
  
   // "title" has MAX_LINE_LENGTH+1 characters, but the last one is reserved
   // for null termination so we only copy MAX_LINE_LENGTH chars.
   title.assign(params.title, params.title + MAX_LINE_LENGTH);
   num_dim = params.num_dim;
   num_nodes = params.num_nodes;
   num_elem = params.num_elem;
   num_elem_blk = params.num_elem_blk;
   num_node_sets = params.num_node_sets;
   num_side_sets = params.num_side_sets;
   num_edge_blk = params.num_edge_blk;
   num_edge = params.num_edge;
  
   this->read_num_time_steps();
  
   ex_err = exII::ex_get_var_param(ex_id, "n", &num_nodal_vars);
   EX_CHECK_ERR(ex_err, "Error reading number of nodal variables.");
  
   ex_err = exII::ex_get_var_param(ex_id, "e", &num_elem_vars);
   EX_CHECK_ERR(ex_err, "Error reading number of elemental variables.");
  
   ex_err = exII::ex_get_var_param(ex_id, "g", &num_global_vars);
   EX_CHECK_ERR(ex_err, "Error reading number of global variables.");
  
   ex_err = exII::ex_get_var_param(ex_id, "s", &num_sideset_vars);
   EX_CHECK_ERR(ex_err, "Error reading number of sideset variables.");
  
   message("Exodus header info retrieved successfully.");
 }

◆ read_nodal_var_values()

void libMesh::ExodusII_IO_Helper::read_nodal_var_values	(	std::string	nodal_var_name,
		int	time_step
	)

inherited

Reads the nodal values for the variable 'nodal_var_name' at the specified time into the 'nodal_var_values' array.

Definition at line 1291 of file exodusII_io_helper.C.

 {
   // Read the nodal variable names from file, so we can see if we have the one we're looking for
   this->read_var_names(NODAL);
  
   // See if we can find the variable we are looking for
   unsigned int var_index = 0;
   bool found = false;
  
   // Do a linear search for nodal_var_name in nodal_var_names
   for (; var_index<nodal_var_names.size(); ++var_index)
     {
       found = (nodal_var_names[var_index] == nodal_var_name);
       if (found)
         break;
     }
  
   if (!found)
     {
       libMesh::err << "Available variables: " << std::endl;
       for (const auto & var_name : nodal_var_names)
         libMesh::err << var_name << std::endl;
  
       libmesh_error_msg("Unable to locate variable named: " << nodal_var_name);
     }
  
   // Allocate enough space to store the nodal variable values
   nodal_var_values.resize(num_nodes);
  
   std::vector<Real> unmapped_nodal_var_values(num_nodes);
  
   // Call the Exodus API to read the nodal variable values
   ex_err = exII::ex_get_nodal_var
     (ex_id,
      time_step,
      var_index+1,
      num_nodes,
      MappedInputVector(unmapped_nodal_var_values, _single_precision).data());
   EX_CHECK_ERR(ex_err, "Error reading nodal variable values!");
  
   for (unsigned i=0; i<static_cast<unsigned>(num_nodes); i++)
     {
       // Use the node_num_map to obtain the ID of this node in the Exodus file,
       // and remember to subtract 1 since libmesh is zero-based and Exodus is 1-based.
       unsigned mapped_node_id = this->node_num_map[i] - 1;
  
       // Store the nodal value in the map.
       nodal_var_values[mapped_node_id] = unmapped_nodal_var_values[i];
     }
 }

References libMesh::ExodusII_IO_Helper::_single_precision, data, libMesh::err, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::NODAL, libMesh::ExodusII_IO_Helper::nodal_var_names, libMesh::ExodusII_IO_Helper::nodal_var_values, libMesh::ExodusII_IO_Helper::node_num_map, libMesh::ExodusII_IO_Helper::num_nodes, and libMesh::ExodusII_IO_Helper::read_var_names().

◆ read_node_num_map()

void libMesh::ExodusII_IO_Helper::read_node_num_map ( )

inherited

Reads the optional node_num_map from the ExodusII mesh file.

Definition at line 660 of file exodusII_io_helper.C.

 {
   node_num_map.resize(num_nodes);
  
   // Note: we cannot use the exII::ex_get_num_map() here because it
   // (apparently) does not behave like ex_get_node_num_map() when
   // there is no node number map in the file: it throws an error
   // instead of returning a default identity array (1,2,3,...).
   ex_err = exII::ex_get_node_num_map
     (ex_id, node_num_map.empty() ? nullptr : node_num_map.data());
  
   EX_CHECK_ERR(ex_err, "Error retrieving nodal number map.");
   message("Nodal numbering map retrieved successfully.");
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] node_num_map[i] = ";
       for (unsigned int i=0; i<static_cast<unsigned int>(std::min(10, num_nodes-1)); ++i)
         libMesh::out << node_num_map[i] << ", ";
       libMesh::out << "... " << node_num_map.back() << std::endl;
     }
 }

References libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::message(), libMesh::ExodusII_IO_Helper::node_num_map, libMesh::ExodusII_IO_Helper::num_nodes, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ read_nodes()

void libMesh::ExodusII_IO_Helper::read_nodes ( )

inherited

Reads the nodal data (x,y,z coordinates) from the ExodusII mesh file.

Definition at line 639 of file exodusII_io_helper.C.

 {
   x.resize(num_nodes);
   y.resize(num_nodes);
   z.resize(num_nodes);
  
   if (num_nodes)
     {
       ex_err = exII::ex_get_coord
         (ex_id,
          MappedInputVector(x, _single_precision).data(),
          MappedInputVector(y, _single_precision).data(),
          MappedInputVector(z, _single_precision).data());
  
       EX_CHECK_ERR(ex_err, "Error retrieving nodal data.");
       message("Nodal data retrieved successfully.");
     }
 }

References libMesh::ExodusII_IO_Helper::_single_precision, data, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::message(), libMesh::ExodusII_IO_Helper::num_nodes, libMesh::ExodusII_IO_Helper::x, libMesh::ExodusII_IO_Helper::y, and libMesh::ExodusII_IO_Helper::z.

◆ read_nodeset()

void libMesh::ExodusII_IO_Helper::read_nodeset ( int id )

inherited

Reads information about nodeset id and inserts it into the global nodeset array at the position offset.

Definition at line 1131 of file exodusII_io_helper.C.

 {
   libmesh_assert_less (id, nodeset_ids.size());
   libmesh_assert_less (id, num_nodes_per_set.size());
   libmesh_assert_less (id, num_node_df_per_set.size());
  
   ex_err = exII::ex_get_set_param(ex_id,
                                   exII::EX_NODE_SET,
                                   nodeset_ids[id],
                                   &num_nodes_per_set[id],
                                   &num_node_df_per_set[id]);
   EX_CHECK_ERR(ex_err, "Error retrieving nodeset parameters.");
   message("Parameters retrieved successfully for nodeset: ", id);
  
   node_list.resize(num_nodes_per_set[id]);
  
   // Don't call ex_get_set unless there are actually nodes there to get.
   // Exodus prints an annoying warning message in DEBUG mode otherwise...
   if (num_nodes_per_set[id] > 0)
     {
       ex_err = exII::ex_get_set(ex_id,
                                 exII::EX_NODE_SET,
                                 nodeset_ids[id],
                                 node_list.data(),
                                 nullptr); // set_extra_list, ignored for node sets
  
       EX_CHECK_ERR(ex_err, "Error retrieving nodeset data.");
       message("Data retrieved successfully for nodeset: ", id);
     }
 }

References libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::message(), libMesh::ExodusII_IO_Helper::node_list, libMesh::ExodusII_IO_Helper::nodeset_ids, libMesh::ExodusII_IO_Helper::num_node_df_per_set, and libMesh::ExodusII_IO_Helper::num_nodes_per_set.

◆ read_nodeset_info()

void libMesh::ExodusII_IO_Helper::read_nodeset_info ( )

inherited

Reads information about all of the nodesets in the ExodusII mesh file.

Definition at line 1057 of file exodusII_io_helper.C.

 {
   nodeset_ids.resize(num_node_sets);
   if (num_node_sets > 0)
     {
       ex_err = exII::ex_get_ids(ex_id,
                                 exII::EX_NODE_SET,
                                 nodeset_ids.data());
       EX_CHECK_ERR(ex_err, "Error retrieving nodeset information.");
       message("All nodeset information retrieved successfully.");
  
       // Resize appropriate data structures -- only do this once outnode the loop
       num_nodes_per_set.resize(num_node_sets);
       num_node_df_per_set.resize(num_node_sets);
     }
  
   char name_buffer[MAX_STR_LENGTH+1];
   for (int i=0; i<num_node_sets; ++i)
     {
       ex_err = exII::ex_get_name(ex_id, exII::EX_NODE_SET,
                                  nodeset_ids[i], name_buffer);
       EX_CHECK_ERR(ex_err, "Error getting node set name.");
       id_to_ns_names[nodeset_ids[i]] = name_buffer;
     }
   message("All node set names retrieved successfully.");
 }

References libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::id_to_ns_names, libMesh::ExodusII_IO_Helper::message(), libMesh::ExodusII_IO_Helper::nodeset_ids, libMesh::ExodusII_IO_Helper::num_node_df_per_set, libMesh::ExodusII_IO_Helper::num_node_sets, and libMesh::ExodusII_IO_Helper::num_nodes_per_set.

◆ read_num_time_steps()

void libMesh::ExodusII_IO_Helper::read_num_time_steps ( )

inherited

Reads the number of timesteps currently stored in the Exodus file and stores it in the num_time_steps variable.

Definition at line 1283 of file exodusII_io_helper.C.

 {
   num_time_steps =
     this->inquire(exII::EX_INQ_TIME, "Error retrieving number of time steps");
 }

References libMesh::ExodusII_IO_Helper::inquire(), and libMesh::ExodusII_IO_Helper::num_time_steps.

Referenced by libMesh::ExodusII_IO_Helper::read_header(), and libMesh::ExodusII_IO_Helper::read_time_steps().

◆ read_qa_records()

void libMesh::ExodusII_IO_Helper::read_qa_records ( )

inherited

Reads the QA records from an ExodusII file.

We can use this to detect when e.g. CUBIT 14+ was used to generate a Mesh file, and work around certain known bugs in that version.

Definition at line 575 of file exodusII_io_helper.C.

 {
   // The QA records are four MAX_STR_LENGTH-byte character strings.
   int num_qa_rec =
     this->inquire(exII::EX_INQ_QA, "Error retrieving number of QA records");
  
   if (verbose)
     libMesh::out << "Found "
                  << num_qa_rec
                  << " QA record(s) in the Exodus file."
                  << std::endl;
  
   if (num_qa_rec > 0)
     {
       // How to dynamically allocate an array of fixed-size char * arrays in C++.
       // http://stackoverflow.com/questions/8529359/creating-a-dynamic-sized-array-of-fixed-sized-int-arrays-in-c
       typedef char * inner_array_t[4];
       inner_array_t * qa_record = new inner_array_t[num_qa_rec];
  
       for (int i=0; i<num_qa_rec; i++)
         for (int j=0; j<4; j++)
           qa_record[i][j] = new char[MAX_STR_LENGTH+1];
  
       ex_err = exII::ex_get_qa (ex_id, qa_record);
       EX_CHECK_ERR(ex_err, "Error reading the QA records.");
  
       // Print the QA records
       if (verbose)
         {
           for (int i=0; i<num_qa_rec; i++)
             {
               libMesh::out << "QA Record: " << i << std::endl;
               for (int j=0; j<4; j++)
                 libMesh::out << qa_record[i][j] << std::endl;
             }
         }
  
  
       // Clean up dynamically-allocated memory
       for (int i=0; i<num_qa_rec; i++)
         for (int j=0; j<4; j++)
           delete [] qa_record[i][j];
  
       delete [] qa_record;
     }
 }

References libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::inquire(), libMesh::out, and libMesh::ExodusII_IO_Helper::verbose.

◆ read_sideset()

void libMesh::ExodusII_IO_Helper::read_sideset	(	int	id,
		int	offset
	)

inherited

Reads information about sideset id and inserts it into the global sideset array at the position offset.

Definition at line 1086 of file exodusII_io_helper.C.

 {
   libmesh_assert_less (id, ss_ids.size());
   libmesh_assert_less (id, num_sides_per_set.size());
   libmesh_assert_less (id, num_df_per_set.size());
   libmesh_assert_less_equal (offset, elem_list.size());
   libmesh_assert_less_equal (offset, side_list.size());
  
   ex_err = exII::ex_get_set_param(ex_id,
                                   exII::EX_SIDE_SET,
                                   ss_ids[id],
                                   &num_sides_per_set[id],
                                   &num_df_per_set[id]);
   EX_CHECK_ERR(ex_err, "Error retrieving sideset parameters.");
   message("Parameters retrieved successfully for sideset: ", id);
  
  
   // It's OK for offset==elem_list.size() as long as num_sides_per_set[id]==0
   // because in that case we don't actually read anything...
 #ifdef DEBUG
   if (static_cast<unsigned int>(offset) == elem_list.size() ||
       static_cast<unsigned int>(offset) == side_list.size() )
     libmesh_assert_equal_to (num_sides_per_set[id], 0);
 #endif
  
  
   // Don't call ex_get_set unless there are actually sides there to get.
   // Exodus prints an annoying warning in DEBUG mode otherwise...
   if (num_sides_per_set[id] > 0)
     {
       ex_err = exII::ex_get_set(ex_id,
                                 exII::EX_SIDE_SET,
                                 ss_ids[id],
                                 &elem_list[offset],
                                 &side_list[offset]);
       EX_CHECK_ERR(ex_err, "Error retrieving sideset data.");
       message("Data retrieved successfully for sideset: ", id);
  
       for (int i=0; i<num_sides_per_set[id]; i++)
         id_list[i+offset] = ss_ids[id];
     }
 }

References libMesh::ExodusII_IO_Helper::elem_list, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::id_list, libMesh::ExodusII_IO_Helper::message(), libMesh::ExodusII_IO_Helper::num_df_per_set, libMesh::ExodusII_IO_Helper::num_sides_per_set, libMesh::ExodusII_IO_Helper::side_list, and libMesh::ExodusII_IO_Helper::ss_ids.

Referenced by libMesh::ExodusII_IO_Helper::write_sideset_data().

◆ read_sideset_data()

void libMesh::ExodusII_IO_Helper::read_sideset_data	(	const MeshBase &	mesh,
		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
	)

inherited

Read sideset variables, if any, into the provided data structures.

Definition at line 2736 of file exodusII_io_helper.C.

 {
   // This reads the sideset variable names into the local
   // sideset_var_names data structure.
   this->read_var_names(SIDESET);
  
   if (num_sideset_vars)
     {
       // Read the sideset data truth table
       std::vector<int> sset_var_tab(num_side_sets * num_sideset_vars);
       ex_err = exII::ex_get_sset_var_tab
         (ex_id,
          num_side_sets,
          num_sideset_vars,
          sset_var_tab.data());
       EX_CHECK_ERR(ex_err, "Error reading sideset variable truth table.");
  
       // Set up/allocate space in incoming data structures.
       var_names = sideset_var_names;
       side_ids.resize(num_sideset_vars);
       bc_vals.resize(num_sideset_vars);
  
       // Read the sideset data.
       //
       // Note: we assume that read_sideset() has already been called
       // for each sideset, so the required values in elem_list and
       // side_list are already present.
       //
       // TODO: As a future optimization, we could read only the values
       // requested by the user by looking at the input parameter
       // var_names and checking whether it already has entries in
       // it. We could do the same thing with the input side_ids
       // container and only read values for requested sidesets.
       int offset=0;
       for (int ss=0; ss<num_side_sets; ++ss)
         {
           offset += (ss > 0 ? num_sides_per_set[ss-1] : 0);
           for (int var=0; var<num_sideset_vars; ++var)
             {
               int is_present = sset_var_tab[num_sideset_vars*ss + var];
  
               if (is_present)
                 {
                   // Record the fact that this variable is defined on this sideset.
                   side_ids[var].insert(ss_ids[ss]);
  
                   // Note: the assumption here is that a previous call
                   // to this->read_sideset_info() has already set the
                   // values of num_sides_per_set, so we just use those values here.
                   std::vector<Real> sset_var_vals(num_sides_per_set[ss]);
                   ex_err = exII::ex_get_sset_var
                     (ex_id,
                      timestep,
                      var + 1, // 1-based sideset variable index!
                      ss_ids[ss],
                      num_sides_per_set[ss],
                      MappedInputVector(sset_var_vals, _single_precision).data());
                   EX_CHECK_ERR(ex_err, "Error reading sideset variable.");
  
                   // Debugging:
                   // libMesh::out << "Variable " << sideset_var_names[var]
                   //              << " is defined on side set " << ss_ids[ss]
                   //              << " and has values: " << std::endl;
                   // for (int i=0; i<num_sides_per_set[ss]; ++i)
                   //   libMesh::out << sset_var_vals[i] << " ";
                   // libMesh::out << std::endl;
  
                   for (int i=0; i<num_sides_per_set[ss]; ++i)
                     {
                       dof_id_type exodus_elem_id = elem_list[i + offset];
                       unsigned int exodus_side_id = side_list[i + offset];
  
                       // FIXME: We should use exodus_elem_num_to_libmesh for this,
                       // but it apparently is never set up, so just
                       // subtract 1 from the Exodus elem id.
                       dof_id_type converted_elem_id = exodus_elem_id - 1;
  
                       // Map Exodus side id to libmesh side id.
                       // Map from Exodus side ids to libmesh side ids.
                       const auto & conv = get_conversion(mesh.elem_ptr(converted_elem_id)->type());
  
                       // Map from Exodus side id to libmesh side id.
                       // Note: the mapping is defined on 0-based indices, so subtract
                       // 1 before doing the mapping.
                       unsigned int converted_side_id = conv.get_side_map(exodus_side_id - 1);
  
                       // Debugging:
                       // libMesh::out << "exodus_elem_id = " << exodus_elem_id
                       //              << "\n"
                       //              << "converted_elem_id = " << converted_elem_id
                       //              << "\n\n"
                       //              << "exodus_side_id = " << exodus_side_id
                       //              << "\n"
                       //              << "converted_side_id = " << converted_side_id
                       //              << std::endl;
  
                       // Make a BCTuple key from the converted information.
                       BoundaryInfo::BCTuple key = std::make_tuple
                         (converted_elem_id,
                          converted_side_id,
                          ss_ids[ss]);
  
                       // Store (elem, side, b_id) tuples in bc_vals[var]
                       bc_vals[var].insert(std::make_pair(key, sset_var_vals[i]));
                     } // end for (i)
                 } // end if (present)
             } // end for (var)
         } // end for (ss)
     } // end if (num_sideset_vars)
 }

◆ read_sideset_info()

void libMesh::ExodusII_IO_Helper::read_sideset_info ( )

inherited

Reads information about all of the sidesets in the ExodusII mesh file.

Definition at line 1022 of file exodusII_io_helper.C.

 {
   ss_ids.resize(num_side_sets);
   if (num_side_sets > 0)
     {
       ex_err = exII::ex_get_ids(ex_id,
                                 exII::EX_SIDE_SET,
                                 ss_ids.data());
       EX_CHECK_ERR(ex_err, "Error retrieving sideset information.");
       message("All sideset information retrieved successfully.");
  
       // Resize appropriate data structures -- only do this once outside the loop
       num_sides_per_set.resize(num_side_sets);
       num_df_per_set.resize(num_side_sets);
  
       // Inquire about the length of the concatenated side sets element list
       num_elem_all_sidesets = inquire(exII::EX_INQ_SS_ELEM_LEN, "Error retrieving length of the concatenated side sets element list!");
  
       elem_list.resize (num_elem_all_sidesets);
       side_list.resize (num_elem_all_sidesets);
       id_list.resize   (num_elem_all_sidesets);
     }
  
   char name_buffer[MAX_STR_LENGTH+1];
   for (int i=0; i<num_side_sets; ++i)
     {
       ex_err = exII::ex_get_name(ex_id, exII::EX_SIDE_SET,
                                  ss_ids[i], name_buffer);
       EX_CHECK_ERR(ex_err, "Error getting side set name.");
       id_to_ss_names[ss_ids[i]] = name_buffer;
     }
   message("All side set names retrieved successfully.");
 }

Referenced by libMesh::ExodusII_IO_Helper::write_sideset_data().

◆ read_time_steps()

void libMesh::ExodusII_IO_Helper::read_time_steps ( )

inherited

Reads and stores the timesteps in the 'time_steps' array.

Definition at line 1266 of file exodusII_io_helper.C.

 {
   // Make sure we have an up-to-date count of the number of time steps in the file.
   this->read_num_time_steps();
  
   if (num_time_steps > 0)
     {
       time_steps.resize(num_time_steps);
       ex_err = exII::ex_get_all_times
         (ex_id,
          MappedInputVector(time_steps, _single_precision).data());
       EX_CHECK_ERR(ex_err, "Error reading timesteps!");
     }
 }

References libMesh::ExodusII_IO_Helper::_single_precision, data, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::num_time_steps, libMesh::ExodusII_IO_Helper::read_num_time_steps(), and libMesh::ExodusII_IO_Helper::time_steps.

◆ read_var_names()

void libMesh::ExodusII_IO_Helper::read_var_names ( ExodusVarType type )

inherited

Definition at line 1344 of file exodusII_io_helper.C.

 {
   switch (type)
     {
     case NODAL:
       this->read_var_names_impl("n", num_nodal_vars, nodal_var_names);
       break;
     case ELEMENTAL:
       this->read_var_names_impl("e", num_elem_vars, elem_var_names);
       break;
     case GLOBAL:
       this->read_var_names_impl("g", num_global_vars, global_var_names);
       break;
     case SIDESET:
       this->read_var_names_impl("s", num_sideset_vars, sideset_var_names);
       break;
     default:
       libmesh_error_msg("Unrecognized ExodusVarType " << type);
     }
 }

Referenced by libMesh::ExodusII_IO_Helper::check_existing_vars(), libMesh::ExodusII_IO_Helper::read_elemental_var_values(), libMesh::ExodusII_IO_Helper::read_nodal_var_values(), and libMesh::ExodusII_IO_Helper::read_sideset_data().

◆ read_var_names_impl()

void libMesh::ExodusII_IO_Helper::read_var_names_impl	(	const char *	var_type,
		int &	count,
		std::vector< std::string > &	result
	)

privateinherited

read_var_names() dispatches to this function.

Definition at line 1367 of file exodusII_io_helper.C.

 {
   // First read and store the number of names we have
   ex_err = exII::ex_get_var_param(ex_id, var_type, &count);
   EX_CHECK_ERR(ex_err, "Error reading number of variables.");
  
   // Do nothing if no variables are detected
   if (count == 0)
     return;
  
   // Second read the actual names and convert them into a format we can use
   NamesData names_table(count, MAX_STR_LENGTH);
  
   ex_err = exII::ex_get_var_names(ex_id,
                                   var_type,
                                   count,
                                   names_table.get_char_star_star()
                                   );
   EX_CHECK_ERR(ex_err, "Error reading variable names!");
  
   if (verbose)
     {
       libMesh::out << "Read the variable(s) from the file:" << std::endl;
       for (int i=0; i<count; i++)
         libMesh::out << names_table.get_char_star(i) << std::endl;
     }
  
   // Allocate enough space for our variable name strings.
   result.resize(count);
  
   // Copy the char buffers into strings.
   for (int i=0; i<count; i++)
     result[i] = names_table.get_char_star(i); // calls string::op=(const char *)
 }

References libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::NamesData::get_char_star(), libMesh::ExodusII_IO_Helper::NamesData::get_char_star_star(), libMesh::out, and libMesh::ExodusII_IO_Helper::verbose.

Referenced by libMesh::ExodusII_IO_Helper::read_var_names().

◆ set_coordinate_offset()

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

inherited

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

Definition at line 3181 of file exodusII_io_helper.C.

 {
   _coordinate_offset = p;
 }

References libMesh::ExodusII_IO_Helper::_coordinate_offset.

◆ use_mesh_dimension_instead_of_spatial_dimension()

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

inherited

Sets the underlying value of the boolean flag _use_mesh_dimension_instead_of_spatial_dimension.

By default, the value of this flag is false.

See the ExodusII_IO class documentation for a detailed description of this flag.

Definition at line 3167 of file exodusII_io_helper.C.

 {
   _use_mesh_dimension_instead_of_spatial_dimension = val;
 }

References libMesh::ExodusII_IO_Helper::_use_mesh_dimension_instead_of_spatial_dimension.

◆ write_as_dimension()

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

inherited

Sets the value of _write_as_dimension.

This directly controls 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 3174 of file exodusII_io_helper.C.

 {
   _write_as_dimension = dim;
 }

References libMesh::ExodusII_IO_Helper::_write_as_dimension, and dim.

◆ write_element_values() [1/2]

void libMesh::Nemesis_IO_Helper::write_element_values	(	const MeshBase &	mesh,
		const EquationSystems &	es,
		const std::vector< std::pair< unsigned int, unsigned int >> &	var_nums,
		int	timestep,
		const std::vector< std::set< subdomain_id_type >> &	vars_active_subdomains
	)

Writes the vector of elemental variable values, one variable and one subdomain at a time.

Definition at line 2596 of file nemesis_io_helper.C.

 {
   // For each variable in names,
   //   For each subdomain in subdomain_map,
   //     If this (subdomain, variable) combination is active
   //       Extract our values into local_soln (localize is a collective)
   //       Write local_soln to file
   for (auto v : index_range(var_nums))
     {
       // Get list of active subdomains for variable v
       const auto & active_subdomains = vars_active_subdomains[v];
  
       for (const int sbd_id_int : global_elem_blk_ids)
         {
           const subdomain_id_type sbd_id =
             cast_int<subdomain_id_type>(sbd_id_int);
           auto it = subdomain_map.find(sbd_id);
           const std::vector<dof_id_type> empty_vec;
           const std::vector<dof_id_type> & elem_ids =
             (it == subdomain_map.end()) ? empty_vec : it->second;
  
           // Possibly skip this (variable, subdomain) combination
           if (active_subdomains.empty() || active_subdomains.count(sbd_id))
             {
               std::vector<numeric_index_type> required_indices;
               required_indices.reserve(elem_ids.size());
  
               const unsigned int sys_num = var_nums[v].first;
  
               for (const auto & id : elem_ids)
                 required_indices.push_back
                   (mesh.elem_ref(id).dof_number
                     (sys_num, var_nums[v].second, 0));
  
               std::vector<Number> local_soln;
               es.get_system(sys_num).current_local_solution->get
                 (required_indices, local_soln);
  
               // It's possible that there's nothing for us to write:
               // we may not be responsible for any elements on the
               // current subdomain.  We did still have to participate
               // in the localize() call above, however, since it is a
               // collective.
               if (local_soln.size())
                 {
                   ex_err = exII::ex_put_elem_var(ex_id,
                                                  timestep,
                                                  static_cast<int>(v+1),
                                                  static_cast<int>(sbd_id),
                                                  static_cast<int>(local_soln.size()),
                                                  local_soln.data());
                   EX_CHECK_ERR(ex_err, "Error writing element values.");
                 }
             }
         }
     } // end loop over vars
  
   ex_err = exII::ex_update(ex_id);
   EX_CHECK_ERR(ex_err, "Error flushing buffers to file.");
 }

References libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::EquationSystems::get_system(), global_elem_blk_ids, libMesh::if(), libMesh::index_range(), mesh, and subdomain_map.

◆ write_element_values() [2/2]

void libMesh::ExodusII_IO_Helper::write_element_values	(	const MeshBase &	mesh,
		const std::vector< Real > &	values,
		int	timestep,
		const std::vector< std::set< subdomain_id_type >> &	vars_active_subdomains
	)

inherited

Writes the vector of values to the element variables.

The 'values' vector is assumed to be in the order: {(u1, u2, u3, ..., uN), (v1, v2, v3, ..., vN), (w1, w2, w3, ..., wN)} where N is the number of elements.

This ordering is produced by calls to ES::build_elemental_solution_vector(). ES::build_discontinuous_solution_vector(), on the other hand, produces an element-major ordering. See the function below for that case.

Definition at line 2854 of file exodusII_io_helper.C.

 {
   if ((_run_only_on_proc0) && (this->processor_id() != 0))
     return;
  
   // Ask the file how many element vars it has, store it in the num_elem_vars variable.
   ex_err = exII::ex_get_var_param(ex_id, "e", &num_elem_vars);
   EX_CHECK_ERR(ex_err, "Error reading number of elemental variables.");
  
   // We will eventually loop over the element blocks (subdomains) and
   // write the data one block at a time. Build a data structure that
   // maps each subdomain to a list of element ids it contains.
   std::map<subdomain_id_type, std::vector<unsigned int>> subdomain_map;
   for (const auto & elem : mesh.active_element_ptr_range())
     subdomain_map[elem->subdomain_id()].push_back(elem->id());
  
   // Use mesh.n_elem() to access into the values vector rather than
   // the number of elements the Exodus writer thinks the mesh has,
   // which may not include inactive elements.
   dof_id_type n_elem = mesh.n_elem();
  
   // Sanity check: we must have an entry in vars_active_subdomains for
   // each variable that we are potentially writing out.
   libmesh_assert_equal_to
     (vars_active_subdomains.size(),
      static_cast<unsigned>(num_elem_vars));
  
   // For each variable, create a 'data' array which holds all the elemental variable
   // values *for a given block* on this processor, then write that data vector to file
   // before moving onto the next block.
   for (unsigned int var_id=0; var_id<static_cast<unsigned>(num_elem_vars); ++var_id)
     {
       // The size of the subdomain map is the number of blocks.
       auto it = subdomain_map.begin();
  
       // Reference to the set of active subdomains for the current variable.
       const auto & active_subdomains
         = vars_active_subdomains[var_id];
  
       for (unsigned int j=0; it!=subdomain_map.end(); ++it, ++j)
         {
           // Skip any variable/subdomain pairs that are inactive.
           // Note that if active_subdomains is empty, it is interpreted
           // as being active on *all* subdomains.
           if (!(active_subdomains.empty() || active_subdomains.count(it->first)))
             continue;
  
           // Get reference to list of elem ids which are in the
           // current subdomain and count, allocate storage to hold
           // data that will be written to file.
           const auto & elem_nums = it->second;
           const unsigned int num_elems_this_block =
             cast_int<unsigned int>(elem_nums.size());
           std::vector<Real> data(num_elems_this_block);
  
           // variable-major ordering is:
           // (u1, u2, u3, ..., uN), (v1, v2, v3, ..., vN), ...
           // where N is the number of elements.
           for (unsigned int k=0; k<num_elems_this_block; ++k)
             data[k] = values[var_id*n_elem + elem_nums[k]];
  
           ex_err = exII::ex_put_elem_var
             (ex_id,
              timestep,
              var_id+1,
              this->get_block_id(j),
              num_elems_this_block,
              MappedOutputVector(data, _single_precision).data());
  
           EX_CHECK_ERR(ex_err, "Error writing element values.");
         }
     }
  
   ex_err = exII::ex_update(ex_id);
   EX_CHECK_ERR(ex_err, "Error flushing buffers to file.");
 }

References libMesh::ExodusII_IO_Helper::_run_only_on_proc0, libMesh::ExodusII_IO_Helper::_single_precision, data, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::get_block_id(), mesh, libMesh::MeshTools::n_elem(), libMesh::ExodusII_IO_Helper::num_elem_vars, and libMesh::ParallelObject::processor_id().

◆ write_element_values_element_major()

void libMesh::ExodusII_IO_Helper::write_element_values_element_major	(	const MeshBase &	mesh,
		const std::vector< Real > &	values,
		int	timestep,
		const std::vector< std::set< subdomain_id_type >> &	vars_active_subdomains,
		const std::vector< std::string > &	derived_var_names,
		const std::map< subdomain_id_type, std::vector< std::string >> &	subdomain_to_var_names
	)

inherited

Same as the function above, but assume the input 'values' vector is in element-major order, i.e.

{(u1,v1,w1), (u2,v2,w2), ... (uN,vN,wN)} This function is called by ExodusII_IO::write_element_data_from_discontinuous_nodal_data() because ES::build_discontinuous_solution_vector() builds the solution vector in this order.

Note: If some variables are subdomain-restricted, then the tuples will be of different lengths for each element, i.e. {(u1,v1,w1), (u2,v2), ... (uN,vN,wN)} if variable w is not active on element 2.

Definition at line 2937 of file exodusII_io_helper.C.

 {
   if ((_run_only_on_proc0) && (this->processor_id() != 0))
     return;
  
   // Ask the file how many element vars it has, store it in the num_elem_vars variable.
   ex_err = exII::ex_get_var_param(ex_id, "e", &num_elem_vars);
   EX_CHECK_ERR(ex_err, "Error reading number of elemental variables.");
  
   // We will eventually loop over the element blocks (subdomains) and
   // write the data one block (subdomain) at a time. Build a data
   // structure that keeps track of how many elements are in each
   // subdomain. This will allow us to reserve space in the data vector
   // we are going to write.
   std::map<subdomain_id_type, unsigned int> subdomain_to_n_elem;
   for (const auto & elem : mesh.active_element_ptr_range())
     subdomain_to_n_elem[elem->subdomain_id()] += 1;
  
   // Sanity check: we must have an entry in vars_active_subdomains for
   // each variable that we are potentially writing out.
   libmesh_assert_equal_to
     (vars_active_subdomains.size(),
      static_cast<unsigned>(num_elem_vars));
  
   // The size of the subdomain map is the number of blocks.
   auto subdomain_to_n_elem_iter = subdomain_to_n_elem.begin();
  
   // Store range of active Elem pointers. We are going to loop over
   // the elements n_vars * n_subdomains times, so let's make sure
   // the predicated iterators aren't slowing us down too much.
   ConstElemRange elem_range
     (mesh.active_elements_begin(),
      mesh.active_elements_end());
  
   for (unsigned int sbd_idx=0;
        subdomain_to_n_elem_iter != subdomain_to_n_elem.end();
        ++subdomain_to_n_elem_iter, ++sbd_idx)
     for (unsigned int var_id=0; var_id<static_cast<unsigned>(num_elem_vars); ++var_id)
       {
         // Reference to the set of active subdomains for the current variable.
         const auto & active_subdomains
           = vars_active_subdomains[var_id];
  
         // If the vars_active_subdomains container passed to this function
         // has an empty entry, it means the variable really is not active on
         // _any_ subdomains, not that it is active on _all_ subdomains. This
         // is just due to the way that we build the vars_active_subdomains
         // container.
         if (!active_subdomains.count(subdomain_to_n_elem_iter->first))
           continue;
  
         // Vector to hold values that will be written to Exodus file.
         std::vector<Real> data;
         data.reserve(subdomain_to_n_elem_iter->second);
  
         unsigned int values_offset = 0;
         for (auto & elem : elem_range)
           {
             // We'll use the Elem's subdomain id in several places below.
             subdomain_id_type sbd_id = elem->subdomain_id();
  
             // Get reference to the list of variable names defining
             // the indexing for the current Elem's subdomain.
             auto subdomain_to_var_names_iter =
               subdomain_to_var_names.find(sbd_id);
  
             // It's possible, but unusual, for there to be an Elem
             // from a subdomain that has no active variables from the
             // set of variables we are currently writing. If that
             // happens, we can just go to the next Elem because we
             // don't need to advance the offset into the values
             // vector, etc.
             if (subdomain_to_var_names_iter == subdomain_to_var_names.end())
               continue;
  
             const auto & var_names_this_sbd
               = subdomain_to_var_names_iter->second;
  
             // Only extract values if Elem is in the current subdomain.
             if (sbd_id == subdomain_to_n_elem_iter->first)
               {
                 // Location of current var_id in the list of all variables on this
                 // subdomain. FIXME: linear search but it's over a typically relatively
                 // short vector of active variable names on this subdomain. We could do
                 // a nested std::map<string,index> instead of a std::vector where the
                 // location of the string is implicitly the index..
                 auto pos =
                   std::find(var_names_this_sbd.begin(),
                             var_names_this_sbd.end(),
                             derived_var_names[var_id]);
  
                 if (pos == var_names_this_sbd.end())
                   libmesh_error_msg("Derived name " << derived_var_names[var_id] << " not found!");
  
                 // Find the current variable's location in the list of all variable
                 // names on the current Elem's subdomain.
                 auto true_index =
                   std::distance(var_names_this_sbd.begin(), pos);
  
                 data.push_back(values[values_offset + true_index]);
               }
  
             // The "true" offset is how much we have to advance the index for each Elem
             // in this subdomain.
             auto true_offset = var_names_this_sbd.size();
  
             // Increment to the next Elem's values
             values_offset += true_offset;
           } // for elem
  
         // Now write 'data' to Exodus file, in single precision if requested.
         if (!data.empty())
           {
             ex_err = exII::ex_put_elem_var
               (ex_id, timestep, var_id+1, this->get_block_id(sbd_idx), data.size(),
                MappedOutputVector(data, _single_precision).data());
  
             EX_CHECK_ERR(ex_err, "Error writing element values.");
           }
       } // for each var_id
  
   ex_err = exII::ex_update(ex_id);
   EX_CHECK_ERR(ex_err, "Error flushing buffers to file.");
 }

References libMesh::ExodusII_IO_Helper::_run_only_on_proc0, libMesh::ExodusII_IO_Helper::_single_precision, data, libMesh::ExodusII_IO_Helper::MappedOutputVector::data(), distance(), libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::get_block_id(), mesh, libMesh::ExodusII_IO_Helper::num_elem_vars, and libMesh::ParallelObject::processor_id().

◆ write_elements()

void libMesh::Nemesis_IO_Helper::write_elements	(	const MeshBase &	mesh,
		bool	use_discontinuous = `false`
	)

overridevirtual

This function is specialized to write the connectivity.

Reimplemented from libMesh::ExodusII_IO_Helper.

Definition at line 2287 of file nemesis_io_helper.C.

 {
   // Only write elements if there are elements blocks available.
   if (this->num_elem_blks_global > 0)
     {
       // Data structure to store element block names that will be used to
       // write the element block names to file.
       NamesData names_table(this->num_elem_blks_global, MAX_STR_LENGTH);
  
       // Loop over all blocks, even if we don't have elements in each block.
       // If we don't have elements we need to write out a 0 for that block...
       for (auto i : IntRange<int>(0, this->num_elem_blks_global))
         {
           // Even if there are no elements for this block on the current
           // processor, we still want to write its name to file, if
           // possible. MeshBase::subdomain_name() will just return an
           // empty string if there is no name associated with the current
           // block.
           names_table.push_back_entry
             (mesh.subdomain_name(cast_int<subdomain_id_type>(this->global_elem_blk_ids[i])));
  
           // Search for the current global block ID in the map
           std::map<int, std::vector<int>>::iterator it =
             this->block_id_to_elem_connectivity.find( this->global_elem_blk_ids[i] );
  
           // If not found, write a zero to file....
           if (it == this->block_id_to_elem_connectivity.end())
             {
               this->ex_err = exII::ex_put_elem_block(this->ex_id,
                                                      this->global_elem_blk_ids[i],
                                                      "Empty",
                                                      0, /* n. elements in this block */
                                                      0, /* n. nodes per element */
                                                      0);  /* number of attributes per element */
  
               EX_CHECK_ERR(this->ex_err, "Error writing element block from Nemesis.");
             }
  
           // Otherwise, write the actual block information and connectivity to file
           else
             {
               subdomain_id_type block =
                 cast_int<subdomain_id_type>(it->first);
               std::vector<int> & this_block_connectivity = it->second;
               std::vector<dof_id_type> & elements_in_this_block = subdomain_map[block];
  
               // Use the first element in this block to get representative information.
               // Note that Exodus assumes all elements in a block are of the same type!
               // We are using that same assumption here!
               const auto & conv =
                 get_conversion(mesh.elem_ref(elements_in_this_block[0]).type());
  
               this->num_nodes_per_elem =
                 mesh.elem_ref(elements_in_this_block[0]).n_nodes();
  
               ex_err = exII::ex_put_elem_block(ex_id,
                                                block,
                                                conv.exodus_elem_type().c_str(),
                                                elements_in_this_block.size(),
                                                num_nodes_per_elem,
                                                0);
               EX_CHECK_ERR(ex_err, "Error writing element block from Nemesis.");
  
               ex_err = exII::ex_put_elem_conn(ex_id,
                                               block,
                                               this_block_connectivity.data());
               EX_CHECK_ERR(ex_err, "Error writing element connectivities from Nemesis.");
             }
         } // end loop over global block IDs
  
       // Only call this once, not in the loop above!
       ex_err = exII::ex_put_elem_num_map(ex_id,
                                          exodus_elem_num_to_libmesh.empty() ? nullptr : exodus_elem_num_to_libmesh.data());
       EX_CHECK_ERR(ex_err, "Error writing element map");
  
       // Write the element block names to file.
       ex_err = exII::ex_put_names(ex_id, exII::EX_ELEM_BLOCK, names_table.get_char_star_star());
       EX_CHECK_ERR(ex_err, "Error writing element block names");
     } // end if (this->num_elem_blks_global > 0)
 }

References block_id_to_elem_connectivity, libMesh::MeshBase::elem_ref(), libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::exodus_elem_num_to_libmesh, libMesh::ExodusII_IO_Helper::NamesData::get_char_star_star(), libMesh::ExodusII_IO_Helper::get_conversion(), global_elem_blk_ids, mesh, libMesh::Elem::n_nodes(), num_elem_blks_global, libMesh::ExodusII_IO_Helper::num_nodes_per_elem, libMesh::ExodusII_IO_Helper::NamesData::push_back_entry(), and subdomain_map.

◆ write_exodus_initialization_info()

void libMesh::Nemesis_IO_Helper::write_exodus_initialization_info	(	const MeshBase &	pmesh,
		const std::string &	title
	)

private

This function writes exodus-specific initialization information.

This information is slightly different when you are working with Nemesis, as it depends on some global information being known.

Definition at line 933 of file nemesis_io_helper.C.

 {
   // This follows the convention of Exodus: we always write out the mesh as LIBMESH_DIM-dimensional,
   // even if it is 2D...
   this->num_dim = LIBMESH_DIM;
  
   this->num_elem = static_cast<unsigned int>(std::distance (pmesh.active_local_elements_begin(),
                                                             pmesh.active_local_elements_end()));
  
   // Exodus will also use *global* number of side and node sets,
   // though it will not write out entries for all of them...
   this->num_side_sets =
     cast_int<int>(this->global_sideset_ids.size());
   this->num_node_sets =
     cast_int<int>(this->global_nodeset_ids.size());
  
   // We need to write the global number of blocks, even though this processor might not have
   // elements in some of them!
   this->num_elem_blk = this->num_elem_blks_global;
  
   ex_err = exII::ex_put_init(ex_id,
                              title_in.c_str(),
                              this->num_dim,
                              this->num_nodes,
                              this->num_elem,
                              this->num_elem_blk,
                              this->num_node_sets,
                              this->num_side_sets);
  
   EX_CHECK_ERR(ex_err, "Error initializing new Nemesis file.");
 }

Referenced by initialize().

◆ write_global_values()

void libMesh::ExodusII_IO_Helper::write_global_values	(	const std::vector< Real > &	values,
		int	timestep
	)

inherited

Writes the vector of global variables.

Definition at line 3131 of file exodusII_io_helper.C.

 {
   if ((_run_only_on_proc0) && (this->processor_id() != 0))
     return;
  
   if (!values.empty())
     {
       ex_err = exII::ex_put_glob_vars
         (ex_id, timestep, num_global_vars,
          MappedOutputVector(values, _single_precision).data());
  
       EX_CHECK_ERR(ex_err, "Error writing global values.");
  
       ex_err = exII::ex_update(ex_id);
       EX_CHECK_ERR(ex_err, "Error flushing buffers to file.");
     }
 }

References libMesh::ExodusII_IO_Helper::_run_only_on_proc0, libMesh::ExodusII_IO_Helper::_single_precision, data, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::num_global_vars, and libMesh::ParallelObject::processor_id().

◆ write_information_records()

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

inherited

Writes the vector of information records.

Definition at line 3092 of file exodusII_io_helper.C.

 {
   if ((_run_only_on_proc0) && (this->processor_id() != 0))
     return;
  
   // There may already be information records in the file (for
   // example, if we're appending) and in that case, according to the
   // Exodus documentation, writing more information records is not
   // supported.
   int num_info = inquire(exII::EX_INQ_INFO, "Error retrieving the number of information records from file!");
   if (num_info > 0)
     {
       libMesh::err << "Warning! The Exodus file already contains information records.\n"
                    << "Exodus does not support writing additional records in this situation."
                    << std::endl;
       return;
     }
  
   int num_records = cast_int<int>(records.size());
  
   if (num_records > 0)
     {
       NamesData info(num_records, MAX_LINE_LENGTH);
  
       // If an entry is longer than MAX_LINE_LENGTH characters it's not an error, we just
       // write the first MAX_LINE_LENGTH characters to the file.
       for (const auto & record : records)
         info.push_back_entry(record);
  
       ex_err = exII::ex_put_info(ex_id, num_records, info.get_char_star_star());
       EX_CHECK_ERR(ex_err, "Error writing global values.");
  
       ex_err = exII::ex_update(ex_id);
       EX_CHECK_ERR(ex_err, "Error flushing buffers to file.");
     }
 }

References libMesh::ExodusII_IO_Helper::_run_only_on_proc0, libMesh::err, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::NamesData::get_char_star_star(), libMesh::ExodusII_IO_Helper::inquire(), libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::NamesData::push_back_entry().

◆ write_nodal_coordinates()

void libMesh::Nemesis_IO_Helper::write_nodal_coordinates	(	const MeshBase &	mesh,
		bool	use_discontinuous = `false`
	)

overridevirtual

This function is specialized from ExodusII_IO_Helper to write only the nodal coordinates stored on the local piece of the Mesh.

Reimplemented from libMesh::ExodusII_IO_Helper.

Definition at line 2228 of file nemesis_io_helper.C.

 {
   auto local_num_nodes = this->exodus_node_num_to_libmesh.size();
  
   x.resize(local_num_nodes);
   y.resize(local_num_nodes);
   z.resize(local_num_nodes);
  
   // Just loop over our list outputting the nodes the way we built the map
   for (auto i : IntRange<int>(0, local_num_nodes))
     {
       const Point & pt = mesh.point(this->exodus_node_num_to_libmesh[i]);
       x[i]=pt(0);
       y[i]=pt(1);
       z[i]=pt(2);
     }
  
   if (local_num_nodes)
     {
       if (_single_precision)
         {
           std::vector<float>
             x_single(x.begin(), x.end()),
             y_single(y.begin(), y.end()),
             z_single(z.begin(), z.end());
  
           ex_err = exII::ex_put_coord(ex_id,
                                       x_single.data(),
                                       y_single.data(),
                                       z_single.data());
         }
       else
         {
           // Call Exodus API to write nodal coordinates...
           ex_err = exII::ex_put_coord(ex_id,
                                       x.data(),
                                       y.data(),
                                       z.data());
         }
       EX_CHECK_ERR(ex_err, "Error writing node coordinates");
  
       // And write the nodal map we created for them
       ex_err = exII::ex_put_node_num_map(ex_id, this->exodus_node_num_to_libmesh.data());
       EX_CHECK_ERR(ex_err, "Error writing node num map");
     }
   else // Does the Exodus API want us to write empty nodal coordinates?
     {
       ex_err = exII::ex_put_coord(ex_id, nullptr, nullptr, nullptr);
       EX_CHECK_ERR(ex_err, "Error writing empty node coordinates");
  
       ex_err = exII::ex_put_node_num_map(ex_id, nullptr);
       EX_CHECK_ERR(ex_err, "Error writing empty node num map");
     }
 }

References libMesh::ExodusII_IO_Helper::_single_precision, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::exodus_node_num_to_libmesh, mesh, libMesh::ExodusII_IO_Helper::x, libMesh::ExodusII_IO_Helper::y, and libMesh::ExodusII_IO_Helper::z.

◆ write_nodal_solution() [1/3]

void libMesh::Nemesis_IO_Helper::write_nodal_solution	(	const EquationSystems &	es,
		const std::vector< std::pair< unsigned int, unsigned int >> &	var_nums,
		int	timestep,
		const std::vector< std::string > &	output_names
	)

Outputs EquationSystems current_local_solution nodal values.

Definition at line 2470 of file nemesis_io_helper.C.

 {
   const MeshBase & mesh = es.get_mesh();
  
   for (auto & var_num : var_nums)
     {
       const unsigned int sys_num = var_num.first;
       const unsigned int var = var_num.second;
  
       const System & sys = es.get_system(sys_num);
       const std::string & name = sys.variable_name(var);
  
       auto pos = std::find(output_names.begin(), output_names.end(), name);
  
       // Skip this name if it's not supposed to be output.
       if (pos == output_names.end())
         continue;
  
       // Compute the (zero-based) index which determines which
       // variable this will be as far as Nemesis is concerned.  This
       // will be used below in the write_nodal_values() call.
       int variable_name_position =
         cast_int<int>(std::distance(output_names.begin(), pos));
  
       // Fill up a std::vector with the dofs for the current variable
       std::vector<numeric_index_type> required_indices(num_nodes);
  
       const FEType type = sys.variable_type(var);
       if (type.family == SCALAR)
         {
           std::vector<numeric_index_type> scalar_indices;
           sys.get_dof_map().SCALAR_dof_indices(scalar_indices, var);
           for (int i=0; i<num_nodes; i++)
             required_indices[i] = scalar_indices[0];
         }
       else
         for (int i=0; i<num_nodes; i++)
           {
             const Node & node = mesh.node_ref(this->exodus_node_num_to_libmesh[i]);
             required_indices[i] = node.dof_number(sys_num, var, 0);
           }
  
       // Get the dof values required to write just our local part of
       // the solution vector.
       std::vector<Number> local_soln;
       sys.current_local_solution->get(required_indices, local_soln);
  
 #ifndef LIBMESH_USE_COMPLEX_NUMBERS
       // Call the ExodusII_IO_Helper function to write the data.
       write_nodal_values(variable_name_position + 1, local_soln, timestep);
 #else
       // We have the local (complex) values. Now extract the real,
       // imaginary, and magnitude values from them.
       std::vector<Real> real_parts(num_nodes);
       std::vector<Real> imag_parts(num_nodes);
       std::vector<Real> magnitudes(num_nodes);
  
       for (int i=0; i<num_nodes; ++i)
         {
           real_parts[i] = local_soln[i].real();
           imag_parts[i] = local_soln[i].imag();
           magnitudes[i] = std::abs(local_soln[i]);
         }
  
       // Write the real, imaginary, and magnitude values to file.
       write_nodal_values(3 * variable_name_position + 1, real_parts, timestep);
       write_nodal_values(3 * variable_name_position + 2, imag_parts, timestep);
       write_nodal_values(3 * variable_name_position + 3, magnitudes, timestep);
 #endif
     }
 }

References std::abs(), libMesh::System::current_local_solution, distance(), libMesh::DofObject::dof_number(), libMesh::ExodusII_IO_Helper::exodus_node_num_to_libmesh, libMesh::FEType::family, libMesh::System::get_dof_map(), libMesh::EquationSystems::get_mesh(), libMesh::EquationSystems::get_system(), mesh, libMesh::Quality::name(), libMesh::ExodusII_IO_Helper::num_nodes, libMesh::SCALAR, libMesh::DofMap::SCALAR_dof_indices(), libMesh::System::variable_name(), libMesh::System::variable_type(), and libMesh::ExodusII_IO_Helper::write_nodal_values().

◆ write_nodal_solution() [2/3]

void libMesh::Nemesis_IO_Helper::write_nodal_solution	(	const NumericVector< Number > &	parallel_soln,
		const std::vector< std::string > &	names,
		int	timestep,
		const std::vector< std::string > &	output_names
	)

Takes a parallel solution vector containing the node-major solution vector for all variables and outputs it to the files.

Parameters

parallel_soln
names	A vector containing the names of all variables in parallel_soln.
timestep	To be passed to the ExodusII_IO_Helper::write_nodal_values() function.
output_names	A vector containing the names of variables in parallel_soln that should actually be written (whitelist).

Note: This version of write_nodal_solution() is called by the parallel version of Nemesis_IO::write_nodal_data(), which is called by MeshOutput::write_equation_systems() for parallel I/O formats like Nemesis. The other version is still available to continue supporting things like NamebasedIO::write_nodal_data(), but this version should be preferred when running in parallel.

Definition at line 2410 of file nemesis_io_helper.C.

 {
   int num_vars = cast_int<int>(names.size());
  
   for (int c=0; c<num_vars; c++)
     {
       // Find the position of names[c] in the output_names vector, if it exists.
       auto pos = std::find(output_names.begin(), output_names.end(), names[c]);
  
       // Skip names[c] if it's not supposed to be output.
       if (pos == output_names.end())
         continue;
  
       // Compute the (zero-based) index which determines which
       // variable this will be as far as Nemesis is concerned.  This
       // will be used below in the write_nodal_values() call.
       int variable_name_position =
         cast_int<int>(std::distance(output_names.begin(), pos));
  
       // Fill up a std::vector with the dofs for the current variable
       std::vector<numeric_index_type> required_indices(num_nodes);
  
       for (int i=0; i<num_nodes; i++)
         required_indices[i] = static_cast<dof_id_type>(this->exodus_node_num_to_libmesh[i]) * num_vars + c;
  
       // Get the dof values required to write just our local part of
       // the solution vector.
       std::vector<Number> local_soln;
       parallel_soln.localize(local_soln, required_indices);
  
 #ifndef LIBMESH_USE_COMPLEX_NUMBERS
       // Call the ExodusII_IO_Helper function to write the data.
       write_nodal_values(variable_name_position + 1, local_soln, timestep);
 #else
       // We have the local (complex) values. Now extract the real,
       // imaginary, and magnitude values from them.
       std::vector<Real> real_parts(num_nodes);
       std::vector<Real> imag_parts(num_nodes);
       std::vector<Real> magnitudes(num_nodes);
  
       for (int i=0; i<num_nodes; ++i)
         {
           real_parts[i] = local_soln[i].real();
           imag_parts[i] = local_soln[i].imag();
           magnitudes[i] = std::abs(local_soln[i]);
         }
  
       // Write the real, imaginary, and magnitude values to file.
       write_nodal_values(3 * variable_name_position + 1, real_parts, timestep);
       write_nodal_values(3 * variable_name_position + 2, imag_parts, timestep);
       write_nodal_values(3 * variable_name_position + 3, magnitudes, timestep);
 #endif
     }
 }

References std::abs(), distance(), libMesh::ExodusII_IO_Helper::exodus_node_num_to_libmesh, libMesh::NumericVector< T >::localize(), libMesh::ExodusII_IO_Helper::num_nodes, and libMesh::ExodusII_IO_Helper::write_nodal_values().

◆ write_nodal_solution() [3/3]

void libMesh::Nemesis_IO_Helper::write_nodal_solution	(	const std::vector< Number > &	values,
		const std::vector< std::string > &	names,
		int	timestep
	)

Takes a solution vector containing the solution for all variables and outputs it to the files.

Definition at line 2372 of file nemesis_io_helper.C.

 {
   int num_vars = cast_int<int>(names.size());
   //int num_values = values.size(); // Not used?
  
   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(num_nodes);
  
       for (int i=0; i<num_nodes; ++i)
         {
           Number value = values[this->exodus_node_num_to_libmesh[i]*num_vars + c];
           real_parts[i] = value.real();
           imag_parts[i] = value.imag();
           magnitudes[i] = std::abs(value);
         }
       write_nodal_values(3*c+1,real_parts,timestep);
       write_nodal_values(3*c+2,imag_parts,timestep);
       write_nodal_values(3*c+3,magnitudes,timestep);
 #else
       std::vector<Number> cur_soln(num_nodes);
  
       // Copy out this variable's solution
       for (int i=0; i<num_nodes; i++)
         cur_soln[i] = values[this->exodus_node_num_to_libmesh[i]*num_vars + c];
  
       write_nodal_values(c+1,cur_soln,timestep);
 #endif
     }
 }

References std::abs(), libMesh::ExodusII_IO_Helper::exodus_node_num_to_libmesh, libMesh::ExodusII_IO_Helper::num_nodes, value, and libMesh::ExodusII_IO_Helper::write_nodal_values().

◆ write_nodal_values()

void libMesh::ExodusII_IO_Helper::write_nodal_values	(	int	var_id,
		const std::vector< Real > &	values,
		int	timestep
	)

inherited

Writes the vector of values to a nodal variable.

Definition at line 3070 of file exodusII_io_helper.C.

 {
   if ((_run_only_on_proc0) && (this->processor_id() != 0))
     return;
  
   if (!values.empty())
     {
       ex_err = exII::ex_put_nodal_var
         (ex_id, timestep, var_id, num_nodes,
          MappedOutputVector(values, _single_precision).data());
  
       EX_CHECK_ERR(ex_err, "Error writing nodal values.");
  
       ex_err = exII::ex_update(ex_id);
       EX_CHECK_ERR(ex_err, "Error flushing buffers to file.");
     }
 }

References libMesh::ExodusII_IO_Helper::_run_only_on_proc0, libMesh::ExodusII_IO_Helper::_single_precision, data, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::num_nodes, and libMesh::ParallelObject::processor_id().

Referenced by write_nodal_solution().

◆ write_nodesets()

void libMesh::Nemesis_IO_Helper::write_nodesets ( const MeshBase & mesh )

overridevirtual

Writes the nodesets for this processor.

Reimplemented from libMesh::ExodusII_IO_Helper.

Definition at line 1933 of file nemesis_io_helper.C.

 {
   // Write the nodesets.  In Nemesis, the idea is to "create space" for the global
   // set of boundary nodesets, but to only write node IDs which are local to the current
   // processor.  This is what is done in Nemesis files created by the "loadbal" script.
  
   // Store a map of vectors for boundary node IDs on this processor.
   // Use a vector of int here so it can be passed directly to Exodus.
   std::map<boundary_id_type, std::vector<int>> local_node_boundary_id_lists;
  
   // FIXME: We should build this list only one time!!  We already built it above, but we
   // did not have the libmesh to exodus node mapping at that time... for now we'll just
   // build it here again, hopefully it's small relative to the size of the entire mesh.
  
   // Build list of (node-id, bc-id) tuples.
   typedef std::tuple<dof_id_type, boundary_id_type> Tuple;
   std::vector<Tuple> bc_tuples = mesh.get_boundary_info().build_node_list();
  
   if (verbose)
     {
       libMesh::out << "[" << this->processor_id() << "] boundary_node_list.size()="
                    << bc_tuples.size() << std::endl;
       libMesh::out << "[" << this->processor_id() << "] (boundary_node_id, boundary_id) = ";
       for (const auto & t : bc_tuples)
         libMesh::out << "(" << std::get<0>(t) << ", " << std::get<1>(t) << ") ";
       libMesh::out << std::endl;
     }
  
   // For each node in the node list, add it to the vector of node IDs for that
   // set for the local processor.  This will be used later when writing Exodus
   // nodesets.
   for (const auto & t : bc_tuples)
     {
       // Don't try to grab a reference to the vector unless the current node is attached
       // to a local element.  Otherwise, another processor will be responsible for writing it in its nodeset.
       std::map<int, int>::iterator it = this->libmesh_node_num_to_exodus.find(std::get<0>(t));
  
       if (it != this->libmesh_node_num_to_exodus.end())
         {
           // Get reference to the vector where this node ID will be inserted.  If it
           // doesn't yet exist, this will create it.
           std::vector<int> & current_id_set = local_node_boundary_id_lists[std::get<1>(t)];
  
           // Push back Exodus-mapped node ID for this set
           // TODO: reserve space in these vectors somehow.
           current_id_set.push_back( it->second );
         }
     }
  
   // See what we got
   if (verbose)
     {
       for (const auto & pr : local_node_boundary_id_lists)
         {
           libMesh::out << "[" << this->processor_id() << "] ID: " << pr.first << ", ";
  
           // Libmesh node ID (Exodus Node ID)
           for (const auto & id : pr.second)
             libMesh::out << id << ", ";
           libMesh::out << std::endl;
         }
     }
  
   // Loop over *global* nodeset IDs, call the Exodus API.  Note that some nodesets may be empty
   // for a given processor.
   if (global_nodeset_ids.size() > 0)
     {
       NamesData names_table(global_nodeset_ids.size(), MAX_STR_LENGTH);
  
       for (const auto & nodeset_id : this->global_nodeset_ids)
         {
           const std::string & current_ns_name =
             mesh.get_boundary_info().get_nodeset_name
               (cast_int<boundary_id_type>(nodeset_id));
  
           // Store this name in a data structure that will be used to
           // write sideset names to file.
           names_table.push_back_entry(current_ns_name);
  
           if (verbose)
             {
               libMesh::out << "[" << this->processor_id()
                            << "] Writing out Exodus nodeset info for ID: " << nodeset_id
                            << ", Name: " << current_ns_name
                            << std::endl;
             }
  
           // Convert current global_nodeset_id into an exodus ID, which can't be zero...
           int exodus_id = nodeset_id;
  
           /*
           // Exodus can't handle zero nodeset IDs (?)  Use max short here since
           // when libmesh reads it back in, it will want to store it as a short...
           if (exodus_id==0)
           exodus_id = std::numeric_limits<short>::max();
           */
  
           // Try to find this boundary ID in the local list we created
           auto it =
             local_node_boundary_id_lists.find (cast_int<boundary_id_type>(nodeset_id));
  
           // No nodes found for this boundary ID on this processor
           if (it == local_node_boundary_id_lists.end())
             {
               if (verbose)
                 libMesh::out << "[" << this->processor_id()
                              << "] No nodeset data for ID: " << nodeset_id
                              << " on this processor." << std::endl;
  
               // Call the Exodus interface to write the parameters of this node set
               this->ex_err = exII::ex_put_node_set_param(this->ex_id,
                                                          exodus_id,
                                                          0, /* No nodes for this ID */
                                                          0  /* No distribution factors */);
               EX_CHECK_ERR(this->ex_err, "Error writing nodeset parameters in Nemesis");
  
             }
           else // Boundary ID *was* found in list
             {
               // Get reference to the vector of node IDs
               std::vector<int> & current_nodeset_ids = it->second;
  
               // Call the Exodus interface to write the parameters of this node set
               this->ex_err = exII::ex_put_node_set_param(this->ex_id,
                                                          exodus_id,
                                                          current_nodeset_ids.size(),
                                                          0  /* No distribution factors */);
  
               EX_CHECK_ERR(this->ex_err, "Error writing nodeset parameters in Nemesis");
  
               // Call Exodus interface to write the actual node IDs for this boundary ID
               this->ex_err = exII::ex_put_node_set(this->ex_id,
                                                    exodus_id,
                                                    current_nodeset_ids.data());
  
               EX_CHECK_ERR(this->ex_err, "Error writing nodesets in Nemesis");
  
             }
         } // end loop over global nodeset IDs
  
       // Write out the nodeset names
       ex_err = exII::ex_put_names(ex_id,
                                   exII::EX_NODE_SET,
                                   names_table.get_char_star_star());
       EX_CHECK_ERR(ex_err, "Error writing nodeset names");
     } // end for loop over global nodeset IDs
 }

References libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, global_nodeset_ids, libMesh::ExodusII_IO_Helper::libmesh_node_num_to_exodus, mesh, libMesh::out, libMesh::ParallelObject::processor_id(), and libMesh::ExodusII_IO_Helper::verbose.

◆ write_sideset_data()

void libMesh::ExodusII_IO_Helper::write_sideset_data	(	const MeshBase &	mesh,
		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
	)

inherited

Write sideset data for the requested timestep.

Definition at line 2587 of file exodusII_io_helper.C.

 {
   if ((_run_only_on_proc0) && (this->processor_id() != 0))
     return;
  
   // Write the sideset variable names to file. This function should
   // only be called once for SIDESET variables, repeated calls to
   // write_var_names overwrites/changes the order of names that were
   // there previously, and will mess up any data that has already been
   // written.
   this->write_var_names(SIDESET, var_names);
  
   // I hope that we are allowed to call read_sideset_info() even
   // though we are in the middle of writing? It seems to work provided
   // that you have already written the mesh itself... read_sideset_info()
   // fills in the following data members:
   // .) num_side_sets
   // .) ss_ids
   this->read_sideset_info();
  
   // Debugging:
   // libMesh::out << "File has " << num_side_sets << " side sets." << std::endl;
  
   // Write "truth" table for sideset variables.  The function
   // exII::ex_put_var_param() must be called before
   // exII::ex_put_sset_var_tab(). For us, this happens during the call
   // to ExodusII_IO_Helper::write_var_names(). sset_var_tab is a logically
   // (num_side_sets x num_sset_var) integer array of 0s and 1s
   // indicating which sidesets a given sideset variable is defined on.
   std::vector<int> sset_var_tab(num_side_sets * var_names.size());
  
   // We now call read_sideset() once per sideset and write any sideset
   // variable values which are defined there.
   int offset=0;
   for (int ss=0; ss<num_side_sets; ++ss)
     {
       // We don't know num_sides_per_set for each set until we call
       // read_sideset(). The values for each sideset are stored (using
       // the offsets) into the 'elem_list' and 'side_list' arrays of
       // this class.
       offset += (ss > 0 ? num_sides_per_set[ss-1] : 0);
       this->read_sideset(ss, offset);
  
       // Debugging:
       // libMesh::out << "Sideset " << ss_ids[ss]
       //              << " has " << num_sides_per_set[ss]
       //              << " entries."
       //              << std::endl;
  
       // For each variable in var_names, write the values for the
       // current sideset, if any.
       for (auto var : index_range(var_names))
         {
           // If this var has no values on this sideset, go to the next one.
           if (!side_ids[var].count(ss_ids[ss]))
             continue;
  
           // Otherwise, fill in this entry of the sideset truth table.
           sset_var_tab[ss*var_names.size() + var] = 1;
  
           // Data vector that will eventually be passed to exII::ex_put_sset_var().
           std::vector<Real> sset_var_vals(num_sides_per_set[ss]);
  
           // Get reference to the BCTuple -> Real map for this variable.
           const auto & data_map = bc_vals[var];
  
           // Loop over elem_list, side_list entries in current sideset.
           for (int i=0; i<num_sides_per_set[ss]; ++i)
             {
               // Get elem_id and side_id from the respective lists that
               // are filled in by calling read_sideset().
               //
               // Note: these are Exodus-specific ids, so we have to convert them
               // to libmesh ids, as that is what will be in the bc_tuples.
               //
               // TODO: we should probably consult the exodus_elem_num_to_libmesh
               // mapping in order to figure out which libmesh element id 'elem_id'
               // actually corresponds to here, instead of just assuming it will be
               // off by one. Unfortunately that data structure does not seem to
               // be used at the moment. If we assume that write_sideset_data() is
               // always called following write(), then this should be a fairly safe
               // assumption...
               dof_id_type elem_id = elem_list[i + offset] - 1;
               unsigned int side_id = side_list[i + offset] - 1;
  
               // Sanity check: make sure that the "off by one"
               // assumption we used above to set 'elem_id' is valid.
               if (libmesh_map_find(libmesh_elem_num_to_exodus, cast_int<int>(elem_id))
                   != elem_list[i + offset])
                 libmesh_error_msg("Error mapping Exodus elem id to libmesh elem id.");
  
               // Map from Exodus side ids to libmesh side ids.
               const auto & conv = get_conversion(mesh.elem_ptr(elem_id)->type());
  
               // Map from Exodus side ids to libmesh side ids.
               unsigned int converted_side_id = conv.get_side_map(side_id);
  
               // Debugging:
               // libMesh::out << "side_id=" << side_id
               //              << ", converted_side_id=" << converted_side_id
               //              << std::endl;
  
               // Construct a key so we can quickly see whether there is any
               // data for this variable in the map.
               BoundaryInfo::BCTuple key = std::make_tuple
                 (elem_id,
                  converted_side_id,
                  ss_ids[ss]);
  
               // Find the data for this (elem,side,id) tuple. Throw an
               // error if not found. Then store value in vector which
               // will be passed to Exodus.
               sset_var_vals[i] = libmesh_map_find(data_map, key);
             } // end for (i)
  
           // As far as I can tell, there is no "concat" version of writing
           // sideset data, you have to call ex_put_sset_var() once per (variable,
           // sideset) pair.
           if (sset_var_vals.size() > 0)
             {
               ex_err = exII::ex_put_sset_var
                 (ex_id,
                  timestep,
                  var + 1, // 1-based variable index of current variable
                  ss_ids[ss],
                  num_sides_per_set[ss],
                  MappedOutputVector(sset_var_vals, _single_precision).data());
               EX_CHECK_ERR(ex_err, "Error writing sideset vars.");
             }
         } // end for (var)
     } // end for (ss)
  
   // Finally, write the sideset truth table.
   ex_err =
     exII::ex_put_sset_var_tab(ex_id,
                               num_side_sets,
                               cast_int<int>(var_names.size()),
                               sset_var_tab.data());
   EX_CHECK_ERR(ex_err, "Error writing sideset var truth table.");
 }

◆ write_sidesets()

void libMesh::Nemesis_IO_Helper::write_sidesets ( const MeshBase & mesh )

overridevirtual

Writes the sidesets for this processor.

Reimplemented from libMesh::ExodusII_IO_Helper.

Definition at line 2084 of file nemesis_io_helper.C.

 {
   // Write the sidesets.  In Nemesis, the idea is to "create space" for the global
   // set of boundary sidesets, but to only write sideset IDs which are local to the current
   // processor.  This is what is done in Nemesis files created by the "loadbal" script.
   // See also: ExodusII_IO_Helper::write_sidesets()...
  
  
   // Store a map of vectors for boundary side IDs on this processor.
   // Use a vector of int here so it can be passed directly to Exodus.
   std::map<boundary_id_type, std::vector<int>> local_elem_boundary_id_lists;
   std::map<boundary_id_type, std::vector<int>> local_elem_boundary_id_side_lists;
  
   // FIXME: We already built this list once, we should reuse that information!
   std::vector<std::tuple<dof_id_type, unsigned short int, boundary_id_type>> bndry_elem_side_id_list =
     mesh.get_boundary_info().build_side_list();
  
   // Integer looping, skipping non-local elements
   for (const auto & t : bndry_elem_side_id_list)
     {
       // Get pointer to current Elem
       const Elem * elem = mesh.elem_ptr(std::get<0>(t));
  
       std::vector<const Elem *> family;
 #ifdef LIBMESH_ENABLE_AMR
       // We need to build up active elements if AMR is enabled and add
       // them to the exodus sidesets instead of the potentially inactive "parent" elements
       // Technically we don't need to "reset" the tree since the vector was just created.
       elem->active_family_tree_by_side(family, std::get<1>(t), /*reset tree=*/false);
 #else
       // If AMR is not even enabled, just push back the element itself
       family.push_back( elem );
 #endif
  
       // Loop over all the elements in the family tree, store their converted IDs
       // and side IDs to the map's vectors.  TODO: Somehow reserve enough space for these
       // push_back's...
       for (const auto & tree_elem : family)
         {
           const dof_id_type f_id = tree_elem->id();
           const Elem & f = mesh.elem_ref(f_id);
  
           // If element is local, process it
           if (f.processor_id() == this->processor_id())
             {
               const auto & conv = get_conversion(f.type());
  
               // Use the libmesh to exodus data structure map to get the proper sideset IDs
               // The data structure contains the "collapsed" contiguous ids.
               //
               // We know the parent element is local, but let's be absolutely sure that all the children have been
               // actually mapped to Exodus IDs before we blindly try to add them...
               local_elem_boundary_id_lists[ std::get<2>(t) ].push_back( libmesh_map_find(libmesh_elem_num_to_exodus, f_id) );
               local_elem_boundary_id_side_lists[ std::get<2>(t) ].push_back(conv.get_inverse_side_map( std::get<1>(t) ));
             }
         }
     }
  
  
   // Loop over *global* sideset IDs, call the Exodus API.  Note that some sidesets may be empty
   // for a given processor.
   if (global_sideset_ids.size() > 0)
     {
       NamesData names_table(global_sideset_ids.size(), MAX_STR_LENGTH);
  
       for (const auto & exodus_id : this->global_sideset_ids)
         {
           const std::string & current_ss_name =
             mesh.get_boundary_info().get_sideset_name
               (cast_int<boundary_id_type>(exodus_id));
  
           // Store this name in a data structure that will be used to
           // write sideset names to file.
           names_table.push_back_entry(current_ss_name);
  
           if (verbose)
             {
               libMesh::out << "[" << this->processor_id()
                            << "] Writing out Exodus sideset info for ID: " << exodus_id
                            << ", Name: " << current_ss_name
                            << std::endl;
             }
  
           // Try to find this boundary ID in the local list we created
           auto it =
             local_elem_boundary_id_lists.find (cast_int<boundary_id_type>(exodus_id));
  
           // No sides found for this boundary ID on this processor
           if (it == local_elem_boundary_id_lists.end())
             {
               if (verbose)
                 libMesh::out << "[" << this->processor_id()
                              << "] No sideset data for ID: " << exodus_id
                              << " on this processor." << std::endl;
  
               // Call the Exodus interface to write the parameters of this side set
               this->ex_err = exII::ex_put_side_set_param(this->ex_id,
                                                          exodus_id,
                                                          0, /* No sides for this ID */
                                                          0  /* No distribution factors */);
               EX_CHECK_ERR(this->ex_err, "Error writing sideset parameters in Nemesis");
  
             }
           else // Boundary ID *was* found in list
             {
               // Get reference to the vector of elem IDs
               std::vector<int> & current_sideset_elem_ids = it->second;
  
               // Get reference to the vector of side IDs
               std::vector<int> & current_sideset_side_ids =
                 libmesh_map_find(local_elem_boundary_id_side_lists,
                                  cast_int<boundary_id_type>(exodus_id));
  
               // Call the Exodus interface to write the parameters of this side set
               this->ex_err = exII::ex_put_side_set_param(this->ex_id,
                                                          exodus_id,
                                                          current_sideset_elem_ids.size(),
                                                          0  /* No distribution factors */);
  
               EX_CHECK_ERR(this->ex_err, "Error writing sideset parameters in Nemesis");
  
               // Call Exodus interface to write the actual side IDs for this boundary ID
               this->ex_err = exII::ex_put_side_set(this->ex_id,
                                                    exodus_id,
                                                    current_sideset_elem_ids.data(),
                                                    current_sideset_side_ids.data());
  
               EX_CHECK_ERR(this->ex_err, "Error writing sidesets in Nemesis");
             }
         } // end for loop over global sideset IDs
  
       // Write sideset names to file.  Some of these may be blank strings
       // if the current processor didn't have all the sideset names for
       // any reason...
       ex_err = exII::ex_put_names(this->ex_id,
                                   exII::EX_SIDE_SET,
                                   names_table.get_char_star_star());
       EX_CHECK_ERR(ex_err, "Error writing sideset names");
  
     } // end if (global_sideset_ids.size() > 0)
 }

References libMesh::Elem::active_family_tree_by_side(), libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::get_conversion(), global_sideset_ids, libMesh::ExodusII_IO_Helper::libmesh_elem_num_to_exodus, mesh, libMesh::out, libMesh::ParallelObject::processor_id(), libMesh::DofObject::processor_id(), libMesh::Elem::type(), and libMesh::ExodusII_IO_Helper::verbose.

◆ write_timestep()

void libMesh::ExodusII_IO_Helper::write_timestep	(	int	timestep,
		Real	time
	)

inherited

Writes the time for the timestep.

Definition at line 2562 of file exodusII_io_helper.C.

 {
   if ((_run_only_on_proc0) && (this->processor_id() != 0))
     return;
  
   if (_single_precision)
     {
       float cast_time = float(time);
       ex_err = exII::ex_put_time(ex_id, timestep, &cast_time);
     }
   else
     {
       double cast_time = double(time);
       ex_err = exII::ex_put_time(ex_id, timestep, &cast_time);
     }
   EX_CHECK_ERR(ex_err, "Error writing timestep.");
  
   ex_err = exII::ex_update(ex_id);
   EX_CHECK_ERR(ex_err, "Error flushing buffers to file.");
 }

References libMesh::ExodusII_IO_Helper::_run_only_on_proc0, libMesh::ExodusII_IO_Helper::_single_precision, libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, and libMesh::ParallelObject::processor_id().

◆ write_var_names()

void libMesh::ExodusII_IO_Helper::write_var_names	(	ExodusVarType	type,
		const std::vector< std::string > &	names
	)

protectedinherited

Wraps calls to exII::ex_put_var_names() and exII::ex_put_var_param().

The enumeration controls whether nodal, elemental, or global variable names are read and which class members are filled in.

Definition at line 1408 of file exodusII_io_helper.C.

 {
   switch (type)
     {
     case NODAL:
       this->write_var_names_impl("n", num_nodal_vars, names);
       break;
     case ELEMENTAL:
       this->write_var_names_impl("e", num_elem_vars, names);
       break;
     case GLOBAL:
       this->write_var_names_impl("g", num_global_vars, names);
       break;
     case SIDESET:
       {
         // Note: calling this function *sets* num_sideset_vars to the
         // number of entries in the 'names' vector, num_sideset_vars
         // does not already need to be set before calling this.
         this->write_var_names_impl("s", num_sideset_vars, names);
         break;
       }
     default:
       libmesh_error_msg("Unrecognized ExodusVarType " << type);
     }
 }

References libMesh::ExodusII_IO_Helper::ELEMENTAL, libMesh::ExodusII_IO_Helper::GLOBAL, libMesh::ExodusII_IO_Helper::NODAL, libMesh::ExodusII_IO_Helper::num_elem_vars, libMesh::ExodusII_IO_Helper::num_global_vars, libMesh::ExodusII_IO_Helper::num_nodal_vars, libMesh::ExodusII_IO_Helper::num_sideset_vars, libMesh::ExodusII_IO_Helper::SIDESET, and libMesh::ExodusII_IO_Helper::write_var_names_impl().

Referenced by libMesh::ExodusII_IO_Helper::initialize_element_variables(), initialize_element_variables(), libMesh::ExodusII_IO_Helper::initialize_global_variables(), libMesh::ExodusII_IO_Helper::initialize_nodal_variables(), and libMesh::ExodusII_IO_Helper::write_sideset_data().

◆ write_var_names_impl()

void libMesh::ExodusII_IO_Helper::write_var_names_impl	(	const char *	var_type,
		int &	count,
		const std::vector< std::string > &	names
	)

privateinherited

write_var_names() dispatches to this function.

Definition at line 1438 of file exodusII_io_helper.C.

 {
   // Update the count variable so that it's available to other parts of the class.
   count = cast_int<int>(names.size());
  
   // Write that number of variables to the file.
   ex_err = exII::ex_put_var_param(ex_id, var_type, count);
   EX_CHECK_ERR(ex_err, "Error setting number of vars.");
  
   if (count > 0)
     {
       NamesData names_table(count, MAX_STR_LENGTH);
  
       // Store the input names in the format required by Exodus.
       for (int i=0; i != count; ++i)
         names_table.push_back_entry(names[i]);
  
       if (verbose)
         {
           libMesh::out << "Writing variable name(s) to file: " << std::endl;
           for (int i=0; i != count; ++i)
             libMesh::out << names_table.get_char_star(i) << std::endl;
         }
  
       ex_err = exII::ex_put_var_names(ex_id,
                                       var_type,
                                       count,
                                       names_table.get_char_star_star()
                                       );
  
       EX_CHECK_ERR(ex_err, "Error writing variable names.");
     }
 }

References libMesh::ExodusII_IO_Helper::ex_err, libMesh::ExodusII_IO_Helper::ex_id, libMesh::ExodusII_IO_Helper::NamesData::get_char_star(), libMesh::ExodusII_IO_Helper::NamesData::get_char_star_star(), libMesh::out, libMesh::ExodusII_IO_Helper::NamesData::push_back_entry(), and libMesh::ExodusII_IO_Helper::verbose.

Referenced by libMesh::ExodusII_IO_Helper::write_var_names().

Member Data Documentation

◆ _communicator

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

protectedinherited

Definition at line 112 of file parallel_object.h.

Referenced by libMesh::EquationSystems::build_parallel_elemental_solution_vector(), libMesh::EquationSystems::build_parallel_solution_vector(), libMesh::ParallelObject::comm(), libMesh::ParallelObject::n_processors(), libMesh::ParallelObject::operator=(), and libMesh::ParallelObject::processor_id().

◆ _coordinate_offset

Point libMesh::ExodusII_IO_Helper::_coordinate_offset

protectedinherited

Definition at line 737 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::set_coordinate_offset(), and libMesh::ExodusII_IO_Helper::write_nodal_coordinates().

◆ _elem_vars_initialized

bool libMesh::ExodusII_IO_Helper::_elem_vars_initialized

protectedinherited

Definition at line 719 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize_element_variables(), and initialize_element_variables().

◆ _global_vars_initialized

bool libMesh::ExodusII_IO_Helper::_global_vars_initialized

protectedinherited

Definition at line 722 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize_global_variables().

◆ _nodal_vars_initialized

bool libMesh::ExodusII_IO_Helper::_nodal_vars_initialized

protectedinherited

Definition at line 725 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize_nodal_variables().

◆ _run_only_on_proc0

bool libMesh::ExodusII_IO_Helper::_run_only_on_proc0

protectedinherited

Definition at line 716 of file exodusII_io_helper.h.

◆ _single_precision

bool libMesh::ExodusII_IO_Helper::_single_precision

protectedinherited

Definition at line 740 of file exodusII_io_helper.h.

◆ _use_mesh_dimension_instead_of_spatial_dimension

bool libMesh::ExodusII_IO_Helper::_use_mesh_dimension_instead_of_spatial_dimension

protectedinherited

Definition at line 730 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize(), and libMesh::ExodusII_IO_Helper::use_mesh_dimension_instead_of_spatial_dimension().

◆ _write_as_dimension

unsigned libMesh::ExodusII_IO_Helper::_write_as_dimension

protectedinherited

Definition at line 734 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize(), and libMesh::ExodusII_IO_Helper::write_as_dimension().

◆ block_id_to_elem_connectivity

std::map<int, std::vector<int> > libMesh::Nemesis_IO_Helper::block_id_to_elem_connectivity

This is the block connectivity, i.e.

for each subdomain (block) there is an element connectivity list. This map associates the block ID to that vector.

Definition at line 443 of file nemesis_io_helper.h.

Referenced by build_element_and_node_maps(), and write_elements().

◆ block_ids

std::vector<int> libMesh::ExodusII_IO_Helper::block_ids

inherited

Definition at line 541 of file exodusII_io_helper.h.

Referenced by build_element_and_node_maps(), libMesh::ExodusII_IO_Helper::get_block_id(), libMesh::ExodusII_IO_Helper::get_block_name(), libMesh::ExodusII_IO_Helper::initialize_element_variables(), libMesh::ExodusII_IO_Helper::read_block_info(), libMesh::ExodusII_IO_Helper::read_elem_in_block(), libMesh::ExodusII_IO_Helper::read_elemental_var_values(), and libMesh::ExodusII_IO_Helper::write_elements().

◆ border_elem_ids

std::set<unsigned> libMesh::Nemesis_IO_Helper::border_elem_ids

private

A set of border elem IDs for this processor.

Definition at line 618 of file nemesis_io_helper.h.

Referenced by compute_element_maps(), and compute_internal_and_border_elems_and_internal_nodes().

◆ border_node_ids

std::set<unsigned> libMesh::Nemesis_IO_Helper::border_node_ids

private

The set which will eventually contain the IDs of "border nodes".

These are nodes that lie on the boundary between one or more processors.

Definition at line 581 of file nemesis_io_helper.h.

Referenced by compute_border_node_ids(), compute_internal_and_border_elems_and_internal_nodes(), and compute_node_maps().

◆ connect

std::vector<int> libMesh::ExodusII_IO_Helper::connect

inherited

Definition at line 547 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_edge_blocks(), libMesh::ExodusII_IO_Helper::read_elem_in_block(), and libMesh::ExodusII_IO_Helper::write_elements().

◆ conversion_map

std::map<ElemType, ExodusII_IO_Helper::Conversion> libMesh::ExodusII_IO_Helper::conversion_map

privateinherited

Associates libMesh ElemTypes with node/face/edge/etc.

mappings of the corresponding Exodus element types.

Definition at line 823 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::get_conversion(), and libMesh::ExodusII_IO_Helper::init_conversion_map().

◆ current_filename

std::string libMesh::ExodusII_IO_Helper::current_filename

inherited

Definition at line 681 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::create(), and libMesh::ExodusII_IO_Helper::open().

◆ edge_block_ids

std::vector<int> libMesh::ExodusII_IO_Helper::edge_block_ids

inherited

Definition at line 544 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_block_info(), and libMesh::ExodusII_IO_Helper::read_edge_blocks().

◆ elem_cmap_elem_cnts

std::vector<int> libMesh::Nemesis_IO_Helper::elem_cmap_elem_cnts

Definition at line 545 of file nemesis_io_helper.h.

Referenced by compute_communication_map_parameters(), get_cmap_params(), get_elem_cmap(), and initialize().

◆ elem_cmap_elem_ids

std::vector<std::vector<int> > libMesh::Nemesis_IO_Helper::elem_cmap_elem_ids

3 vectors of vectors for storing element communication IDs for this processor.

There will be num_elem_cmaps rows, row i will have elem_cmap_elem_cnts[i] entries. To be used with Nemesis::ne_get_elem_cmap().

Definition at line 565 of file nemesis_io_helper.h.

Referenced by compute_elem_communication_maps(), get_elem_cmap(), and initialize().

◆ elem_cmap_ids

std::vector<int> libMesh::Nemesis_IO_Helper::elem_cmap_ids

Definition at line 544 of file nemesis_io_helper.h.

Referenced by compute_communication_map_parameters(), compute_elem_communication_maps(), get_cmap_params(), get_elem_cmap(), initialize(), and put_elem_cmap().

◆ elem_cmap_proc_ids

std::vector<std::vector<int> > libMesh::Nemesis_IO_Helper::elem_cmap_proc_ids

Definition at line 567 of file nemesis_io_helper.h.

Referenced by compute_elem_communication_maps(), get_elem_cmap(), and initialize().

◆ elem_cmap_side_ids

std::vector<std::vector<int> > libMesh::Nemesis_IO_Helper::elem_cmap_side_ids

Definition at line 566 of file nemesis_io_helper.h.

Referenced by compute_elem_communication_maps(), get_elem_cmap(), and initialize().

◆ elem_list

std::vector<int> libMesh::ExodusII_IO_Helper::elem_list

inherited

Definition at line 584 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_sideset(), libMesh::ExodusII_IO_Helper::read_sideset_data(), libMesh::ExodusII_IO_Helper::read_sideset_info(), and libMesh::ExodusII_IO_Helper::write_sideset_data().

◆ elem_mapb

std::vector<int> libMesh::Nemesis_IO_Helper::elem_mapb

Vector which stores border element IDs.

Will have length num_border_elems. To be used with Nemesis::ne_get_elem_map().

Definition at line 510 of file nemesis_io_helper.h.

Referenced by compute_element_maps(), get_elem_map(), and initialize().

◆ elem_mapi

std::vector<int> libMesh::Nemesis_IO_Helper::elem_mapi

Vector which stores internal element IDs.

Will have length num_internal_elems. To be used with Nemesis::ne_get_elem_map().

Definition at line 503 of file nemesis_io_helper.h.

Referenced by compute_element_maps(), get_elem_map(), and initialize().

◆ elem_num_map

std::vector<int> libMesh::ExodusII_IO_Helper::elem_num_map

inherited

Definition at line 599 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_elem_num_map(), libMesh::ExodusII_IO_Helper::read_elemental_var_values(), and libMesh::ExodusII_IO_Helper::write_elements().

◆ elem_type

std::vector<char> libMesh::ExodusII_IO_Helper::elem_type

inherited

Definition at line 614 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::ExodusII_IO_Helper(), libMesh::ExodusII_IO_Helper::get_elem_type(), libMesh::ExodusII_IO_Helper::read_edge_blocks(), and libMesh::ExodusII_IO_Helper::read_elem_in_block().

◆ elem_var_names

std::vector<std::string> libMesh::ExodusII_IO_Helper::elem_var_names

inherited

Definition at line 645 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize_element_variables(), initialize_element_variables(), libMesh::ExodusII_IO_Helper::read_elemental_var_values(), and libMesh::ExodusII_IO_Helper::read_var_names().

◆ elem_var_values

std::vector<Real> libMesh::ExodusII_IO_Helper::elem_var_values

inherited

Definition at line 648 of file exodusII_io_helper.h.

◆ element_equivalence_map

std::map<std::string, ElemType> libMesh::ExodusII_IO_Helper::element_equivalence_map

privateinherited

Defines equivalence classes of Exodus element types that map to libmesh ElemTypes.

Definition at line 816 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::get_conversion(), and libMesh::ExodusII_IO_Helper::init_element_equivalence_map().

◆ ex_err

int libMesh::ExodusII_IO_Helper::ex_err

inherited

Definition at line 495 of file exodusII_io_helper.h.

◆ ex_id

int libMesh::ExodusII_IO_Helper::ex_id

inherited

Definition at line 492 of file exodusII_io_helper.h.

◆ exodus_elem_num_to_libmesh

std::vector<int> libMesh::ExodusII_IO_Helper::exodus_elem_num_to_libmesh

inherited

Definition at line 619 of file exodusII_io_helper.h.

Referenced by build_element_and_node_maps(), and write_elements().

◆ exodus_node_num_to_libmesh

std::vector<int> libMesh::ExodusII_IO_Helper::exodus_node_num_to_libmesh

inherited

Definition at line 624 of file exodusII_io_helper.h.

Referenced by build_element_and_node_maps(), write_nodal_coordinates(), and write_nodal_solution().

◆ ftype

char libMesh::Nemesis_IO_Helper::ftype

The type of file to be written.

Either 's', for a scalar load-balance file, or 'p' for a parallel file. To be used with Nemesis::ne_get_init_info().

Definition at line 395 of file nemesis_io_helper.h.

Referenced by get_init_info().

◆ global_elem_blk_cnts

std::vector<int> libMesh::Nemesis_IO_Helper::global_elem_blk_cnts

Definition at line 427 of file nemesis_io_helper.h.

Referenced by compute_num_global_elem_blocks(), get_eb_info_global(), and initialize().

◆ global_elem_blk_ids

std::vector<int> libMesh::Nemesis_IO_Helper::global_elem_blk_ids

Read the global element block IDs and counts.

These vectors will eventually have num_elem_blks_global entries. To be used with Nemesis::ne_get_eb_info_global().

Definition at line 426 of file nemesis_io_helper.h.

Referenced by compute_num_global_elem_blocks(), get_eb_info_global(), initialize(), initialize_element_variables(), write_element_values(), and write_elements().

◆ global_nodeset_ids

std::vector<int> libMesh::Nemesis_IO_Helper::global_nodeset_ids

Containers for reading global nodeset information.

One vector entry per nodeset. Each vector will eventually have num_node_sets_global entries, and will be used in calls to Nemesis::ne_get_ns_param_global().

It's an error to call ne_get_ns_param_global when num_node_sets_global==0

Definition at line 416 of file nemesis_io_helper.h.

Referenced by compute_num_global_nodesets(), get_ns_param_global(), initialize(), put_ns_param_global(), write_exodus_initialization_info(), and write_nodesets().

◆ global_sideset_ids

std::vector<int> libMesh::Nemesis_IO_Helper::global_sideset_ids

Containers for reading global sideset (boundary conditions) information.

Each vector will eventually have num_side_sets_global entries, and be used in calls to Nemesis::ne_get_ss_param_global().

It's an error to call ne_get_ss_param_global when num_side_sets_global==0

Definition at line 404 of file nemesis_io_helper.h.

Referenced by compute_num_global_sidesets(), get_ss_param_global(), initialize(), put_ss_param_global(), write_exodus_initialization_info(), and write_sidesets().

◆ global_var_names

std::vector<std::string> libMesh::ExodusII_IO_Helper::global_var_names

inherited

Definition at line 651 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize_global_variables(), and libMesh::ExodusII_IO_Helper::read_var_names().

◆ id_list

std::vector<int> libMesh::ExodusII_IO_Helper::id_list

inherited

Definition at line 593 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_sideset(), and libMesh::ExodusII_IO_Helper::read_sideset_info().

◆ id_to_block_names

std::map<int, std::string> libMesh::ExodusII_IO_Helper::id_to_block_names

inherited

Definition at line 657 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::get_block_name(), and libMesh::ExodusII_IO_Helper::read_block_info().

◆ id_to_edge_block_names

std::map<int, std::string> libMesh::ExodusII_IO_Helper::id_to_edge_block_names

inherited

Definition at line 658 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_block_info(), and libMesh::ExodusII_IO_Helper::read_edge_blocks().

◆ id_to_ns_names

std::map<int, std::string> libMesh::ExodusII_IO_Helper::id_to_ns_names

inherited

Definition at line 660 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::get_node_set_name(), libMesh::ExodusII_IO_Helper::read_all_nodesets(), and libMesh::ExodusII_IO_Helper::read_nodeset_info().

◆ id_to_ss_names

std::map<int, std::string> libMesh::ExodusII_IO_Helper::id_to_ss_names

inherited

Definition at line 659 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::get_side_set_name(), and libMesh::ExodusII_IO_Helper::read_sideset_info().

◆ internal_elem_ids

std::set<unsigned> libMesh::Nemesis_IO_Helper::internal_elem_ids

private

A set of internal elem IDs for this processor.

Definition at line 613 of file nemesis_io_helper.h.

Referenced by compute_element_maps(), and compute_internal_and_border_elems_and_internal_nodes().

◆ internal_node_ids

std::set<unsigned> libMesh::Nemesis_IO_Helper::internal_node_ids

private

A set of internal node IDs for this processor.

Definition at line 608 of file nemesis_io_helper.h.

Referenced by compute_internal_and_border_elems_and_internal_nodes(), and compute_node_maps().

◆ libmesh_elem_num_to_exodus

std::map<int, int> libMesh::ExodusII_IO_Helper::libmesh_elem_num_to_exodus

inherited

Definition at line 618 of file exodusII_io_helper.h.

Referenced by build_element_and_node_maps(), compute_elem_communication_maps(), compute_element_maps(), libMesh::ExodusII_IO_Helper::write_elements(), libMesh::ExodusII_IO_Helper::write_sideset_data(), write_sidesets(), and libMesh::ExodusII_IO_Helper::write_sidesets().

◆ libmesh_node_num_to_exodus

std::map<int, int> libMesh::ExodusII_IO_Helper::libmesh_node_num_to_exodus

inherited

Definition at line 623 of file exodusII_io_helper.h.

Referenced by build_element_and_node_maps(), compute_node_communication_maps(), compute_node_maps(), libMesh::ExodusII_IO_Helper::write_elements(), libMesh::ExodusII_IO_Helper::write_nodal_coordinates(), and write_nodesets().

◆ local_subdomain_counts

std::map<subdomain_id_type, unsigned> libMesh::Nemesis_IO_Helper::local_subdomain_counts

private

This map keeps track of the number of elements in each subdomain (block) for this processor.

Definition at line 575 of file nemesis_io_helper.h.

Referenced by build_element_and_node_maps(), and compute_num_global_elem_blocks().

◆ nemesis_err_flag

int libMesh::Nemesis_IO_Helper::nemesis_err_flag

Member data.

All (?) Nemesis functions return an int. If it's negative that signals an error! Internally, we use the ExodusII_IO_Helper::check_err() function to check for errors.

Definition at line 365 of file nemesis_io_helper.h.

Referenced by get_cmap_params(), get_eb_info_global(), get_elem_cmap(), get_elem_map(), get_init_global(), get_init_info(), get_loadbal_param(), get_node_cmap(), get_node_map(), get_ns_param_global(), get_ss_param_global(), put_cmap_params(), put_eb_info_global(), put_elem_cmap(), put_elem_map(), put_init_global(), put_init_info(), put_loadbal_param(), put_n_coord(), put_node_cmap(), put_node_map(), put_ns_param_global(), and put_ss_param_global().

◆ nodal_var_names

std::vector<std::string> libMesh::ExodusII_IO_Helper::nodal_var_names

inherited

Definition at line 636 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize_nodal_variables(), libMesh::ExodusII_IO_Helper::read_nodal_var_values(), and libMesh::ExodusII_IO_Helper::read_var_names().

◆ nodal_var_values

std::vector<Real> libMesh::ExodusII_IO_Helper::nodal_var_values

inherited

Definition at line 639 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_nodal_var_values().

◆ node_cmap_ids

std::vector<int> libMesh::Nemesis_IO_Helper::node_cmap_ids

Vectors for storing the communication map parameters.

Each will eventually have length num_node_cmaps OR num_elem_cmaps as appropriate. For use with Nemesis::ne_get_cmap_params().

Definition at line 542 of file nemesis_io_helper.h.

Referenced by compute_communication_map_parameters(), compute_node_communication_maps(), get_cmap_params(), get_node_cmap(), initialize(), and put_node_cmap().

◆ node_cmap_node_cnts

std::vector<int> libMesh::Nemesis_IO_Helper::node_cmap_node_cnts

Definition at line 543 of file nemesis_io_helper.h.

Referenced by compute_communication_map_parameters(), get_cmap_params(), get_node_cmap(), and initialize().

◆ node_cmap_node_ids

std::vector<std::vector<int> > libMesh::Nemesis_IO_Helper::node_cmap_node_ids

2 vectors of vectors for storing the node communication IDs for this processor.

There will be num_node_cmaps rows, row i will have node_cmap_node_cnts[i] entries. To be used with Nemesis::ne_get_node_cmap().

Remark: node_cmap_proc_ids is a vector, all entries of which are = node_cmap_ids[i] Not sure what the point of that is...

Definition at line 556 of file nemesis_io_helper.h.

Referenced by compute_node_communication_maps(), get_node_cmap(), and initialize().

◆ node_cmap_proc_ids

std::vector<std::vector<int> > libMesh::Nemesis_IO_Helper::node_cmap_proc_ids

Definition at line 557 of file nemesis_io_helper.h.

Referenced by compute_node_communication_maps(), get_node_cmap(), and initialize().

◆ node_list

std::vector<int> libMesh::ExodusII_IO_Helper::node_list

inherited

Definition at line 590 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_nodeset().

◆ node_mapb

std::vector<int> libMesh::Nemesis_IO_Helper::node_mapb

Vector which stores border node IDs.

Will have length num_border_nodes. To be used with Nemesis::ne_get_node_map().

Definition at line 526 of file nemesis_io_helper.h.

Referenced by compute_node_maps(), get_node_map(), and initialize().

◆ node_mape

std::vector<int> libMesh::Nemesis_IO_Helper::node_mape

Vector which stores external node IDs.

Will have length num_external_nodes. To be used with Nemesis::ne_get_node_map().

Definition at line 533 of file nemesis_io_helper.h.

Referenced by compute_node_maps(), get_node_map(), and initialize().

◆ node_mapi

std::vector<int> libMesh::Nemesis_IO_Helper::node_mapi

Vector which stores internal node IDs.

Will have length num_internal_nodes. To be used with Nemesis::ne_get_node_map().

Definition at line 519 of file nemesis_io_helper.h.

Referenced by compute_node_maps(), get_node_map(), and initialize().

◆ node_num_map

std::vector<int> libMesh::ExodusII_IO_Helper::node_num_map

inherited

Definition at line 596 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_edge_blocks(), libMesh::ExodusII_IO_Helper::read_nodal_var_values(), libMesh::ExodusII_IO_Helper::read_node_num_map(), and libMesh::ExodusII_IO_Helper::write_nodal_coordinates().

◆ node_sets_dist_fact

std::vector<Real> libMesh::ExodusII_IO_Helper::node_sets_dist_fact

inherited

Definition at line 581 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_all_nodesets().

◆ node_sets_dist_index

std::vector<int> libMesh::ExodusII_IO_Helper::node_sets_dist_index

inherited

Definition at line 573 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_all_nodesets().

◆ node_sets_node_index

std::vector<int> libMesh::ExodusII_IO_Helper::node_sets_node_index

inherited

Definition at line 569 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_all_nodesets(), and libMesh::ExodusII_IO_Helper::write_nodesets().

◆ node_sets_node_list

std::vector<int> libMesh::ExodusII_IO_Helper::node_sets_node_list

inherited

Definition at line 577 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_all_nodesets(), and libMesh::ExodusII_IO_Helper::write_nodesets().

◆ nodes_attached_to_local_elems

std::set<int> libMesh::Nemesis_IO_Helper::nodes_attached_to_local_elems

libMesh numbered node ids attached to local elems.

Definition at line 432 of file nemesis_io_helper.h.

Referenced by build_element_and_node_maps(), and compute_internal_and_border_elems_and_internal_nodes().

◆ nodeset_ids

std::vector<int> libMesh::ExodusII_IO_Helper::nodeset_ids

inherited

Definition at line 553 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::get_node_set_id(), libMesh::ExodusII_IO_Helper::get_node_set_name(), libMesh::ExodusII_IO_Helper::read_all_nodesets(), libMesh::ExodusII_IO_Helper::read_nodeset(), libMesh::ExodusII_IO_Helper::read_nodeset_info(), and libMesh::ExodusII_IO_Helper::write_nodesets().

◆ num_attr

int libMesh::ExodusII_IO_Helper::num_attr

inherited

Definition at line 535 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_elem_in_block().

◆ num_border_elems

int libMesh::Nemesis_IO_Helper::num_border_elems

The number of border FEM elements.

Elements local to this processor but whose FEM nodes reside on other processors as well. To be used with the Nemesis::ne_get_loadbal_param() routine.

Definition at line 480 of file nemesis_io_helper.h.

Referenced by compute_internal_and_border_elems_and_internal_nodes(), get_elem_map(), get_loadbal_param(), and initialize().

◆ num_border_nodes

int libMesh::Nemesis_IO_Helper::num_border_nodes

The number of FEM nodes local to a processor but residing in an element which also has FEM nodes on other processors.

To be used with the Nemesis::ne_get_loadbal_param() routine.

Definition at line 460 of file nemesis_io_helper.h.

Referenced by compute_border_node_ids(), compute_internal_and_border_elems_and_internal_nodes(), get_loadbal_param(), get_node_map(), and initialize().

◆ num_df_per_set

std::vector<int> libMesh::ExodusII_IO_Helper::num_df_per_set

inherited

Definition at line 562 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_sideset(), and libMesh::ExodusII_IO_Helper::read_sideset_info().

◆ num_dim

int libMesh::ExodusII_IO_Helper::num_dim

inherited

Definition at line 498 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize(), libMesh::ExodusII_IO_Helper::print_header(), libMesh::ExodusII_IO_Helper::read_header(), and write_exodus_initialization_info().

◆ num_edge

int libMesh::ExodusII_IO_Helper::num_edge

inherited

Definition at line 516 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize(), and libMesh::ExodusII_IO_Helper::read_header().

◆ num_edge_blk

int libMesh::ExodusII_IO_Helper::num_edge_blk

inherited

Definition at line 520 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize(), libMesh::ExodusII_IO_Helper::read_block_info(), libMesh::ExodusII_IO_Helper::read_edge_blocks(), libMesh::ExodusII_IO_Helper::read_header(), and libMesh::ExodusII_IO_Helper::write_elements().

◆ num_elem

int libMesh::ExodusII_IO_Helper::num_elem

inherited

Definition at line 510 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize(), libMesh::ExodusII_IO_Helper::print_header(), libMesh::ExodusII_IO_Helper::read_elem_num_map(), libMesh::ExodusII_IO_Helper::read_header(), and write_exodus_initialization_info().

◆ num_elem_all_sidesets

int libMesh::ExodusII_IO_Helper::num_elem_all_sidesets

inherited

Definition at line 538 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_sideset_info().

◆ num_elem_blk

int libMesh::ExodusII_IO_Helper::num_elem_blk

inherited

Definition at line 513 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize(), libMesh::ExodusII_IO_Helper::initialize_element_variables(), libMesh::ExodusII_IO_Helper::print_header(), libMesh::ExodusII_IO_Helper::read_block_info(), libMesh::ExodusII_IO_Helper::read_elemental_var_values(), libMesh::ExodusII_IO_Helper::read_header(), libMesh::ExodusII_IO_Helper::write_elements(), and write_exodus_initialization_info().

◆ num_elem_blks_global

int libMesh::Nemesis_IO_Helper::num_elem_blks_global

Definition at line 373 of file nemesis_io_helper.h.

Referenced by compute_num_global_elem_blocks(), get_eb_info_global(), get_init_global(), initialize(), write_elements(), and write_exodus_initialization_info().

◆ num_elem_cmaps

int libMesh::Nemesis_IO_Helper::num_elem_cmaps

The number of elemental communication maps for this processor.

(One per neighboring proc?) To be used with the Nemesis::ne_get_loadbal_param() routine.

Definition at line 494 of file nemesis_io_helper.h.

Referenced by compute_communication_map_parameters(), compute_elem_communication_maps(), compute_internal_and_border_elems_and_internal_nodes(), get_cmap_params(), get_elem_cmap(), get_loadbal_param(), and initialize().

◆ num_elem_this_blk

int libMesh::ExodusII_IO_Helper::num_elem_this_blk

inherited

Definition at line 529 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_elem_in_block(), and libMesh::ExodusII_IO_Helper::read_elemental_var_values().

◆ num_elem_vars

int libMesh::ExodusII_IO_Helper::num_elem_vars

inherited

Definition at line 642 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize_element_variables(), initialize_element_variables(), libMesh::ExodusII_IO_Helper::read_header(), libMesh::ExodusII_IO_Helper::read_var_names(), libMesh::ExodusII_IO_Helper::write_element_values(), libMesh::ExodusII_IO_Helper::write_element_values_element_major(), and libMesh::ExodusII_IO_Helper::write_var_names().

◆ num_elems_global

int libMesh::Nemesis_IO_Helper::num_elems_global

Definition at line 372 of file nemesis_io_helper.h.

Referenced by get_init_global().

◆ num_external_nodes

int libMesh::Nemesis_IO_Helper::num_external_nodes

The number of FEM nodes that reside on another processor but whose element partially resides on the current processor.

To be used with the Nemesis::ne_get_loadbal_param() routine.

Definition at line 467 of file nemesis_io_helper.h.

Referenced by get_loadbal_param(), get_node_map(), and initialize().

◆ num_global_node_counts

std::vector<int> libMesh::Nemesis_IO_Helper::num_global_node_counts

Definition at line 417 of file nemesis_io_helper.h.

Referenced by compute_num_global_nodesets(), get_ns_param_global(), and initialize().

◆ num_global_node_df_counts

std::vector<int> libMesh::Nemesis_IO_Helper::num_global_node_df_counts

Definition at line 418 of file nemesis_io_helper.h.

Referenced by get_ns_param_global(), and initialize().

◆ num_global_side_counts

std::vector<int> libMesh::Nemesis_IO_Helper::num_global_side_counts

Definition at line 405 of file nemesis_io_helper.h.

Referenced by compute_num_global_sidesets(), get_ss_param_global(), and initialize().

◆ num_global_side_df_counts

std::vector<int> libMesh::Nemesis_IO_Helper::num_global_side_df_counts

Definition at line 406 of file nemesis_io_helper.h.

Referenced by get_ss_param_global(), and initialize().

◆ num_global_vars

int libMesh::ExodusII_IO_Helper::num_global_vars

inherited

Definition at line 501 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize_global_variables(), libMesh::ExodusII_IO_Helper::read_global_values(), libMesh::ExodusII_IO_Helper::read_header(), libMesh::ExodusII_IO_Helper::read_var_names(), libMesh::ExodusII_IO_Helper::write_global_values(), and libMesh::ExodusII_IO_Helper::write_var_names().

◆ num_internal_elems

int libMesh::Nemesis_IO_Helper::num_internal_elems

The number of internal FEM elements.

Elements local to this processor. To be used with the Nemesis::ne_get_loadbal_param() routine.

Definition at line 473 of file nemesis_io_helper.h.

Referenced by compute_internal_and_border_elems_and_internal_nodes(), get_elem_map(), get_loadbal_param(), and initialize().

◆ num_internal_nodes

int libMesh::Nemesis_IO_Helper::num_internal_nodes

To be used with the Nemesis::ne_get_loadbal_param() routine.

The number of FEM nodes contained in FEM elements wholly owned by the current processor. To be used with the Nemesis::ne_get_loadbal_param() routine.

Definition at line 453 of file nemesis_io_helper.h.

Referenced by compute_internal_and_border_elems_and_internal_nodes(), get_loadbal_param(), get_node_map(), and initialize().

◆ num_nodal_vars

int libMesh::ExodusII_IO_Helper::num_nodal_vars

inherited

Definition at line 633 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize_nodal_variables(), libMesh::ExodusII_IO_Helper::read_header(), libMesh::ExodusII_IO_Helper::read_var_names(), and libMesh::ExodusII_IO_Helper::write_var_names().

◆ num_node_cmaps

int libMesh::Nemesis_IO_Helper::num_node_cmaps

The number of nodal communication maps for this processor.

(One per neighboring proc?) To be used with the Nemesis::ne_get_loadbal_param() routine.

Definition at line 487 of file nemesis_io_helper.h.

Referenced by compute_border_node_ids(), compute_communication_map_parameters(), compute_node_communication_maps(), get_cmap_params(), get_loadbal_param(), get_node_cmap(), and initialize().

◆ num_node_df_per_set

std::vector<int> libMesh::ExodusII_IO_Helper::num_node_df_per_set

inherited

Definition at line 565 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_all_nodesets(), libMesh::ExodusII_IO_Helper::read_nodeset(), libMesh::ExodusII_IO_Helper::read_nodeset_info(), and libMesh::ExodusII_IO_Helper::write_nodesets().

◆ num_node_sets

int libMesh::ExodusII_IO_Helper::num_node_sets

inherited

Definition at line 523 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize(), libMesh::ExodusII_IO_Helper::print_header(), libMesh::ExodusII_IO_Helper::read_all_nodesets(), libMesh::ExodusII_IO_Helper::read_header(), libMesh::ExodusII_IO_Helper::read_nodeset_info(), and write_exodus_initialization_info().

◆ num_node_sets_global

int libMesh::Nemesis_IO_Helper::num_node_sets_global

Definition at line 374 of file nemesis_io_helper.h.

Referenced by compute_num_global_nodesets(), get_init_global(), get_ns_param_global(), and initialize().

◆ num_nodes

int libMesh::ExodusII_IO_Helper::num_nodes

inherited

Definition at line 507 of file exodusII_io_helper.h.

◆ num_nodes_global

int libMesh::Nemesis_IO_Helper::num_nodes_global

Global initial information.

The names are self-explanatory for the most part. Used with Nemesis::ne_get_init_global().

Definition at line 371 of file nemesis_io_helper.h.

Referenced by get_init_global().

◆ num_nodes_per_elem

int libMesh::ExodusII_IO_Helper::num_nodes_per_elem

inherited

Definition at line 532 of file exodusII_io_helper.h.

Referenced by build_element_and_node_maps(), libMesh::ExodusII_IO_Helper::read_elem_in_block(), write_elements(), and libMesh::ExodusII_IO_Helper::write_elements().

◆ num_nodes_per_set

std::vector<int> libMesh::ExodusII_IO_Helper::num_nodes_per_set

inherited

Definition at line 559 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_all_nodesets(), libMesh::ExodusII_IO_Helper::read_nodeset(), libMesh::ExodusII_IO_Helper::read_nodeset_info(), and libMesh::ExodusII_IO_Helper::write_nodesets().

◆ num_proc

int libMesh::Nemesis_IO_Helper::num_proc

The number of processors for which the NEMESIS I file was created.

To be used with Nemesis::ne_get_init_info().

Definition at line 381 of file nemesis_io_helper.h.

Referenced by get_init_info().

◆ num_proc_in_file

int libMesh::Nemesis_IO_Helper::num_proc_in_file

The number of processors for which the NEMESIS I file stores information.

This is generally equal to 1 (1 CPU/file) at least for the splitting Derek gave us. To be used with Nemesis::ne_get_init_info().

Definition at line 388 of file nemesis_io_helper.h.

Referenced by get_init_info().

◆ num_side_sets

int libMesh::ExodusII_IO_Helper::num_side_sets

inherited

Definition at line 526 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::initialize(), libMesh::ExodusII_IO_Helper::print_header(), libMesh::ExodusII_IO_Helper::read_header(), libMesh::ExodusII_IO_Helper::read_sideset_data(), libMesh::ExodusII_IO_Helper::read_sideset_info(), write_exodus_initialization_info(), and libMesh::ExodusII_IO_Helper::write_sideset_data().

◆ num_side_sets_global

int libMesh::Nemesis_IO_Helper::num_side_sets_global

Definition at line 375 of file nemesis_io_helper.h.

Referenced by compute_num_global_sidesets(), get_init_global(), get_ss_param_global(), and initialize().

◆ num_sides_per_set

std::vector<int> libMesh::ExodusII_IO_Helper::num_sides_per_set

inherited

Definition at line 556 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_sideset(), libMesh::ExodusII_IO_Helper::read_sideset_data(), libMesh::ExodusII_IO_Helper::read_sideset_info(), and libMesh::ExodusII_IO_Helper::write_sideset_data().

◆ num_sideset_vars

int libMesh::ExodusII_IO_Helper::num_sideset_vars

inherited

Definition at line 504 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_header(), libMesh::ExodusII_IO_Helper::read_sideset_data(), libMesh::ExodusII_IO_Helper::read_var_names(), and libMesh::ExodusII_IO_Helper::write_var_names().

◆ num_time_steps

int libMesh::ExodusII_IO_Helper::num_time_steps

inherited

Definition at line 627 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_num_time_steps(), and libMesh::ExodusII_IO_Helper::read_time_steps().

◆ opened_for_reading

bool libMesh::ExodusII_IO_Helper::opened_for_reading

inherited

Definition at line 671 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::close(), and libMesh::ExodusII_IO_Helper::open().

◆ opened_for_writing

bool libMesh::ExodusII_IO_Helper::opened_for_writing

inherited

Definition at line 667 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::close(), libMesh::ExodusII_IO_Helper::create(), create(), libMesh::ExodusII_IO_Helper::open(), and ~Nemesis_IO_Helper().

◆ proc_border_elem_sets

std::map<unsigned, std::set<std::pair<unsigned,unsigned> > > libMesh::Nemesis_IO_Helper::proc_border_elem_sets

private

Map between processor ID and (element,side) pairs bordering that processor ID.

Definition at line 598 of file nemesis_io_helper.h.

Referenced by compute_communication_map_parameters(), compute_elem_communication_maps(), and compute_internal_and_border_elems_and_internal_nodes().

◆ proc_nodes_touched_intersections

std::map<unsigned, std::set<unsigned> > libMesh::Nemesis_IO_Helper::proc_nodes_touched_intersections

private

Another map to store sets of intersections with each other processor (other than ourself, of course).

A node which appears in one of these vectors belongs to element owned by at least this processor and one other.

Definition at line 588 of file nemesis_io_helper.h.

Referenced by compute_border_node_ids(), compute_communication_map_parameters(), and compute_node_communication_maps().

◆ side_list

std::vector<int> libMesh::ExodusII_IO_Helper::side_list

inherited

Definition at line 587 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_sideset(), libMesh::ExodusII_IO_Helper::read_sideset_data(), libMesh::ExodusII_IO_Helper::read_sideset_info(), and libMesh::ExodusII_IO_Helper::write_sideset_data().

◆ sideset_var_names

std::vector<std::string> libMesh::ExodusII_IO_Helper::sideset_var_names

inherited

Definition at line 654 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_sideset_data(), and libMesh::ExodusII_IO_Helper::read_var_names().

◆ ss_ids

std::vector<int> libMesh::ExodusII_IO_Helper::ss_ids

inherited

Definition at line 550 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::get_side_set_id(), libMesh::ExodusII_IO_Helper::get_side_set_name(), libMesh::ExodusII_IO_Helper::read_sideset(), libMesh::ExodusII_IO_Helper::read_sideset_data(), libMesh::ExodusII_IO_Helper::read_sideset_info(), and libMesh::ExodusII_IO_Helper::write_sideset_data().

◆ subdomain_map

std::map<subdomain_id_type, std::vector<dof_id_type> > libMesh::Nemesis_IO_Helper::subdomain_map

Map of subdomains to element numbers.

Definition at line 437 of file nemesis_io_helper.h.

Referenced by build_element_and_node_maps(), write_element_values(), and write_elements().

◆ time_steps

std::vector<Real> libMesh::ExodusII_IO_Helper::time_steps

inherited

Definition at line 630 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::read_time_steps().

◆ title

std::vector<char> libMesh::ExodusII_IO_Helper::title

inherited

Definition at line 611 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::ExodusII_IO_Helper(), libMesh::ExodusII_IO_Helper::print_header(), and libMesh::ExodusII_IO_Helper::read_header().

◆ verbose

bool libMesh::ExodusII_IO_Helper::verbose

inherited

Definition at line 663 of file exodusII_io_helper.h.

◆ x

std::vector<Real> libMesh::ExodusII_IO_Helper::x

inherited

Definition at line 602 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::print_nodes(), libMesh::ExodusII_IO_Helper::read_nodes(), write_nodal_coordinates(), and libMesh::ExodusII_IO_Helper::write_nodal_coordinates().

◆ y

std::vector<Real> libMesh::ExodusII_IO_Helper::y

inherited

Definition at line 605 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::print_nodes(), libMesh::ExodusII_IO_Helper::read_nodes(), write_nodal_coordinates(), and libMesh::ExodusII_IO_Helper::write_nodal_coordinates().

◆ z

std::vector<Real> libMesh::ExodusII_IO_Helper::z

inherited

Definition at line 608 of file exodusII_io_helper.h.

Referenced by libMesh::ExodusII_IO_Helper::print_nodes(), libMesh::ExodusII_IO_Helper::read_nodes(), write_nodal_coordinates(), and libMesh::ExodusII_IO_Helper::write_nodal_coordinates().

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

include/mesh/nemesis_io_helper.h
src/mesh/nemesis_io_helper.C

Public Types

Public Member Functions

Public Attributes

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ proc_border_elem_sets_iterator

◆ proc_nodes_touched_iterator

Member Enumeration Documentation

◆ ExodusVarType

Constructor & Destructor Documentation

◆ Nemesis_IO_Helper()

◆ ~Nemesis_IO_Helper()

Member Function Documentation

◆ build_element_and_node_maps()

◆ check_existing_vars()

◆ close()

◆ comm()

◆ compute_border_node_ids()

◆ compute_communication_map_parameters()

◆ compute_elem_communication_maps()

◆ compute_element_maps()

◆ compute_internal_and_border_elems_and_internal_nodes()

◆ compute_node_communication_maps()

◆ compute_node_maps()

◆ compute_num_global_elem_blocks()

◆ compute_num_global_nodesets()

◆ compute_num_global_sidesets()

◆ construct_nemesis_filename()

◆ create()

◆ get_block_id()

◆ get_block_name()

◆ get_cmap_params()

◆ get_complex_names()

◆ get_complex_subdomain_to_var_names()

◆ get_complex_vars_active_subdomains()

◆ get_conversion() [1/2]

◆ get_conversion() [2/2]

◆ get_eb_info_global()

◆ get_elem_cmap()

◆ get_elem_map()

◆ get_elem_type()

◆ get_init_global()

◆ get_init_info()

◆ get_loadbal_param()

◆ get_node_cmap()

◆ get_node_map()

◆ get_node_set_id()

◆ get_node_set_name()

◆ get_ns_param_global()

◆ get_side_set_id()

◆ get_side_set_name()

◆ get_ss_param_global()

◆ init_conversion_map()

◆ init_element_equivalence_map()

◆ initialize()

◆ initialize_element_variables()

◆ initialize_global_variables()

◆ initialize_nodal_variables()

◆ inquire()

◆ message() [1/2]

◆ message() [2/2]

◆ n_processors()

◆ open()

◆ print_header()

◆ print_nodes()

◆ processor_id()

◆ put_cmap_params()

◆ put_eb_info_global()

◆ put_elem_cmap()

◆ put_elem_map()

◆ put_init_global()

◆ put_init_info()

◆ put_loadbal_param()

◆ put_n_coord()

◆ put_node_cmap()