在浮点中实现pow()函数的最有效方法

Question

我试图实现我自己的pow()和sqrt()函数版本,因为我的自定义库没有pow()/ sqrt()浮点支持.

有人可以帮忙吗？

Answer 1

是的,Sun可以(我现在认为是Oracle):

fdlibm,"可自由分发的数学库",具有sqrt和pow,以及许多其他数学函数.

然而,它们是相当高科技的实现,当然,没有什么是这样的"最有效"的实现.你的源代码把它完成之后,或者是你真的没有那么多寻找pow和sqrt,但实际上寻找浮点算法编程的教育？

Answer 2

当然 - 如果你有指数和自然的日志功能,这很容易.

既然如此y = x^n,你可以采取双方的自然对数:

ln(y) = n*ln(x)

然后采取双方的指数为您提供您想要的:

y = exp(n*ln(x))

如果你想要更好的东西,我所知道的最好的地方是Abramowitz和Stegun.