Numérotation des lignes :

apat3
C APAT3     SOURCE    PV        22/04/15    13:20:05     11344                SUBROUTINE APAT3(X,XELTRI,IELEM,IARTRI,INOTRI,TLI)CC**********************************************************************CC   OBJET  : CE SOUS-PROGRAMME RECHERCHE L'APPARTENANCE D'UN POINTC   -----    A UN TRI3CC   ARGUMENTS:C   ---------CC   ENTREE  :  IDIM   = DIMENSION DE L'ESPACEC              X      = TABLEAU DES COORDONNEES DU POINTC              XELTRI = TABLEAU DES COORDONNEES DU TRI3CC   SORTIE  :  IELEM  = 0 SI LE POINT N'APPARTIENT PAS A L'ELEMENTC                     = 1 SI LE POINT APPARTIENT A L'ELEMENTC              IARTRI = 0 SI LE POINT N'APPARTIENT PAS A UNE ARETEC                     = N NUMERO DE L'ARETE SINONC              INOTRI = 0 SI LE POINT N'APPARTIENT PAS A UN NOEUDC                     = N NUMERO DU NOEUD SINONCC***********************************************************************C      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8 (A-H,O-Z)      DIMENSION X(3),XELTRI(2,1),IART(3)      DIMENSION TLI(3)C      DATA IART/2,3,1/C                  IART numero local de l arete opose au noeud IC      IDIM=2      IELEM=0      IARTRI=0      INOTRI=0      NOEU=3C     EPS=1.D-5      EPS=1.D-10      UN=1.D0+EPSCC***  CALCUL DES COORDONNEES BARYCENTRIQUESC      CALL COBAR3(X,XELTRI,TLI)CC      IF(TLI(1).LT.-EPS.OR.TLI(1).GT.UN)GO TO 9999      IF(TLI(2).LT.-EPS.OR.TLI(2).GT.UN)GO TO 9999      IF(TLI(3).LT.-EPS.OR.TLI(3).GT.UN)GO TO 9999C           LE POINT APPARTIENT A L ELEMENT      IELEM=1 C***  RECHERCHE DE L'APPARTENANCE A UNE ARETEC      DO 50 I=1,NOEU      IARTRI=IART(I)      VAL=ABS(TLI(I))      IF(VAL.LE.EPS) GO TO 60   50 CONTINUE      IARTRI=0      GO TO 9999CC***  RECHERCHE DE L'APPARTENANCE A UN NOEUD DE L'ARETEC   60 CONTINUE      DO 70 I=1,3      INOTRI=I      VAL=ABS(TLI(INOTRI)-1.D0)      IF(VAL.LE.EPS) GO TO 9999   70 CONTINUE      INOTRI=0 9999 RETURN      END   

© Cast3M 2003 - Tous droits réservés.
Mentions légales