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}