[論文レビュー] Variational Principles for Minkowski Type Problems, Discrete Optimal Transport, and Discrete Monge-Ampere Equations

Xianfeng Gu, Feng Luo+2 more|arXiv (Cornell University)|Feb 22, 2013

Geometric Analysis and Curvature Flows参考文献 7被引用数 33

ひとこと要約

本稿は、支持高さ上の凸最適化問題として解を定式化することにより、離散的最適輸送、ミンコフスキー型問題、離散的モンジュ＝アンペール方程式の有限次元変分原理を確立する。アレクサンドロフの定理の変分的証明を、臨界点が平行移動を除いて一意な解を与える強い凸エネルギー関数を用いて行い、ヘッシアン行列に対するニュートン法を用いることで効率的な計算が可能である。

ABSTRACT

In this paper, we develop several related finite dimensional variational principles for discrete optimal transport (DOT), Minkowski type problems for convex polytopes and discrete Monge-Ampere equation (DMAE). A link between the discrete optimal transport, discrete Monge-Ampere equation and the power diagram in computational geometry is established.

研究の動機と目的

• 所与の法線ベクトルと面積をもつ凸多面体の存在および一意性に関するアレクサンドロフの定理の有限次元的変分的証明を提供すること。
• 計算幾何学における離散的最適輸送、離散的モンジュ＝アンペール方程式、およびパワー図の間の変分的枠組みを確立すること。
• 凸最適化とニュートン法を用いて、有限な像を持つアレクサンドロフ写像を効率的に計算するアルゴリズムを開発すること。
• エネルギー関数の強い凸性を用いて、離散的ヘッシアン写像の無限小剛性を証明すること。

提案手法

• 離散的最適輸送問題の解を、支持高さ空間上でのエネルギー関数の臨界点として定式化する。
• エネルギー関数を、上位集合の体積と目的面積との差分の積分として定義し、高さ変数に関する線形項を補正する。
• 体積勾配に関連する微分形式の閉じ性を用いて、適切に定義されたポテンシャル関数の存在を保証する。
• エネルギー関数のヘッシアン行列を導出し、それが正定値かつ対角優勢であることを示し、強い凸性を保証する。
• エネルギー関数に対してニュートン法を適用し、最適高さベクトルの効率的数値計算を実現する。
• 区分的線形凸関数とそのレジェンドル変換の双対性を活用し、離散的ヘッシアンをボロノイ型細胞の体積として解釈する。

実験結果

リサーチクエスチョン

• RQ1 所与の面積と法線をもつ凸多面体の存在および一意性に関するアレクサンドロフの定理に対して、変分的証明を構築できるか？
• RQ2 離散的最適輸送を、幾何的意味を持つ有限次元凸最適化問題としてどのように定式化できるか？
• RQ3 離散的モンジュ＝アンペール方程式の文脈において、エネルギー関数のヘッシアン行列の幾何的解釈は何か？
• RQ4 高さベクトルから目的体積への勾配写像は局所微分同相写像か？また、どのような条件下で全単射となるか？
• RQ5 強い凸エネルギー関数に対するニュートン法を用いて、離散的モンジュ＝アンペール方程式の解を効率的に計算できるか？

主な発見

• エネルギー関数 E(h) はアフィン部分空間 H₀ 上で強く凸であるため、平行移動を除いて一意な最大化点を持つ。
• E(h) のヘッシアン行列は正定値であり、これによりニュートン法の局所収束性が保証され、勾配写像が局所微分同相写像であることが示される。
• 解ベクトル h は E(h) の一意な最大化点であり、対応する区分的線形関数 u_h は二次コストを最小化する離散的最適輸送写像を実現する。
• 勾配写像 ∇E: H₀ → A は、目的体積の正の象限への全単射であるため、解の存在および一意性が保証される。
• 各点 p_i における双対関数の離散的ヘッシアンは、対応する細胞 W_i(h) の体積に等しく、離散的モンジュ＝アンペール方程式の幾何的解釈が得られる。
• 本手法により、強い凸性のおかげで収束が保証されるため、ニュートン法を用いたアレクサンドロフ写像の効率的数値計算が可能となる。