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.

Macros
#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)

Macro Definition Documentation

CHECKIERR

#define CHECKIERR	(		rc,
			message
	)

Value:

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

Definition at line 14 of file imoab_coupler_utils.hpp.

POP_TIMER

#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.

PUSH_TIMER

#define PUSH_TIMER

(

operation

)

Value:

 { \
 timer_ops = timer.time_since_birth(); \
 opName = operation; \
 }

Definition at line 21 of file imoab_coupler_utils.hpp.

Function Documentation

create_group_and_comm()

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.

{
 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;
}

References CHECKIERR, groupTasks, ierr, and MPI_COMM_WORLD.

Referenced by main().

create_joint_comm_group()

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.

{
 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;
}

References CHECKIERR, ierr, and MPI_COMM_WORLD.

Referenced by main().

setup_component_coupler_meshes()

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.

{
 int ierr = 0;
 if( *cmpcomm != MPI_COMM_NULL )
 {
 // load first mesh
 ierr = iMOAB_LoadMesh( cmpId, filename.c_str(), readopts.c_str(), &nghlay );
 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 = cmpcouTag;
 ierr = iMOAB_FreeSenderBuffers( cmpId, &context_id );
 CHECKIERR( ierr, "cannot free buffers used to send atm mesh" )
 }
 return 0;
}

References CHECKIERR, filename, ierr, iMOAB_LoadMesh(), nghlay, and readopts().

Referenced by main().