p_dmax_array.F

Go to the documentation of this file.

!                   ***********************
                    SUBROUTINE p_dmax_array
!                   ***********************
!
     &(mypart, n, values)
!
!***********************************************************************
! PARALLEL   V6P3                                   21/08/2010
!***********************************************************************
!
!brief    MAXIMUM VALUE FROM ALL THE PROCESSORS.
!
!history  Y AUDOUIN (LNHE)
!+        09/03/2018
!+
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| MYPART         |-->| SEND BUFFER
!| N              |-->| SIZE OF BUFFER
!| VALUES         |-->| RECEIVED BUFFER
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_parallel
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN) :: N
      DOUBLE PRECISION, INTENT(IN) :: MYPART(n)
      DOUBLE PRECISION, INTENT(INOUT) :: VALUES(n)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
#if defined HAVE_MPI
      INTEGER IER
!
!-----------------------------------------------------------------------
!

#  ifdef COMPAD
      values=0.d0
      CALL ampi_allreduce(mypart,values,n,ampi_type,ampi_op_max,
     &                    comm,ier)
#  else
      CALL mpi_allreduce(mypart,values,n,mpi_double_precision,mpi_max,
     &                   comm,ier)
#  endif
!
      IF(ier.NE.0) THEN
        WRITE(lu,*) 'P_DMAX_ARRAY: ERROR IN MPI_ALLREDUCE'
        WRITE(lu,*) 'MPI ERROR ',ier
        CALL plante(1)
        stop
      ENDIF
      RETURN
#else
!
      values=mypart
!
#endif
!
!-----------------------------------------------------------------------
!
      END