rtadet
C RTADET SOURCE CB215821 17/11/30 21:17:12 9639 > NBE,COORD,SPH, > IADET,NBADET,NADMAX,ZERO,iarr) C ********************************************************************** C OBJET : RECHERCHE DES TRIANGLES A DETRUIRE (IE "NON-DELAUNAY") C A L'AJOUT D'UN POINT C EN ENTREE : C XYZPT : COORDONNEES DU POINT AJOUTE C IDIMC : DIMENSION DE L'ESPACE C ITRNOE,NBNMAX,ITRTRI,NBCMAX,NBE,COORD : LA TRIANGULATION C SPH : LES SPHERES CIRCONSCRITES AUX TRIANGLES C NBADETMNAX : TAILLE DU TABLEAU IADET C ZERO : PRECISION DU TEST "POINT DANS SPHERE" C EN SORTIE : C IADET : TABLEAU DES TRIANGLES "NON-DELAUNAY" C NBADET : NOMBRE DE TRIANGLES " " " " C ********************************************************************** IMPLICIT INTEGER(I-N) REAL*8 XYZPT(*) INTEGER IDIMC INTEGER ITRNOE(*),NBNMAX,ITRTRI(*),NBCMAX INTEGER NBE,IADET(*),NBADET,NADMAX,iarr C C C --- POUR LE DEBUG --- C C COMMON /DEBUG/ ITRACE, ITEST, ierrOR, IMESS C INTEGER ITRACE, ITEST, ierrOR C CHARACTER*256 IMESS C --------------------------------------------------------------------- C --- POUR LES STATS --- C C COMMON /STATS/ ICARD(100) C INTEGER ICARD C --------------------------------------------------------------------- C --- VARIABLES INTERNES --- INTEGER I,NT,IPTDS,IPTDSC,IPTDS2,ITRACE REAL*8 SPHC(4) INTEGER SPPOIN EXTERNAL SPPOIN C ITRACE = 0 DO i = 1, 4 SPHC(i) = 0.D0 ENDDO iarr = 0 NBADET = 0 DO 30 I=1,NBE C ----- ON PREND LE DERNIER NOEUD --- NT =ITRNOE((I-1)*NBNMAX+1+IDIMC) IF ( NT.EQ. 0 ) GO TO 30 C IF( IPTDS.EQ.1 )THEN IADET(1) = I NBADET = 1 > COORD,SPH,IADET,NBADET,NADMAX,ZERO,iarr) RETURN ENDIF C 30 CONTINUE RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales