[論文レビュー] Computing height persistence and homology generators in R3 efficiently

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

Topological and Geometric Data Analysis被引用数 2

ひとこと要約

本稿では、R³に線形に埋め込まれた単体的2複体における高さの持続的性質とホモロジー生成子（H₁およびH₂）を計算するO(n log n)時間のアルゴリズムを提示する。これは、ジグザグ持続的性質、リーブ・グラフの性質、および効率的な幾何データ構造を活用することで、既存の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.

研究の動機と目的

• 一般の単体的2複体がR³に埋め込まれた場合に、高さの持続的性質とホモロジー生成子を効率的に計算するという未解決問題に取り組む。
• グラフや曲面に対しては既知の結果があるものの、このケースにおけるO(n²)未満または近い時間計算量のアルゴリズムの欠如を克服する。
• R³に埋め込まれた複体におけるホモロジー生成子の計算において、行列乗算より優れた時間計算量を達成する。
• 表面やグラフといった特殊ケースにとどまらず、一般の2複体に対しても効率的な位相的計算の適用範囲を拡大する。

提案手法

• 最近のジグザグ持続的性質の進展を活用し、2複体上の高さ関数の持続的性質をモデル化および計算する。
• リーブ・グラフに関する新しい幾何的観察を適用し、位相的不変量の計算を簡素化および高速化する。
• 効率的な幾何データ構造を用いて、持続的性質計算中に位相的特徴を維持および照会する。
• ホモロジー生成子の計算問題を、O(n log n)時間で処理可能な一連の持続的ホモロジー演算に還元する。
• 分割統治戦略と高さ関数に沿ったスイープライン技術を組み合わせ、動的な連結性およびサイクル情報の維持を図る。
• 高さ関数の持続的性質を、ホモロジカル構造を保持しつつ効率的な更新を可能にする一連のジグザグフィルトレーションに還元する。

実験結果

リサーチクエスチョン

• RQ1 R³に埋め込まれた2複体における高さの持続的性質は、サブ・クアドレティック時間で計算可能か？
• RQ2 H₁およびH₂の基底を、出力サイズkを用いてO(n log n + k)時間で計算することは可能か？
• RQ3 R³に埋め込まれた複体におけるホモロジー計算で、行列乗算の計算ボトルネックを回避できるか？
• RQ4 リーブ・グラフの性質をどのように活用すれば、3次元複体における位相的持続的性質計算を高速化できるか？

主な発見

• 本稿では、R³に埋め込まれた2複体における高さの持続的性質をO(n log n)時間で計算でき、行列乗算によるO(n^ω)の境界を改善した。
• H₁(K)およびH₂(K)の基底を、出力サイズkを用いてO(n log n + k)時間で計算可能であり、これは従来の手法に比べ顕著な改善である。
• ジグザグ持続的性質の活用により、高さ関数が変化する際の位相的特徴の追跡が効率的に行えるようになり、計算が高速化された。
• 幾何データ構造が、スイープ処理中に動的な位相的情報を効果的に維持できることを示した。これにより、改善された時間計算量が実現された。
• 密行列演算への依存を減らしたことで、理論的および実用的な利点が示され、大規模な複体では計算が非現実的になることが避けられた。
• 表面やグラフに限らず、一般の2複体に対しても、効率的に計算可能な位相的不変量の範囲がR³において拡大された。