快速开平方取倒数的算法

最新推荐文章于 2024-07-18 22:34:36 发布

rtoax

最新推荐文章于 2024-07-18 22:34:36 发布

阅读量759

点赞数

分类专栏：【算法与数学】

本文链接：https://blog.csdn.net/Rong_Toa/article/details/101191445

版权

【算法与数学】专栏收录该内容

112 篇文章 9 订阅

订阅专栏

/* 快速平方根取倒数 */
float __rt_Q_rsqrt(float number)
{
#define _MAGIC_NUM1  0x5f3759df
#define _MAGIC_NUM2  0x5f375a86

    long i;
    float x2, y;
    const float threehalfs = 1.5F;
    x2  = number * 0.5F;
    y   = number;
    i   = * ( long * ) &y; 
    // evil floating point bit level hacking
    i   = _MAGIC_NUM2 - ( i >> 1 ); // what the fuck?
    y   = * ( float * ) &i;
    y   = y * ( threehalfs - ( x2 * y * y ) ); 
    // 1st iteration
    // y   = y * ( threehalfs - ( x2 * y * y ) );
    // 2nd iteration, this can be removed
#ifndef Q3_VM
#ifdef __linux__
    rt_assert( !t_isnan(y) );
    // bk010122 - FPE?
#endif   
#endif    
    return y;
}

https://mp.weixin.qq.com/s/MX6b1AZwabKjiHz2TpQQ9w