A routine to split a subcommunicator into sub-groups based on each procs specific colour "col". The sub communicator's handle is passed back in new_comm
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | col |
Processors colour |
||
type(comm_type), | intent(out) | :: | new_comm |
The new sub communicator's handle |
||
integer, | intent(in) | :: | sub |
subroutine split_nokey_to_commtype_sub (col,new_comm,sub)
implicit none
integer, intent(in) :: col !< Processors colour
type(comm_type), intent(out) :: new_comm !< The new sub communicator's handle
integer, intent(in) :: sub
# ifdef MPI
integer :: ierror
# endif
integer :: comm_id, nproc, iproc
!Split the comm group, note we use a constant key of 0 across all procs
!meaning that the rank order is the same in the old and new communicators
# ifdef MPI
!$OMP MASTER
call time_message(.false., time_mp_other, ' MPI Overheads')
!$OMP END MASTER
call mpi_comm_split(sub,col,0,comm_id,ierror)
!$OMP MASTER
call time_message(.false., time_mp_other, ' MPI Overheads')
!$OMP END MASTER
!If the resulting communicator is the null communicator then
!just return the default type with the appropriate id.
if(comm_id == mp_comm_null) then
new_comm%id = comm_id
return
end if
# else
UNUSED_DUMMY(col); UNUSED_DUMMY(sub)
comm_id = mp_comm_null
# endif
new_comm%id=comm_id
call nproc_comm(comm_id,nproc)
new_comm%nproc=nproc
call rank_comm(comm_id,iproc)
new_comm%iproc=iproc
new_comm%proc0=iproc.eq.0
end subroutine split_nokey_to_commtype_sub