GetEntities.cpp File Reference

#include "moab/Core.hpp"
#include "moab/Range.hpp"
#include "moab/CN.hpp"
#include <iostream>
#include <memory>
#include <string>
#include <string_view>

Include dependency graph for GetEntities.cpp:

Go to the source code of this file.

Namespaces
	anonymous_namespace{GetEntities.cpp}

Functions
void	print_usage (const char *program_name)
int	main (int argc, char **argv)

Variables
const char *const	anonymous_namespace{GetEntities.cpp}::DEFAULT_MESH_FILE = "hex01.vtk"
const int	anonymous_namespace{GetEntities.cpp}::DIMENSION = 3

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

Examples

GetEntities.cpp.

Definition at line 44 of file GetEntities.cpp.

{
    try
    {
        // Parse command line arguments
        if( argc > 1 && ( std::string( argv[1] ) == "-h" || std::string( argv[1] ) == "--help" ) )
        {
            print_usage( argv[0] );
            return 0;
        }
 
        // Get mesh file path
        const std::string mesh_file =
            ( argc > 1 ) ? argv[1] : std::string( MESH_DIR ) + "/" + std::string( DEFAULT_MESH_FILE );
 
        // Initialize MOAB
        auto moab = std::make_unique< Core >();
        if( !moab )
        {
            std::cerr << "Error: Failed to create MOAB instance\n";
            return 1;
        }
 
        // Load the mesh
        MB_CHK_SET_ERR( moab->load_mesh( mesh_file.c_str() ), "Failed to load mesh file: " << mesh_file );
 
        // Get all entities in the mesh
        Range entities;
        MB_CHK_SET_ERR( moab->get_entities_by_handle( 0, entities ), "Failed to retrieve entities from mesh" );
 
        // Process each entity
        for( const auto& entity : entities )
        {
            const auto entity_type = moab->type_from_handle( entity );
            const auto entity_id   = moab->id_from_handle( entity );
 
            if( entity_type == moab::MBVERTEX )
            {
                // For vertices, get adjacent elements
                Range adjacent_elements;
                MB_CHK_SET_ERR( moab->get_adjacencies( &entity, 1, DIMENSION, false, adjacent_elements ),
                                "Failed to get adjacencies for vertex " << entity_id );
 
                std::cout << "Vertex " << entity_id << " is connected to " << adjacent_elements.size() << " elements\n";
                if( !adjacent_elements.empty() )
                {
                    std::cout << "  Adjacent elements: ";
                    for( const auto& elem : adjacent_elements )
                    {
                        std::cout << moab->id_from_handle( elem ) << " ";
                    }
                    std::cout << "\n";
                }
            }
            else if( entity_type < moab::MBENTITYSET )
            {
                // For elements, get connectivity
                const EntityHandle* connectivity = nullptr;
                int num_vertices                 = 0;
 
                MB_CHK_SET_ERR( moab->get_connectivity( entity, connectivity, num_vertices ),
                                "Failed to get connectivity for " << moab::CN::EntityTypeName( entity_type ) << " "
                                                                  << entity_id );
 
                std::cout << moab::CN::EntityTypeName( entity_type ) << " " << entity_id << " has " << num_vertices
                          << " vertices: ";
 
                for( int i = 0; i < num_vertices; ++i )
                {
                    std::cout << moab->id_from_handle( connectivity[i] );
                    if( i < num_vertices - 1 ) std::cout << ", ";
                }
                std::cout << "\n";
            }
        }
 
        return 0;
    }
    catch( const std::exception& e )
    {
        std::cerr << "Error: " << e.what() << "\n";
        return 1;
    }
    catch( ... )
    {
        std::cerr << "Error: Unknown exception occurred\n";
        return 1;
    }
}

References anonymous_namespace{GetEntities.cpp}::DEFAULT_MESH_FILE, anonymous_namespace{GetEntities.cpp}::DIMENSION, moab::Range::empty(), moab::CN::EntityTypeName(), MB_CHK_SET_ERR, MBENTITYSET, MBVERTEX, MESH_DIR, print_usage(), and moab::Range::size().

print_usage()

void print_usage

(

const char *

program_name

)

Examples

GetEntities.cpp.

Definition at line 38 of file GetEntities.cpp.

{
    std::cout << "Usage: " << program_name << " [meshfile]\n"
              << "  meshfile - Path to the mesh file (default: " << MESH_DIR << "/" << DEFAULT_MESH_FILE << ")\n";
}

References anonymous_namespace{GetEntities.cpp}::DEFAULT_MESH_FILE, and MESH_DIR.

Referenced by main().