* fichier : test_inter.dgibi
* Teste l'opérateur INTER
* en 2D :
'OPTION' 'DIME' 2 'ELEM' 'QUA4' ;
p1 = 0 0 ;
p2 = 0 1 ;
p3 = 1 1 ;
p4 = 1 0 ;
d1 = 'DROIT' p1 10 p2 ;
* intersection partielle
e1 = 'POINT' ('COOR' 2 d1) 'EGSU' .4 ;
e2 = 'POINT' ('COOR' 2 d1) 'EGIN' .6 ;
a = 'INTERSECTION' e1 e2 ;
* inclusion
e3 = 'CHANGER' d1 'POI1' ;
a = 'INTERSECTION' e1 e3 ;
* intersection avec un maillage vide
vide = 'DIFF' d1 d1 ;
a = 'INTERSECTION' d1 vide 'NOVERIF' ;
a = 'INTERSECTION' vide d1 'NOVERIF' ;
* intersection vide
e1 = 'POINT' ('COOR' 2 d1) 'EGSU' .6 ;
e2 = 'POINT' ('COOR' 2 d1) 'EGIN' .4 ;
a = 'INTERSECTION' e1 e2 'NOVERIF' ;
* produisant un message d'erreur :
*'OPTION' 'ERREUR' 'IGNORE' ;
*a = 'INTERSECTION' e1 e2 ;
*ee = valeu erma ;mess ' erreur maxi ' ee;
*'OPTION' 'ERREUR' 'NORMAL' ;
* en 3D :
'OPTION' 'DIME' 3 ;
p1 = 0 0 0 ;
p2 = 0 1 0 ;
p3 = 1 1 0 ;
p4 = 1 0 0 ;
p5 = 0 0 1 ;
p6 = 0 1 1 ;
p7 = 1 1 1 ;
p8 = 1 0 1 ;
a = 'MANU' 'CUB8' p1 p2 p3 p4 p5 p6 p7 p8;
b = ('MANU' 'CUB8' p1 p2 p3 p4 p5 p6 p7 p8) 'PLUS' (1 0 0) ;
d = ('CHANGER' a 'TET4') 'PLUS' (0 0 1);
e = ('CHANGER' b 'TET4') 'PLUS' (0 0 1);
'ELIMINATION' (d 'ET' e) 0.001 ;
* multizones croisées :
* sans conservation du pointeur de maillage
h = 'INTERSECTION' d (d 'ET' e) ;
* avec, car maillage inclus
h = 'INTERSECTION' d (a 'ET' d) ;
h = 'INTERSECTION' (a 'ET' d) d ;
* 2 maillages complexes
h = 'INTERSECTION' (d 'ET' e 'ET' a) (a 'ET' b 'ET' d) ;
h = 'INTERSECTION' (d 'ET' a) (a 'ET' (CHANGER POI1 e) 'ET' d) ;
'TEMPS' 'SGAC' 'IMPR';
'FIN' ;
© Cast3M 2003 - All rights reserved.
Disclaimer