Gathers rank-3 data from all processes and all processes receive the same data
Note
see mp_allgather_z12
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
complex(kind=real_8), | intent(in), | CONTIGUOUS | :: | msgout(:,:,:) |
Rank-3 data to send |
|
complex(kind=real_8), | intent(out), | CONTIGUOUS | :: | msgin(:,:,:,:) | ||
type(mp_comm_type), | intent(in) | :: | gid |
SUBROUTINE mp_allgather_z34(msgout, msgin, gid) !! Gathers rank-3 data from all processes and all processes receive the !! same data !! @note see mp_allgather_z12 COMPLEX(kind=real_8), CONTIGUOUS, INTENT(IN) :: msgout(:, :, :) !! Rank-3 data to send COMPLEX(kind=real_8), CONTIGUOUS, INTENT(OUT) :: msgin(:, :, :, :) TYPE(mp_comm_type), INTENT(IN) :: gid CHARACTER(LEN=*), PARAMETER :: routineN = 'mp_allgather_z34' 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_ALLGATHER(msgout, scount, MPI_DOUBLE_COMPLEX, & msgin, rcount, MPI_DOUBLE_COMPLEX, & gid%handle, ierr) IF (ierr /= 0) CALL mp_stop(ierr, "mpi_allgather @ "//routineN) #else MARK_USED(gid) msgin(:, :, :, 1) = msgout(:, :, :) #endif CALL timestop(handle) END SUBROUTINE mp_allgather_z34