#include "moab/Core.hpp"
#include "moab/Interface.hpp"
#include "moab/Range.hpp"
#include "MBTagConventions.hpp"
#include <iostream>

Include dependency graph for SetsNTags.cpp:

Macros
#define	MESH_DIR "."

Functions
int	main (int argc, char **argv)

Variables
string	test_file_name = string( MESH_DIR ) + string( "/hex01.vtk" )

const char *	tag_nms [] = { MATERIAL_SET_TAG_NAME, DIRICHLET_SET_TAG_NAME, NEUMANN_SET_TAG_NAME }

Macro Definition Documentation

◆ MESH_DIR

#define MESH_DIR "."

Examples: SetsNTags.cpp.

Definition at line 20 of file SetsNTags.cpp.

Function Documentation

◆ main()

int main	(	int	argc,
		char **	argv
	)

Examples: SetsNTags.cpp.

Definition at line 28 of file SetsNTags.cpp.

 {
  // Get the material set tag handle
  Tag mtag;
  ErrorCode rval;
  Range sets, set_ents;
  
  // Get MOAB instance
  Interface* mb = new( std::nothrow ) Core;
  if( NULL == mb ) return 1;
  
  // Need option handling here for input filename
  if( argc > 1 )
  {
  // User has input a mesh file
  test_file_name = argv[1];
  }
  
  // Load a file
  rval = mb->load_file( test_file_name.c_str() );MB_CHK_ERR( rval );
  
  // Loop over set types
  for( int i = 0; i < 3; i++ )
  {
  // Get the tag handle for this tag name; tag should already exist (it was created during
  // file read)
  rval = mb->tag_get_handle( tag_nms[i], 1, MB_TYPE_INTEGER, mtag );MB_CHK_ERR( rval );
  
  // Get all the sets having that tag (with any value for that tag)
  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 the entities and printing them
  Range::iterator set_it;
  for( set_it = sets.begin(); set_it != sets.end(); ++set_it )
  {
  // Get the id for this set
  int set_id;
  rval = mb->tag_get_data( mtag, &( *set_it ), 1, &set_id );MB_CHK_ERR( rval );
  
  // Get the entities in the set, recursively
  rval = mb->get_entities_by_handle( *set_it, set_ents, true );MB_CHK_ERR( rval );
  
  cout << tag_nms[i] << " " << set_id << " has " << set_ents.size() << " entities:" << endl;
  set_ents.print( " " );
  set_ents.clear();
  }
  }
  
  delete mb;
  
  return 0;
 }