p_isum_array.F

Go to the documentation of this file.

!                   ***********************
                    SUBROUTINE p_isum_array
!                   ***********************
!
     &(send_buffer, recv_buffer, ncount, ierr)
!
!***********************************************************************
! PARALLEL   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    SUM VALUES FROM ALL PROCESSES AND DISTRIBUTES THE RESULT BACK TO ALL PROCESSES
!
!history  AUDOUIN YOANN
!+        25/04/2013
!+
!+   CREATED
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| SEND_BUFFER       |-->| SEND BUFFER.
!| RECV_BUFFER       |<--| RECV BUFFER.
!| NCOUNT            |-->| SIZE IF THE BUFFERS
!| IERR              |<--| RETURN VALUE OF THE MPI_FUNCTION
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_parallel
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)  :: NCOUNT
      INTEGER, INTENT(IN)  :: SEND_BUFFER(ncount)
      INTEGER, INTENT(OUT) :: RECV_BUFFER(ncount)
      INTEGER, INTENT(OUT) :: IERR
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
!
!-----------------------------------------------------------------------
!
      ierr = 0
      recv_buffer = 0
#if defined HAVE_MPI
      CALL mpi_allreduce(send_buffer,recv_buffer,ncount,mpi_integer,
     &                   mpi_sum,comm,ierr)
!
      IF(ierr.NE.0) THEN
        WRITE(lu,*) 'P_ISUM_ARRAY: ERROR IN MPI_ALLREDUCE'
        WRITE(lu,*) 'MPI ERROR ',ierr
        CALL plante(1)
        stop
      ENDIF
#else
      WRITE(lu,*)'CALL OF P_ISUM_ARRAY IN ITS VOID VERSION'
!
!-----------------------------------------------------------------------
!
#endif
!
!-----------------------------------------------------------------------
!
      RETURN
      END