Télécharger la plaquette Cast3M


fr en

Télécharger arete3.eso

Retour à la liste

Numérotation des lignes : 

arete3

  1. C ARETE3    SOURCE    PV        09/01/08    21:15:16     6254
  2.       SUBROUTINE ARETE3 (TS,ISEG,TN,INOR,CANGLE,NBSEG)
  3.  
  4. *       APPELE PAR ARETE
  5. *          ELIMINE LES SEG2 PAR COMPARAISON D'ANGLE
  6. *
  7. *-------------------------------------------------------------------
  8. *
  9. *   VARIABLES PRINCIPALES :
  10. *   ---------------------
  11. *
  12. *      TSEG : SEGMENT 2 DIMENSIONS  CONTENANT
  13. *
  14. *     | A TRACER (= 0  NON,   =1 OUI)    |
  15. *     | COULEUR                          |
  16. *     | N› NOEUD MIN                     |
  17. *     | N› NOEUD MAX                     |
  18. *     | ENTREES DANS TSEG                |
  19. *
  20. *      TNOR : SEGMENT 2 DIMENSIONS CONTENANT LES NORMALES A COMPARER
  21. *             POUR CHAQUE FACE ELEMENTAIRE
  22. *
  23. *      ISEG : TAILLE REELLE DU SEGMENT TSEG
  24. *      INOR : TAILLE REELLE DU SEGMENT TNOR
  25. *
  26. *--------------------------------------------------------------------
  27. *      AUTEUR : J.BRUN (JUIN 90)
  28.  
  29.       IMPLICIT INTEGER(I-N)
  30.       IMPLICIT REAL*8 (A-H,O-Y)
  31.       IMPLICIT LOGICAL (Z)
  32.       SEGMENT TN
  33.          REAL*8 TNOR (N1,3)
  34.       ENDSEGMENT
  35.       SEGMENT TS
  36.          INTEGER TSEG (M1,M2)
  37.       ENDSEGMENT
  38.  
  39.       NBSEG=0
  40.       DO 1 I=1,ISEG
  41.  
  42. *     DEMANDE TRACER SEGMENT NON PARTAGE ENTRE 2 ELEMENTS
  43.            IF (TSEG(I,6).EQ.0) THEN
  44.                 NBSEG=NBSEG+1
  45.                 TSEG(I,1)=1
  46.            ELSE
  47.                 J=5
  48.                 ZFINI=.FALSE.
  49. *     RECHERCHE LE NB DE NORMALE ASSOCIEE A UN SEGMENT
  50. 2               J=J+1
  51.                 ZFINI=(TSEG(I,J+1).EQ.0)
  52.                 IF (.NOT.((J.EQ.TSEG(/2)).OR.(ZFINI))) GOTO 2
  53. *     COMPARE LES NORMALES 2 A 2 ET S'ARRETE S'IL DOIT TRACER
  54.  
  55.  
  56.                 DO 3 K=5,J
  57.  
  58.                      DO 4 KK=K+1,J
  59.  
  60.                 CBIS=TNOR(TSEG(I,K),1)*TNOR(TSEG(I,KK),1)
  61.      .              +TNOR(TSEG(I,K),2)*TNOR(TSEG(I,KK),2)
  62.      .              +TNOR(TSEG(I,K),3)*TNOR(TSEG(I,KK),3)
  63.                 IF (CBIS.LT.0D0) CBIS=-CBIS
  64.                     IF (CBIS.LT.CANGLE) THEN
  65.                               TSEG(I,1)=1
  66.                               NBSEG=NBSEG+1
  67.                               GOTO 1
  68.                      ENDIF
  69. 4                    CONTINUE
  70. 3               CONTINUE
  71.            ENDIF
  72. 1     CONTINUE
  73.       RETURN
  74.       END
  75.  
  76.  
  77.  
  78.  

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