[論文レビュー] Unmixing the mixed volume computation

Tianran Chen|arXiv (Cornell University)|Mar 5, 2017

Polynomial and algebraic computation参考文献 59被引用数 7

ひとこと要約

本稿では、$ n $ 個の凸多面体の混合体積が、それらの凸包の正規化体積の $ n! $ 倍に等しくなる十分条件を確立し、複雑な混合体積の計算を単一の多面体体積問題に変換する。主な結果により、元の多面体よりも凸包の幾何構造が単純な場合に顕著な計算コスト削減が可能となり、特に多面体ホモトピー法による多項式方程式系の解法において利点をもたらす。

ABSTRACT

Computing mixed volume of convex polytopes is an important problem in computational algebraic geometry. This paper establishes sufficient conditions under which the mixed volume of several convex polytopes exactly equals the normalized volume of the convex hull of their union. Under these conditions the problem of computing mixed volume of several polytopes can be transformed into a volume computation problem for a single polytope in the same dimension. We demonstrate through problems from real world applications that substantial reduction in computational costs can be achieved via this transformation in situations where the convex hull of the union of the polytopes has less complex geometry than the original polytopes. We also discuss the important implications of this result in the polyhedral homotopy method for solving polynomial systems.

研究の動機と目的

• 複数の凸多面体の混合体積が $ n! $ 倍のそれらの凸包の体積に等しくなる十分な幾何的条件を特定すること。
• 混合体積の計算複雑度を、単一の多面体体積問題に変換することで低減すること。
• 多面体ホモトピー法を用いた多項式方程式系の解法などの応用分野における実用的性能向上を可能にすること。
• 凸包体積計算の過程で自動的に検証可能な条件を提供すること。

提案手法

• 本稿では、定理1.1を導入し、$ \operatorname{conv}(\tilde{S}) $ のすべての真の正次元の面が各 $ S_i $ と交わるならば、$ \operatorname{MV}(\operatorname{conv}S_1,\dots,\operatorname{conv}S_n) = n!\operatorname{vol}_n(\operatorname{conv}(\tilde{S})) $ が成り立つと述べる。
• 定理1.2は、より弱い条件を許容することで一般化を図る：ある $ S_i $ との面の交わり、あるいは特定の次元制約を満たす低次元部分空間への面の射影が許可される。
• 混合体積の単調性と凸包の性質を活用し、混合体積を正規化体積の計算に等価に変換する。
• これらの条件は、凸包体積計算の副産物として検証可能であり、既存のアルゴリズムへの統合が可能である。
• 本手法は、潮流方程式、 Kuramoto モデル、テンソル固有値問題といった実世界の応用事例を通じて検証されている。
• 実験的検証を支援するため、正則分割と正規化体積を効率的に計算できる C++ ソフトウェアライブラリ libtropicana が開発された。

実験結果

リサーチクエスチョン

• RQ1 どの幾何的条件下で $ n $ 個の凸多面体の混合体積が、それらの凸包の体積の $ n! $ 倍に等しくなるか？
• RQ2 多面体の凸包の幾何構造が元の多面体よりも単純な場合、混合体積の計算を単一の体積計算に還元できるか？
• RQ3 この等価性のための条件は、体積計算の過程でどのように効率的に検証できるか？
• RQ4 多項式方程式系の解法などの実用的応用において、この変換がもたらす計算上の利点は何か？
• RQ5 この手法は、より一般的な代数幾何の文脈におけるニュートン・オクヌコフ体に対しても拡張可能か？

主な発見

• 定理1.1により、凸包のすべての正次元の面がすべての $ S_i $ と交わる場合、$ Q_1, \dots, Q_n $ の混合体積は $ n! \cdot \operatorname{vol}_n(\operatorname{conv}(Q_1 \cup \cdots \cup Q_n)) $ に等しくなる。
• 定理1.2におけるより一般的な条件では、一部の面が特定の $ S_i $ と交わらない場合でも、特定の交わりまたは射影制約を満たしていれば等価性が保たれる。
• $ \mathbb{R}^2 $ 内の2つの多角形の例では、混合体積は正確に $ 2 \cdot \operatorname{vol}_2(\operatorname{conv}(Q_1 \cup Q_2)) = 32 $ に等しくなり、簡略化による32単位の体積的利得が確認された。
• 凸包の頂点数が少なく、または構造が単純な場合、計算コストの低減が達成される。
• 本手法は半混合系にも適用可能であり、PHoM や Hom4PS-2.0 などの多面体ホモトピー・ソルバーへの統合が可能である。
• ソフトウェア libtropicana により、正則分割と正規化体積の効率的計算が可能となり、理論的結果の実験的検証が支援された。