p_dsum.F

Go to the documentation of this file.

!                   ********************************
                    DOUBLE PRECISION FUNCTION p_dsum
!                   ********************************
!
     &(mypart)
!
!***********************************************************************
! PARALLEL   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    SUM OF VALUES 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_dsum=0.d0
      CALL ampi_allreduce(mypart,p_dsum,1,ampi_type,
     &                    ampi_op_sum,comm,ier)
#  else
      CALL mpi_allreduce(mypart,p_dsum,1,mpi_double_precision,
     &                   mpi_sum,comm,ier)
#  endif
!
      IF(ier.NE.0) THEN
        WRITE(lu,*) 'P_DSUM: ERROR IN MPI_ALLREDUCE'
        WRITE(lu,*) 'MPI ERROR: ',ier
        CALL plante(1)
        stop
      ENDIF
      RETURN
#else
      WRITE(lu,*) 'CALL OF P_DSUM IN ITS VOID VERSION'
!
      p_dsum=mypart
!
#endif
!
!-----------------------------------------------------------------------
!
      RETURN
      END