[論文レビュー] A Generalization of the Persistent Laplacian to Simplicial Maps

Aziz Burak Gülen, Facundo Mémoli+2 more|arXiv (Cornell University)|Jan 1, 2023

Graph theory and applications被引用数 2

ひとこと要約

この論文は、包含写像に限らない重み付き単体写像へ、シュール補元を用いて上方向および下方向の永続ラプラシアンを定義することで、永続ラプラシアンを一般化する。永続ラプラシアンの核の次元が永続ベッチ数に等しいことを証明し、その計算に向けたアルゴリズムを提示するとともに、写像の合成における本質的固有値の単調性を確立し、包含に基づくフィルトレーションを超えたスペクトル永続性を拡張する。

ABSTRACT

The (combinatorial) graph Laplacian is a fundamental object in the analysis of, and optimization on, graphs. Via a topological view, this operator can be extended to a simplicial complex K and therefore offers a way to perform "signal processing" on p-(co)chains of K. Recently, the concept of persistent Laplacian was proposed and studied for a pair of simplicial complexes K ↪ L connected by an inclusion relation, further broadening the use of Laplace-based operators. In this paper, we significantly expand the scope of the persistent Laplacian by generalizing it to a pair of weighted simplicial complexes connected by a weight preserving simplicial map f: K → L. Such a simplicial map setting arises frequently, e.g., when relating a coarsened simplicial representation with an original representation, or the case when the two simplicial complexes are spanned by different point sets, i.e. cases in which it does not hold that K ⊂ L. However, the simplicial map setting is much more challenging than the inclusion setting since the underlying algebraic structure is much more complicated. We present a natural generalization of the persistent Laplacian to the simplicial setting. To shed insight on the structure behind it, as well as to develop an algorithm to compute it, we exploit the relationship between the persistent Laplacian and the Schur complement of a matrix. A critical step is to view the Schur complement as a functorial way of restricting a self-adjoint positive semi-definite operator to a given subspace. As a consequence of this relation, we prove that the qth persistent Betti number of the simplicial map f: K → L equals the nullity of the qth persistent Laplacian Δ_q^{K,L}. We then propose an algorithm for finding the matrix representation of Δ_q^{K,L} which in turn yields a fundamentally different algorithm for computing the qth persistent Betti number of a simplicial map. Finally, we study the persistent Laplacian on simplicial towers under weight-preserving simplicial maps and establish monotonicity results for their eigenvalues.

研究の動機と目的

• 重みを保存する単体複体間の一般化された単体写像へ、包含に基づくフィルトレーションを超えて永続ラプラシアンを拡張すること。
• K ⊆ L が不要なより広範な単体写像の文脈において、永続ラプラシアンの理論的基盤を確立すること。
• シュール補元技術を用いて、永続ラプラシアン行列表現を計算するためのアルゴリズムを開発すること。
• 任意の単体写像において、永続ラプラシアンの核の次元が永続ベッチ数に等しいことを証明すること。
• 特に上方向永続ラプラシアンの本質的スペクトルに関して、単体写像の合成における固有値の単調性を分析すること。

提案手法

• q番目の永続ラプラシアン ∆K,Lq を、fqの像へのL上の組合せ的ラプラシアンのシュール制限として定義し、上方向および下方向成分を用いる。
• カテゴリー理論とローエン順序を用いて、自己随伴非負定値作用素の部分空間へのファンクター的制限としてシュール補元を形式化する。
• fqの像に制限された∆Lq,upおよび∆Lq,downのシュール補元を計算することで、∆K,Lqの行列表現を構築する。
• 下方向永続ラプラシアンの核を用いて、上方向成分の変換行列を構築し、完全な行列再構成を可能にする。
• 必然的に生じるゼロ固有値を除いたコア作用素としての「本質的上方向永続ラプラシアン」の概念を導入する。
• 本質的上方向永続ラプラシアンのスペクトルが、重みを保存する単体写像の合成に関して単調であることを証明する。

実験結果

リサーチクエスチョン

• RQ1 包含に基づくフィルトレーションを超えて、単体複体間の任意の重みを保存する単体写像へ、永続ラプラシアンを意味的に一般化できるか？
• RQ2 包含の場合と同様に、単体写像の文脈においても、永続ラプラシアンの核の次元が永続ベッチ数に等しいか？
• RQ3 シュール補元が、この一般化された文脈における永続ラプラシアンの定義および計算にファンクター的ツールとして用いられるか？
• RQ4 単体写像の合成に関して、永続ラプラシアンの固有値は単調であるか。もしそうならば、どのような条件下で成立するか？
• RQ5 必然的ゼロ固有値は上方向永続ラプラシアンにどのような役割を果たすか。また、本質的スペクトルに注目することで、単調性を回復できるか？

主な発見

• q番目の永続ラプラシアン ∆K,Lq の核の次元は、単体写像 f:K→L におけるq番目の永続ベッチ数に等しく、包含の場合の古典的結果を一般化する。
• 永続ラプラシアンは、L上の組合せ的ラプラシアンをfqの像に制限したシュール補元を用いて計算可能であり、時間計算量 O((nKq)³ + (nLq)³ + nLq+1) の新しいアルゴリズムを提供する。
• このアルゴリズムは永続ラプラシアン行列を計算するだけでなく、同時に同じ時間計算量で永続ベッチ数も算出でき、従来の手法とは代替可能な手法を提供する。
• 上方向永続ラプラシアンには、核拡張に起因する必然的ゼロ固有値が含まれるが、これらを除去することで本質的上方向永続ラプラシアンが得られる。
• 本質的上方向永続ラプラシアンは単調性を満たす：重みを保存する単体写像の合成に関して、λK,M,essq,up,k ≥ λL,M,essq,up,k が成り立つ。
• 下方向永続ラプラシアンの単調性は全射写像において成り立ち、重みを保存する仮定のもとで合成を行うと、下方向永続ラプラシアンの固有値に対しても λK,M,q,down,k ≥ λL,M,q,down,k が成り立つ。