Cubical Complexes

ArrayToPureCubicalComplex
ArrayToPureCubicalComplexA,n)

Inputs an integer array A of dimension d and an integer n. It returns a d-dimensional pure cubical complex corresponding to the black/white "image" determined by the threshold n and the values of the entries in A. (Integers below the threshold correspond to a black pixel, and higher integers correspond to a white pixel.)
PureCubicalComplex
PureCubicalComplexA,n)

Inputs a binary array A of dimension d. It returns the corresponding d-dimensional pure cubical complex.
FramedPureCubicalComplex
FramedPureCubicalComplex(M)

Inputs a pure cubical complex M and returns the pure cubical complex with a border of zeros attached the each face of the boundary array M!.boundaryArray. This function just adds a bit of space for performing operations such as thickenings to M.
RandomCubeOfPureCubicalComplex
RandomCubeOfPureCubicalComplex(M)

Inputs a pure cubical complex M and returns a pure cubical complex R with precisely the same dimensions as M. The complex R consist of one cube selected at random from M.

PureCubicalComplexIntersection
PureCubicalComplexIntersection(S,T)

Inputs two pure cubical complexes with common dimension and array size. It returns the intersection of the two complexes. (An entry in the binary array of the intersection has value 1 if and only if the corresponding entries in the binary arrays of S and T both have value 1.)
PureCubicalComplexUnion
PureCubicalComplexUnion(S,T)

Inputs two pure cubical complexes with common dimension and array size. It returns the union of the two complexes. (An entry in the binary array of the union has value 1 if and only if at least one of the corresponding entries in the binary arrays of S and T has value 1.)

PureCubicalComplexDifference
PureCubicalComplexDifference(S,T)

Inputs two pure cubical complexes with common dimension and array size. It returns the difference S-T. (An entry in the binary array of the difference has value 1 if and only if the corresponding entry in the binary array of S is 1 and the corresponding entry in the binary array of T is 0.)

ReadImageAsPureCubicalComplex
ReadImageAsPureCubicalComplex("file.png",n)

Reads an image file ("file.png", "file.eps", "file.bmp" etc) and an integer n between 0 and 765. It returns a 2-dimensional pure cubical complex based on the black/white version of the image determined by the threshold n.
ReadLinkImageAsPureCubicalComplex
ReadLinkImageAsPureCubicalComplex("file.png")
ReadLinkImageAsPureCubicalComplex("file.png",n)

Reads an image file ("file.png", "file.eps", "file.bmp" etc) containing a knot or link diagram, and optionally a positive integer n. The integer n should be a little larger than the line thickness in the link diagram, and if not provided then n is set equal to 10.

The function tries to output the corresponding knot or link as a 3-dimensional pure cubical complex. Ideally the link diagram should be produced with line thickness 6 in Xfig, and the under-crossing spaces should not be too large or too small or too near one another. The function does not always succeed: it applies several checks, and if one of these checks fails then the function returns "fail".
ReadImageSequenceAsPureCubicalComplex
ReadImageSequenceAsPureCubicalComplex("directory",n)

Reads the name of a directory containing a sequence of image files (ordered alphanumerically), and an integer n between 0 and 765. It returns a 3-dimensional pure cubical complex based on the black/white version of the images determined by the threshold n.
Size(T)

This returns the number of non-zero entries in the binary array of the cubical complex, or pure cubical complex T.

Dimension(T)

This returns the dimension of the cubical complex, or pure cubical complex T.

WritePureCubicalComplexAsImage
WritePureCubicalComplexAsImage(T,"filename","ext")

Inputs a 2-dimensional pure cubical complex T, and a filename followed by its extension (e.g. "myfile" followed by "png"). A black/white image is saved to the file.

ViewPureCubicalComplex
ViewPureCubicalComplex(T)
ViewPureCubicalComplex(T,"mozilla")

Inputs a 2-dimensional pure cubical complex T, and optionally a command such as "mozilla" for viewing image files. A black/white image is displayed.
Homology
Homology(T,n)
Homology(T)

Inputs a pure cubical complex, or cubical complex, or simplicial complex T and a non-negative integer n. It returns the n-th integral homology of T as a list of torsion integers. If no value of n is input then the list of all homologies of T in dimensions 0 to Dimension(T) is returned .
Bettinumbers
Bettinumbers(T,n)
Bettinumbers(T)

Inputs a pure cubical complex, or cubical complex, simplicial complex or chain complex T and a non-negative integer n. The rank of the n-th rational homology group H_n(T,\mathbb Q) is returned. If no value for n is input then the list of Betti numbers in dimensions 0 to Dimension(T) is returned .
DirectProductOfPureCubicalComplexes
DirectProductOfPureCubicalComplexes(M,N)

Inputs two pure cubical complexes M,N and returns their direct product D as a pure cubical complex. The dimension of D is the sum of the dimensions of M and N.
SuspensionOfPureCubicalComplex
SuspensionOfPureCubicalComplex(M)

Inputs a pure cubical complex M and returns a pure cubical complex with the homotopy type of the suspension of M.
EulerCharacteristic(T)

Inputs a pure cubical complex, or cubical complex, or simplicial complex T and returns its Euler characteristic.
PathComponentOfPureCubicalComplex
PathComponentOfPureCubicalComplex(T,n)

Inputs a pure cubical complex T and an integer n in the rane 1, ..., Bettinumbers(T)[1] . It returns the n-th path component of T as a pure cubical complex. The value n=0 is also allowed, in which case the number of path components is returned.
ChainComplex(T)

Inputs a pure cubical complex, or cubical complex, or simplicial complex T and returns the (often very large) cellular chain complex of T.
ChainComplexOfPair(T,S)

Inputs a pure cubical complex or cubical complex T and subcomplex S. It returns the quotient C(T)/C(S) of cellular chain complexes.
ExcisedPureCubicalPair(T,S)

Inputs a pure cubical complex T and subcomplex S. It returns the pair [T\setminus intS, S\setminus intS]) of pure cubical complexes where intS is the pure cubical complex obtained from S by removing its boundary.
ChainInclusionOfPureCubicalPair
ChainInclusionOfPureCubicalPair(S,T)

Inputs a pure cubical complex T and subcomplex S. It returns the chain inclusion C(S) \rightarrow C(T) of cellular chain complexes.
ChainMapOfPureCubicalPairs
ChainMapOfPureCubicalPairs(M,S,N,T)

Inputs a pure cubical complex N and subcomplexes M, T and S in T. It returns the chain map C(M/S) \rightarrow C(N/T) of quotient cellular chain complexes.
ContractPureCubicalComplex
ContractPureCubicalComplex(T)

Inputs a pure cubical complex T of dimension d and removes d-dimensional cells from T without changing the homotopy type of T. When the function has been applied, no further d-cells can be removed from T without changing its homotopy type. This function modifies T.

ContractedComplex
ContractedComplex(T)

Inputs a pure cubical complex T and returns a structural copy of the complex obtained from T by applying the function ContractPureCubicalComplex(T).

ZigZagContractedPureCubicalComplex
ZigZagContractedPureCubicalComplex(T)

Inputs a pure cubical complex T and returns a homotopy equivalent pure cubical complex S. The aim is for S to involve fewer cells than T and certainly to involve no more cells than T.

ContractCubicalComplex
ContractCubicalComplex(T)

Inputs a cubical complex T and removes cells without changing the homotopy type of T. It changes T. In particular, it adds the components T.vectors and T.rewrite of a discrete vector field.

At present this function only works for cubical complexes of dimension 2 or 3.
DVFReducedCubicalComplex
DVFReducedCubicalComplex(T)

Inputs a cubical complex T and returns a non-regular cubical complex R by constructing a discrete vector field. The vector field is designed to minimize the number of critical cells in R at the cost of allowing cell attaching maps that are not homeomorphisms on boundaries.

At present this function works only for 2- and 3-dimensional cubical complexes.

The function ChainComplex(R) can be used to obtain the cellular chain complex of R.

SkeletonOfCubicalComplex
SkeletonOfCubicalComplex(T,n)

Inputs a cubical complex, or pure cubical complex T and positive integer n. It returns the n-skeleton of T as a cubical complex.

ContractibleSubomplexOfPureCubicalComplex
ContractibleSubomplexOfPureCubicalComplex(T)

Inputs a pure cubical complex T and returns a maximal contractible pure cubical subcomplex.
AcyclicSubomplexOfPureCubicalComplex
AcyclicSubomplexOfPureCubicalComplex(T)

Inputs a pure cubical complex T and returns a (not necessarily connected) pure cubical subcomplex having trivial homology in all degrees greater than 0.
HomotopyEquivalentMaximalPureCubicalSubcomplex
HomotopyEquivalentMaximalPureCubicalSubcomplex(T,S)

Inputs a pure cubical complex T together with a pure cubical subcomplex S. It returns a pure cubical subcomplex H of T which contains S and is maximal with respect to the property that it is homotopy equivalent to S.
HomotopyEquivalentMinimalPureCubicalSubcomplex
HomotopyEquivalentMinimalPureCubicalSubcomplex(T,S)

Inputs a pure cubical complex T together with a pure cubical subcomplex S. It returns a pure cubical subcomplex H of T which contains S and is minimal with respect to the property that it is homotopy equivalent to T.
BoundaryOfPureCubicalComplex
BoundaryOfPureCubicalComplex(T)

Inputs a pure cubical complex T and returns its boundary as a pure cubical complex. The boundary consists of all cubes which have one or more facets that lie in just the one cube.

SingularitiesOfPureCubicalComplex
SingularitiesOfPureCubicalComplex(T,radius,tolerance)

Inputs a pure cubical complex T together with a positive integer "radius" and an integer "tolerance" in the range 1..100. It returns the pure cubical subcomplex of those cells in the boundary where the boundary is not differentiable. (The method for deciding differentiability at a point is crude/discrete, prone to errors and depends on the radius and tolerance.)
ThickenedPureCubicalComplex
ThickenedPureCubicalComplex(T)

Inputs a pure cubical complex T and returns a pure cubical complex S. If a euclidean cube is in T then this cube and all its neighbouring cubes are included in S.

CropPureCubicalComplex
CropPureCubicalComplex(T)

Inputs a pure cubical complex T and returns a pure cubical complex S obtained from T by removing any "zero boundary sheets" of the binary array. Thus S and T are isometric as euclidean spaces but there may be fewer zero entries in the binary array for S.
BoundingPureCubicalComplex
BoundingPureCubicalComplex(T)

Inputs a pure cubical complex T and returns a contractible pure cubical complex S containing T.
MorseFiltration
MorseFiltration(M,i,t,bool)
MorseFiltration(M,i,t)

Inputs a pure cubical complex M of dimension d, an integer i between 1 and d, a positive integer t and a boolean value True or False. The function returns a list [M_1, M_2, ..., M_t] of pure cubical complexes with M_k a subcomplex of M_{k+1}. The list is constructed by setting all slices of M perpendicular to the i-th axis equal to zero if they meet the ith axis at a sufficiently high coordinate (if bool=True) or sufficiently low coordinate (if bool=False).

If the variable bool is not specified then it is assumed to have the value True.
ComplementOfPureCubicalComplex
ComplementOfPureCubicalComplex(T)

Inputs a pure cubical complex T and returns a pure cubical complex S. A euclidean cube is in S precisely when the cube is not in T.

PureCubicalComplexToTextFile
PureCubicalComplexToTextFile(file,M)

Inputs a pure cubical complex M and a string containing the address of a file. A representation of this complex is written to the file in a format that can be read by the CAPD (Computer Assisted Proofs in Dynamics) software developed by Marian Mrozek and others.

ThickeningFiltration
ThickeningFiltration(M,n)
ThickeningFiltration(M,n,k)

Inputs a pure cubical complex M and a positive integer n. It returns a filtered pure cubical complex constructed frim n thickenings of M. If a positive integer k is supplied as an optional third argument, then each step of the filtration is obtained from a k-fold thickening.
Dendrogram
Dendrogram(M)

Inputs a filtered pure cubical complex M and returns data that specifies the dendrogram (or phylogenetic tree) describing how path components are born and then merge during the filtration.

DendrogramDisplay
DendrogramDisplay(M)

Inputs a filtered pure cubical complex M, or alternatively inputs the out from the command Dendrogram(M), and then uses GraphViz software to display the path component dendrogram of M.
DendrogramToPersistenceMat
DendrogramToPersistenceMat(D)

Inputs the output of the function Dendrogram(M) and returns the corresponding degree 0 Betti bar code.
ReadImageAsFilteredPureCubicalComplex
ReadImageAsFilteredPureCubicalComplex(file,n)

Inputs a string containing the path to an image file, together with a positive integer n. It returns a filtered pure cubical complex of filtration length n.
ComplementOfFilteredPureCubicalComplex
ComplementOfFilteredPureCubicalComplex(M)

Inputs a filtered pure cubical complex M and returns the complement as a filtered pure cubical complex.
PersistentHomologyOfFilteredPureCubicalComplex
PersistentHomologyOfFilteredPureCubicalComplex(M,n)

Inputs a filtered pure cubical complex M and a non-negative integer n. It returns the degree n persistent homology of M with rational coefficients.