p_dmax.F

Go to the documentation of this file.

!                   ********************************
                    DOUBLE PRECISION FUNCTION p_dmax
!                   ********************************
!
     &(mypart)
!
!***********************************************************************
! PARALLEL   V6P3                                   21/08/2010
!***********************************************************************
!
!brief    MAXIMUM VALUE FROM ALL THE PROCESSORS.
!
!history  RAINER JOHANNI (SGI MUNICH)
!+        **/10/1999
!+
!+   ADAPTED FOR MPI
!
!history  J.A. JANKOWSKI (BAW KARLSRUHE)
!+        28/12/1999
!+
!+   RELEASE 5.0 MODIFIED
!
!history  J-M HERVOUET (LNHE)
!+        23/06/2008
!+        V5P9
!+
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| MYPART         |-->| SEND BUFFER
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_parallel
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION, INTENT(IN) :: MYPART
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
#if defined HAVE_MPI
      INTEGER IER
!
!-----------------------------------------------------------------------
!

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