spmd_fl_sum.F File Reference

#include "implicit_f.inc"
#include "spmd.inc"
#include "com01_c.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	spmd_fl_sum (lsum, len, lsumt)

Function/Subroutine Documentation

spmd_fl_sum()

subroutine spmd_fl_sum	(		lsum,
		integer	len,
			lsumt )

Definition at line 36 of file spmd_fl_sum.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
        USE spmd_comm_world_mod, ONLY : spmd_comm_world
#include      "implicit_f.inc"
C-----------------------------------------------------------------
C   M e s s a g e   P a s s i n g
C-----------------------------------------------
#include "spmd.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER LEN
      my_real
     .        lsum(*), lsumt(*)
C-----------------------------------------------
C   L o c a l  V a r i a b l e s
C-----------------------------------------------
#ifdef MPI
      INTEGER MSGOFF, MSGOFF2,ITAG, REQ(2), TSTAT(MPI_STATUS_SIZE,2), IERR,
     .        I, II, REQ1(NSPMD-1), J,
     .        STAT1(MPI_STATUS_SIZE,NSPMD-1)
      my_real
     .        lsumc(len)
C
      DATA msgoff/186/
      DATA msgoff2/187/
C
      IF (ispmd/=0) THEN
         itag=msgoff
         CALL mpi_isend(lsum, len, real, it_spmd(1), itag,
     .                  spmd_comm_world, req(1), ierr)
C
         itag=msgoff2
         CALL mpi_irecv(lsumt, len, real, it_spmd(1), itag,
     .                  spmd_comm_world, req(2), ierr)
C
         CALL mpi_waitall(2, req, tstat, ierr)
      ELSE
         DO i=1,len
            lsumt(i)=lsum(i)
         ENDDO
C
         DO i=1,nspmd-1
            ii=i+1
            itag=msgoff
            CALL mpi_irecv(lsumc, len, real, it_spmd(ii), itag,
     .                     spmd_comm_world, req1(1), ierr)
            CALL mpi_wait(req1(1), stat1(1,1), ierr)
C
            DO j=1,len
               lsumt(j)=lsumt(j)+lsumc(j)
            ENDDO
         ENDDO
C
         DO i=1,nspmd-1
            ii=i+1
            itag=msgoff2
            CALL mpi_isend(lsumt, len, real, it_spmd(ii), itag,
     .                     spmd_comm_world, req1(i), ierr)
         ENDDO
C
         CALL mpi_waitall(nspmd-1, req1, stat1, ierr)
      ENDIF
C
#endif
      RETURN