各种快速SQRT

最新推荐文章于 2022-04-17 08:41:52 发布

潘李亮

最新推荐文章于 2022-04-17 08:41:52 发布

阅读量2.8k

点赞数

文章标签： float

本文链接：https://blog.csdn.net/Nhsoft/article/details/146096

版权

float mysqrt(float x)
{
        union
       {
                    int intPart;
                          float floatPart;
         } convertor;
          union
         {
                        int intPart;
                             float floatPart;
          } convertor2;

          convertor.floatPart = x;
          convertor2.floatPart = x;
          convertor.intPart = 0x1FBCF800 + (convertor.intPart >> 1);
          convertor2.intPart = 0x5f3759df - (convertor2.intPart >> 1);
          return 0.5f*(convertor.floatPart + (x * convertor2.floatPart));
}

float InvSqrt_Lomont(float x)
{
         float xhalf = 0.5f*x;
         int i = *(int*)&x;
         i = 0x5f375a86 - (i>>1);
         x = *(float*)&i;
         x = x*(1.5f-xhalf*x*x);
         // x = x*(1.5f-xhalf*x*x); // added precision
         return x;
} // InvSqrt_Lomont

float InvSqrt (float x)
{
    float xhalf = 0.5f*x;
    int i = *(int*)&x;
    i = 0x5f3759df - (i >> 1);
    x = *(float*)&i;
    x = x*(1.5f - xhalf*x*x);
    return x;
}