gap> A:=[[[0,0,0],[0,0,0],[0,0,0]],
>        [[1,1,1],[1,0,1],[1,1,1]],
>        [[0,0,0],[0,0,0],[0,0,0]]];;
gap> M:=PureCubicalComplex(A);
Pure cubical complex of dimension 3.

gap> P:=PurePermutahedralComplex(A);
Pure permutahedral complex of dimension 3.

gap> Display(M);
gap> Display(P);
