 Read On file "maillage.msh"
 Nb of Vertex 73  Nb of Triangles 116
 End of read: area = 1
verification matelement 
0.000271369	3.98265e-05	1.82216e-05	-0.000102957	-7.93678e-05	-7.28863e-05	0	0	0	0	0	0	0	0.0042517	-0.0042517	
3.98265e-05	0.000271222	-2.26637e-06	-7.10881e-06	-0.000131969	-8.98763e-05	0	0	0	0	0	0	0	0.00192744	0.00192744	
1.82216e-05	-2.26637e-06	0.00024221	-9.14158e-05	-2.76923e-05	-0.00010281	0	0	0	0	0	0	0	-0.00232426	-0.00232426	
-0.000102957	-7.10881e-06	-9.14158e-05	0.000802754	0.000242526	0.000276179	0	0	0	0	0	0	0	-0.000581066	0.00406746	
-7.93678e-05	-0.000131969	-2.76923e-05	0.000242526	0.000818016	0.000329023	0	0	0	0	0	0	0	-0.0042517	-0.0042517	
-7.28863e-05	-8.98763e-05	-0.00010281	0.000276179	0.000329023	0.000817869	0	0	0	0	0	0	0	-0.00192744	-0.00192744	
0	0	0	0	0	0	0.000271369	3.98265e-05	1.82216e-05	-0.000102957	-7.93678e-05	-7.28863e-05	-0.00119048	-0.00119048	0	
0	0	0	0	0	0	3.98265e-05	0.000271222	-2.26637e-06	-7.10881e-06	-0.000131969	-8.98763e-05	-0.00396825	-0.00396825	0	
0	0	0	0	0	0	1.82216e-05	-2.26637e-06	0.00024221	-9.14158e-05	-2.76923e-05	-0.00010281	0.00277778	-0.00277778	0	
0	0	0	0	0	0	-0.000102957	-7.10881e-06	-9.14158e-05	0.000802754	0.000242526	0.000276179	-0.00277778	-0.00277778	0	
0	0	0	0	0	0	-7.93678e-05	-0.000131969	-2.76923e-05	0.000242526	0.000818016	0.000329023	0.00119048	0.00119048	0	
0	0	0	0	0	0	-7.28863e-05	-8.98763e-05	-0.00010281	0.000276179	0.000329023	0.000817869	-0.00396825	0.00396825	0	
0	0	0	0	0	0	-0.00119048	-0.00396825	0.00277778	-0.00277778	0.00119048	-0.00396825	-1e-07	0	0	
0.0042517	0.00192744	-0.00232426	-0.000581066	-0.0042517	-0.00192744	-0.00119048	-0.00396825	-0.00277778	-0.00277778	0.00119048	0.00396825	0	-1e-07	0	
-0.0042517	0.00192744	-0.00232426	0.00406746	-0.0042517	-0.00192744	0	0	0	0	0	0	0	0	-1e-07	
fin de verification
boucle 1 : 1
coordonne du point 0 : 0.142857
b[] 4e+100
boucle 1 : 4
coordonne du point 0 : 0.285714
b[] 4e+100
boucle 1 : 9
coordonne du point 0 : 0.428571
b[] 4e+100
boucle 1 : 12
coordonne du point 0 : 0.714286
b[] 4e+100
boucle 1 : 13
coordonne du point 0 : 0.571429
b[] 4e+100
boucle 1 : 23
coordonne du point 0 : 0.857143
b[] 4e+100
 nb coef = 5859 c max = 1e+100
valeur de status =0
valeur de status =1
valeur de status =1
-- FreeFem++ v  2.200000 (date Thu Sep 20 21:50:05 CEST 2007)
 Load: lg_fem lg_mesh eigenvalue  UMFPACK 
    1 : //border a(t=1,0){x=0;y=1+t;label=1;};
    2 : //border b(t=0,1){x=5*t;y=1;label=2;};
    3 : //border c(t=1,0){x=5 ;y=t;label=2;};
    4 : //border d(t=0,1){x=5+15*t;y=0;label=2;};
    5 : //border e(t=0,1){x=20 ;y=2*t;label=3;};
    6 : //border h(t=1,0){x=20*t;y=2;label=4;};
    7 : //mesh Th = buildmesh( a(7) + b(40) + c(10) + d(150) + e(5) + h(100));
    8 : 
    9 : mesh Th = readmesh("maillage.msh");
   10 : 
   11 : fespace Xh(Th,P2);		//discretisation pour la vitesse
   12 : fespace Mh(Th,P1);		//discretisation pour la pression
   13 : 
   14 : 
   15 : plot(Th,wait=1);
   16 : 
   17 : Xh ll;
   18 : 
   19 : ifstream ff("freefemU1.sol"); 
   20 : 
   21 : ff>>ll[]; 
   22 : 
   23 : 
   24 : 
   25 : plot(ll,wait=1,cmm="u1",value=1, ps="sol_cpp(U1).ps");
   26 : 
   27 : Xh ll2;
   28 : 
   29 : ifstream ff1("freefemU2.sol"); 
   30 : 
   31 : ff1>>ll2[]; 
   32 : 
   33 : 
   34 : 
   35 : plot(ll2,wait=1,cmm=" u2 ",value=1, ps="sol_cpp2(U2).ps");
   36 : 
   37 : 
   38 : 
   39 : Mh ll3;
   40 : ifstream ff2("freefemP.sol"); 
   41 : 
   42 : ff2>>ll3[]; 
   43 : 
   44 : 
   45 : 
   46 : 
   47 : plot(ll3,wait=1,cmm=" p ",value=1, ps="sol_cpp3(P).ps");
   48 : plot([ll,ll2],wait=1,cmm=" [u1,u2] ",value=1, ps="sol_cpp4(u1,u2).ps");
   49 : plot(coef=0.2,cmm=" [u1,u2] et p ",value=1,ll3,[ll,ll2], ps="sol_val.ps");
   50 :  sizestack + 1024 =1704  ( 680 )

  --  mesh:  Nb of Triangles =    116, Nb of Vertices 73
 Nb Of Nodes = 261
 Nb of DF = 261
 Nb Of Nodes = 73
 Nb of DF = 73
 TrueColor 
times: compile 0s, execution 0.04s
 CodeAlloc : nb ptr  2011,  size :134880
Bien: On a fini Normalement
the end
