load "msh3" load "medit" load "UMFPACK64" load "iovtk" load"gmsh" lockOrientation = false; // Global parameters verbosity = 0; real cpu=clock(); int N=1; mesh3 Th=readmesh3("Th.mesh"); mesh3 Th2=readmesh3("CC_simple_geo_2.mesh"); mesh3 Th3 = Th + Th2; lockOrientation = true; fespace Ph(Th, P2); fespace Ph2(Th2, P2); Ph condair = 0.024; Ph2 condpu = 0.29; fespace Ph3(Th3, P2,periodic=[[3,y,z],[5,y,z],[4,y,z],[6,y,z],[26,y,x],[600,y,x],[2,y,x],[1,y,x],[10,y,x],[34,y,x],[30,y,x],[14,y,x],[18,y,x],[22,y,x],[28,y,z],[21,y,z],[24,y,z],[33,y,z],[77,y,z],[8,y,z],[13,y,z],[17,y,z],[27,x,z],[32,x,z],[19,x,z],[25,x,z],[100,x,z],[15,x,z],[9,x,z],[11,x,z]]); Ph3 cond= condair+condpu; mesh3 renumbered = Th3; fespace Vh(renumbered, P2,periodic=[[3,y,z],[5,y,z],[4,y,z],[6,y,z],[26,y,x],[600,y,x],[2,y,x],[1,y,x],[10,y,x],[34,y,x],[30,y,x],[14,y,x],[18,y,x],[22,y,x],[28,y,z],[21,y,z],[24,y,z],[33,y,z],[77,y,z],[8,y,z],[13,y,z],[17,y,z],[27,x,z],[32,x,z],[19,x,z],[25,x,z],[100,x,z],[15,x,z],[9,x,z],[11,x,z]]); Vh u = x + y + z, wh1,wh2,wh3,vh,coeffh;