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