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
 Range::iterator set_it;
 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;
}