forints.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "param_c.inc"
#include "vect01_c.inc"
#include "scr17_c.inc"
#include "task_c.inc"
#include "parit_c.inc"
#include "scr23_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	forints (pm, geo, x, a, ar, v, vr, ms, in, w, elbuf, val2, veul, fv, stifn, stifr, fsky, tf, bufmat, partsav, fani, fsav, skew, dt2t, iads, iparg, npc, ixs, neltst, ityptst, ipart, iparts, itab, fskyi, bufgeo, kxx, ixx, isky, ipartx, gresav, grth, igrth, elbuf_tab)

Function/Subroutine Documentation

forints()

subroutine forints	(		pm,
			geo,
			x,
			a,
			ar,
			v,
			vr,
			ms,
			in,
			w,
			elbuf,
			val2,
			veul,
			fv,
			stifn,
			stifr,
			fsky,
			tf,
			bufmat,
			partsav,
			fani,
			fsav,
			skew,
			dt2t,
		integer, dimension(8,*)	iads,
		integer, dimension(nparg,*)	iparg,
		integer, dimension(*)	npc,
		integer, dimension(nixs,*)	ixs,
		integer	neltst,
		integer	ityptst,
		integer, dimension(lipart1,*)	ipart,
		integer, dimension(*)	iparts,
		integer, dimension(*)	itab,
			fskyi,
			bufgeo,
		integer, dimension(nixx,*)	kxx,
			ixx,
		integer, dimension(*)	isky,
		integer, dimension(*)	ipartx,
			gresav,
		integer, dimension(*)	grth,
		integer, dimension(*)	igrth,
		type(elbuf_struct_), dimension(ngroup), target	elbuf_tab )

Definition at line 36 of file forints.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE initbuf_mod
      USE elbufdef_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "param_c.inc"
#include      "vect01_c.inc"
#include      "scr17_c.inc"
#include      "task_c.inc"
#include      "parit_c.inc"
#include      "scr23_c.inc"
C-----------------------------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXS(NIXS,*),IPART(LIPART1,*) ,IPARTS(*) ,
     .   NPC(*), IPARG(NPARG,*),IADS(8,*),NELTST,ITYPTST,
     .   KXX(NIXX,*),I XX(*), ITAB(*),
     .   IPARTX(*), ISKY(*),GRTH(*),IGRTH(*)
C     REAL
      my_real
     .   x(3,*)    ,v(3,*)   ,vr(3,*),
     .   ms(*)   ,in(*)   ,pm(npropm,*),skew(lskew,*),geo(npropg,*),
     .   bufmat(*) ,w(3,*)    ,veul(*),tf(*)  ,
     .   fsav(nthvki,*)       ,elbuf(*) ,
     .   fv(*)     ,a(3,*)    ,ar(3,*)     ,val2(*)  ,
     .   fani(3,*) ,partsav(*),stifn(*) ,stifr(*),
     .   fsky(*) ,
     .   dt2t,fskyi(lskyi,nfskyi) ,bufgeo(*),gresav(*)
C
      TYPE(ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,N, NG, NVC, MLW, JFT, JLT ,K ,
     .   K1, K2, KAD,IAD2,NF1,IPRI,NGLOC, NELEM, OFFSET,
     .   NSG, NEL,ISOLNOD
      INTEGER IADUIX,IADUX ,IADUV ,IADUVR,IADUMS,
     .        IADUIN,IADUSM,IADUSR,IADUMV,IADURV,
     .        IADUF ,IADUTQ
C-----------------------------------------------
C
      DO 240 ng=1,ngroup
          IF(iparg(8,ng)==1.OR.iparg(5,ng)/=100)GOTO 240
          IF (iddw>0) CALL startimeg(ng)
          DO nelem = 1,iparg(2,ng),nvsiz
C---
C
            offset = nelem - 1
C
           CALL initbuf(iparg    ,ng      ,
     2        mtn     ,nel     ,nft     ,iad     ,ity     ,
     3        npt     ,jale    ,ismstr  ,jeul    ,jtur    ,
     4        jthe    ,jlag    ,jmult   ,jhbe    ,jivf    ,
     5        nvaux   ,jpor    ,jcvt    ,jclose  ,jplasol ,
     6        irep    ,iint    ,igtyp   ,israt   ,isrot   ,
     7        icsen   ,isorth  ,isorthg ,ifailure,jsms    )
C
            nsg   =iparg(10,ng)
            nvc   =iparg(19,ng)
            jhbe  =iparg(23,ng)
            igre  =iparg(51,ng)
C
            lft=1
            llt=min(nvsiz,nel-nelem+1)
C           Multi-purpose element Groups are constituted,
C           but Multi-purpose element is not supposed to be vectorized ..
             CALL xforc3(
     1   nft,          nel,          geo,          pm,
     2   itab,         kxx,          ixx,          x,
     3   v,            vr,           a,            ar,
     4   elbuf(iad),   npc,          tf,           skew,
     5   dt2t,         neltst,       ityptst,      stifn,
     6   stifr,        ms,           in,           fskyi,
     7   isky,         partsav,      ipartx,       bufmat,
     8   bufgeo,       gresav,       grth,         igrth,
     9   elbuf_tab(ng),igre)
            IF (iddw>0) CALL stoptimeg(ng)
         END DO
 240  CONTINUE
C----------------------------------
      RETURN