[論文レビュー] A JKO splitting scheme for Kantorovich-Fisher-Rao gradient flows

Thomas Gallouët, Léonard Monsaingeon|arXiv (Cornell University)|Feb 14, 2016

Geometry and complex manifolds参考文献 32被引用数 23

ひとこと要約

本稿では、ケンタロビッチ＝ファイラー＝ラオ（KFR）距離における勾配流のための新規なJ.K.O.型分割スキームを導入し、流れを逐次的なワサービングト（拡散）およびファイラー＝ラオ（反応）ステップに分解する。この手法は、反応-対流-拡散PDEの弱解への収束を証明し、KFR距離のMKおよびFR成分への直交分解を活用してエネルギー散逸不等式（EDI）を確立する。

ABSTRACT

In this article we set up a splitting variant of the JKO scheme in order to handle gradient flows with respect to the Kantorovich-Fisher-Rao metric, recently introduced and defined on the space of positive Radon measure with varying masses. We perform successively a time step for the quadratic Wasserstein/Monge-Kantorovich distance, and then for the Hellinger/Fisher-Rao distance. Exploiting some inf-convolution structure of the metric we show convergence of the whole process for the standard class of energy functionals under suitable compactness assumptions, and investigate in details the case of internal energies. The interest is double: On the one hand we prove existence of weak solutions for a certain class of reaction-advection-diffusion equations, and on the other hand this process is constructive and well adapted to available numerical solvers.

研究の動機と目的

• 測度における質量の変動を許容するケンタロビッチ＝ファイラー＝ラオ（KFR）距離における勾配流のための構成的かつ数値的に取り扱いやすいスキームの開発を目的とする。
• KFR距離におけるジョルダン＝キンダーレーラー＝オットー（JKO）スキームの時間分割変種の収束を、ワサービングト（MK）およびファイラー＝ラオ（FR）成分に分離することで確立すること。
• 基礎となる距離空間に最小限の幾何的仮定しか課さない条件下で、反応-対流-拡散PDEの弱解の存在を、提案された分割スキームを用いて証明すること。
• KFR距離の全流れにおける厳密なエネルギー散逸不等式（EDI）を回復するため、連続的なMKおよびFRステップによるエネルギー減少の正確な近似が可能であることを示すこと。
• 古典的なJKO推定（エネルギー単調性、質量制御、BVバインディング）が分割フレームワーク下でも維持されることを示し、安定性および収束性を保証すること。

提案手法

• 時間分割JKOスキームを提案し、連続する時間ステップでワサービングト（MK）距離とファイラー＝ラオ（FR）距離の最小化を交互に実行することで、全KFR勾配流を近似する。
• KFR距離の形式的リーマン構造を活用し、その無限小ノルムが直交和として分解されることを示す：$\|\operatorname{grad}_{\mathtt{KFR}}\mathcal{F}(\rho)\|^{2} = \|\operatorname{grad}_{\mathtt{MK}}\mathcal{F}(\rho)\|^{2} + \|\operatorname{grad}_{\mathtt{FR}}\mathcal{F}(\rho)\|^{2}$、これにより分割の正当化がなされる。
• KFR距離のインフィミウム畳み込み構造を用いて、逐次的最小化プロセスを厳密に正当化し、離散エネルギー推定を導出する。
• 標準JKOスキーム（MK用）の既知の結果を適用し、変数変換$s = \sqrt{\rho}$を用いてFRステップを凸ヒルバート問題に再定式化する。
• 密度列に対して一様な$L^1 \cap L^\infty$バインディングを確立し、密度の強$\mathrm{L}^1$収束および勾配項の弱収束を証明する。
• 両ステップからの推定を組み合わせ、離散エネルギー散逸不等式（EDI）を導出し、$\tau \to 0$の極限でそれが成立することを示す。

実験結果

リサーチクエスチョン

• RQ1 ワサービングトおよびファイラー＝ラオ成分に分離することで、KFR勾配流のための分割スキームを構築可能であり、弱解への収束を保つことができるか？
• RQ2 距離の形式的直交分解が成り立つ場合、MKおよびFR距離への逐次的最小化が、全KFR勾長流を適切に近似できるか？
• RQ3 古典的なJKO推定（エネルギー単調性、全二乗距離バインディング、質量制御）が分割フレームワーク下でも保持されるか？
• RQ4 エネルギー汎関数$\mathcal{F}$にどのような条件下で、離散EDIが得られ、それが連続版に収束するか？
• RQ5 提案されたスキームは数値的に実行可能であり、特にMKステップのモンジュ＝アンペールソルバーおよびFRステップの凸最適化ソルバーと互換性があるか？

主な発見

• 適切なコンパクト性および凸性仮定の下で、提案された分割スキームは反応-対流-拡散PDE $\partial_t\rho = \operatorname{div}(\rho \nabla(U'(\rho) + \Psi + K*\rho)) - \rho(U'(\rho) + \Psi + K*\rho)$ の弱解に収束することが示された。
• エネルギー散逸不等式 $\mathcal{F}(\rho(t_2)) + \int_{t_1}^{t_2} \left( \int_\Omega |\nabla U'(\rho)|^2 \, d\rho \, dt + \int_\Omega |U'(\rho)|^2 \, d\rho \, dt \right) \leq \mathcal{F}(\rho(t_1))$ が $\tau \to 0$ の極限で回復され、距離勾配流構造が確認された。
• 各時間ステップで離散エネルギー推定 $\mathcal{F}(\rho^{n+1}) + \tau \left( \int_\Omega |\nabla U'(\rho^{n+1/2})|^2 \, d\rho^{n+1/2} + \int_\Omega |U'(\rho^{n+1})|^2 \, d\rho^{n+1} \right) \leq \mathcal{F}(\rho^n)$ が成立し、安定性が保証された。
• 密度 $\rho^\tau(t) \to \rho(t)$ の強$\mathrm{L}^1$収束および $\tilde{\rho}^\tau \nabla U'(\tilde{\rho}^\tau) \rightharpoonup \rho \nabla U'(\rho)$ の弱収束が確立され、エネルギー不等式における極限操作が可能となった。
• この手法は構成的であり、既存の数値ソルバーと互換性がある：MKステップでは標準的なモントゥ＝アンペールソルバーが使用可能で、FRステップは $s = \sqrt{\rho}$ 変換により凸最適化問題に還元される。
• 構造的条件 $\rho U''(\rho) + U'(\rho)/2 \geq 0$ が、FR距離における $U$ の測地的凸性と同値であることが示され、MK-凸性と併せて全KFR-EDIを保証する。