Datasets examples

* fichier :  test_norm_env.dgibi

 GRAPH=VRAI ;
 GRAPH=FAUX ;
OPTI ECHO 0 ;

opti dime 2 elem tri6 ;
p1=1 0; p2=0 0 ;
l1= p1 d 1 p2 ;                                                         
m1= l1 trans 1 (0 1) ;                                                  
opti dime 2 elem qua8 ;
m2=(cote 3 m1) trans 1 (0 1) ;
m1 = m1 et m2 ;
Si GRAPH ;trace m1 ;Finsi;

Mm1= chan m1 QUAF ;
ERROR = 0 ;
*                        ** 2D ***

******************** LINE *******************************
 DISCR='LINE' ;
$m1=mode Mm1 'NAVIER_STOKES' DISCR ;
an= doma $m1  'NORMALEV' ;
an1= vect an 0.1 ux uy jaune;
Si GRAPH ; trace an1 (DOMA $m1 'MAILLAGE') ; Finsi;
*list an ;

'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
DB1 = 'DBIT' ('MODE' ('CONT' Mm1) 'NAVIER_STOKES' DISCR) an 'IMPR' ;
DB2 = 4.82842712 ;
err1=1.E-8 ;
Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;

******************** QUAF *******************************
 DISCR='QUAF' ;
$m1=mode Mm1 'NAVIER_STOKES' DISCR ;
an= doma $m1  'NORMALEV' ;
an1= vect an 0.1 ux uy jaune;
Si GRAPH ; trace an1 (DOMA $m1 'MAILLAGE') ; Finsi;
*list an ;

'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
DB1 = 'DBIT' ('MODE' ('CONT' Mm1) 'NAVIER_STOKES' DISCR) an 'IMPR' ;
DB2 = 5.60947571 ;
err1=1.E-8 ;
Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;

******************** MACRO ******************************
 DISCR='MACRO';
$m1=mode Mm1 'NAVIER_STOKES' DISCR ;
an= doma $m1  'NORMALEV' ;
an1= vect an 0.1 ux uy jaune;
Si GRAPH ; trace an1 (DOMA $m1 'MAILLAGE') ; Finsi;
*list an ;

'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
DB1 = 'DBIT' ('MODE' ('CONT' Mm1) 'NAVIER_STOKES' DISCR) an 'IMPR' ;
DB2 = 5.41421356 ;
err1=1.E-8 ;
Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;


**********************************
*                        ** 3D ***
**********************************

opti dime 3 elem cu20 ;
v1= m1 volu 1 trans (0 0  1) ;
Si GRAPH; trace v1 ; Finsi;
Mv1= chan QUAF v1 ;                                                     

******************** LINE *******************************
 DISCR='LINE' ;
$v1= mode Mv1 'NAVIER_STOKES' DISCR ;
na= doma $v1 'NORMALEV' ;
na1= vect na 0.1 ux uy uz jaune ;
Si GRAPH ; trace na1 (DOMA $v1 'MAILLAGE') ; Finsi;
*list an ;

'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
DB1 = 'DBIT' ('MODE' ('DOMA' $v1 'ENVELOPP') 'NAVIER_STOKES' DISCR)
        na 'IMPR' ;
DB2 = 6.32391111 ;
err1=1.E-8 ;
Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;

******************** QUAF *******************************
 DISCR='QUAF' ;
$v1= mode Mv1 'NAVIER_STOKES' DISCR ;
na= doma $v1 'NORMALEV' ;
na1= vect na 0.1 ux uy uz jaune ;
Si GRAPH ; trace na1 (DOMA $v1 'MAILLAGE') ; Finsi;
*list an ;

'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
DB1 = 'DBIT' ('MODE' ('DOMA' $v1 'ENVELOPP') 'NAVIER_STOKES' DISCR)
        na 'IMPR' ;
DB2 = 8.62820114 ;
err1=1.E-8 ;
Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;

******************** MACRO ******************************
 DISCR='MACRO' ;
$v1= mode Mv1 'NAVIER_STOKES' DISCR ;
na= doma $v1 'NORMALEV' ;
na1= vect na 0.1 ux uy uz jaune ;
Si GRAPH ; trace na1 (DOMA $v1 'MAILLAGE') ; Finsi;
*list an ;

'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
DB1 = 'DBIT' ('MODE' ('DOMA' $v1 'ENVELOPP') 'NAVIER_STOKES' DISCR)
        na 'IMPR' ;
DB2 = 7.87130633 ;
err1=1.E-8 ;
Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;

'SI' (ERROR > 0) ;
'MESS' 'Il y a ' ERROR ' erreur(s)' ;
'ERRE' 5 ;
'SINON' ;
'MESS' 'Tout s est bien passé !' ;
'FINSI' ;


FIN ;