persistence_distance_matrix.py

import gtda.homology as hl

# represent data as a matrix of pairwise distances
distance_matrix = ...

# define topological features to track
homology_dimensions = [0, 1, 2]

# define simplicial complex to construct
persistence = hl.VietorisRipsPersistence(
    metric="precomputed", homology_dimensions=homology_dimensions
)

# calculate persistence diagram
persistence_diagram = persistence.fit_transform(distance_matrix)