[论文解读] Computing height persistence and homology generators in R3 efficiently

Tamal K. Dey|arXiv (Cornell University)|Jan 6, 2019

Topological and Geometric Data Analysis被引用 2

一句话总结

本文提出了一种时间复杂度为 O(n log n) 的算法，用于计算在 R³ 中线性嵌入的单纯复形的高程持久性以及同调生成元（H₁ 和 H₂），通过利用交错持久性、Reeb 图的性质以及高效的几何数据结构，显著优于此前已知的 O(n^ω) 时间复杂度上界。

ABSTRACT

Recently it has been shown that computing the dimension of the first homology group H1(K) of a simplicial 2-complex K embedded linearly in R4 is as hard as computing the rank of a sparse 0 − 1 matrix. This puts a major roadblock to computing persistence and a homology basis (generators) for complexes embedded in R4 and beyond in less than quadratic or even near-quadratic time. But, what about dimension three? It is known that when K is a graph or a surface with n simplices linearly embedded in R3, the persistence for piecewise linear functions on K can be computed in O(n log n) time and a set of generators of total size k can be computed in O(n + k) time. However, the question for general simplicial complexes K linearly embedded in R3 is not completely settled. No algorithm with a complexity better than that of the matrix multiplication is known for this important case. We show that the persistence for height functions on such complexes, hence called height persistence, can be computed in O(n log n) time. This allows us to compute a basis (generators) of Hi(K), i = 1, 2, in O(n log n + k) time where k is the size of the output. This improves significantly the current best bound of O(nω), ω being the exponent of matrix multiplication. We achieve these improved bounds by leveraging recent results on zigzag persistence in computational topology, new observations about Reeb graphs, and some efficient geometric data structures.

研究动机与目标

• 解决在 R³ 中嵌入的一般单纯复形的高程持久性与同调生成元高效计算这一开放问题。
• 克服尽管在图和曲面上已有相关结果，但该情形下仍缺乏时间复杂度低于二次方或接近二次方的算法这一局限。
• 实现计算 R³ 嵌入复形中同调生成元的时间复杂度优于矩阵乘法。
• 将高效拓扑计算的适用范围扩展至一般 2-复形，而不仅限于曲面或图等特殊情况。

提出的方法

• 利用交错持久性领域的最新进展，对 2-复形上的高程函数建模并计算持久性。
• 应用关于 Reeb 图的新几何观察，以简化并加速拓扑不变量的计算。
• 采用高效的几何数据结构，在持久性计算过程中维护并查询拓扑特征。
• 将同调生成元的计算问题转化为一系列可在 O(n log n) 时间内处理的持久同调操作序列。
• 结合分治策略与基于高程函数的扫描线技术，以动态维护连通性与环路信息。
• 通过将高程函数持久性问题约化为一系列保持同调结构的同时支持高效更新的交错滤子序列。

实验结果

研究问题

• RQ1 能否在亚二次时间内计算 R³ 中嵌入的 2-复形的高程持久性？
• RQ2 是否可能以 O(n log n + k) 时间计算 H₁ 和 H₂ 的基，其中 k 为输出大小？
• RQ3 能否在 R³ 嵌入复形的同调计算中避免矩阵乘法的计算瓶颈？
• RQ4 如何利用 Reeb 图的性质来加速 3D 复形中拓扑持久性的计算？

主要发现

• 本文在 R³ 中嵌入的 2-复形上实现了 O(n log n) 的时间复杂度来计算高程持久性，优于基于矩阵乘法的 O(n^ω) 上界。
• H₁(K) 和 H₂(K) 的基可在 O(n log n + k) 时间内计算，其中 k 为输出大小，相比先前方法有显著改进。
• 交错持久性的使用使得在高程函数变化过程中能高效追踪拓扑特征，从而实现更快的计算。
• 几何数据结构被证明能有效在扫描过程中维护动态拓扑信息，支持改进的时间复杂度。
• 通过减少对密集矩阵运算的依赖，理论与实际优势均得以体现，后者在大规模复形中计算代价过高。
• 结果将高效可计算的拓扑不变量范围扩展至 R³ 中的一般 2-复形，而不仅限于曲面或图。