Gathers data from all processes to one
Data length All data (msg) is equal-sized
MPI mapping mpi_gather
Note
see mp_gather_z
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
complex(kind=real_8), | intent(in), | CONTIGUOUS | :: | msg(:) |
Datum to send to root |
|
complex(kind=real_8), | intent(out), | CONTIGUOUS | :: | msg_gather(:) | ||
integer, | intent(in) | :: | root | |||
type(mp_comm_type), | intent(in) | :: | gid |
SUBROUTINE mp_gather_zv(msg, msg_gather, root, gid) !! Gathers data from all processes to one !! !! Data length !! All data (msg) is equal-sized !! !! MPI mapping !! mpi_gather !! @note see mp_gather_z COMPLEX(kind=real_8), CONTIGUOUS, INTENT(IN) :: msg(:) !! Datum to send to root COMPLEX(kind=real_8), CONTIGUOUS, INTENT(OUT) :: msg_gather(:) INTEGER, INTENT(IN) :: root TYPE(mp_comm_type), INTENT(IN) :: gid CHARACTER(LEN=*), PARAMETER :: routineN = 'mp_gather_zv' INTEGER :: handle, ierr, msglen ierr = 0 CALL timeset(routineN, handle) msglen = SIZE(msg) #if defined(__parallel) CALL mpi_gather(msg, msglen, MPI_DOUBLE_COMPLEX, msg_gather, & msglen, MPI_DOUBLE_COMPLEX, root, gid%handle, ierr) IF (ierr /= 0) CALL mp_stop(ierr, "mpi_gather @ "//routineN) CALL add_perf(perf_id=4, msg_size=msglen*(2*real_8_size)) #else MARK_USED(root) MARK_USED(gid) msg_gather = msg #endif CALL timestop(handle) END SUBROUTINE mp_gather_zv