rmatforp.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/.
!||====================================================================
!||    rmatforp          ../engine/source/materials/mat/mat013/rmatforp.F
!||--- called by ------------------------------------------------------
!||    resol             ../engine/source/engine/resol.F
!||--- calls      -----------------------------------------------------
!||    rmatpon           ../engine/source/materials/mat/mat013/rmatpon.F
!||    spmd_exch_a_rm6   ../engine/source/mpi/kinematic_conditions/spmd_exch_a_rm6.F
!||    startime          ../engine/source/system/timer_mod.F90
!||    stoptime          ../engine/source/system/timer_mod.F90
!||--- uses       -----------------------------------------------------
!||    timer_mod         ../engine/source/system/timer_mod.F90
!||====================================================================
      SUBROUTINE rmatforp(TIMERS,
     1   A       ,AR       ,X        ,VR     ,IN     ,
     2   STIFN   ,STIFR    ,IRBY     ,LNRBY  ,RBY    , 
     3   ICODRBY ,WEIGHT   ,MS       ,V      ,FR_RBM ,
     4   IAD_RBY ,ARBY    ,VRBY     ,ARRBY   ,VRRBY  ,
     5   KIND    ,RBYM6   )
        USE timer_mod
C----6---------------------------------------------------------------7---------8
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      "com01_c.inc"
#include      "com04_c.inc"
#include      "timeri_c.inc"
C-----------------------------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE(timer_), INTENT(INOUT) :: TIMERS
      INTEGER FR_RBM(*),IAD_RBY(*),
     .        ICODRBY(2,*),WEIGHT(*),LNRBY(*),IRBY(NIRBYM,*),
     .        KIND(NRBYM)
C     REAL
      my_real
     .   rby(nfrbym,*) ,a(3,*) ,ar(3,*) ,x(3,*) ,vr(3,*),
     .   in(*) ,stifn(*),stifr(*),ms(*),v(3,*) , 
     .   vrrby(3,*),arby(3,*),arrby(3,*),
     .   vrby(3,*)
      DOUBLE PRECISION RBYM6(6,6,NRBYM)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  K,M
 
C-------------------------------------
C Traitement rigid material  
C-------------------------------------
 
!$OMP SINGLE
 
      k=1
      DO m=1,nrbym
        kind(m) = k
        k  = k  + irby(2,m)
      ENDDO
 
!$OMP END SINGLE
 
!$OMP DO SCHEDULE(DYNAMIC,1)
        DO m=1,nrbym
          k = kind(m)
          CALL rmatpon(
     1      a           ,ar       ,x      ,rby(1,m),lnrby(k) ,
     2      irby(1,m)   ,stifn    ,stifr  ,weight  ,irby(2,m),
     3      rbym6(1,1,m),icodrby  ,arby   ,vrby    ,arrby    ,
     3      vrrby    ,1        )
        ENDDO
!$OMP END DO
 
        IF (nspmd > 1) THEN
!$OMP SINGLE
          IF (imon>0) CALL startime(timers,11)
C
C Comm non multi-thread
C
          CALL spmd_exch_a_rm6(
     1         iad_rby,fr_rbm,iad_rby(nspmd+1),rbym6 )
          IF (imon>0) CALL stoptime(timers,11)
!$OMP END SINGLE
 
        END IF
C
C Traitement fin parith/on multi-thread
C
 
!$OMP DO SCHEDULE(DYNAMIC,1)
        DO m=1,nrbym
          k = kind(m)
          CALL rmatpon(
     1      a           ,ar      ,x     ,rby(1,m),lnrby(k) ,
     2      irby(1,m)   ,stifn   ,stifr ,weight  ,irby(2,m),
     4      rbym6(1,1,m),icodrby ,arby  ,vrby    ,arrby    ,
     3      vrrby       ,2       )
        ENDDO
!$OMP END DO
C
      RETURN
      END