如何计算mpf_class的epsilon



我试图通过使用GMP(GNU多精度算术库)来改进c++程序。

有一个计数器变量,它将通过double类型的计算结果周期性地增加。添加发生在ISR(中断服务例程)内部,因此应该很快。

到目前为止,计数器变量的类型是double,这会导致问题。一些设备(运行多年)达到了每次添加都会被吸收的程度,因此再也不会发生任何事情。其他设备已经开始丢弃较小的添加。

正因为如此,现在是时候通过使用mpf_class而不是double来改进它了。但是为了正确设置mpf_class的大小,我需要知道;(ε)。

我的问题是:如何计算mpf_class对象的ε?

这个计算的一个例子(c++)将是很好的。

GMP中的mpf类型不提供精确的基-2表示。所请求的精度被底层数组元素中的位数四舍五入,然后添加至少一个附加元素。如果在典型的32位系统上要求精度为53,则实际值的计算值在65到96位之间。(从技术上讲,在这种情况下,mpf类型使用基数2^32算术)。

我建议您改用MPFR(带舍入的多精度浮点)库。它使用GMP库进行底层计算,但提供了更可预测的行为。GMP文件中的注释建议使用MPFR。

最新更新