load "msh3" load "iovtk" load "gmsh" load "MUMPS_seq" /* LABELS FRONT 1 BACK 2 LEFT 3 RIGHT 4 TOP 5 BOTTOM 6 INTERFACE 7 */ mesh3 cellunit=readmesh3("Unitcell.mesh"); savemesh(cellunit,"total.mesh"); plot(cellunit); func perio=[[1,x,z],[3,x,z],[2,y,z],[4,y,z],[5,x,y],[6,x,y]]; real Ce0=1.73e4; real k1=0.003; real k2=0.002; real n=0.2; real D=1.41e-14; real l=0.5; real Free=1; real Symmetry=2; fespace Vh(cellunit,P1, periodic=perio); fespace Ph(cellunit,P0);