假设我们有一组浮点数,尾数为“m”,指数为“e”。假设我们想要逼近一个函数“f”。
从理论上我们知道,通常使用“范围缩减函数”,然后从这样的函数中我们得出全局函数值。
例如 let x = (sx,ex,mx) (sign exp and mantissa) then... log2(x) = ex + log2(1.mx) 所以基本上范围缩小函数是“log2(1.mx)” .
目前我已经实现了倒数、平方根、log2 和 exp2,最近我开始使用三角函数。但是,如果给定一个全局误差界限(尤其是 ulp 误差),我可能会为范围缩减函数推导出一个误差界限,是否有一些关于此类问题的研究?说到 log2(x)(例如),我想说...
“好的,我想要带有 k ulp 错误的 log2(x),为了实现这一点,考虑到我们的浮点系统,我们需要用 p ulp 错误来近似 log2(1.mx)”
请记住,正如我所说,我们知道我们正在使用浮点数,但格式是通用的,所以它可能是经典的 F32,但即使例如 e=10,m = 8 end 等等。
我实际上找不到任何显示此类研究的参考资料。我拥有的参考资料(即穆勒书)没有以这种方式处理该主题,因此我正在寻找某种论文或类似的东西。你知道有什么参考吗?
我也试图自己推导出这样的约束,但这并不容易......