TestErrorHandlingPar.cpp

Go to the documentation of this file.

/** @example TestErrorHandlingPar.cpp \n
 * Description: This example tests MOAB's trace back error handler in parallel.\n
 *
 * <b>To run</b>: mpiexec -np <n> ./TestErrorHandlingPar <test_case_num(1 to 2)> \n
 */

#include "moab/Core.hpp"
#ifdef MOAB_HAVE_MPI
#include "moab_mpi.h"
#endif

#include <iostream>

using namespace moab;
using namespace std;

// In this test case, a global fatal error MB_NOT_IMPLEMENTED is returned by MOAB
// Note, it is printed by root processor 0 only
ErrorCode TestErrorHandlingPar_1()
{
    Core moab;
    Interface& mb = moab;

    string opts = ";;";
#ifdef MOAB_HAVE_MPI
    // Use parallel options
    opts += "PARALLEL=READ_PART;PARTITION_METHOD=SQIJ";
#endif

    // Load a CAM-FV file and read a variable on edges (not supported yet)
    string test_file = string( MESH_DIR ) + string( "/io/fv3x46x72.t.3.nc" );
    opts += ";VARIABLE=US";
    ErrorCode rval = mb.load_file( test_file.c_str(), NULL, opts.c_str() );MB_CHK_ERR( rval );

    return MB_SUCCESS;
}

// In this test case, a per-processor relevant error MB_FAILURE is returned by MOAB
// Note, it is printed by all processors
ErrorCode TestErrorHandlingPar_2()
{
    Core moab;
    Interface& mb = moab;

    string opts = ";;";
#ifdef MOAB_HAVE_MPI
    // Use parallel options
    opts += "PARALLEL=READ_PART;PARTITION_METHOD=UNKNOWN";
#endif

    // Load a CAM-FV file with an unknown partition method specified
    string test_file = string( MESH_DIR ) + string( "/io/fv3x46x72.t.3.nc" );
    opts += ";VARIABLE=T";
    ErrorCode rval = mb.load_file( test_file.c_str(), NULL, opts.c_str() );MB_CHK_ERR( rval );

    return MB_SUCCESS;
}

int main( int argc, char** argv )
{
    if( argc < 2 )
    {
        cout << "Usage: " << argv[0] << " <test_case_num(1 to 2>" << endl;
        return 0;
    }

#ifdef MOAB_HAVE_MPI
    MPI_Init( &argc, &argv );
#endif

    // Initialize error handler, optional for this example (using moab instances)
    MBErrorHandler_Init();

    ErrorCode rval = MB_SUCCESS;

    int test_case_num = atoi( argv[1] );
    switch( test_case_num )
    {
        case 1:
            rval = TestErrorHandlingPar_1();MB_CHK_ERR( rval );
            break;
        case 2:
            rval = TestErrorHandlingPar_2();MB_CHK_ERR( rval );
            break;
        default:
            break;
    }

    // Finalize error handler, optional for this example (using moab instances)
    MBErrorHandler_Finalize();

#ifdef MOAB_HAVE_MPI
    MPI_Finalize();
#endif

    return 0;
}