Persistent Cohomology

## Classes | |

class | Gudhi::persistent_cohomology::Field_Zp |

Structure representing the coefficient field . More... | |

class | Gudhi::persistent_cohomology::Multi_field |

Structure representing coefficients in a set of finite fields simultaneously using the chinese remainder theorem. More... | |

class | Gudhi::persistent_cohomology::Persistent_cohomology< FilteredComplex, CoefficientField > |

Computes the persistent cohomology of a filtered complex. More... | |

Computation of persistent cohomology using the algorithm of [9] and [10] and the Compressed Annotation Matrix implementation of [6]

The theory of homology consists in attaching to a topological space a sequence of (homology) groups, capturing global topological features like connected components, holes, cavities, etc. Persistent homology studies the evolution – birth, life and death – of these features when the topological space is changing. Consequently, the theory is essentially composed of three elements: topological spaces, their homology groups and an evolution scheme.

Topological spaces are represented by simplicial complexes. Let be a set of *vertices*. A *simplex* is a subset of vertices . A *simplicial complex* on is a collection of simplices , , such that . The dimension of is its number of elements minus 1. A *filtration* of a simplicial complex is a function satisfying whenever .

We define the concept FilteredComplex which enumerates the requirements for a class to represent a filtered complex from which persistent homology may be computed. We use the vocabulary of simplicial complexes, but the concept is valid for any type of cell complex. The main requirements are the definition of:

- type
`Indexing_tag`

, which is a model of the concept`IndexingTag`

, describing the nature of the indexing scheme, - type Simplex_handle to manipulate simplices,
- method
`int dimension(Simplex_handle)`

returning the dimension of a simplex, - type and method
`Boundary_simplex_range boundary_simplex_range(Simplex_handle)`

that returns a range giving access to the codimension 1 subsimplices of the input simplex, as-well-as the coefficients in the definition of the operator . The iterators have value type`Simplex_handle`

, - type and method
`Filtration_simplex_range filtration_simplex_range ()`

that returns a range giving access to all the simplices of the complex read in the order assigned by the indexing scheme, - type and method
`Filtration_value filtration (Simplex_handle)`

that returns the value of the filtration on the simplex represented by the handle.

For a ring , the group of *n-chains*, denoted , of is the group of formal sums of n-simplices with coefficients. The *boundary operator* is a linear operator such that , where means is omitted from the list. The chain groups form a sequence:

of finitely many groups and homomorphisms , indexed by the dimension . The boundary operators satisfy the property for every and we define the homology groups:

We refer to [17] for an introduction to homology theory and to [11] for an introduction to persistent homology.

"Changing" a simplicial complex consists in applying a simplicial map. An *indexing scheme* is a directed graph together with a traversal order, such that two consecutive nodes in the graph are connected by an arrow (either forward or backward). The nodes represent simplicial complexes and the directed edges simplicial maps.

From the computational point of view, there are two types of indexing schemes of interest in persistent homology: *linear* ones in persistent homology [21] , and *zigzag* ones in zigzag persistent homology [7]. These indexing schemes have a natural left-to-right traversal order, and we describe them with ranges and iterators. In the current release of the Gudhi library, only the linear case is implemented.

In the following, we consider the case where the indexing scheme is induced by a filtration. Ordering the simplices by increasing filtration values (breaking ties so as a simplex appears after its subsimplices of same filtration value) provides an indexing scheme.

We provide several example files: run these examples with -h for details on their use, and read the README file.

- Persistent_cohomology/rips_persistence.cpp computes the Rips complex of a point cloud and outputs its persistence diagram. $> ./rips_persistence ../../data/points/tore3D_1307.off -r 0.25 -m 0.5 -d 3 -p 3The complex contains 177838 simplicesand has dimension 33 0 0 inf3 1 0.0983494 inf3 1 0.104347 inf3 2 0.138335 inf

- Persistent_cohomology/rips_multifield_persistence.cpp computes the Rips complex of a point cloud and outputs its persistence diagram with a family of field coefficients.

- Persistent_cohomology/rips_distance_matrix_persistence.cpp computes the Rips complex of a distance matrix and outputs its persistence diagram.

- Persistent_cohomology/alpha_complex_3d_persistence.cpp computes the persistent homology with coefficients of the alpha complex on points sampling from an OFF file. $> ./alpha_complex_3d_persistence ../../data/points/tore3D_300.off 2 0.45Simplex_tree dim: 32 0 0 inf2 1 0.0682162 1.00012 1 0.0934117 1.000032 2 0.56444 1.03938

- Persistent_cohomology/exact_alpha_complex_3d_persistence.cpp computes the persistent homology with coefficients of the alpha complex on points sampling from an OFF file. Here, as CGAL computes the exact values, it is slower, but it is necessary when points are on a grid for instance. $> ./exact_alpha_complex_3d_persistence ../../data/points/sphere3D_pts_on_grid.off 2 0.1Simplex_tree dim: 32 0 0 inf2 2 0.0002 0.2028

- Persistent_cohomology/weighted_alpha_complex_3d_persistence.cpp computes the persistent homology with coefficients of the weighted alpha complex on points sampling from an OFF file and a weights file. $> ./weighted_alpha_complex_3d_persistence ../../data/points/tore3D_300.off../../data/points/tore3D_300.weights 2 0.45Simplex_tree dim: 32 -0 0 inf2 1 0.0682162 1.00012 1 0.0934117 1.000032 2 0.56444 1.03938

- Persistent_cohomology/alpha_complex_persistence.cpp computes the persistent homology with coefficients of the alpha complex on points sampling from an OFF file. $> ./alpha_complex_persistence -r 32 -p 2 -m 0.45 ../../data/points/tore3D_300.offAlpha complex is of dimension 3 - 9273 simplices - 300 vertices.Simplex_tree dim: 32 0 0 inf2 1 0.0682162 1.00012 1 0.0934117 1.000032 2 0.56444 1.03938

- Persistent_cohomology/periodic_alpha_complex_3d_persistence.cpp computes the persistent homology with coefficients of the periodic alpha complex on points sampling from an OFF file. $> ./periodic_alpha_complex_3d_persistence ../../data/points/grid_10_10_10_in_0_1.off 3 1.0Periodic Delaunay computed.Simplex_tree dim: 33 0 0 inf3 1 0.0025 inf3 1 0.0025 inf3 1 0.0025 inf3 2 0.005 inf3 2 0.005 inf3 2 0.005 inf3 3 0.0075 inf

- Persistent_cohomology/plain_homology.cpp computes the plain homology of a simple simplicial complex without filtration values.

- Copyright
- GNU General Public License v3.

GUDHI Version 2.0.0 - C++ library for Topological Data Analysis (TDA) and Higher Dimensional Geometry Understanding. | Generated on Wed Apr 19 2017 22:26:16 for GUDHI by 1.8.11 |