Math.Pow(等等)如何实际工作

Question

所以我谷歌搜索了很长时间,我几乎找不到任何东西.我从这个url找到了一些关于可能实现Math.Pow的信息,但它们是不准确的,例如这段代码

public static double PowerA(double a, double b)
{
    int tmp = (int)(BitConverter.DoubleToInt64Bits(a) >> 32);
    int tmp2 = (int)(b * (tmp - 1072632447) + 1072632447);
    return BitConverter.Int64BitsToDouble(((long)tmp2) << 32);
}
static void Main(string[] args)
{
    double x = 12.53, y = 16.45;
    Console.WriteLine(Math.Pow(x, y));
    Console.WriteLine(PowerA(x, y));
}

提供输出:

1,15158266266297E+18
8,9966384455562E+17

如此不准确......

我当时认为它的作用类似于一系列,但我不确定.

Answer 1

通常通过以下公式评估pow:

x^y = exp2(y*log2(x))

功能exp2(x),log2(x)直接在FPU中实现.如果你想实现bignums,那么它们也可以由基本运算符使用sqrt-powers的预计算表进行评估,如:

2^1/2, 2^1/4, 2^1/8, 2^1/16, 2^1/32 ...

加快这个过程