[論文レビュー] Chain complex reduction via fast digraph traversal

Leon Lampret|arXiv (Cornell University)|Mar 2, 2019

Topological and Geometric Data Analysis参考文献 73被引用数 4

ひとこと要約

この論文は、代数的モース理論を活用してコストの高い行列演算を回避する、高速な有向グラフ走査に基づくチェーン複体の簡約化のための新規アルゴリズムを導入する。双方向有向グラフ表現上の非巡回マッチングを計算することで、満たしの増加とメモリ使用量を著しく削減し、500万×500万のスパース行列ですら効率的なホモロジー計算を可能にする。ホモトピー型を保持するとともに、明示的なホモロジー生成子と誘導写像を返す。

ABSTRACT

Reducing a chain complex (whilst preserving its homotopy-type) using algebraic Morse theory gives the same end-result as Gaussian elimination, but AMT does it only on certain rows/columns and with several pivots (in all matrices simultaneously). Crucially, instead of doing costly row/column operations on a sparse matrix, it computes traversals of a bipartite digraph. This significantly reduces the running time and memory load (smaller fill-in and coefficient growth of the matrices). However, computing with AMT requires the construction of a valid set of pivots (called a Morse matching). We discover a family of Morse matchings on any chain complex of free modules of finite rank. We show that every acyclic matching is a subset of some member of our family, so all maximal Morse matchings are of this type. Both the input and output of AMT are chain complexes, so the procedure can be used iteratively. When working over a field or a local PID, this process ends in a chain complex with zero matrices, which produces homology. However, even over more general rings, the process often reveals homology, or at least reduces the complex so much that other algorithms can finish the job. Moreover, it also returns homotopy equivalences to the reduced complexes, which reveal the generators of the homology and the induced maps $H_\ast(φ)$. We design a new algorithm for reducing a chain complex and implement it. We test that it outperforms other CASs. As a special case, given a sparse matrix over any field, the algorithm offers a new way of computing the rank and a sparse basis of the kernel (or null space), cokernel (or quotient space, or complementary subspace), image, preimage, sum and intersection subspace. It outperforms built-in algorithms in other CASs.

研究の動機と目的

• ホモトピー型を保持しつつ、メモリと時間の両面で効率的なチェーン複体の簡約化手法を開発すること。
• スパース行列に対する従来の行・列操作の計算ボトルネックを、グラフベースの走査に置き換えること。
• 一般の環上でも、ホモロジー、ホモロジー生成子、および誘導写像 $ H_*(\varphi) $ を計算できる実用的なアルゴリズムを提供すること。
• 特に体およびPID上での反復的簡約を可能とし、ゼロ行列のみを含む複体に到達するようにすること。
• Mathematicaにおけるアルゴリズムの実装とベンチマークを実施し、既存のコンピュータ代数システム（CAS）を上回ることを示すこと。

提案手法

• チェーン複体から双方向有向グラフを構築し、非ゼロ成分を出発モジュールと到達モジュールの間の有向辺として表現する。
• この有向グラフの高速走査により、非巡回性と最大簡約性を保証するモースマッチング（ピボットの集合）を計算する。
• 代数的モース理論を用いて、スパース行列に対する高価なガウスの消去法を避ける形で、行と列を同時に削除する。
• すべての最大非巡回マッチングを含む、普遍的なモースマッチングの族を用いることで、完全性と最適性を保証する。
• 反復的な簡約プロセスはホモトピー同値を返し、ホモロジー生成子と誘導写像を抽出するのに必要な構造を保持する。
• Mathematicaにおける実装では、最適化されたデータ構造とグラフ走査を用いて、満たしの増加と係数の成長を最小限に抑える。

実験結果

リサーチクエスチョン

• RQ1 双方向有向グラフ表現上のグラフ走査に置き換えることで、チェーン複体の簡約が高速化可能か？
• RQ2 任意の有限ランク自由加群のチェーン複体に対して、すべての最大非巡回マッチングを含む普遍的なモースマッチングの族が存在するか？
• RQ3 この手法は、大規模スパースチェーン複体において、従来のCASアルゴリズムと比較して、著しくメモリ使用量と実行時間を削減可能か？
• RQ4 標準的手法が係数の成長により失敗する場合でも、このアルゴリズムが一般の環、特に $\mathbb{Z}$ 上でホモロジーを計算できる範囲はどの程度か？
• RQ5 行列が簡約後に密行列化された場合でも、このアルゴリズムは、トロイティーや誘導写像を含む完全なホモロジー情報を、簡約された複体から回復可能か？

主な発見

• 8,3-Kneserグラフの独立性複体から得られる500万×500万の行列を、10分未満で最小限のメモリで簡約し、SageMathを上回る性能を示した。
• 対称半群代数 $ A_3 $ のホッフシュイルド複体に関して、$ \mathbb{Z}_2 $ 上でホモロジー $ H_* $ を34分、14GBのRAMで計算したが、SageMathは1日後にクラッシュした。
• 外積スターリング＝ライスナー代数 $ A_3 $ の場合、$ \mathbb{Z}_2 $ 上でホモロジー $ H_* $ を2時間43分、46GBのRAMで計算したが、SageMathは4日後に停止した。
• 有理数体 $ \mathbb{Q} $ 上では、ヘイゼンベルクリー代数 $ \mathfrak{g}_{13} $ のホモロジー $ H_* $ を1時間47分、18GBのRAMで計算したが、SageMathは8日後にクラッシュした。
• リーリー代数 $ \mathfrak{g}_5 $ のチェヴァリの複体において、最大の行列を520万×520万から非ゼロ要素が11,760個のスパース行列に簡約し、すべての絶対値が28以下に制限された。
• テストされたすべての複体に対して、$ \mathbb{Z} $ 上でトロイティーや生成子を含む完全なホモロジー情報を返し、標準的なCASが失敗する状況でも $ H_* $ の計算が可能であることを示した。