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

Description: Get entities and report non-vertex entity connectivity and vertex adjacencies.
then delete edges, and write result To run: ./DeleteEdges [meshfile] [outfile]

/** @example GetEntities.cpp
* Description: Get entities and report non-vertex entity connectivity and vertex adjacencies.\n
* This example shows how to get connectivity and adjacencies.\n
*
* To run: ./GetEntities [meshfile]\n
* (default values can run if users don't specify a mesh file)
*/
#include "moab/Core.hpp"
#include "moab/Range.hpp"
#include "moab/CN.hpp"
#include <iostream>
using namespace moab;
using namespace std;
#ifndef MESH_DIR
#define MESH_DIR "."
#endif
string test_file_name = string( MESH_DIR ) + string( "/hex01.vtk" );
int main( int argc, char** argv )
{
if( argc > 1 )
{
// User has input a mesh file
test_file_name = argv[1];
}
// Instantiate & load a mesh from a file
Core* mb = new( std::nothrow ) Core;
if( NULL == mb ) return 1;
ErrorCode rval = mb->load_mesh( test_file_name.c_str() );MB_CHK_ERR( rval );
Range ents;
// Get all entities in the database
rval = mb->get_entities_by_handle( 0, ents );MB_CHK_ERR( rval );
for( Range::iterator it = ents.begin(); it != ents.end(); ++it )
{
if( MBVERTEX == mb->type_from_handle( *it ) )
{
Range adjs;
rval = mb->get_adjacencies( &( *it ), 1, 3, false, adjs );MB_CHK_ERR( rval );
cout << "Vertex " << mb->id_from_handle( *it ) << " adjacencies:" << endl;
adjs.print();
}
else if( mb->type_from_handle( *it ) < MBENTITYSET )
{
const EntityHandle* connect;
int num_connect;
rval = mb->get_connectivity( *it, connect, num_connect );MB_CHK_ERR( rval );
cout << CN::EntityTypeName( mb->type_from_handle( *it ) ) << " " << mb->id_from_handle( *it )
<< " vertex connectivity is: ";
for( int i = 0; i < num_connect; i++ )
cout << mb->id_from_handle( connect[i] ) << " ";
cout << endl;
}
}
delete mb;
return 0;
}