exact ca peut marcher comme ca aussi mais les racines carrées sont pas mal lentes non plus.

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
    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
    
    #ifndef Q3_VM
    #ifdef __linux__
        assert( !isnan(y) ); // bk010122 - FPE?
    #endif
    #endif
    return y;
}

Fonction de racine carrée utilisée par John Carmack dans Quake 3 (license GPL attention), beaucoup plus rapide que la fonction native. (Oui oui chelou hein ca utilise le nombre 0x5f3759df me demandez pas à quoi il correspond )