Télécharger la plaquette Cast3M


fr en

Télécharger sens.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  sens.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. *       FICHIER DGIBI POUR TESTER L'OPERATEURR SENS
  6. *       option b)
  7. *
  8. *
  9. *       Arnaud de Gayffier
  10. *
  11. opti dime 2 elem seg2 ;
  12. densite 1. ;
  13. *       on cree deux cercles concentriques orientés 
  14. *       dans des sens différents.
  15. pc0 = 0. 0. ;
  16. p1 = 10. 0. ;
  17. p2 = -10. 0.01 ;
  18. c1 = c 10 p1 pc0 p2;
  19. pc1 = 0. 0.02 ;
  20. c2 = c 10 p2 pc1 p1;
  21. *
  22. p3 = 5. 0. ;
  23. p4 = -5. 0.01 ;
  24. c3 = c 10 p4 pc0 p3 ;
  25. c4 = c 10 p3 pc1 p4 ;
  26. *
  27. *       on maille l'interieur 
  28. su1 = surf (c1 et c2 et c3 et c4 ) 'PLAN';
  29. *       on extrait le contour
  30. cont1 = cont su1 ;
  31. *       on extrait les composantes connexes du contours
  32. tab1 = ccon cont1 ;
  33. *       on extrait les orientations du contours
  34. tab2 = 'SENS' tab1 ;
  35. *
  36. *       controle de validite
  37. si ( (nbelem ( diff tab1.1 c1 )) NEG  31) ;
  38. * tab1.1 est le cercle exterieur
  39.   si ( (tab2.1) 'NEG' 1 ) ;
  40.    erreur 5 ;
  41.   finsi ;
  42.   mess 'OK' ;
  43. finsi ;
  44. *
  45. si ( (nbelem ( diff tab1.1 c3 )) NEG  30) ;
  46. * tab1.1 est le cercle interieur
  47.   si ( (tab2.1) 'NEG' -1 ) ;
  48.    erreur 5 ;
  49.   finsi ;
  50.   mess 'OK' ;
  51. finsi ;
  52. *****
  53. si ( (nbelem ( diff tab1.2 c1 )) NEG  30) ;
  54. * tab1.2 est le cercle exterieur
  55.   si ( (tab2.2) 'NEG' 1 ) ;
  56.    erreur 5 ;
  57.   finsi ;
  58.   mess 'OK' ;
  59. finsi ;
  60. *
  61. si ( (nbelem ( diff tab1.2  c3 )) NEG  30) ;
  62. * tab1.2 est le cercle interieur
  63.   si ( (tab2.2) 'NEG' -1 ) ;
  64.    erreur 5 ;
  65.   finsi ;
  66.   mess 'OK' ;
  67. finsi ;
  68. fin ;
  69.  
  70.  
  71.  
  72.  

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