inter_minmax_node.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "com04_c.inc"
#include "task_c.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	inter_minmax_node (itask, box_limit, x)

Function/Subroutine Documentation

inter_minmax_node()

subroutine inter_minmax_node	(	integer, intent(in)	itask,
		intent(inout)	box_limit,
		intent(in)	x )

Definition at line 34 of file inter_minmax_node.F.

!$COMMENT
!       INTER_MINMAX_NODE description :
!       compute the minimal & maximal positions      
!
!       INTER_MINMAX_NODE organization :
!       loop over the nodes & min+max computation
!$ENDCOMMENT
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
        USE intbufdef_mod  
        USE inter_struct_mod
        USE message_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER, INTENT(in) :: ITASK
        my_real, dimension(6), intent(inout) :: box_limit !< upper & lower bound of the box
        my_real, DIMENSION(3,NUMNOD), INTENT(in) :: x            !   position
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER :: FIRST, LAST
        INTEGER :: I,J
        my_real, DIMENSION(3) :: max_limit,min_limit
!   ----------------------------------------
 
        max_limit(1:3) = -ep30
        min_limit(1:3) = ep30
        first = 1+itask*numnod/nthread
        last = (itask+1)*numnod/nthread
        !   ------------------
        !   find the min & max position 
        DO i=first,last
            DO j=1,3
                max_limit(j) = max(max_limit(j),x(j,i))
                min_limit(j) = min(min_limit(j),x(j,i))
            ENDDO            
        ENDDO
        !   ------------------
 
        !   ------------------
        !   omp reduction 
#include "lockon.inc"
        box_limit(1) = max(box_limit(1),max_limit(1))
        box_limit(2) = max(box_limit(2),max_limit(2))
        box_limit(3) = max(box_limit(3),max_limit(3))
 
        box_limit(4) = min(box_limit(4),min_limit(1))
        box_limit(5) = min(box_limit(5),min_limit(2))
        box_limit(6) = min(box_limit(6),min_limit(3))
#include "lockoff.inc"
        !   ------------------
        CALL my_barrier()
 
        RETURN