spmd_reduce_db.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	spmd_reduce_db (value, res, my_size, my_rank, my_operation)

Function/Subroutine Documentation

spmd_reduce_db()

subroutine spmd_reduce_db	(	real(kind=8), dimension(my_size), intent(in)	value,
		real(kind=8), dimension(my_size), intent(out)	res,
		integer, intent(in)	my_size,
		integer, intent(in)	my_rank,
		character(len=4), intent(in)	my_operation )

Definition at line 31 of file spmd_reduce_db.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   D u m m y   A r g u m e n t s
C-----------------------------------------------
        CHARACTER(len=4), INTENT(in) :: MY_OPERATION
        INTEGER, INTENT(in) :: MY_SIZE,MY_RANK
        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 reduce
        CALL mpi_reduce(VALUE,res,my_size,mpi_double_precision,mpi_operation,my_rank,spmd_comm_world,error)
!       ------------------
#else
        res = 0
#endif
        RETURN