[論文レビュー] Differentiation and regularity of semi-discrete optimal transport with respect to the parameters of the discrete measure

Frédéric de Gournay, Jonas Kahn+1 more|arXiv (Cornell University)|Mar 2, 2018

Advanced Mathematical Modeling in Engineering参考文献 27被引用数 21

ひとこと要約

本稿では、最適輸送双対性および等高線公式の道具を用いて、ディラックデルタ質量の位置に関して、半離散的最適輸送の2階微分可能性を満たす十分条件を確立する。背景測度にやや弱い正則性仮定を課すと、最適輸送コストが2回連続的に微分可能であることが証明され、ステイプリングやブルーノイズ生成などの応用分野における効率的な2階最適化を可能にする。

ABSTRACT

This paper aims at determining under which conditions the semi-discrete optimal transport is twice differentiable with respect to the parameters of the discrete measure and exhibits numerical applications. The discussion focuses on minimal conditions on the background measure to ensure differentiability. We provide numerical illustrations in stippling and blue noise problems.

研究の動機と目的

• 背景測度に関する最小限の条件を同定し、ディラック質量の位置に関して半離散的最適輸送コストが2階微分可能であることを保証すること。
• 従来のボロノイ図エネルギーの微分可能性に関する結果を、より高次元および正則性が低い背景測度へと拡張すること。
• ステイプリングやブルーノイズ生成などの数値的応用において2階最適化法を可能にすること。
• ヘッセ行列がラグレアセルの間を情報伝播する役割をどのように果たすかを明らかにすること、特に局所的依存性とグローバル依存性の問題において。
• 非滑らか状態における背景測度νのパrameterに関して最適輸送を微分する理論的基盤を提供すること。

提案手法

• カントロビッチ双対性を用いて半離散的最適輸送問題を定式化し、ラグランジュ乗数φの関数として双対目的関数を凹関数として表現する。
• ラグレアタイルセーションフレームワークを用いて、最適輸送セルをc(zi,x)−c(zj,x)の下位集合の交差として表現し、幾何的解析を可能にする。
• 等高線公式を適用して、双対目的関数のφおよびzに関する微分を計算し、タイルセーションの構造を活用する。
• 境界測度の正則性および一点xにおけるセルの交差数（#M⁻¹({x}））の解析を通じて、2階微分可能性を確立する。
• 背景測度νに関する条件を導出し、輸送コストのヘッセ行列が存在し連続であることを示す。
• ステイプリングおよびブルーノイズ問題において、ニュートン法とリドウのアルゴリズムを比較し、数値的に手法の有効性を検証する。

実験結果

リサーチクエスチョン

• RQ1 背景測度νに関して、半離散的最適輸送コストがディラック質量の位置に関して2階微分可能であるための最小限の条件は何か？
• RQ2 最適輸送コストのヘッセ行列は、ラグレアタイルセレーションにおける接続性および情報伝播をどのように表現するか？
• RQ3 ステイプリングのような問題では、局所的変更がゆっくりと伝播するが、なぜニュートン法がリドウのアルゴリズムを上回るのか？
• RQ4 背景測度の正則性とボロノイ図エネルギーの微分可能性の関係は何か？
• RQ5 曲線に台を持つような非滑らか背景測度に対しても、2階法を効果的に適用できるか？

主な発見

• 背景測度νにやや弱い正則性仮定を課すと、ディラック質量の位置に関して半離散的最適輸送コストが2回連続的に微分可能であることが保証される。これは、従来の仮定よりも低い正則性を許容する。
• 背景測度νが有界な密度を持つ、または低次元多様体に台を持つ場合、ヘッセ行列が適切に定義され計算可能である。ただし、#M⁻¹({x})が有界であることが条件となる。
• 数値実験では、ニュートン法がリドウのアルゴリズムよりもステイプリング問題において著しく高速に収束することが示され、ヘッセ行列のおかげで領域全体に即座に情報が伝播するからである。
• ブルーノイズ問題では、ステイプリング問題よりもニュートン法の吸引域が広く、すべてのセルにグローバル依存性があるため、2階最適化に適していることが示された。
• ボロノイ図エネルギーも同様の条件下で2階微分可能であり、従来の結果を高次元および正則性が低い測度へと一般化した。
• ステイプリング問題では算術精度の制限が性能を制限する一方、ブルーノイズ問題はグローバルな性質と高速収束のおかげに安定的かつスケーラブルである。