Find on MPI rank 'toRank' MAX/MIN/SUM of 2D array (assumed-shape) across all MPI ranks.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(inout) | :: | arr(:,:) | |||
| character(len=3), | intent(in) | :: | parOP | |||
| integer, | intent(in) | :: | toRank |
SUBROUTINE par_Reduce_array2D(arr,parOP,toRank) ! MODULES USE MODgvec_Globals, ONLY : wp,myRank IMPLICIT NONE !-------------------------------------------------------------------------------------------------------------------------------- ! INPUT VARIABLES REAL(wp), INTENT(INOUT) :: arr(:,:) CHARACTER(LEN=3), INTENT(IN) :: parOP INTEGER, INTENT(IN) :: toRank ! =0 by default !-------------------------------------------------------------------------------------------------------------------------------- ! LOCAL VARIABLES # if MPI MPI_op_TYPE :: mpiOP INTEGER :: ierr INTEGER :: sz !================================================================================================================================ ! BODY sz=SIZE(arr) SELECT CASE(parOP) CASE('MAX') mpiOP=MPI_MAX CASE('MIN') mpiOP=MPI_MIN CASE('SUM') mpiOP=MPI_SUM END SELECT IF (myRank.EQ.toRank) THEN CALL MPI_Reduce(MPI_IN_PLACE, arr, sz, dType, mpiOP, toRank, worldComm, ierr) ELSE CALL MPI_Reduce(arr, arr, sz, dType, mpiOP, toRank, worldComm, ierr) END IF # endif END SUBROUTINE par_Reduce_array2D