[论文解读] 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 \cdot \operatorname{vol}_2(\operatorname{conv}(Q_1 \cup Q_2)) = 32 $，表明简化带来了 32 单位的体积增益。
• 当凸包的顶点数更少或结构更简单时，该变换可显著降低计算成本。
• 该方法适用于半混合系统，并可集成到 PHoM 和 Hom4PS-2.0 等多面体同伦求解器中。
• 软件 libtropicana 支持高效计算正规剖分和归一化体积，为理论结果的实验验证提供了支持。