[论文解读] Perturbed affine arithmetic for invariant computation in numerical program analysis

Éric Goubault, Sylvie Putot|ArXiv.org|Jul 18, 2008

Numerical Methods and Algorithms参考文献 18被引用 12

一句话总结

本论文提出扰动仿射算术作为一种新颖的抽象域，用于数值程序分析，结合了序理论与拓扑结构（具体为有序巴拿赫空间），以实现收敛的不动点迭代。它证明了任意阶稳定线性递归滤波器的有限精确边界，并通过改进的并集算子保留变量间依赖关系，为非线性情况提供了准确性的有力证据。

ABSTRACT

We completely describe a new domain for abstract interpretation of numerical programs. Fixpoint iteration in this domain is proved to converge to finite precise invariants for (at least) the class of stable linear recursive filters of any order. Good evidence shows it behaves well also for some non-linear schemes. The result, and the structure of the domain, rely on an interesting interplay between order and topology.

研究动机与目标

• 开发一种稳健的抽象域，用于计算涉及浮点数与实数运算的数值程序中的精确不变量。
• 解决在安全关键系统（如航空与控制系统）中高效约束变量值的挑战。
• 通过引入兼具序结构与拓扑结构的抽象域，克服标准区间算术与仿射算术的局限性，以实现更好的收敛性与精度。
• 证明在所提出的抽象域中，不动点迭代可收敛至任意阶稳定线性递归滤波器的有限精确边界。

提出的方法

• 本文构建了一个广义仿射形式抽象域，使其构成一个有序巴拿赫空间，从而支持序理论与拓扑推理的结合。
• 提出一种扰动并集算子，将不动点迭代的效果建模为对仿射形式的扰动，从而保留变量间的依赖关系。
• 该抽象域被证明是‘几乎’有界完全格，具有有界的下确界与上确界，利用了序结构与范数拓扑之间的相互作用。
• 开发了一种克莱尼迭代技术以计算不动点，并证明了其在稳定线性递归滤波器中的收敛性。
• 该方法通过误差项（εi）的符号表示不确定性，并应用平移与并集操作在程序状态间传播边界。
• 通过允许相关变量间的联合并集并保持仿射关系（如 y = 2x - 2），扩展了标准仿射算术，从而在循环与条件语句中提升精度。

实验结果

研究问题

• RQ1 能否通过结合序结构与拓扑结构的数值抽象域，确保任意阶线性递归滤波器的不动点迭代收敛？
• RQ2 在抽象解释过程中，如何保留程序变量间的依赖关系以提升精度？
• RQ3 扰动仿射算术在非线性方案中相较于标准区间算术与仿射算术能有多大程度的性能提升？
• RQ4 该抽象域的结构能否支持更高阶近似（如泰勒模型）以提升精度？
• RQ5 不动点迭代的收敛行为是否可推广至非线性系统或其他赋范空间？

主要发现

• 所提出的抽象域确保了所有稳定线性递归滤波器（任意阶）的不动点迭代收敛至有限精确边界。
• 该抽象域为有序巴拿赫空间，且‘几乎’是Riesz空间，其有界完全性通过拓扑与序理论对偶性得以建立。
• 并集算子被重新定义为对仿射形式的扰动，从而实现依赖关系的精确传播，并优于朴素并集。
• 以线性递归滤波器为例，该方法在三次迭代后计算出紧致边界：γ(ˆxn³) = [-0.32, 0.8]。
• 在条件程序中，该方法能保持仿射关系（如 y = 2x - 2）在不同分支间的传递，实现精确的联合抽象。
• 通过将误差项限制为离散值（如 ε₂ = ±1），该框架支持析取分析，实现在条件分支中的精确推理。