Problème avec <limits>

Bonjour,

J'ai un problème dont je ne parviens pas à trouver l'origine.
Quand j'exécute le code:

Code:

---

1 2 3 4 5 6 7 8 9 10

#include <limits> #include <iostream>   int main() { double test_val=std::numeric_limits<double>::max()/10.; std::cout<<"diff 1: "<<test_val-test_val<<"\n"; std::cout<<"diff 2 : "<<std::numeric_limits<double>::max()/10.-std::numeric_limits<double>::max()/10.<<"\n"; }

---

j'obtiens:


diff 1: 0
diff 2 : 4.99447e+290

Sincèrement je ne sais pas trop quoi penser ...
Est que quelqu'un pourrait compiler et exécuter ce petit programme et me dire ce qu'il
obtient.
L'idéal serait une personne qui utilise gnu gcc c++ version 4.0.1 sur Mandriva 2006.
Merci.
Amicalement,
Michel

Ma conf :

mandriva 2007.1
gcc version 4.1.2 20070302 (prerelease) (4.1.2-1mdv2007.1)

résultat :

Code:

---

1 2 3

diff 1: 0 diff 2 : 4.99447e+290

---

Hope this will help !

Merci beaucoup !

Ca ressemble bien à un bug ...

Je ne sais pas, j'ai "bêtement" copié/collé ton code sans chercher à comprendre.

Qu'est-ce que tu essayes de mettre en évidence ?

Les deux valeurs "diff 1" et "diff 2" devraient être nulles.
C'est supposé faire la même différence.

Cela me semble être un comportement normal, les opérations sur les flottants ce n'est pas du calcul symbolique, c'est du calcul numérique avec une précision finie.

C'est du au meme problème que:

Code:

---

std::numeric_limits<double>::max()-10 == std::numeric_limits<double>::max()

---

Les flottant sont une approximation assez proche mais qu'une approximation des reels:

en faisant:

std::numeric_limits<double>::max()/10.-std::numeric_limits<double>::max()/10,

Le programme te retourne le un resultat qui "approxime" zero, de la meme façon que parfois x/10-x/10 = 1.E-71 ou du genre.

En fait, il faut comprendre que 4.99447e+290, c'est le plus petit "pas" pour un flottant aussi grand, c'est à dire que la valeur de la mantisse est multiplié par ce nombre.

Je te propose d'aller jeter un coup d'oeil sur le net sur la norme IEEE 754, tu comprendras bien mieux que j'explique :aie:

Oui mais pourquoi les deux facons de calculer cette différence ne me retournent pas le même résultat ?

test_val-test_val doit etre optimisé par le compilateur qui se dit que "var moins var vaut zero", enfin j'imagine ;)