Télécharger tria.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : tria.dgib
  2. ************************************************************************
  3. ************************************************************************
  4. ************************************************************************
  5. * tria.dgibi : CAS TEST de l'operateur TRIA *
  6. * Triangulation de Delaunay d'un maillage de points (type POI1). *
  7. * On teste plusieurs cas : en 1D, 2D et 3D, avec, a chaque fois : *
  8. * - un nuage de points distribues aleatoirement ; *
  9. * - un reseau regulier (test de robustesse de l'algorithme). *
  10. ************************************************************************
  11.  
  12. ** INDICATEUR DE TRACE
  13. ITRAC = FAUX ;
  14. OPTI 'ECHO' 1 'TRAC' 'X' ;
  15.  
  16.  
  17. ************************************************************************
  18. * EN DIMENSION 1 *
  19. ************************************************************************
  20.  
  21. OPTI 'DIME' 1 'ELEM' 'SEG2' ;
  22.  
  23. ** POINTS ALEATOIRES
  24. N1 = 10 ;
  25. LX = BRUI 'BLAN' 'UNIF' 0.5 0.5 N1 ;
  26. REPE B1 N1 ;
  27. SI (&B1 EGA 1) ;
  28. MAIL1 = POIN (EXTR LX &B1) ;
  29. SINO ;
  30. MAIL1 = MAIL1 ET (POIN (EXTR LX &B1)) ;
  31. FINS ;
  32. FIN B1 ;
  33. MAIL1 = MAIL1 COUL 'ROUG' ;
  34. MAIL2 = TRIA MAIL1 ;
  35. SI ITRAC ;
  36. TRAC (MAIL2 ET MAIL1) 'TITR' 'Points aleatoires (1D)' ;
  37. FINS ;
  38.  
  39.  
  40.  
  41.  
  42. ************************************************************************
  43. * EN DIMENSION 2 *
  44. ************************************************************************
  45.  
  46. OPTI 'DIME' 2 'ELEM' 'TRI3' ;
  47.  
  48. ** POINTS ALEATOIRES DANS LE PLAN
  49. N1 = 8 ;
  50. LX = BRUI 'BLAN' 'UNIF' 0.5 0.5 N1 ;
  51. LY = BRUI 'BLAN' 'UNIF' 0.5 0.5 N1 ;
  52. REPE B1 N1 ;
  53. SI (&B1 EGA 1) ;
  54. MAIL1 = ((EXTR LX &B1) (EXTR LY &B1)) ;
  55. SINO ;
  56. MAIL1 = MAIL1 ET ((EXTR LX &B1) (EXTR LY &B1)) ;
  57. FINS ;
  58. FIN B1 ;
  59. MAIL1 = MAIL1 COUL 'ROUG' ;
  60. MAIL2 = TRIA MAIL1 ;
  61. MAIL21 = TRIA MAIL1 'CONV' ;
  62. MAIL22 = TRIA MAIL1 0.2 ;
  63. MAIL23 = TRIA MAIL1 0.2 'CONV' ;
  64. MAIL3 = CONT MAIL2 ;
  65. MES2 = MESU MAIL2 ;
  66. MES3 = MESU MAIL3 'SURF' ;
  67. SI (('ABS' (MES2 - MES3)) '>' 1.E-13) ;
  68. MESS ' ****** Erreur dans TRIA :'
  69. MESS ' Les mesures des deux triangulations ne sont pas identiques' ;
  70. ERRE 5 ;
  71. FINS ;
  72. SI ITRAC ;
  73. TRAC (MAIL2 ET MAIL1) 'TITR' 'Points aleatoires (2D)' ;
  74. TRAC (MAIL21 ET MAIL1) 'TITR' 'Points aleatoires (2D) + CONV' ;
  75. TRAC (MAIL22 ET MAIL1) 'TITR' 'Points aleatoires (2D) + 0.2' ;
  76. TRAC (MAIL23 ET MAIL1) 'TITR' 'Points aleatoires (2D) + 0.2 + CONV' ;
  77. FINS ;
  78.  
  79.  
  80.  
  81. ** POINTS REPARTIS SUR UN RESEAU CARRE REGULIER
  82. N1 = 3 ;
  83. L1 = D N1 (0. 0.) (1. 0.) ;
  84. MAIL1 = (CHAN 'POI1' (L1 TRAN N1 (0. 1.))) COUL 'ROUG' ;
  85. MAIL2 = TRIA MAIL1 ;
  86. MAIL21 = TRIA MAIL1 'CONV' ;
  87. MAIL22 = TRIA MAIL1 0.2 ;
  88. MAIL23 = TRIA MAIL1 0.2 'CONV' ;
  89. MAIL3 = CONT MAIL2 ;
  90. MES2 = MESU MAIL2 ;
  91. MES3 = MESU MAIL3 'SURF' ;
  92. SI (('ABS' (MES2 - MES3)) '>' 1.E-13) ;
  93. MESS ' ****** Erreur dans TRIA :'
  94. MESS ' Les mesures des deux triangulations ne sont pas identiques' ;
  95. ERRE 5 ;
  96. FINS ;
  97. SI ITRAC ;
  98. TRAC (MAIL2 ET MAIL1) 'TITR' 'Points reguliers (2D)' ;
  99. TRAC (MAIL21 ET MAIL1) 'TITR' 'Points reguliers (2D) + CONV' ;
  100. TRAC (MAIL22 ET MAIL1) 'TITR' 'Points reguliers (2D) + 0.2' ;
  101. TRAC (MAIL23 ET MAIL1) 'TITR' 'Points reguliers (2D) + 0.2 + CONV' ;
  102. FINS ;
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110. ** POINTS ALEATOIRES MAIS ALIGNES
  111. N1 = 5 ;
  112. LX = BRUI 'BLAN' 'UNIF' 0.5 0.5 N1 ;
  113. REPE B1 (DIME LX) ;
  114. SI (&B1 EGA 1) ;
  115. MAIL1 = ((EXTR LX &B1) 0.) ;
  116. SINO ;
  117. MAIL1 = MAIL1 ET ((EXTR LX &B1) 0.) ;
  118. FINS ;
  119. FIN B1 ;
  120. DEPL MAIL1 'PROJ' 'CYLI' (0. 1.) 'DROI' (0. 2.) (2. 3.) ;
  121. MAIL1 = MAIL1 COUL 'ROUG' ;
  122. MAIL2 = TRIA MAIL1 ;
  123. MAIL21 = TRIA MAIL1 'CONV' ;
  124. MAIL22 = TRIA MAIL1 0.05 ;
  125. MAIL23 = TRIA MAIL1 0.05 'CONV' ;
  126. SI ITRAC ;
  127. TRAC (MAIL2 ET MAIL1) 'ELEM' 'TITR' 'Points aleatoires alignes (2D)' ;
  128. TRAC (MAIL21 ET MAIL1) 'ELEM' 'TITR'
  129. 'Points aleatoires alignes (2D) + CONV' ;
  130. TRAC (MAIL22 ET MAIL1) 'ELEM' 'TITR'
  131. 'Points aleatoires alignes (2D) + 0.05' ;
  132. TRAC (MAIL23 ET MAIL1) 'ELEM' 'TITR'
  133. 'Points aleatoires alignes (2D) + 0.05 + CONV' ;
  134. FINS ;
  135.  
  136.  
  137.  
  138.  
  139. ************************************************************************
  140. * EN DIMENSION 3 *
  141. ************************************************************************
  142. OPTI 'DIME' 3 'ELEM' 'TET4' ;
  143.  
  144. R1 = @REPERE (-0.5 -0.5 -0.5) (PROG 1. 1. 1.) VRAI 'VERT' ;
  145.  
  146. ** POINTS ALEATOIRES DANS L'ESPACE
  147. N1 = 30 ;
  148. LX = BRUI 'BLAN' 'UNIF' 0.5 0.5 N1 ;
  149. LY = BRUI 'BLAN' 'UNIF' 0.5 0.5 N1 ;
  150. LZ = BRUI 'BLAN' 'UNIF' 0.5 0.5 N1 ;
  151. REPE B1 N1 ;
  152. SI (&B1 EGA 1) ;
  153. MAIL1 = ((EXTR LX &B1) (EXTR LY &B1) (EXTR LZ &B1)) ;
  154. SINO ;
  155. MAIL1 = MAIL1 ET ((EXTR LX &B1) (EXTR LY &B1) (EXTR LZ &B1)) ;
  156. FINS ;
  157. FIN B1 ;
  158. MAIL1 = MAIL1 COUL 'ROUG' ;
  159. MAIL2 = TRIA MAIL1 ;
  160. MAIL21 = TRIA MAIL1 'CONV' ;
  161. MAIL22 = TRIA MAIL1 0.4 ;
  162. MAIL23 = TRIA MAIL1 0.4 'CONV' ;
  163. MAIL3 = ENVE MAIL2 ;
  164. MES2 = MESU MAIL2 ;
  165. MES3 = MESU MAIL3 'VOLU' ;
  166. SI (('ABS' (MES2 - MES3)) '>' 1.E-13) ;
  167. MESS ' **** Erreur dans TRIA:' ;
  168. MESS ' Les mesures des deux triangulations ne sont pas identiques' ;
  169. ERRE 5 ;
  170. FINS ;
  171. SI ITRAC ;
  172. TRAC (MAIL2 ET MAIL1 ET R1) 'TITR' 'Points aleatoires (3D)' ;
  173. TRAC (MAIL21 ET MAIL1 ET R1) 'TITR' 'Points aleatoires (3D) + CONV' ;
  174. TRAC (MAIL22 ET MAIL1 ET R1) 'TITR' 'Points aleatoires (3D) + 0.4' ;
  175. TRAC (MAIL23 ET MAIL1 ET R1) 'TITR'
  176. 'Points aleatoires (3D) + 0.4 + CONV' ;
  177. FINS ;
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184. ** POINTS REPARTIS SUR UN RESEAU CUBIQUE REGULIER
  185. N1 = 3 ;
  186. L1 = D N1 (0. 0. 0.) (1. 0. 0.) ;
  187. S1 = L1 TRAN N1 (0. 1. 0.) ;
  188. MAIL1 = (CHAN 'POI1' (S1 VOLU 'TRAN' N1 (0. 0. 1.))) COUL 'ROUG' ;
  189. MAIL2 = TRIA MAIL1 ;
  190. MAIL21 = TRIA MAIL1 'CONV' ;
  191. MAIL22 = TRIA MAIL1 0.2 ;
  192. MAIL23 = TRIA MAIL1 0.2 'CONV' ;
  193. MAIL3 = ENVE MAIL2 ;
  194. MES2 = MESU MAIL2 ;
  195. MES3 = MESU MAIL3 'VOLU' ;
  196. SI (('ABS' (MES2 - MES3)) '>' 1.E-13) ;
  197. MESS ' ****** Erreur dans TRIA :'
  198. MESS ' Les mesures des deux triangulations ne sont pas identiques' ;
  199. ERRE 5 ;
  200. FINS ;
  201. SI ITRAC ;
  202. TRAC (MAIL2 ET MAIL1 ET R1) 'TITR' 'Points reguliers (3D)' ;
  203. TRAC (MAIL21 ET MAIL1 ET R1) 'TITR' 'Points reguliers (3D) + CONV' ;
  204. TRAC (MAIL22 ET MAIL1 ET R1) 'TITR' 'Points reguliers (3D) + 0.2' ;
  205. TRAC (MAIL23 ET MAIL1 ET R1) 'TITR'
  206. 'Points reguliers (3D) + 0.2 + CONV' ;
  207. FINS ;
  208.  
  209.  
  210.  
  211.  
  212.  
  213. ** POINTS ALEATOIRES MAIS COPLANAIRES
  214. N1 = 20 ;
  215. LX = BRUI 'BLAN' 'UNIF' 0.5 0.5 N1 ;
  216. LY = BRUI 'BLAN' 'UNIF' 0.5 0.5 N1 ;
  217. REPE B1 N1 ;
  218. SI (&B1 EGA 1) ;
  219. MAIL1 = ((EXTR LX &B1) (EXTR LY &B1) 0.) ;
  220. SINO ;
  221. MAIL1 = MAIL1 ET ((EXTR LX &B1) (EXTR LY &B1) 0.) ;
  222. FINS ;
  223. FIN B1 ;
  224. DEPL MAIL1 'PROJ' 'CYLI' (0. 0. 1.) 'PLAN'
  225. (0. 0. 1.) (0. 1. 1.5) (1. 0. 0.) ;
  226. MAIL1 = MAIL1 COUL 'ROUG' ;
  227. MAIL2 = TRIA MAIL1 ;
  228. MAIL21 = TRIA MAIL1 'CONV' ;
  229. MAIL22 = TRIA MAIL1 0.2 ;
  230. MAIL23 = TRIA MAIL1 0.2 'CONV' ;
  231. SI ITRAC ;
  232. TRAC (MAIL2 ET MAIL1 ET R1) 'TITR'
  233. 'Points aleatoires coplanaires (3D)' ;
  234. TRAC (MAIL21 ET MAIL1 ET R1) 'TITR'
  235. 'Points aleatoires coplanaires (3D) + CONV' ;
  236. TRAC (MAIL22 ET MAIL1 ET R1) 'TITR'
  237. 'Points aleatoires coplanaires (3D) + 0.2' ;
  238. TRAC (MAIL23 ET MAIL1 ET R1) 'TITR'
  239. 'Points aleatoires coplanaires (3D) + 0.2 + CONV' ;
  240. FINS ;
  241.  
  242.  
  243. ** POINTS REPARTIS SUR UN RESEAU CARRE COPLANAIRES
  244. N1 = 4 ;
  245. L1 = D N1 (0. 0. 0.) (1. 0. 0.) ;
  246. MAIL1 = L1 TRAN N1 (0. 1. 0.) ;
  247. DEPL MAIL1 'PROJ' 'CYLI' (0. 0. 1.) 'PLAN'
  248. (0. 0. 0.3) (0. 1. 0.5) (1. 0. 0.) ;
  249. MAIL1 = (CHAN 'POI1' MAIL1) COUL 'ROUG' ;
  250. MAIL2 = TRIA MAIL1 ;
  251. MAIL21 = TRIA MAIL1 'CONV' ;
  252. MAIL22 = TRIA MAIL1 0.1 ;
  253. MAIL23 = TRIA MAIL1 0.1 'CONV' ;
  254. SI ITRAC ;
  255. TRAC (MAIL2 ET MAIL1 ET R1) 'TITR'
  256. 'Points reguliers coplanaires (3D)' ;
  257. TRAC (MAIL21 ET MAIL1 ET R1) 'TITR'
  258. 'Points reguliers coplanaires (3D) + CONV' ;
  259. TRAC (MAIL22 ET MAIL1 ET R1) 'TITR'
  260. 'Points reguliers coplanaires (3D) + 0.1' ;
  261. TRAC (MAIL23 ET MAIL1 ET R1) 'TITR'
  262. 'Points reguliers coplanaires (3D) + 0.1 + CONV' ;
  263. FINS ;
  264.  
  265.  
  266. ** POINTS ALEATOIRES MAIS ALIGNES
  267. N1 = 4 ;
  268. LX = BRUI 'BLAN' 'UNIF' 0.5 0.5 N1 ;
  269. REPE B1 N1 ;
  270. SI (&B1 EGA 1) ;
  271. MAIL1 = ((EXTR LX &B1) 0. 0.) ;
  272. SINO ;
  273. MAIL1 = MAIL1 ET ((EXTR LX &B1) 0. 0.) ;
  274. FINS ;
  275. FIN B1 ;
  276. DEPL MAIL1 'PROJ' 'CYLI' (0. 1. 0.) 'PLAN'
  277. (0. 0. 0.) (0. 0. 1.) (1. 2. 0.) ;
  278. DEPL MAIL1 'PROJ' 'CYLI' (0. 0. 1.) 'PLAN'
  279. (0. 0. 0.5) (1. 0. 1.) (0. 1. 1.) ;
  280. MAIL1 = MAIL1 COUL 'ROUG' ;
  281. MAIL2 = TRIA MAIL1 ;
  282. MAIL21 = TRIA MAIL1 'CONV' ;
  283. MAIL22 = TRIA MAIL1 0.2 ;
  284. MAIL23 = TRIA MAIL1 0.2 'CONV' ;
  285. SI ITRAC ;
  286. TRAC (MAIL2 ET MAIL1 ET R1) 'ELEM' 'TITR'
  287. 'Points aleatoires coplanaires (3D)' ;
  288. TRAC (MAIL21 ET MAIL1 ET R1) 'ELEM' 'TITR'
  289. 'Points aleatoires coplanaires (3D) + CONV' ;
  290. TRAC (MAIL22 ET MAIL1 ET R1) 'ELEM' 'TITR'
  291. 'Points aleatoires coplanaires (3D) + 0.2' ;
  292. TRAC (MAIL23 ET MAIL1 ET R1) 'ELEM' 'TITR'
  293. 'Points aleatoires coplanaires (3D) + 0.2 + CONV' ;
  294. FINS ;
  295.  
  296.  
  297.  
  298.  
  299.  
  300.  
  301.  
  302.  
  303. ************************************************************************
  304. * MAILLAGE A PARTIR D'UN CONTOUR/SURFACE (EN DIMENSION 2) *
  305. * ENVELOPPE/VOLUME (EN DIMENSION 3) *
  306. ************************************************************************
  307. OPTI 'DIME' 2 'ELEM' 'TRI3' ;
  308.  
  309. ** TRIANGULATION A PARTIR D'UN CONTOUR FERME EN U
  310. X1 = 0.05 ;
  311. N1 = 15 ;
  312. N2 = 4 ;
  313. L1 = D N1 (0. 0.) (1. 0.) ;
  314. L2 = D N1 L1 (1. 1.) ;
  315. L3 = D N2 L2 ((0.5+X1) 1.) ;
  316. L4 = D N1 L3 ((0.5+X1) 0.1) ;
  317. L5 = D N2 L4 ((0.5-X1) 0.1) ;
  318. L6 = D N1 L5 ((0.5-X1) 1.) ;
  319. L7 = D N2 L6 (0. 1.) ;
  320. L8 = D N1 L7 (0. 0.) ;
  321. ELIM L8 1.E-9 ;
  322. MAIL1 = TRIA L8 0.2 ;
  323. SI ITRAC ;
  324. TRAC MAIL1 'TITR' 'Remplissage a partir du contour' ;
  325. FINS ;
  326. S8 = SURF L8 ;
  327. MAIL1 = TRIA S8 0.2 ;
  328. SI ITRAC ;
  329. TRAC MAIL1 'TITR' 'Re-maillage a partir d une surface' ;
  330. FINS ;
  331.  
  332. OPTI 'DIME' 3 'ELEM' 'TET4' ;
  333. V8 = CHAN 'TET4' (MAIL1 VOLU 'TRAN' 4 (0. 0. 0.5)) ;
  334. ENV8 = ENVE V8 ;
  335. MAIL1 = TRIA ENV8 0.2 ;
  336. SI ITRAC ;
  337. TRAC 'FACE' MAIL1 'TITR' 'Maillage a partir de l enveloppe' ;
  338. FINS ;
  339. MAIL1 = TRIA V8 0.2 ;
  340. SI ITRAC ;
  341. TRAC 'FACE' MAIL1 'TITR' 'Re-maillage a partir d un volume' ;
  342. FINS ;
  343.  
  344.  
  345. FIN ;
  346.  
  347.  
  348.  
  349.  
  350.  
  351.  
  352.  
  353.  

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