reduce_mmx.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	reduce_mmx (rin, rinout, len, type)

Function/Subroutine Documentation

reduce_mmx()

subroutine reduce_mmx	(		rin,
			rinout,
		integer	len,
		integer	type )

Definition at line 27 of file reduce_mmx.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER LEN, TYPE
      my_real
     .        rin(*), rinout(*)
C-----------------------------------------------
C   L o c a l  V a r i a b l e s
C-----------------------------------------------
#ifdef MPI
      INTEGER I, L, I_LEN
      parameter(i_len = 46)
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
      l = 1
      DO i=1,len
C X MIN/MAX
        rinout(l)   = max(rinout(l),rin(l))
        rinout(l+1) = max(rinout(l+1),rin(l+1))
        rinout(l+2) = max(rinout(l+2),rin(l+2))
        rinout(l+3) = min(rinout(l+3),rin(l+3))
        rinout(l+4) = min(rinout(l+4),rin(l+4))
        rinout(l+5) = min(rinout(l+5),rin(l+5))
 
        rinout(l+6) = max(rinout(l+6),rin(l+6))
        rinout(l+7) = max(rinout(l+7),rin(l+7))
        rinout(l+8) = max(rinout(l+8),rin(l+8))
 
        rinout(l+ 9) = min(rinout(l+9),rin(l+9))
        rinout(l+10) = min(rinout(l+10),rin(l+10))
        rinout(l+11) = min(rinout(l+11),rin(l+11))
 
        rinout(l+12) = max(rinout(l+12),rin(l+12))
        rinout(l+13) = max(rinout(l+13),rin(l+13))
        rinout(l+14) = max(rinout(l+14),rin(l+14))
 
        rinout(l+15) = min(rinout(l+15),rin(l+15))
        rinout(l+16) = min(rinout(l+16),rin(l+16))
        rinout(l+17) = min(rinout(l+17),rin(l+17))
 
        l = l + 18
C
        rinout(l)   = max(rinout(l),rin(l))
        rinout(l+1) = max(rinout(l+1),rin(l+1))
        rinout(l+2) = max(rinout(l+2),rin(l+2))
        rinout(l+3) = min(rinout(l+3),rin(l+3))
        rinout(l+4) = min(rinout(l+4),rin(l+4))
        rinout(l+5) = min(rinout(l+5),rin(l+5))
        rinout(l+6) = max(rinout(l+6),rin(l+6))
        rinout(l+7) = max(rinout(l+7),rin(l+7))
        rinout(l+8) = max(rinout(l+8),rin(l+8))
        rinout(l+9) = min(rinout(l+9),rin(l+9))
        rinout(l+10) = min(rinout(l+10),rin(l+10))
        rinout(l+11) = min(rinout(l+11),rin(l+11))
 
        l = l + 12
C V MIN/MAX
        rinout(l)   = max(rinout(l),rin(l))
        rinout(l+1) = max(rinout(l+1),rin(l+1))
        rinout(l+2) = max(rinout(l+2),rin(l+2))
        rinout(l+3) = min(rinout(l+3),rin(l+3))
        rinout(l+4) = min(rinout(l+4),rin(l+4))
        rinout(l+5) = min(rinout(l+5),rin(l+5))
        l = l + 6
C
        rinout(l)   = max(rinout(l),rin(l))
        rinout(l+1) = max(rinout(l+1),rin(l+1))
        rinout(l+2) = max(rinout(l+2),rin(l+2))
        rinout(l+3) = min(rinout(l+3),rin(l+3))
        rinout(l+4) = min(rinout(l+4),rin(l+4))
        rinout(l+5) = min(rinout(l+5),rin(l+5))
        l = l + 6
C NEWFRONT / SIZE
        rinout(l) = rinout(l)+rin(l)
        l = l + 1
C TZINF
        rinout(l) = min(rinout(l),rin(l))              
        l = l + 1
C TZINF
        rinout(l) = max(rinout(l),rin(l))              
        l = l + 1
C T25 main gap changes with thickness change
        rinout(l) = max(rinout(l),rin(l))              
        l = l + 1
      END DO
C
#endif
      RETURN