UniformRefinement.cpp

Go to the documentation of this file.

/*This unit test is for the uniform refinement capability based on AHF datastructures*/
#include <iostream>
#include "moab/Core.hpp"
#include "moab/NestedRefine.hpp"

using namespace moab;

int main( int argc, char* argv[] )
{
    Core mb;
    Interface* mbImpl = &mb;
    ErrorCode error;

    if( argc == 1 )
    {
        std::cerr << "Usage: " << argv[0] << " [filename]" << std::endl;
        return 1;
    }

    const char* filename = argv[1];
    error                = mbImpl->load_file( filename );MB_CHK_ERR( error );

    NestedRefine uref( &mb );

    // Usage: The level degrees array controls number of refinemetns and
    // the degree of refinement at each level.
    // Example: int level_degrees[4] = {2,3,2,3};
    std::vector< int > level_degrees;
    if( argc > 2 )
    {
        level_degrees.resize( argc - 2 );
        for( int i = 2; i < argc; ++i )
            level_degrees[i - 2] = atoi( argv[i] );
    }
    else
    {
        level_degrees.resize( 2 );
    }
    int num_levels = static_cast< int >( level_degrees.size() );

    std::cout << "Starting hierarchy generation" << std::endl;
    std::vector< EntityHandle > set;
    error = uref.generate_mesh_hierarchy( num_levels, level_degrees.data(), set );MB_CHK_ERR( error );
    std::cout << "Finished hierarchy generation" << std::endl;

    std::stringstream file;
    file << "mesh_hierarchy.h5m";
    error = mbImpl->write_file( file.str().c_str() );MB_CHK_ERR( error );
    return 0;
}