reduce_mmx.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
C
!||====================================================================
!||    reduce_mmx   ../engine/source/mpi/interfaces/reduce_mmx.F
!||====================================================================
      SUBROUTINE reduce_mmx(RIN,RINOUT,LEN,TYPE)
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
      END