我遇到过的神奇的难以解释的算法

一. 求刚好大于某个数的2的幂的数.

　　即 Round up to the next highest power of 2.

　　需求来源: 在开发图形引擎的时候,一般的纹理对象的长和宽必须是2的幂. 比如128,256,1024这样的数.

     但是我们的纹理不一定是这样的,因此有时候需要进行纹理的补充.

　　如果有人能知道该算法原理,求告知求地址

　　http://graphics.stanford.edu/~seander/bithacks.html

　　devised by Sean Eron Aderson

　　http://graphics.stanford.edu/~seander/

　　方法1:

int nextPowerOf2(int n) { 
    n -= 1; 
    n |= n >> 1; 
    n |= n >> 2; 
    n |= n >> 4; 
    n |= n >> 8; 
    n |= n >> 16; 
    return n + 1; 
} 

　　方法2:

　　　　下面的函数仅在un为2的幂时返回0;

　　　　

uint32_t is2n(uint32_t un)
{
    return un&(un-1);
}

　　　　下面的函数返回不大于un的2的最大幂；

uint32_t max2n(uint32_t un)
{
    uint32_t mi = is2n(un);
    return mi?max2n(mi):un;
}

二. 大名鼎鼎的平方根倒数速算法

　　http://zh.wikipedia.org/wiki/%E5%B9%B3%E6%96%B9%E6%A0%B9%E5%80%92%E6%95%B0%E9%80%9F%E7%AE%97%E6%B3%95

float Q_rsqrt( float number )
{
    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（对浮点数的邪恶位级hack）
    i  = 0x5f3759df - ( 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（第二次迭代，可以删除）
 
    return y;
}

 

未完待续.......