Find MAX/MIN/SUM of 2D array (assumed-shape) across all MPI ranks and bradcast result back to all MPI ranks.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(inout) | :: | arr(:,:) | |||
| character(len=3), | intent(in) | :: | parOP |
SUBROUTINE par_AllReduce_array2D(arr,parOP) ! MODULES USE MODgvec_Globals, ONLY : wp IMPLICIT NONE !-------------------------------------------------------------------------------------------------------------------------------- ! INPUT VARIABLES REAL(wp), INTENT(INOUT) :: arr(:,:) CHARACTER(LEN=3), INTENT(IN) :: parOP !-------------------------------------------------------------------------------------------------------------------------------- ! 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 CALL MPI_AllReduce(MPI_IN_PLACE, arr, sz, dType, mpiOP, worldComm, ierr) # endif END SUBROUTINE par_AllReduce_array2D