gap> K:=ClosedSurface(1);;
gap> facets:=MaximalSimplicesOfSimplicialComplex(K);;
gap> LoadPackage("simpcomp");;
gap> scK:=SCFromFacets(S);;
gap> scY:=SCCartesianProduct(scK,scK);;
gap> Y:=MaximalSimplicesToSimplicialComplex(SCFacetsEx(scY));
Simplicial complex of dimension 4.

gap> cup:=CupProduct(Y);;
gap> cup(2,2,[1,0,0,0,0,0],[0,1,0,0,0,0]);
[ 0 ]
