* fichier :  test_lapn.dgibi
'OPTI' 'DIME' 3 'ELEM' 'QUA4' ;

*** AX
AX11 =  1./3. ; AX12 = -1./3. ; AX13 = -1./6. ; AX14 =  1./6. ;
AX21 = -1./3. ; AX22 =  1./3. ; AX23 =  1./6. ; AX24 = -1./6. ;
AX31 = -1./6. ; AX32 =  1./6. ; AX33 =  1./3. ; AX34 = -1./3. ;
AX41 =  1./6. ; AX42 = -1./6. ; AX43 = -1./3. ; AX44 =  1./3. ;

*** AY
AY11 =  1./3. ; AY12 =  1./6. ; AY13 = -1./6. ; AY14 = -1./3. ;
AY21 =  1./6. ; AY22 =  1./3. ; AY23 = -1./3. ; AY24 = -1./6. ;
AY31 = -1./6. ; AY32 = -1./3. ; AY33 =  1./3. ; AY34 =  1./6. ;
AY41 = -1./3. ; AY42 = -1./6. ; AY43 =  1./6. ; AY44 =  1./3. ;

AS11 = (0.5*AX11) + (2.0*AY11) ; AS12 = (0.5*AX12) + (2.0*AY12) ;
AS13 = (0.5*AX13) + (2.0*AY13) ; AS14 = (0.5*AX14) + (2.0*AY14) ;
AS21 = (0.5*AX21) + (2.0*AY21) ; AS22 = (0.5*AX22) + (2.0*AY22) ;
AS23 = (0.5*AX23) + (2.0*AY23) ; AS24 = (0.5*AX24) + (2.0*AY24) ;
AS31 = (0.5*AX31) + (2.0*AY31) ; AS32 = (0.5*AX32) + (2.0*AY32) ;
AS33 = (0.5*AX33) + (2.0*AY33) ; AS34 = (0.5*AX34) + (2.0*AY34) ;
AS41 = (0.5*AX41) + (2.0*AY41) ; AS42 = (0.5*AX42) + (2.0*AY42) ;
AS43 = (0.5*AX43) + (2.0*AY43) ; AS44 = (0.5*AX44) + (2.0*AY44) ;
LAS1 = 'PROG' AS11 AS12 AS13 AS14 ;
LAS2 = 'PROG' AS21 AS22 AS23 AS24 ;
LAS3 = 'PROG' AS31 AS32 AS33 AS34 ;
LAS4 = 'PROG' AS41 AS42 AS43 AS44 ;

P1 = 0.0  0.0 0.0 ;
P2 = 0.0  2.0 0.0 ;
P3 = 0.0  2.0 1.0 ;
P4 = 0.0  0.0 1.0 ;

MODT = 'MANU' 'QUA4' P1 P2 P3 P4 ;
$MODT = 'MODE' ('CHAN' MODT 'QUAF') 'NAVIER_STOKES' 'LINE' ;


TTT = 'EQEX' 'NITER' 1 'OMEGA' 1. 'ITMA' 0
      'OPTI' 'EF' 'IMPL'
      'ZONE' $MODT 'OPER' 'LAPN' 1.0 'INCO' 'TN' ;
TTT.'INCO' = 'TABLE' 'INCO' ;
TNS  = 'KCHT' $MODT 'SCAL' 'SOMMET' 1.0 ;
TTT.'INCO'.'TN'  = TNS ;
B1 MatA = 'LAPN' TTT.'1LAPN' ;

C1 = 'KOPS' MatA 'MULT' ('MANU' 'CHPO' P1 1 'TN' 1.0 'NATU' 'DISCRET');
C11 = 'EXTR' C1 'TN' P1 ;
C21 = 'EXTR' C1 'TN' P2 ;
C31 = 'EXTR' C1 'TN' P3 ;
C41 = 'EXTR' C1 'TN' P4 ;

C2 = 'KOPS' MatA 'MULT' ('MANU' 'CHPO' P2 1 'TN' 1.0 'NATU' 'DISCRET');
C12 = 'EXTR' C2 'TN' P1 ;
C22 = 'EXTR' C2 'TN' P2 ;
C32 = 'EXTR' C2 'TN' P3 ;
C42 = 'EXTR' C2 'TN' P4 ;

C3 = 'KOPS' MatA 'MULT' ('MANU' 'CHPO' P3 1 'TN' 1.0 'NATU' 'DISCRET');
C13 = 'EXTR' C3 'TN' P1 ;
C23 = 'EXTR' C3 'TN' P2 ;
C33 = 'EXTR' C3 'TN' P3 ;
C43 = 'EXTR' C3 'TN' P4 ;

C4 = 'KOPS' MatA 'MULT' ('MANU' 'CHPO' P4 1 'TN' 1.0 'NATU' 'DISCRET');
C14 = 'EXTR' C4 'TN' P1 ;
C24 = 'EXTR' C4 'TN' P2 ;
C34 = 'EXTR' C4 'TN' P3 ;
C44 = 'EXTR' C4 'TN' P4 ;

LBS1 = 'PROG' C11 C12 C13 C14 ;
LBS2 = 'PROG' C21 C22 C23 C24 ;
LBS3 = 'PROG' C31 C32 C33 C34 ;
LBS4 = 'PROG' C41 C42 C43 C44 ;

LERT = (LAS1 ET LAS2 ET LAS3 ET LAS4) - (LBS1 ET LBS2 ET LBS3 ET LBS4) ;

'MESS' 'Maximum erreur :' ('MAXI' LERT 'ABS') ;

'SI' (('MAXI' LERT 'ABS') > 1.E-14) ;
'MESS' 'Il y a des problèmes sur la matrice du laplacien' ;
'ERRE' 5 ;
'FINSI' ;

'FIN' ;