* fichier : tria.dgibi ************************************************************************ ************************************************************************ ************************************************************************ * tria.dgibi : CAS TEST de l'operateur TRIA * * Triangulation de Delaunay d'un maillage de points (type POI1). * * On teste plusieurs cas : en 1D, 2D et 3D, avec, a chaque fois : * * - un nuage de points distribues aleatoirement ; * * - un reseau regulier (test de robustesse de l'algorithme). * * - ajout de test pour les quadratiques (SGounand 2024/10/07) * ************************************************************************ ** INDICATEUR DE TRACE ITRAC = FAUX ; ************************************************************************ * EN DIMENSION 1 * ************************************************************************ ** POINTS ALEATOIRES N1 = 10 ; REPE B1 N1 ; SI (&B1 EGA 1) ; SINO ; FINS ; FIN B1 ; SI ITRAC ; FINS ; ************************************************************************ * EN DIMENSION 2 * ************************************************************************ ** POINTS ALEATOIRES DANS LE PLAN N1 = 8 ; REPE B1 N1 ; SI (&B1 EGA 1) ; SINO ; FINS ; FIN B1 ; SI (('ABS' (MES2 - MES3)) '>' 1.E-13) ; MESS ' Les mesures des deux triangulations ne sont pas identiques' ; FINS ; SI ITRAC ; FINS ; ** POINTS REPARTIS SUR UN RESEAU CARRE REGULIER N1 = 3 ; L1 = D N1 (0. 0.) (1. 0.) ; SI (('ABS' (MES2 - MES3)) '>' 1.E-13) ; MESS ' Les mesures des deux triangulations ne sont pas identiques' ; FINS ; SI ITRAC ; FINS ; ** POINTS ALEATOIRES MAIS ALIGNES N1 = 5 ; SI (&B1 EGA 1) ; SINO ; FINS ; FIN B1 ; SI ITRAC ; 'Points aleatoires alignes (2D) + 0.05' ; FINS ; ************************************************************************ * EN DIMENSION 3 * ************************************************************************ ** POINTS ALEATOIRES DANS L'ESPACE N1 = 30 ; REPE B1 N1 ; SI (&B1 EGA 1) ; SINO ; FINS ; FIN B1 ; SI (('ABS' (MES2 - MES3)) '>' 1.E-13) ; MESS ' Les mesures des deux triangulations ne sont pas identiques' ; FINS ; SI ITRAC ; FINS ; ** POINTS REPARTIS SUR UN RESEAU CUBIQUE REGULIER N1 = 3 ; L1 = D N1 (0. 0. 0.) (1. 0. 0.) ; SI (('ABS' (MES2 - MES3)) '>' 1.E-13) ; MESS ' Les mesures des deux triangulations ne sont pas identiques' ; FINS ; SI ITRAC ; FINS ; ** POINTS ALEATOIRES MAIS COPLANAIRES N1 = 20 ; REPE B1 N1 ; SI (&B1 EGA 1) ; SINO ; FINS ; FIN B1 ; (0. 0. 1.) (0. 1. 1.5) (1. 0. 0.) ; SI ITRAC ; 'Points aleatoires coplanaires (3D)' ; 'Points aleatoires coplanaires (3D) + 0.2' ; FINS ; ** POINTS REPARTIS SUR UN RESEAU CARRE COPLANAIRES N1 = 4 ; L1 = D N1 (0. 0. 0.) (1. 0. 0.) ; (0. 0. 0.3) (0. 1. 0.5) (1. 0. 0.) ; SI ITRAC ; 'Points reguliers coplanaires (3D)' ; 'Points reguliers coplanaires (3D) + 0.1' ; FINS ; ** POINTS ALEATOIRES MAIS ALIGNES N1 = 4 ; REPE B1 N1 ; SI (&B1 EGA 1) ; SINO ; FINS ; FIN B1 ; (0. 0. 0.) (0. 0. 1.) (1. 2. 0.) ; (0. 0. 0.5) (1. 0. 1.) (0. 1. 1.) ; SI ITRAC ; 'Points aleatoires coplanaires (3D)' ; 'Points aleatoires coplanaires (3D) + 0.2' ; FINS ; ************************************************************************ * MAILLAGE A PARTIR D'UN CONTOUR/SURFACE (EN DIMENSION 2) * * ENVELOPPE/VOLUME (EN DIMENSION 3) * ************************************************************************ ** TRIANGULATION A PARTIR D'UN CONTOUR FERME EN U X1 = 0.05 ; N1 = 15 ; N2 = 4 ; L1 = D N1 (0. 0.) (1. 0.) ; L2 = D N1 L1 (1. 1.) ; L3 = D N2 L2 ((0.5+X1) 1.) ; L4 = D N1 L3 ((0.5+X1) 0.1) ; L5 = D N2 L4 ((0.5-X1) 0.1) ; L6 = D N1 L5 ((0.5-X1) 1.) ; L7 = D N2 L6 (0. 1.) ; L8 = D N1 L7 (0. 0.) ; SI ITRAC ; FINS ; SI ITRAC ; FINS ; ** TRIANGULATION A PARTIR D'UN CONTOUR FERME EN U X1 = 0.05 ; N1 = 15 ; N2 = 4 ; L1 = D N1 (0. 0.) (1. 0.) ; L2 = D N1 L1 (1. 1.) ; L3 = D N2 L2 ((0.5+X1) 1.) ; L4 = D N1 L3 ((0.5+X1) 0.2) ; L5 = C N2 'PASS' L4 (0.5 ('-' 0.2 X1)) ((0.5-X1) 0.2) ; L6 = D N1 L5 ((0.5-X1) 1.) ; L7 = D N2 L6 (0. 1.) ; L8 = D N1 L7 (0. 0.) ; 'TRAC' L8 ; SI ITRAC ; FINS ; SI ITRAC ; FINS ; SI ITRAC ; FINS ; SI ITRAC ; FINS ; SI ITRAC ; FINS ; SI ITRAC ; FINS ; FIN ;
© Cast3M 2003 - Tous droits réservés.
Mentions légales