[论文解读] Quadratic Zonotopes - An Extension of Zonotopes to Quadratic Arithmetics

Assalé Adjé, Pierre-Loïc Garoche+1 more|arXiv (Cornell University)|Nov 1, 2015

Numerical Methods and Algorithms被引用 4

一句话总结

本文提出二次zonotope，作为仿射zonotope的扩展，通过引入二次误差项（ǫiǫj）以提高浮点计算静态分析的精度。通过利用半定规划（SDP）进行区间投影，该方法在非线性函数（如乘法和超越函数）上获得比区间或仿射zonotope更紧的边界，同时保持误差项数量的二次复杂度。

ABSTRACT

Affine forms are a common way to represent convex sets of $\mathbb{R}$ using a base of error terms $\epsilon \in [-1, 1]^m$. Quadratic forms are an extension of affine forms enabling the use of quadratic error terms $\epsilon_i \epsilon_j$. In static analysis, the zonotope domain, a relational abstract domain based on affine forms has been used in a wide set of settings, e.g. set-based simulation for hybrid systems, or floating point analysis, providing relational abstraction of functions with a cost linear in the number of errors terms. In this paper, we propose a quadratic version of zonotopes. We also present a new algorithm based on semi-definite programming to project a quadratic zonotope, and therefore quadratic forms, to intervals. All presented material has been implemented and applied on representative examples.

研究动机与目标

• 将zonotopic抽象域扩展至处理二次算术，以在非线性浮点计算中获得更高精度。
• 开发一种基于半定规划（SDP）的投影算法，用于计算二次型的更紧区间上界近似。
• 在包含超越函数和迭代函数的代表性基准上，评估新域相对于标准抽象域（如区间和仿射zonotope）的性能。
• 证明二次zonotope在精度上优于仿射和区间抽象，同时保持可管理的计算复杂度。
• 探索二次zonotope作为非凸、非对称代数域的潜力，适用于涉及负误差项乘积的性质建模。

提出的方法

• 提出基于二次型的二次zonotope域，通过引入ǫiǫj项以及额外的误差符号ǫ⁺、ǫ⁻、ǫ±，扩展了仿射形式。
• 将二次型的语义定义为误差空间Cm在该形式下的像，通过在定义域上取最小/最大值计算区间边界。
• 引入一种安全的上界近似方法（基于绝对值和符号感知边界的MT方法），并提出一种更紧的替代方法，利用半定规划（SDP）求解器。
• 设计了具有误差项数量二次复杂度的算术运算（加法、取反、标量乘法和乘法）。
• 实现从区间到二次形式的反向抽象映射，支持集成到抽象解释框架中。
• 将该域集成到APRON库中，并使用SDP求解器（CSDP、Mosek）进行精确的区间投影，支持可选启用以管理性能开销。

实验结果

研究问题

• RQ1 对于乘法和超越运算等非线性函数，二次zonotope能否显著优于仿射zonotope或区间，提供更紧的区间上界近似？
• RQ2 尽管计算成本更高，使用半定规划（SDP）进行区间投影是否能带来相对于标准MT上界近似方法的实质性精度提升？
• RQ3 在包含循环和条件分支的基准测试中，二次zonotope域与T1P域（带线性约束的仿射zonotope）相比，在精度和性能上表现如何？
• RQ4 二次zonotope的几何语义是否为非凸且非对称？该性质是否能更好地建模某些程序行为，如符号相关的误差传播？
• RQ5 二次zonotope能否作为有效的二阶泰勒抽象？其如何与分支定界等精度增强技术集成？

主要发现

• 对于输入x ∈ [-1, 1]的反正切函数，二次zonotope的边界为[-1.002866, 1.002866]，显著优于区间和仿射zonotope（均为[-1.919149, 1.919149]），也优于仿射约束zonotope（[-1.349407, 1.349407]）。
• 在x ∈ [-10, 10]的同一函数上，二次zonotope达到边界[-1.597501, 1.591769]，优于区间和仿射zonotope，并接近T1P域的精度。
• 在A ∈ [16, 20]的1/√A的Householder迭代中，二次zonotope展现出更优的收敛性与更低的全局误差，优于区间、仿射zonotope，甚至优于T1P域，尤其在更高展开次数时优势明显。
• 基于SDP的投影方法提供的上界近似比MT方法更紧，尽管计算成本更高，且可选择性启用以平衡精度与性能。
• 尽管基于SDP的语义计算成本较高，该方法在实践中扩展性良好，尤其适用于多项式与非线性计算。
• 结果表明，二次zonotope在非线性运算中精度高于仿射zonotope和区间，其非凸、非对称的语义可能支持涉及负误差项乘积的复杂性质建模。