1load "msh3"
2
3int upper = 1;
4int others = 2;
5int inner = 3;
6int n = 10;
7
8border D01(t=0, 1) {x=0; y=-1+t;z=3; }
9border D02(t=0, 1){x=1.5-1.5*t; y=-1; z=3;label=upper;}
10border D03(t=0, 1){x=1.5; y=-t; z=3;label=upper;}
11border D04(t=0, 1){x=1+0.5*t; y=0; z=3;label=others;}
12border D05(t=0, 1){x=0.5+0.5*t; y=0; z=3;label=others;}
13border D06(t=0, 1){x=0.5*t; y=0; z=3;label=others;}
14border D11(t=0, 1){x=0.5; y=-0.5*t; z=3;label=inner;}
15border D12(t=0, 1){x=0.5+0.5*t; y=-0.5; z=3;label=inner;}
16border D13(t=0, 1){x=1; y=-0.5+0.5*t; z=3;label=inner;}
17
18plot(D01(-n) + D02(-n) + D03(-n) + D04(-n) + D05(-n)
19   + D06(-n) + D11(n) + D12(n) + D13(n), wait=true);
20
21meshL ThL=buildmeshL(D01(-n) + D02(-n) + D03(-n) + D04(-n) + D05(-n)
22    + D06(-n) + D11(n) + D12(n) + D13(n));
23
24savemesh(ThL,"toto.mesh");
25
26//helice conique
27border E1(t=0, 10.*pi){x=(1.)*t*cos(t); y=-(1.)*t*sin(t); z=t;}
28meshL ThL2=buildmeshL(E1(1000));
29plot(ThL2);
30
31mesh Th=square(10,10);
32int[int] ll=[4];
33meshL ThL3 = extract(Th,[x+2,y*5],refedge=ll);