Discussion :

[deadalnix]

Bonjour, c'est encore moi avec mes problèmes tordus Je suis en pleine période je martyrise mon pauvre gcc.

Au menu du jour, la fonction qui calcule l'inverse d'une racine carré a la sauce quake 3. L'original est en C :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
    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
    return y;

Voila ce que donne l'original.

Mais ici on est en C++ et ces cast mode C ne font pas propres ! j'ai donc essayé différentes chose, mais impossible de faire gober ca a mon gcc sans qu'il ne se fache :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

i = reinterpret_cast<long>(y);   // evil floating point bit level hacking

Donne une erreur : invalide cas from type float to type long int.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

i = *reinterpret_cast<long*>(&y);   // evil floating point bit level hacking

Donne quand a lui attention : déréférencement du pointeur type-punned brisera les strictes d'aliases . Il est possible de compiler quand même mais le code ne fonctionne pas.

Google ne m'aide pas beaucoups sur ce coup la. je crois qu'il ne m'aime aps ces temps-ci . . .

A noter que sans l'option de compilation -O3, le second code marche.

[deadalnix]

J'ai trouvé la solution en fouinant sur cette page : http://gcc.gnu.org/onlinedocs/gcc/Op...timize-Options

Il faut utiliser une union du type :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
union converter {
   long i;
   float y;
};

[Laurent Gomila]

Les unions ça fonctionne ne pratique mais pas en théorie, donc à éviter. Ceci-dit je pense qu'il n'y a que ça qui fonctionne avec les options d'optimisations poussées.

A part ça, les hacks de ce genre c'est bien mais c'est plus trop utile. Si tu finis ton jeu et lances un profiler, tu verras que ce n'est pas là que les pertes de performances se trouvent. Alors ne te prends pas trop la tête avec ce genre de code.

[deadalnix]

A part ça, les hacks de ce genre c'est bien mais c'est plus trop utile. Si tu finis ton jeu et lances un profiler, tu verras que ce n'est pas là que les pertes de performances se trouvent. Alors ne te prends pas trop la tête avec ce genre de code.

Je sais bien, mais je suis dans une période ou je teste un peu tout et n'importe quoi pour voir comment ça réagit. Je ne compte pas recoder intégralement la lib mathématique

[NairodDorian]

ça c'est le genre de code C totalement inutile et illisible mais optimisé.
deux instructions asm et c'est réglé avec des perfs supérieur.

[deadalnix]

Je ne penses pas que carmack ou tarolli puissent être traités de branlo faisant du mauvais code

[yan]

Salut.
Attention normalement, le chiffre magique ne fonctionne que sur 32 bits et
long et float doit être de même taille.

Si ca t'intéresse, voici un benchmark
ktd.club.fr/programmation/fichiers/FastMath.pdf