Pow实现双倍

我正在开发一个将用于运动控制的代码,我遇到了powfunction的问题。 我使用VS2010作为IDE。

这是我的问题:我有:

double p = 100.0000; double d = 1000.0000; t1 = pow((p/(8.0000*d),1.00/4.000); 

在评估最后一个函数时,我得不到更好的近似值。 我得到一个正确的7位十进制数字,结果数字都是垃圾。 我猜测pow函数只将任何输入变量转换为float和过程与计算。

  1. 我对吗?
  2. 如果是这样,是否有任何代码我可以“重新启动”以获得更好的精度?

编辑:解决了。

毕竟,我遇到了由OGRE 3D框架正在使用的Direct3D引起的FPU配置位问题。

如果使用OGRE,在配置GUI上,只需设置“浮点模式=一致”。

如果使用原始Direct3D,在调用CreateDevice时,请确保将“D3DCREATE_FPU_PRESERVE”标志传递给它。

原帖:

您可能正在使用将FPU的默认精度更改为单精度的libray。 然后,所有浮点运算,即使是双精度运算,也将实际作为单精度运算执行。

作为测试,您可以尝试调用_controlfp(_CW_DEFAULT,0xfffff); (您需要包括)在执行计算之前,看看是否得到了正确的结果。 这会将浮点控制字重置为默认值。 请注意,它也会重置其他设置,这可能会导致问题。

一个改变浮点精度的公共库是Direct3D 9(也可能是其他版本):默认情况下,它在创建设备时将FPU更改为单精度。 如果使用它,请在创建设备时指定标志D3DCREATE_FPU_PRESERVE以防止其更改FPU精度。

您可能正在使用将FPU的默认精度更改为单精度的libray。 然后,所有浮点运算,即使在double s上,实际上也将作为单精度运算执行。

作为测试,您可以尝试调用_controlfp( _CW_DEFAULT, 0xfffff ); (您需要包含 )才能执行计算以查看是否获得了正确的结果。 这会将浮点控制字重置为默认值。 请注意,它也会重置其他设置,这可能会导致问题。

一个改变浮点精度的公共库是Direct3D 9(也可能是其他版本):默认情况下,它在创建设备时将FPU更改为单精度。 如果使用它, D3DCREATE_FPU_PRESERVE在创建设备时指定标志D3DCREATE_FPU_PRESERVE以防止其更改FPU精度。

你是如何确定你只获得7位精度的? 您是否正在打印t1并指定正确的输出格式? 在我的机器上,使用VS2010,代码如下:

 int main() { double p = 100.0000; double d = 1000.0000; double t1 = pow(p/(8.0000*d),1.00/4.000); printf("t1=%.15f\n", t1); // C std::cout << "t1=" << std::setprecision(15) << t1 << '\n'; // C++ } 

生成此输出:

 t1=0.334370152488211 t1=0.334370152488211