Element-wise sum of a rank-1 array on all processes.
Note
see mp_sum_z
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
complex(kind=real_8), | intent(inout), | CONTIGUOUS | :: | msg(:) |
Vector to sum and result |
|
type(mp_comm_type), | intent(in) | :: | gid | |||
type(mp_request_type), | intent(inout) | :: | request |
SUBROUTINE mp_isum_zv(msg, gid, request) !! Element-wise sum of a rank-1 array on all processes. !! @note see mp_sum_z COMPLEX(kind=real_8), CONTIGUOUS, INTENT(INOUT) :: msg(:) !! Vector to sum and result TYPE(mp_comm_type), INTENT(IN) :: gid TYPE(mp_request_type), INTENT(INOUT) :: request CHARACTER(LEN=*), PARAMETER :: routineN = 'mp_isum_zv' INTEGER :: handle, ierr #if defined(__parallel) INTEGER :: msglen #endif ierr = 0 CALL timeset(routineN, handle) #if defined(__parallel) msglen = SIZE(msg) IF (msglen > 0) THEN CALL mpi_iallreduce(MPI_IN_PLACE, msg, msglen, MPI_DOUBLE_COMPLEX, MPI_SUM, gid%handle, request%handle, ierr) IF (ierr /= 0) CALL mp_stop(ierr, "mpi_iallreduce @ "//routineN) ELSE request = mp_request_null END IF CALL add_perf(perf_id=23, msg_size=msglen*(2*real_8_size)) #else MARK_USED(msg) MARK_USED(gid) request = mp_request_null #endif CALL timestop(handle) END SUBROUTINE mp_isum_zv