Mesh Oriented datABase  (version 5.5.1)
An array-based unstructured mesh library
TestExodusII.cpp

This example demonstrates how to retrieve material, dirichlet and neumann sets from an ExodusII file.
Sets in MOAB contain entities and have a tag on them to give meaning to the entities. Tag names: MATERIAL_SET", "DIRICHLET_SET" and "NEUMANN_SET" are reserved and are associated with their corresponding entity sets. Sets are traversed to find out the type number of entities contained in each set.
Steps in this example :

  1. Instantiate MOAB
  2. Get input mesh file name and load it.
  3. Loop over the three sets: material, dirichlet and neumann
    1. Get TagHandle and EntitySet(corresponding to the TagHandle)
    2. Loop thru all the EntitySet's
      1. Get the set id and entities in this set
  4. Destroy the MOAB instance

To compile: make TestExodusII
To run:

  1. TestExodusII <mesh-file>
  2. TestExodusII (This uses the default <mesh-file>: MeshFiles/unittest/mbtest2.g)
/** @example TestExodusII.cpp
* This example demonstrates how to retrieve material, dirichlet and neumann sets
* from an ExodusII file. \n
* Sets in MOAB contain entities and have a tag on them to give meaning to the entities.
* Tag names: MATERIAL_SET", "DIRICHLET_SET" and "NEUMANN_SET" are reserved and
* are associated with their corresponding entity sets.
* Sets are traversed to find out the type number of entities contained in each set. \n
*
* <b>Steps in this example </b>:
* -# Instantiate MOAB
* -# Get input mesh file name and load it.
* -# Loop over the three sets: material, dirichlet and neumann
* -# Get TagHandle and EntitySet(corresponding to the TagHandle)
* -# Loop thru all the EntitySet's
* -# Get the set id and entities in this set
* -# Destroy the MOAB instance
*
*
* <b> To compile: </b>
* make TestExodusII \n
*
* <b> To run: </b>
* -# TestExodusII <mesh-file> \n
* -# TestExodusII (This uses the default <mesh-file>: MeshFiles/unittest/mbtest2.g)
*/
#include <iostream>
// Include header for MOAB instance and range
#include "moab/Core.hpp"
using namespace moab;
using namespace std;
string test_file_name = string( MESH_DIR ) + string( "/mbtest2.g" );
int main( int argc, char** argv )
{
#ifdef MOAB_HAVE_NETCDF
// Get MOAB instance
Interface* mb = new( std::nothrow ) Core;
if( NULL == mb ) return 1;
// Get the material set tag handle
Tag mtag;
ErrorCode rval;
const char* tag_nms[] = { "MATERIAL_SET", "DIRICHLET_SET", "NEUMANN_SET" };
Range sets, set_ents;
// Load a file
if( argc == 1 )
{
cout << "Running default case, loading " << test_file_name << endl;
cout << "Usage: " << argv[0] << " <filename>\n" << endl;
rval = mb->load_file( test_file_name.c_str() );MB_CHK_ERR( rval );
}
else
{
rval = mb->load_file( argv[argc - 1] );MB_CHK_ERR( rval );
cout << "Loaded mesh file: " << argv[argc - 1] << endl;
}
// Loop over set types
for( int i = 0; i < 3; i++ )
{
rval = mb->tag_get_handle( tag_nms[i], 1, MB_TYPE_INTEGER, mtag );MB_CHK_ERR( rval );
// Get all the sets of that type in the mesh
sets.clear();
rval = mb->get_entities_by_type_and_tag( 0, MBENTITYSET, &mtag, NULL, 1, sets );MB_CHK_ERR( rval );
// Iterate over each set, getting entities
for( set_it = sets.begin(); set_it != sets.end(); ++set_it )
{
EntityHandle this_set = *set_it;
// Get the id for this set
int set_id;
rval = mb->tag_get_data( mtag, &this_set, 1, &set_id );MB_CHK_ERR( rval );
// Get the entities in the set, recursively
rval = mb->get_entities_by_handle( this_set, set_ents, true );MB_CHK_ERR( rval );
cout << tag_nms[i] << " " << set_id << " has " << set_ents.size() << " entities:" << endl;
// Print the entities contained in this set
set_ents.print( " " );
set_ents.clear();
}
}
delete mb;
#else
cout << " This test needs moab configured with netcdf \n";
#endif
return 0;
}