[論文レビュー] Numerical Solution of the Optimal Transportation Problem via viscosity solutions for the Monge-Ampere equation

Jean‐David Benamou, Brittany D. Froese+1 more|arXiv (Cornell University)|Aug 23, 2012

Nonlinear Partial Differential Equations被引用数 5

ひとこと要約

本稿では、コンpakト有限差分スキームを用いて、最適輸送（OT）境界条件を伴う楕円型モンジュ・アンぺール方程式を解く、高速でニュートン法に基づく数値的手法を提示する。この手法は、ラプラス方程式を複数回解くのと同等の計算効率を達成しており、非凸領域や非連結領域を含む、特異解を有する問題に対しても安定した性能を示す。

ABSTRACT

A numerical method for the solution of the elliptic Monge-Ampere Partial Differential Equation, with boundary conditions corresponding to the Optimal Transportation (OT) problem is presented. A local representation of the OT boundary conditions is combined with a finite difference scheme for the Monge-Ampere equation. Newton's method is implemented leading to a fast solver, comparable to solving the Laplace equation on the same grid several times. Theoretical justification for the method is given by a convergence proof in the companion paper (Benamou et al., 2012). In this paper, the algorithm is modified to a simpler compact stencil implementation and details of the implementation are given. Solutions are computed with densities supported on non-convex and disconnected domains. Computational examples demonstrate robust performance on singular solutions and fast computational times.

研究の動機と目的

• 最適輸送問題に現れるモンジュ・アンぺール方程式に対する効率的な数値スケーラーの開発を目的とする。
• OT境界条件を正確に捉えることができるコンパクトな有限差分スキームの実装を目的とする。
• ニュートン法を用いることで高速収束を達成し、ラプラス方程式を複数回解くのと同等の効率性を実現することを目的とする。
• 特異解を有する非凸領域および非連結領域に対しても対応可能な手法の拡張を目的とする。
• 理論的収束保証を伴う、実装可能で実用的なアルゴリズムの提供を目的とする。

提案手法

• 最適輸送境界条件の局所的表現を、モンジュ・アンぺール方程式の有限差分離散化に統合する。
• 構造化グリッド上での精度と安定性を向上させるために、コンパクトステンシルの有限差分スキームを採用する。
• 離散化から生じる非線形連立一次方程式に対してニュートン法を適用し、高速収束を実現する。
• 密度が非凸領域や非連結領域に支持される場合に対応できるように、手法を適応させる。
• 境界条件を慎重に取り扱い、可視化解構造を保持するようにアルゴリズムを実装する。
• 同伴証明（Benamou et al., 2012）により理論的収束が裏付けられ、信頼性が保証される。

実験結果

リサーチクエスチョン

• RQ1 OT境界条件を伴うモンジュ・アンぺール方程式に対して、高速かつ安定な数値スキームを開発可能か？
• RQ2 コンパクトステンシル実装は、標準的なステンシルと比較して、精度と効率性においてどのように異なるか？
• RQ3 ニュートン法に基づくスケーラーは、特異解や複雑な領域幾何構造に対しても高い性能を維持できるか？
• RQ4 ラプラス方程式を解く場合と比較して、計算コストはどの程度か？
• RQ5 密度が非凸領域または非連結領域に支持される場合でも、手法は安定性を保つのか？

主な発見

• スケーラーは、同じグリッド上でラプラス方程式を複数回解くのと同等の計算速度を達成する。
• 特異性を有する解に対しても、困難な領域設定下でも安定した性能を示す。
• 非凸領域および非連結領域に支持される密度に対しても、解が成功裏に計算可能である。
• コンパクトステンシル実装は、標準ステンシルと比較して安定性と精度が向上している。
• 不規則な幾何構造を有するテストケースを含め、多様なケースにおいても高速収束を維持する。
• 同伴証明により理論的収束が裏付けられ、数値的手法の妥当性が検証される。