S:=SimplicialComplex([[1,2],[2,3],[3,1]]);;
S:=RegularCWComplex(S);;
T:=DirectProduct(S,S,S,S);
Print(Cohomology(T,1),"\n");
Print(Cohomology(T,3),"\n");
Print(Cohomology(T,4),"\n");
cup:=CupProduct(T);
Print(cup(3,1,[0,0,1,0],[0,0,1,0]),"\n");
