Télécharger transb.eso

Retour à la liste

Numérotation des lignes :

  1. C TRANSB SOURCE PV 06/04/11 21:16:17 5396
  2. C SOUS PROGRAMME DE VERIFICATION QU'UN SEGMENT POTENTIEL (LL1 LL2)
  3. C EST DU BON COTE DU CONTOUR (LL3 LL1 LL4)
  4. C
  5. SUBROUTINE TRANSB(X,IRECL,LL1,LL2,LL3,LL4)
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8 (A-H,O-Z)
  8. DIMENSION X(3,1)
  9. IRECL=1
  10. X41=X(1,LL1)-X(1,LL4)
  11. Y41=X(2,LL1)-X(2,LL4)
  12. X31=X(1,LL1)-X(1,LL3)
  13. Y31=X(2,LL1)-X(2,LL3)
  14. UV=X41*Y31-X31*Y41
  15. X21=X(1,LL1)-X(1,LL2)
  16. Y21=X(2,LL1)-X(2,LL2)
  17. VX=X41*Y21-Y41*X21
  18. UX=X21*Y31-Y21*X31
  19. IF (UV.GE.0.D0.AND.VX.GT.0.D0.AND.UX.GT.0.D0) IRECL=0
  20. IF (UV.LT.0.D0.AND.(VX.GT.0.D0.OR.UX.GT.0.D0)) IRECL=0
  21. RETURN
  22. END
  23.  
  24.  
  25.  

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