[論文レビュー] A viscosity framework for computing Pogorelov solutions of the Monge-Ampere equation

Jean‐David Benamou, Brittany D. Froese|arXiv (Cornell University)|Jul 4, 2014

Geometric Analysis and Curvature Flows参考文献 15被引用数 5

ひとこと要約

本稿は、原子測度を伴う最適輸送問題におけるMonge-Ampère方程式のPogorelov解を計算するための粘性枠組みを導入する。ハイブリッド粘性-Aleksandrov形式化と片側方向微分を用いることで、一貫性があり単調な有限差分離散化を可能にした。従来の粘性ソルバーが失敗するDirac測度の状況でも、安定した数値収束を示した。

ABSTRACT

We consider the Monge-Kantorovich optimal transportation problem between two measures, one of which is a weighted sum of Diracs. This problem is traditionally solved using expensive geometric methods. It can also be reformulated as an elliptic partial differential equation known as the Monge-Ampere equation. However, existing numerical methods for this non-linear PDE require the measures to have finite density. We introduce a new formulation that couples the viscosity and Aleksandrov solution definitions and show that it is equivalent to the original problem. Moreover, we describe a local reformulation of the subgradient measure at the Diracs, which makes use of one-sided directional derivatives. This leads to a consistent, monotone discretisation of the equation. Computational results demonstrate the correctness of this scheme when methods designed for conventional viscosity solutions fail.

研究の動機と目的

• 測度の一方が重み付きDirac質量の和である場合、密度が存在しないため標準的な粘性法が失敗するという数値的課題に対処する。
• 粘性解とAleksandrov（一般化された）解の間のギャップを埋めるために、一貫性があり単調なままのハイブリッド粘性-Aleksandrov形式化を導入する。
• Dirac点における部分勾配測度の局所的特徴付けを、片側方向微分に基づいて行い、Monge-Ampère作用素の正確な離散化を可能にする。
• 一貫性と単調性を両立する有限差分スキームを構築し、Pogorelov意味での一般化解への収束を保証する。

提案手法

• 本稿は、原子測度を扱えるように、粘性解理論とAleksandrovの一般化解概念を組み合わせたハイブリッド粘性-Aleksandrov形式化を提案する。
• 片側方向微分を用いたDirac点における部分勾配測度の局所的特徴付けを導入し、半離散的最適輸送問題における細胞体積の計算を可能にする。
• 非一様グリッドと広いスティンシルを用いた中央差分による2階方向微分の近似により、一貫性があり単調な有限差分離散化を構築する。
• 境界条件のHamilton-Jacobi方程式は、領域周辺の層で上流スキームにより近似され、符号付き距離関数はLegendre-Fenchel変換を用いて表現される。
• 平均ゼロ制約は、一貫性と単調性を維持する補正項を組み込んだ修正された離散Monge-Ampère作用素により強制される。
• 非定数密度の場合は、極座標における密度の原始関数を事前に計算し、それらを用いてDirac点における局所的測度を計算することで、拡張される。

実験結果

リサーチクエスチョン

• RQ1 測度の一方がDirac質量の和である場合、密度が存在しないため標準的な粘性法が失敗するが、その状況でも粘性に基づく数値枠組みを構築可能か？
• RQ2 方向微分を用いた局所的特徴付けにより、Dirac点における部分勾配測度をどのように特徴付けることができるか？これにより、正確な離散化が可能になるか？
• RQ3 一般化（Pogorelov）解のMonge-Ampère方程式に対して、一貫性と単調性を両立する有限差分スキームを設計可能か？
• RQ4 従来の手法（例：Pogorelovのアルゴリズム）と比較して、数値的に精度と収束性の点で本手法はどのように異なるか？
• RQ5 本枠組みは、非一様かつ非定数密度を伴う源測度へも拡張可能か？

主な発見

• 提案された粘性-Aleksandrov形式化は、Pogorelovの意味で元の最適輸送問題と数学的に同等であり、解の正しさを保証する。
• 有限差分スキームは一貫性があり単調であるため、適切な条件下で一意の一般化解への収束が保証される。
• 数値的結果から、従来の粘性ソルバーが失敗する状況でも、本手法はPogorelov解を正しく計算できることを示した。特にDirac測度が存在する場合に顕著である。
• Dirac点における片側方向微分の使用により、細胞の幾何的再構築を必要とせず、輸送細胞の体積を正確に計算可能である。
• 極座標における密度の原始関数を事前に計算することで、非定数密度に対しても本手法は安定に保たれる。
• 一般密度および背景測度への拡張は、混合Aleksandrov-粘性解枠組みにより達成され、一貫性と単調性が維持されたまま離散化が保たれる。