An output-sensitive algorithm for persistent homology

  • Authors:
  • Chao Chen;Michael Kerber

  • Affiliations:
  • Institute of Science and Technology (IST) Austria, Klosterneuburg, Austria;Institute of Science and Technology (IST) Austria, Klosterneuburg, Austria

  • Venue:
  • Computational Geometry: Theory and Applications
  • Year:
  • 2013

Quantified Score

Hi-index 0.00

Visualization

Abstract

In this paper, we present the first output-sensitive algorithm to compute the persistence diagram of a filtered simplicial complex. For any @C0, it returns only those homology classes with persistence at least @C. Instead of the classical reduction via column operations, our algorithm performs rank computations on submatrices of the boundary matrix. For an arbitrary constant @d@?(0,1), the running time is O(C"("1"-"@d")"@CR"d(n)logn), where C"("1"-"@d")"@C is the number of homology classes with persistence at least (1-@d)@C, n is the total number of simplices in the complex, d its dimension, and R"d(n) is the complexity of computing the rank of an nxn matrix with O(dn) nonzero entries. Depending on the choice of the rank algorithm, this yields a deterministic O(C"("1"-"@d")"@Cn^2^.^3^7^6) algorithm, an O(C"("1"-"@d")"@Cn^2^.^2^8) Las-Vegas algorithm, or an O(C"("1"-"@d")"@Cn^2^+^@e) Monte-Carlo algorithm for an arbitrary @e0. The space complexity of the Monte-Carlo version is bounded by O(dn)=O(nlogn).