[論文レビュー] Efficient Approximation of the Matching Distance for 2-parameter persistence

Michael Kerber, Arnur Nigmetov|arXiv (Cornell University)|Dec 12, 2019

Topological and Geometric Data Analysis被引用数 8

ひとこと要約

本稿では、2パラメータの恒続的ホモロジーにおけるマッチング距離の近似のための幾何学的な再解釈と効率的なアルゴリズムを提示する。四分木細分化フレームワーク内での適応的局所変動バウンドを用いることで、グローバルな変動バウンドに代えて、よりタイトな長方形固有のバウンドを導入し、従来手法に比べ最大15倍の高速化を達成した。HERAライブラリに公開された実装により、このタスクに対する最初の効率的なオープンソースツールが提供された。

ABSTRACT

The matching distance is a computationally tractable topological measure to compare multi-filtered simplicial complexes. We design efficient algorithms for approximating the matching distance of two bi-filtered complexes to any desired precision $\\epsilon>0$. Our approach is based on a quad-tree refinement strategy introduced by Biasotti et al., but we recast their approach entirely in geometric terms. This point of view leads to several novel observations resulting in a practically faster algorithm. We demonstrate this speed-up by experimental comparison and provide our code in a public repository which provides the first efficient publicly available implementation of the matching distance.

研究の動機と目的

• 2パラメータの恒続的ホモロジーにおけるマッチング距離を、マルチフィルター・データを比較するための重要な位相的測度として、高速で実用的なアルゴリズムで近似すること。
• 正確な計算の計算不能性を克服し、高精度で効率的な近似手法を提供すること。
• バイアソッティらの四分木ベースの近似アルゴリズムを、幾何的再解釈と適応的バウンドを用いて改善すること。
• 広範なトポロジカル・データ解析への採用を可能にする、公開可能な高パフォーマンス実装を提供すること。

提案手法

• マッチング距離近似問題を完全に幾何的用語に再定式化し、2次元パラメータ空間内の有界長方形 R でパラメータ化された直線を定義する。
• R 内の各点 p に対して、p でパラメータ化される直線上に制限されたパーシステンス図形のボトルネック距離 f(p) を割り当て、マッチング距離は sup_{p∈R} f(p) として定義される。
• 従来の研究で用いられたグローバルバウンドに代えて、各部分長方形 S 内での f(p) 変動を推定する適応的局所変動バウンドを導入する。
• BFS順序での四分木分割と分枝限定法を用い、ある長方形の上界が現在の下界に対して十分に小さい場合に停止する。
• 優先度キューに基づく細分化戦略を実装し、次に細分化する長方形を、最大の上界を持つものから動的に選択する。
• 各長方形ごとによりタイトな局所バウンドを計算することでパフォーマンスを最適化し、不要な細分化を減らし、収束を加速する。

実験結果

リサーチクエスチョン

• RQ1 マッチング距離の2パラメータのパーシステントホモロジーに対する近似を、アルゴリズムを幾何学的に再解釈することで、より効率的に行えるか？
• RQ2 適応的局所変動バウンドは、グローバルバウンドに比べ、必要な細分化回数をどの程度削減できるか？
• RQ3 グローバルバウンドを1つだけ用いるのではなく、局所的に導出されたバウンドを用いることで、実際のパフォーマンスにどの程度の向上が得られるか？
• RQ4 既存の手法を上回る性能を発揮する、実用的で公開可能なマッチング距離実装を開発できるか？
• RQ5 このアルゴリズムは実世界の2フィルトレーションにおいてどの程度スケーリング可能か？また、正確な計算と比較してどうなるか？

主な発見

• 提案手法は、入力のバイフィルトレーションに応じて、グローバルバウンド法に比べ3〜8倍（一部のケースでは最大15倍）の高速化を達成した。
• 適応的局所変動バウンドの使用により、各長方形における上界が顕著に小さくなり、必要な細分化回数が削減された。
• HERAライブラリへの実装により、2パラメータの恒続的ホモロジーにおけるマッチング距離を計算する最初の公開可能で効率的なツールが提供された。
• 実験結果から、相対誤差が時間経過とともに滑らかに減少し、予測可能な時間枠内で 1+ε 近似を保証することがわかった。
• 局所線形バウンドのバージョンは、1長方形あたりの計算コストが高めだが、並列処理が容易であり、大規模なインスタンスでは他の手法を上回ると予想される。