i16main.F File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	i16main (nin, ipari, intbuf_tab, x, v, a, itask, igrnod, eminx, wat, ms, iadll, lll, jll, sll, xll, n_mul_mx, ixs, ixs16, ixs20, nkmax, ixs10, comntag, igrbric)

Function/Subroutine Documentation

i16main()

subroutine i16main	(	integer	nin,
		integer, dimension(npari,ninter)	ipari,
		type(intbuf_struct_), dimension(*)	intbuf_tab,
			x,
			v,
			a,
		integer	itask,
		type (group_), dimension(ngrnod)	igrnod,
			eminx,
			wat,
			ms,
		integer, dimension(*)	iadll,
		integer, dimension(*)	lll,
		integer, dimension(*)	jll,
		integer, dimension(*)	sll,
			xll,
		integer	n_mul_mx,
		integer, dimension(nixs,*)	ixs,
		integer, dimension(8,*)	ixs16,
		integer, dimension(12,*)	ixs20,
		integer	nkmax,
		integer, dimension(6,*)	ixs10,
		integer, dimension(*)	comntag,
		type (group_), dimension(ngrbric)	igrbric )

Definition at line 40 of file i16main.F.

C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE intbufdef_mod
      USE groupdef_mod
      use element_mod , only : nixs
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C
C-------------------------------------------------------------------------------
C     NOM    DIMENSION               DESCRIPTION                       E/S
C-------------------------------------------------------------------------------
C
C    NIN    1                        INTERFACE NUMBER                   E
C                  
C   IPARI   NPARI,NINTER             PARAMETRES D'INTERFACE             E
C                  
C     X     3,NUMNOD                 COORDONNEES                        E
C                  
C     V     3,NUMNOD                 VELOCITIES                           E
C                  
C    EMINX  6*NME<6*NUMELS           MIN MAX DE CHAQUE ELEMENT        TMP_GLOBAL
C                  
C                  
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      "com04_c.inc"
#include      "com08_c.inc"
#include      "param_c.inc"
#include      "task_c.inc"
 
      COMMON /i16tmp/size
      my_real
     .   SIZE
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NIN,ITASK ,N_MUL_MX ,NKMAX   
      INTEGER IPARI(NPARI,NINTER),
     .        IADLL(*) ,LLL(*) ,JLL(*) ,SLL(*) ,
     .        IXS(NIXS,*)  ,IXS16(8,*)  ,IXS20(12,*) ,IXS10(6,*),
     .        COMNTAG(*)
C     REAL
      my_real 
     .   x(3,*), v(3,*), a(3,*), ms(*),eminx(*),wat(*),
     .   xll(*)
 
      TYPE(INTBUF_STRUCT_) INTBUF_TAB(*)
      TYPE (GROUP_)  , DIMENSION(NGRNOD) :: IGRNOD
      TYPE (GROUP_)  , DIMENSION(NGRBRIC) :: IGRBRIC
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NSN,NME,
     .        NME_T,ESH_T,IGN,IGE,MULTIMP,NOINT,MX_CAND,itri 
C     REAL
      my_real
     .   startt, stopt,xx,xy,xz,tzinf,minbox
      SAVE itri
      data itri/0/
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      noint       = ipari(15,nin)
      multimp     = ipari(23,nin)
      ign         = ipari(36,nin)
      ige         = ipari(34,nin)
      nsn   =igrnod(ign)%NENTITY
      nme   =igrbric(ige)%NENTITY
C
      mx_cand=multimp*nsn
      startt = intbuf_tab(nin)%VARIABLES(3)
      stopt  = intbuf_tab(nin)%VARIABLES(11)
      IF(startt>tt) RETURN
      IF(tt>stopt)  RETURN
C -------------------------------------------------------------
       nme_t=nme/nthread
       esh_t = itask*nme_t
       IF(itask==nthread-1)nme_t=nme-(nthread-1)*(nme/nthread)
       SIZE   = zero
       intbuf_tab(nin)%XSLVG(1) = -ep30
       intbuf_tab(nin)%XSLVG(2) = -ep30
       intbuf_tab(nin)%XSLVG(3) = -ep30
       intbuf_tab(nin)%XSLVG(4) =  ep30
       intbuf_tab(nin)%XSLVG(5) =  ep30
       intbuf_tab(nin)%XSLVG(6) =  ep30
       intbuf_tab(nin)%XMSRG(1) = -ep30
       intbuf_tab(nin)%XMSRG(2) = -ep30
       intbuf_tab(nin)%XMSRG(3) = -ep30
       intbuf_tab(nin)%XMSRG(4) =  ep30
       intbuf_tab(nin)%XMSRG(5) =  ep30
       intbuf_tab(nin)%XMSRG(6) =  ep30       
C -------------------------------------------------------------
       CALL my_barrier
C=======================================================================
C      Test to know whether to retriate + initialization
C=======================================================================
       CALL i16crit(
     1   x    ,igrnod(ign)%ENTITY,igrbric(ige)%ENTITY ,nsn  ,eminx      ,
     2   nme    ,itask       ,intbuf_tab(nin)%XSAV ,ixs  ,ixs16      ,
     3   ixs20    ,ixs10       ,v       ,a  ,intbuf_tab(nin)%XMSRG,
     4   intbuf_tab(nin)%XSLVG )
 
C -------------------------------------------------------------
       CALL my_barrier
C -------------------------------------------------------------
C error in the estimation formula
       tzinf  = intbuf_tab(nin)%VARIABLES(4) * SIZE / nme / 6
       minbox = intbuf_tab(nin)%VARIABLES(5) * SIZE / nme / 6
       xx = max(intbuf_tab(nin)%XSLVG(1)-intbuf_tab(nin)%XMSRG(4),
     *          intbuf_tab(nin)%XMSRG(1)-intbuf_tab(nin)%XSLVG(4))
       xy = max(intbuf_tab(nin)%XSLVG(2)-intbuf_tab(nin)%XMSRG(5),
     *          intbuf_tab(nin)%XMSRG(2)-intbuf_tab(nin)%XSLVG(5))
       xz = max(intbuf_tab(nin)%XSLVG(3)-intbuf_tab(nin)%XMSRG(6),
     *          intbuf_tab(nin)%XMSRG(3)-intbuf_tab(nin)%XSLVG(6))
C -------------------------------------------------------------
       IF(xx**2 +  xy**2 + xz**2 >= tzinf**2) THEN
C=======================================================================
C        TRI
C=======================================================================
          ipari(35,nin) = 0
C -------------------------------------------------------------
         CALL my_barrier
C -------------------------------------------------------------
           CALL i16buce(
     1   igrnod(ign)%ENTITY,ixs   ,ixs16   ,ixs20  ,igrbric(ige)%ENTITY ,
     2   nme_t      ,wat   ,nsn    ,intbuf_tab(nin)%CAND_E,intbuf_tab(nin)%CAND_N,
     3   noint      ,ipari(35,nin) ,tzinf   ,minbox    ,eminx      ,
     4   intbuf_tab(nin)%XSAV,itask   ,x    ,v    ,a        ,
     5   mx_cand    ,ixs10   ,esh_t)
C -------------------------------------------------------------
         CALL my_barrier
C -------------------------------------------------------------
       ENDIF
C=======================================================================
C      CONTACT
C=======================================================================
       CALL i8lagm(
     1   x   ,v         ,lll    ,jll     ,sll        ,
     2   xll   ,intbuf_tab(nin)%CAND_N,intbuf_tab(nin)%CAND_E,ipari(35,nin),ixs  ,
     3   iadll   ,eminx        ,igrnod(ign)%ENTITY,igrbric(ige)%ENTITY,
     4   n_mul_mx     ,itask   ,a         ,ipari(30,nin),
     5   nin   ,nkmax        ,comntag   )
       CALL i10lagm(
     1   x       ,v            ,lll       ,jll           ,sll         ,
     2   xll     ,intbuf_tab(nin)%CAND_N,intbuf_tab(nin)%CAND_E,ipari(35,nin),ixs       ,
     3   ixs10   ,iadll        ,eminx     ,igrnod(ign)%ENTITY,igrbric(ige)%ENTITY,
     4   n_mul_mx     ,itask     ,a            ,ipari(30,nin),
     5   nin     ,nkmax        ,comntag   )
       CALL i16lagm(
     1   x   ,v         ,lll    ,jll     ,sll        ,
     2   xll   ,intbuf_tab(nin)%CAND_N,intbuf_tab(nin)%CAND_E,ipari(35,nin),ixs ,
     3   ixs16   ,iadll        ,eminx    ,igrnod(ign)%ENTITY,igrbric(ige)%ENTITY,
     4   n_mul_mx     ,itask   ,a         ,ipari(30,nin),
     5   nin   ,nkmax        ,comntag   )
       CALL i20lagm(
     1   x   ,v         ,lll    ,jll     ,sll        ,
     2   xll   ,intbuf_tab(nin)%CAND_N,intbuf_tab(nin)%CAND_E,ipari(35,nin),ixs       ,
     3   ixs20   ,iadll        ,eminx    ,igrnod(ign)%ENTITY,igrbric(ige)%ENTITY,
     4   n_mul_mx     ,itask   ,a         ,ipari(30,nin),
     5   nin   ,nkmax        ,comntag   )
!
      CALL my_barrier
!
      RETURN