load "msh3" mesh3 Th3=cube(2,2,2); mesh3 Thg=trunc(Th3,1,split=4);// to have correct graphique verbosity = 4; fespace Gh(Thg,P1); fespace Eh(Th3,P0edge); Th3 = buildBdMesh(Th3); meshS ThS = Th3.Gamma; //plot(ThS,wait=1); // get edges of the 2D surface mesh ThS = buildBdMesh(ThS); meshL ThL=ThS.Gamma; plot(ThL,wait=1); cout << " ThS :"<< ThS.nt << " "<< ThS.mesure <<" "<< ThS.bordermesure << endl; cout << " ThL :"<< ThL.nt << " "<< ThL.mesure <<" "<< ThS.bordermesure << endl; fespace Ph(ThL,P0); Ph un =1; Eh U3=un; Gh GU3=U3; //plot(GU3,wait=1);