MOAB: Mesh Oriented datABase  (version 5.2.1)
imoab_coupler_utils.hpp File Reference
+ This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Defines

#define CHECKIERR(rc, message)
#define PUSH_TIMER(operation)
#define POP_TIMER(localcomm, localrank)

Functions

int create_group_and_comm (int start, int end, MPI_Group worldGroup, MPI_Group *group, MPI_Comm *comm)
int create_joint_comm_group (MPI_Group agroup, MPI_Group bgroup, MPI_Group *abgroup, MPI_Comm *abcomm)
int setup_component_coupler_meshes (iMOAB_AppID cmpId, int cmpTag, iMOAB_AppID cplCmpId, int cmpcouTag, MPI_Comm *cmpcomm, MPI_Group *cmpPEGroup, MPI_Comm *coucomm, MPI_Group *cplPEGroup, MPI_Comm *cmpcoucomm, std::string &filename, std::string &readopts, int nghlay, int repartitioner_scheme)

Define Documentation

#define CHECKIERR (   rc,
  message 
)
Value:
if( 0 != rc )                                      \
    {                                                  \
        printf( "%s. ErrorCode = %d\n", message, rc ); \
        return 1;                                      \
    }

Definition at line 15 of file imoab_coupler_utils.hpp.

Referenced by create_group_and_comm(), create_joint_comm_group(), and setup_component_coupler_meshes().

#define POP_TIMER (   localcomm,
  localrank 
)
Value:
{                                                                                             \
        double locElapsed = timer.time_since_birth() - timer_ops, minElapsed = 0, maxElapsed = 0; \
        MPI_Reduce( &locElapsed, &maxElapsed, 1, MPI_DOUBLE, MPI_MAX, 0, localcomm );             \
        MPI_Reduce( &locElapsed, &minElapsed, 1, MPI_DOUBLE, MPI_MIN, 0, localcomm );             \
        if( !localrank )                                                                          \
            std::cout << "[LOG] Time taken to " << opName.c_str() << ": max = " << maxElapsed     \
                      << ", avg = " << ( maxElapsed + minElapsed ) / 2 << "\n";                   \
        opName.clear();                                                                           \
    }

Definition at line 26 of file imoab_coupler_utils.hpp.

Referenced by main().

#define PUSH_TIMER (   operation)
Value:
{                                         \
        timer_ops = timer.time_since_birth(); \
        opName    = operation;                \
    }

Definition at line 21 of file imoab_coupler_utils.hpp.

Referenced by main().


Function Documentation

int create_group_and_comm ( int  start,
int  end,
MPI_Group  worldGroup,
MPI_Group *  group,
MPI_Comm comm 
)

Definition at line 40 of file imoab_coupler_utils.hpp.

References CHECKIERR, groupTasks, ierr, and MPI_COMM_WORLD.

Referenced by main().

{
    std::vector< int > groupTasks;
    groupTasks.resize( end - start + 1, 0 );
    for( int i = start; i <= end; i++ )
        groupTasks[i - start] = i;

    int ierr = MPI_Group_incl( worldGroup, end - start + 1, &groupTasks[0], group );
    CHECKIERR( ierr, "Cannot create group" )

    ierr = MPI_Comm_create( MPI_COMM_WORLD, *group, comm );
    CHECKIERR( ierr, "Cannot create comm" )

    return 0;
}
int create_joint_comm_group ( MPI_Group  agroup,
MPI_Group  bgroup,
MPI_Group *  abgroup,
MPI_Comm abcomm 
)

Definition at line 56 of file imoab_coupler_utils.hpp.

References CHECKIERR, ierr, and MPI_COMM_WORLD.

Referenced by main().

{
    int ierr = MPI_Group_union( agroup, bgroup, abgroup );
    CHECKIERR( ierr, "Cannot create joint union group" )

    ierr = MPI_Comm_create( MPI_COMM_WORLD, *abgroup, abcomm );
    CHECKIERR( ierr, "Cannot create joint communicator from union group" )

    return 0;
}
int setup_component_coupler_meshes ( iMOAB_AppID  cmpId,
int  cmpTag,
iMOAB_AppID  cplCmpId,
int  cmpcouTag,
MPI_Comm cmpcomm,
MPI_Group *  cmpPEGroup,
MPI_Comm coucomm,
MPI_Group *  cplPEGroup,
MPI_Comm cmpcoucomm,
std::string &  filename,
std::string &  readopts,
int  nghlay,
int  repartitioner_scheme 
)

Definition at line 67 of file imoab_coupler_utils.hpp.

References CHECKIERR, ierr, iMOAB_FreeSenderBuffers, iMOAB_LoadMesh, iMOAB_ReceiveMesh, iMOAB_SendMesh, and nghlay.

Referenced by main().

{
    int ierr = 0;
    if( *cmpcomm != MPI_COMM_NULL )
    {
        // load first mesh
        ierr = iMOAB_LoadMesh( cmpId, filename.c_str(), readopts.c_str(), &nghlay, filename.length(),
                               readopts.length() );
        CHECKIERR( ierr, "Cannot load component mesh" )

        // then send mesh to coupler pes
        ierr = iMOAB_SendMesh( cmpId, cmpcoucomm, cplPEGroup, &cmpcouTag,
                               &repartitioner_scheme );  // send to  coupler pes
        CHECKIERR( ierr, "cannot send elements" )
    }
    // now, receive mesh, on coupler communicator; first mesh 1, atm
    if( *coucomm != MPI_COMM_NULL )
    {

        ierr = iMOAB_ReceiveMesh( cplCmpId, cmpcoucomm, cmpPEGroup,
                                  &cmpTag );  // receive from component
        CHECKIERR( ierr, "cannot receive elements on coupler app" )
    }

    // we can now free the sender buffers
    if( *cmpcomm != MPI_COMM_NULL )
    {
        int context_id = -1;
        ierr = iMOAB_FreeSenderBuffers( cmpId, &context_id );
        CHECKIERR( ierr, "cannot free buffers used to send atm mesh" )
    }
    return 0;
}
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines