Sends and receives vector data
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
complex(kind=real_4), | intent(in), | CONTIGUOUS | :: | msgin(:) |
Data to send |
|
integer, | intent(in) | :: | dest |
Process to send data to |
||
complex(kind=real_4), | intent(out), | CONTIGUOUS | :: | msgout(:) |
Received data |
|
integer, | intent(in) | :: | source |
Process from which to receive |
||
type(mp_comm_type), | intent(in) | :: | comm |
Message passing environment identifier |
SUBROUTINE mp_sendrecv_cv(msgin, dest, msgout, source, comm) !! Sends and receives vector data COMPLEX(kind=real_4), CONTIGUOUS, INTENT(IN) :: msgin(:) !! Data to send INTEGER, INTENT(IN) :: dest !! Process to send data to COMPLEX(kind=real_4), CONTIGUOUS, INTENT(OUT) :: msgout(:) !! Received data INTEGER, INTENT(IN) :: source !! Process from which to receive TYPE(mp_comm_type), INTENT(IN) :: comm !! Message passing environment identifier CHARACTER(LEN=*), PARAMETER :: routineN = 'mp_sendrecv_cv' INTEGER :: handle, ierr #if defined(__parallel) INTEGER :: msglen_in, msglen_out, & recv_tag, send_tag #endif ierr = 0 CALL timeset(routineN, handle) #if defined(__parallel) msglen_in = SIZE(msgin) msglen_out = SIZE(msgout) send_tag = 0 ! cannot think of something better here, this might be dangerous recv_tag = 0 ! cannot think of something better here, this might be dangerous CALL mpi_sendrecv(msgin, msglen_in, MPI_COMPLEX, dest, send_tag, msgout, & msglen_out, MPI_COMPLEX, source, recv_tag, comm%handle, MPI_STATUS_IGNORE, ierr) IF (ierr /= 0) CALL mp_stop(ierr, "mpi_sendrecv @ "//routineN) CALL add_perf(perf_id=7, & msg_size=(msglen_in + msglen_out)*(2*real_4_size)/2) #else MARK_USED(dest) MARK_USED(source) MARK_USED(comm) msgout = msgin #endif CALL timestop(handle) END SUBROUTINE mp_sendrecv_cv