i17crit.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "com04_c.inc"
#include "task_c.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i17crit (x, v, a, nelem, eminx, nme, itask, xsav, ixs, ixs16, ixs20, xmsr_g, size_t_g)

Function/Subroutine Documentation

i17crit()

subroutine i17crit	(		x,
			v,
			a,
		integer, dimension(*)	nelem,
			eminx,
		integer	nme,
		integer	itask,
			xsav,
		integer, dimension(nixs,*)	ixs,
		integer, dimension(8,*)	ixs16,
		integer, dimension(12,*)	ixs20,
			xmsr_g,
			size_t_g )

Definition at line 33 of file i17crit.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ITASK,NME,
     .        NELEM(*),IXS(NIXS,*),IXS16(8,*),IXS20(12,*)
      my_real
     .   x(3,*),v(3,*),a(3,*),xsav(3,*),eminx(6,*),
     .   xmsr_g(6),size_t_g
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NSNF,NMEF,NSNL,NMEL,I, J, II, K,I16,I20,LFT16,LLT16,
     .        LFT20,LLT20,INDEX16(MVSIZ),INDEX20(MVSIZ)
      my_real
     .   xmsr(6), size_t
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
      nmef = 1 + itask*nme / nthread
      nmel = (itask+1)*nme / nthread
C--------------------------------------------------------------
C     0- CALCUL DU CRITERE POUR SAVOIR SI ON DOIT TRIER OU NON
C--------------------------------------------------------------
      xmsr(1) = -ep30
      xmsr(2) = -ep30
      xmsr(3) = -ep30
      xmsr(4) =  ep30
      xmsr(5) =  ep30
      xmsr(6) =  ep30    
      size_t  =  zero
C------------------------------------
C     CALCUL DES BORNES DES ELEMENTS
C------------------------------------
      DO i=nmef,nmel
        eminx(1,i) = ep30
        eminx(2,i) = ep30
        eminx(3,i) = ep30
        eminx(4,i) = -ep30
        eminx(5,i) = -ep30
        eminx(6,i) = -ep30      
      ENDDO
C
      lft16=1
      llt16=0
      lft20=1
      llt20=0
      DO i=nmef,nmel
          i20=nelem(i)-numels8-numels10
          i16=i20-numels20
          IF(i16>=1.AND.i16<=numels16)THEN
            llt16=llt16+1
            index16(llt16)=i
            IF(llt16==mvsiz-1)THEN
              CALL i16box(
     1            lft16,llt16  ,nelem,eminx,nmef ,nmel  ,
     2            x    ,v      ,a    ,ixs  ,ixs16,size_t,
     3            xmsr ,index16,xsav  )
              llt16=0
            ENDIF
          ELSEIF(i20>=1.AND.i20<=numels20)THEN
            llt20=llt20+1
            index20(llt20)=i
            IF(llt20==mvsiz-1)THEN
              CALL i20box(
     1            lft20,llt20  ,nelem,eminx,nmef ,nmel  ,
     2            x    ,v      ,a    ,ixs  ,ixs20,size_t,
     3            xmsr ,index20,xsav  )
              llt20=0
            ENDIF
          ENDIF
      END DO
      IF(llt16>0)CALL i16box(
     1            lft16,llt16  ,nelem,eminx,nmef ,nmel  ,
     2            x    ,v      ,a    ,ixs  ,ixs16,size_t,
     3            xmsr ,index16,xsav  )
      IF(llt20>0)CALL i20box(
     1            lft20,llt20  ,nelem,eminx,nmef ,nmel  ,
     2            x    ,v      ,a    ,ixs  ,ixs20,size_t,
     3            xmsr ,index20,xsav  )
C
#include "lockon.inc"
      xmsr_g(1)=max(xmsr_g(1),xmsr(1))
      xmsr_g(2)=max(xmsr_g(2),xmsr(2))
      xmsr_g(3)=max(xmsr_g(3),xmsr(3))
      xmsr_g(4)=min(xmsr_g(4),xmsr(4))
      xmsr_g(5)=min(xmsr_g(5),xmsr(5))
      xmsr_g(6)=min(xmsr_g(6),xmsr(6))
      size_t_g = size_t_g + size_t
#include "lockoff.inc"
C
      RETURN