1include "cube.idp"
2load "medit"
3int[int, int] LL = [[2,2], [1,1], [1,1]];
4real[int, int] BB = [[0,1], [0,1], [0,1]];
5int[int] CC = [4, 4, 4];
6load "splitmesh4"
7mesh3 Th3=Cube(CC, BB, LL);
8mesh3 Th4=splitmesh4(Th3);
9assert(Th4.nt == Th3.nt * 4);
10cout << Th4.nt << " == 4 x " << Th3.nt << endl;
11assert(Th4.nbe == Th3.nbe);
12cout << Th4.nbe << " == " << Th3.nbe << endl;
13plot(Th3, cmm="Th3", wait=1);
14plot(Th4, cmm="Th4", wait=1);
15fespace Vh(Th4, P1);
16Vh xyz = x + y + z;
17medit("splitmesh4",Th4,xyz);
18