spmd_allreduce_db_comm.F File Reference

#include "implicit_f.inc"
#include "spmd.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	spmd_allreduce_db_comm (value, res, my_size, my_operation, my_comm)

Function/Subroutine Documentation

spmd_allreduce_db_comm()

subroutine spmd_allreduce_db_comm	(	real(kind=8), dimension(my_size), intent(in)	value,
		real(kind=8), dimension(my_size), intent(out)	res,
		integer, intent(in)	my_size,
		character(len=*), intent(in)	my_operation,
		integer, intent(in)	my_comm )

Definition at line 27 of file spmd_allreduce_db_comm.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#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   D u m m y   A r g u m e n t s
C-----------------------------------------------
        CHARACTER(len=*), INTENT(in) :: MY_OPERATION
        INTEGER, INTENT(in) :: MY_SIZE,MY_COMM
        REAL(kind=8), dimension(my_size), INTENT(in) :: VALUE
        REAL(kind=8), dimension(my_size), INTENT(out) :: res
 
#ifdef MPI
C-----------------------------------------------
C   L o c a l  V a r i a b l e s
C-----------------------------------------------
        INTEGER :: error
        INTEGER :: MPI_OPERATION
!       ----------------------------------------
        
        IF(my_operation(1:3)=="SUM") THEN
                mpi_operation=mpi_sum
        ELSEIF(my_operation(1:4)=="PROD") THEN
                mpi_operation=mpi_prod
        ELSEIF(my_operation(1:3)=="MAX") THEN
                mpi_operation=mpi_max
        ELSEIF(my_operation(1:3)=="MIN") THEN
                mpi_operation=mpi_min               
        ENDIF
 
!       ------------------
!       mpi comm allreduce
        CALL mpi_allreduce(VALUE,res,my_size,mpi_double_precision,mpi_operation,my_comm,error)
!       ------------------
#else
        res = 0
#endif
        RETURN