[论文解读] 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

一句话总结

本文在较弱的正则性假设下，通过最优传输对偶性和等高线公式等工具，建立了半离散最优传输关于狄拉克质量位置的二阶可微性的充分条件。证明了在背景测度满足温和正则性条件时，最优传输代价函数是二阶连续可微的，从而在点阵图（stippling）和蓝噪声生成等应用中实现了高效的二阶优化。

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.

研究动机与目标

• 建立确保半离散最优传输代价关于狄拉克质量位置二阶可微性的背景测度的最小条件。
• 将现有关于Voronoi图能量可微性的结果推广至高维空间和正则性更低的背景测度。
• 在点阵图和蓝噪声生成等数值应用中，实现二阶优化方法。
• 阐明海森矩阵在传播信息跨越Laguerre单元中的作用，特别是在具有局部依赖与全局依赖的问题中。
• 为在非光滑设置下对背景测度ν的参数微分最优传输提供理论基础。

提出的方法

• 利用Kantorovitch对偶性表述半离散最优传输问题，将对偶目标函数表示为拉格朗日乘子φ的拟凹函数。
• 采用Laguerre镶嵌框架，将最优传输单元表示为c(zi,x)−c(zj,x)的下水平集的交集，以支持几何分析。
• 应用等高线公式计算对偶目标函数关于φ和z的导数，利用镶嵌结构的特性。
• 通过分析边界测度的正则性以及在一点处相交的单元数量（计数因子#M⁻¹({x})），建立二阶可微性。
• 推导出背景测度ν满足的条件，使得传输代价的海森矩阵存在且连续。
• 通过牛顿法进行数值验证，并与Lloyd算法在点阵图和蓝噪声问题上进行比较。

实验结果

研究问题

• RQ1 在背景测度ν的何种最小条件下，半离散最优传输代价关于狄拉克质量位置是二阶可微的？
• RQ2 最优传输代价的海森矩阵如何编码Laguerre镶嵌中单元间的连通性与信息传播？
• RQ3 为何牛顿法在点阵图等局部变化传播缓慢的问题中优于Lloyd算法？
• RQ4 背景测度ν的正则性与Voronoi图能量可微性之间有何关系？
• RQ5 对于支持在曲线上的非光滑背景测度（如此类测度），二阶方法是否仍可有效应用？

主要发现

• 在背景测度ν满足较弱正则性假设（包括低于以往假设的正则性）时，半离散最优传输代价关于狄拉克质量位置是二阶连续可微的。
• 当背景测度ν具有有界密度或支持在低维流形上，且计数因子#M⁻¹({x})有界时，传输代价的海森矩阵是良好定义且可计算的。
• 数值实验表明，由于海森矩阵能瞬时传播全域信息，牛顿法在点阵图问题中收敛速度显著快于Lloyd算法，尤其在新增点时优势明显。
• 在蓝噪声问题中，牛顿法的吸引域大于点阵图问题，表明由于其对所有单元具有全局依赖性，该问题更适用于二阶优化。
• 在相同条件下，Voronoi图的能量也是二阶可微的，将先前结果推广至高维空间和正则性更低的测度。
• 该方法的性能在点阵图问题中受限于算术精度，而蓝噪声问题因具有全局特性与更快收敛性，保持稳定且可扩展。