问下开方函数

shaoxiaoning

是下面这个 计算快？ 还是用牛顿迭代法算法快啊。

1. 
   
2. #include   <math.h>   
   
3. #include   <stdio.h>   
   
4.    
   
5.   void   main(   void   )   
   
6.   {   
   
7.         double   x   =   8.0,   y   =   1.0/3.0,   z;   
   
8.    
   
9.         z   =   pow(   x,   y   );   
   
10.         printf(   "%.1f   to   the   power   of   %.1f   is   %.1f\n",   x,   y,   z   );   
    
11.   }   
    
12. 
    

复制代码

tarkey

明显Newton Raphson's method要快,贴段代码,以前算布的时候用过的:

1. 
   
2. float SquareRoot(float x)
   
3. {
   
4.         DWORD index;
   
5.         float r;
   
6.         DWORD *dptr = (DWORD *)&r;
   
7.        
   
8.         *(DWORD *)&r = (((EXPONENT_BIAS << MANTISSA_BITS) + EXP_OF(x)) >> 1) & EXPONENT_MASK;
   
9. 
   
10.         index = ((*(DWORD *)&x) >> (MANTISSA_BITS - SQRT_TABLE_ENTRY_BITS + 1))
    
11.                         & (SQRT_TABLE_NUM_ENTRIES - 1);
    
12.         *dptr |= SqrtTab[index];
    
13. 
    
14.         // One iteration of Newton Raphson's method will double the accuracy
    
15.         r = (3.0f - (r * r) / x) * r * 0.5f;
    
16. 
    
17.         return r;
    
18. }
    

复制代码

不过用之前要初始化一次表,全部代码见附件:

shaoxiaoning

多谢大侠， 不过 猛然一看不知所措。

tarkey

原理就是把浮点数拆成sign位,mantissa和exponent位,然后根据这些信息插值

shaoxiaoning

此外 如果用 预计算生成一个表，用查表法 和 牛顿迭代比较呢？ 据说 在一定条件下用
Hash 查表 可以一次就找出来，小弟没试过，希望 介绍下。

tarkey

查完表之后要插值的..

kevindophi

quake3原代码里有一个carmack实现的开方函数。。。