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);