Gathers rank-3 data from all processes and all processes receive the same data
see mp_allgather_c12
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
complex(kind=real_4), | intent(in), | CONTIGUOUS | :: | msgout(:,:,:) |
Rank-3 data to send |
|
complex(kind=real_4), | intent(out), | CONTIGUOUS | :: | msgin(:,:,:) | ||
integer, | intent(in) | :: | gid | |||
integer, | intent(out) | :: | request |
SUBROUTINE mp_iallgather_c33(msgout, msgin, gid, request)
!! Gathers rank-3 data from all processes and all processes receive the
!! same data
!! @note see mp_allgather_c12
COMPLEX(kind=real_4), CONTIGUOUS, INTENT(IN) :: msgout(:, :, :)
!! Rank-3 data to send
COMPLEX(kind=real_4), CONTIGUOUS, INTENT(OUT) :: msgin(:, :, :)
INTEGER, INTENT(IN) :: gid
INTEGER, INTENT(OUT) :: request
CHARACTER(LEN=*), PARAMETER :: routineN = 'mp_iallgather_c33'
INTEGER :: handle, ierr
#if defined(__parallel)
INTEGER :: rcount, scount
#endif
ierr = 0
CALL timeset(routineN, handle)
#if defined(__parallel)
scount = SIZE(msgout(:, :, :))
rcount = scount
CALL MPI_IALLGATHER(msgout, scount, MPI_COMPLEX, &
msgin, rcount, MPI_COMPLEX, &
gid, request, ierr)
IF (ierr /= 0) CALL mp_stop(ierr, "mpi_iallgather @ "//routineN)
#else
MARK_USED(gid)
msgin(:, :, :) = msgout(:, :, :)
request = mp_request_null
#endif
CALL timestop(handle)
END SUBROUTINE mp_iallgather_c33