數值權衡（-0.5 *的std ::日誌（X））

Question

我所遇到的一些舊代碼，計算

double y = 1/std::sqrt(x); 

使用：

constexpr double base16 = 16.0; 
double log_base16 = std::log(base16); 
double y = std::pow(base16, -0.5 * std::log(x)/log_base16); 

基本上是：

double y = std::exp(-0.5 * std::log(x)); 

有什麼理由對於數值的好處（如準確性或更有可能避免下溢/過流程）之間的方法？原作者可能會這樣想。

Answer 1

原代碼被認爲是非常調皮的確，特別是在現代C++標準和IEEE754浮點：由IEEE標準要求

的std ::開方是準確的。 [原文如此]

此外，std::pow具有沒有這樣的要求。

因此我會試圖將它改寫爲1/std::sqrt(x)，當然是測試。

參考：http://en.cppreference.com/w/cpp/numeric/math/sqrt

Answer 2

我看不到實現的sqrt（）戰俘而言什麼特別好的原因（）和log（）。這可能是sqrt（）的實現中存在一個bug，這被用作解決方法。

一般來說，我預計這個實現會更慢，更不精確。