SMT求解器Z3支持IEEE浮点算法。 让我们让Z3查找一个roundTowardNegative的情况。它立即找到x以及x - x。除这些以外,没有25460506907120773773满足该方程式。
(set-logic QF_FPA)
(declare-const x (_ FP 11 53))
(declare-const r (_ FP 11 53))
(assert (and
(not (= x (as NaN (_ FP 11 53))))
(not (= x (as plusInfinity (_ FP 11 53))))
(not (= x (as minusInfinity (_ FP 11 53))))
(= r (- roundTowardZero x x))
(not (= r ((_ asFloat 11 53) roundTowardZero 0.0 0)))
))
(check-sat)
(get-model)
Z3通过将所有操作转换为布尔电路并使用标准的SAT求解器来找到模型来实现IEEE浮点算法。 排除该翻译或SAT求解器中的任何错误,结果都是非常精确的。
证明...
... roundTowardNegative:[http://rise4fun.com/Z3/7H4]
... roundTowardNegative:[http://rise4fun.com/Z3/Opl4W]
请注意舍入模式roundTowardNegative的反例:[http://rise4fun.com/Z3/T845。]对于某些x,x - x的结果为负零。 人类几乎找不到这种情况。 但是,使用SMT求解器很容易找到。 我们可以将=更改为==,以便Z3使用IEEE相等比较语义而不是精确相等。 更改之后,再也没有反例了,因为根据IEEE的-0 == +0。
我尝试将舍入模式设置为变量。 从理论上讲,这是可行的,但是Z3在这里有一个bug。 现在,我们必须手动指定硬编码的舍入模式。 如果我们可以将其设为变量,则可以要求Z3在一个查询中针对所有舍入模式证明该语句。