UniformRefinement.cpp File Reference

#include "moab/Core.hpp"
#include "moab/NestedRefine.hpp"
#include <iostream>
#include <vector>
#include <string>
#include <cstdlib>
#include <stdexcept>

Include dependency graph for UniformRefinement.cpp:

Go to the source code of this file.

Namespaces
	anonymous_namespace{UniformRefinement.cpp}

Functions
void	anonymous_namespace{UniformRefinement.cpp}::print_usage (const char *program_name)
std::vector< int >	anonymous_namespace{UniformRefinement.cpp}::parse_refinement_levels (int argc, char *argv[])
int	main (int argc, char *argv[])

Variables
constexpr const char *	anonymous_namespace{UniformRefinement.cpp}::DEFAULT_OUTPUT_FILE = "mesh_hierarchy.h5m"
constexpr int	anonymous_namespace{UniformRefinement.cpp}::DEFAULT_REFINEMENT_LEVELS = 2

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

Examples

UniformRefinement.cpp.

Definition at line 80 of file UniformRefinement.cpp.

{
    try
    {
        // Parse command line arguments
        if( argc < 2 || std::string( argv[1] ) == "--help" || std::string( argv[1] ) == "-h" )
        {
            print_usage( argv[0] );
            return argc < 2 ? 1 : 0;
        }
 
        const std::string input_file = argv[1];
        auto level_degrees           = parse_refinement_levels( argc, argv );
        const int num_levels         = static_cast< int >( level_degrees.size() );
 
        std::cout << "Input file: " << input_file << "\n"
                  << "Refinement levels: " << num_levels << "\n"
                  << "Level degrees: ";
 
        for( const auto& degree : level_degrees )
        {
            std::cout << degree << " ";
        }
        std::cout << "\n";
 
        // Initialize MOAB and load the mesh
        auto moab = std::make_unique< Core >();
        if( !moab )
        {
            throw std::runtime_error( "Failed to create MOAB instance" );
        }
 
        // Load the input mesh file
        MB_CHK_SET_ERR_CONT( moab->load_file( input_file.c_str() ), "Failed to load mesh file: " << input_file );
 
        // Set up the refinement
        NestedRefine refiner( moab.get() );
 
        std::vector< EntityHandle > mesh_sets;
        std::cout << "Starting mesh hierarchy generation for " << num_levels << " levels...\n";
        // Perform the refinement
        MB_CHK_SET_ERR( refiner.generate_mesh_hierarchy( num_levels, level_degrees.data(), mesh_sets ),
                        "Failed to generate mesh hierarchy" );
        std::cout << "Successfully generated mesh hierarchy\n";
 
        // Write the refined mesh to file
        MB_CHK_SET_ERR( moab->write_file( DEFAULT_OUTPUT_FILE ),
                        "Failed to write output file: " << DEFAULT_OUTPUT_FILE );
 
        std::cout << "Output written to: " << DEFAULT_OUTPUT_FILE << "\n";
        return 0;
    }
    catch( const std::exception& e )
    {
        std::cerr << "\nError: " << e.what() << "\n";
        return 1;
    }
    catch( ... )
    {
        std::cerr << "\nError: Unknown exception occurred\n";
        return 1;
    }
}

References anonymous_namespace{UniformRefinement.cpp}::DEFAULT_OUTPUT_FILE, moab::NestedRefine::generate_mesh_hierarchy(), input_file, MB_CHK_SET_ERR, MB_CHK_SET_ERR_CONT, anonymous_namespace{UniformRefinement.cpp}::parse_refinement_levels(), and anonymous_namespace{UniformRefinement.cpp}::print_usage().