1 2verbosity=0; 3real x0=0.0, x1=1.5; 4real y0=0.0, y1=2.0; 5real xx0=0.5, xx1=1.; 6real yy0=0.5, yy1=1.; 7int acc,ref, i; 8real sH; 9ref=10; 10 11border a( t= x0,x1){ x= t ; y= y0; label=1; }; 12border b( t= y0,y1){ x= x0; y= t ; label=2; }; 13border c( t= x0,x1){ x= t ; y= y1; label=3; }; 14border d( t= y0,y1){ x= x1; y= t ; label=4; }; 15border aa( t= xx0,xx1){ x= t ; y= yy0; label=1; }; 16border bb( t= yy0,yy1){ x= xx0; y= t ; label=2; }; 17border cc( t= xx0,xx1){ x= t ; y= yy1; label=3; }; 18border dd( t= yy0,yy1){ x= xx1; y= t ; label=4; }; 19 20//cout << "enter accuracy = "; 21//cin >> acc; 22acc=0; 23// Creation des maillages grossier et fin 24 25//mesh TH=square(ref,ref,[x0+(x1-x0)*x,y0+(y1-y0)*y]); 26 27mesh TH = buildmesh( a(10) + b(-10) + c(-10) + d(+10) 28 +aa(-10) + bb(+10) + cc(+10) + dd(-10) ); 29plot(TH,wait=1); 30end; 31 32{ // hole 33real pi=4.0*atan(1); 34border a(t=0,2*pi){ x=cos(t); y=sin(t);label=1;}; 35border b(t=0,2*pi){ x=0.3+0.3*cos(t); y=0.3*sin(t);label=2;}; 36mesh Thwithouthole= buildmesh(a(50)+b(+30)); 37mesh Thwithhole = buildmesh(a(50)+b(-30)); 38plot(Thwithouthole,wait=1,ps="Thwithouthole.eps"); 39plot(Thwithhole,wait=1,ps="Thwithhole.eps"); 40}