rmatforp.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "timeri_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	rmatforp (timers, a, ar, x, vr, in, stifn, stifr, irby, lnrby, rby, icodrby, weight, ms, v, fr_rbm, iad_rby, arby, vrby, arrby, vrrby, kind, rbym6)

Function/Subroutine Documentation

rmatforp()

subroutine rmatforp	(	type(timer_), intent(inout)	timers,
			a,
			ar,
			x,
			vr,
			in,
			stifn,
			stifr,
		integer, dimension(nirbym,*)	irby,
		integer, dimension(*)	lnrby,
			rby,
		integer, dimension(2,*)	icodrby,
		integer, dimension(*)	weight,
			ms,
			v,
		integer, dimension(*)	fr_rbm,
		integer, dimension(*)	iad_rby,
			arby,
			vrby,
			arrby,
			vrrby,
		integer, dimension(nrbym)	kind,
		double precision, dimension(6,6,nrbym)	rbym6 )

Definition at line 35 of file rmatforp.F.

        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