Download test_inter.dgibi

Back to the list

* 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