Discussion :

[zaim]

Bonjour à vous tous;
Alors voila un code qui me fait des migraines:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
#include <iostream.h>
#include <iomanip.h>
int main()
{
long double a,b,c,d;
a = 1;
b = 1e-09;
c = a - b;
d = 10*( a - b );
cout << " a = "<< a ;
cout <<"\n b = "<< b ;
cout <<"\n c = a - b = "<< c ;
cout <<"\n d = 10*( a - b )= "<<d<<"\n";
}

mon problème est qu'il me donne a - b = 1 et 10*(a - b) = 10 alors qu'il devrait me donner 0.999999999 et 9.99999999.
J'ai même essayé b = 1 - 1e-09 puis evaluer 1/(a-b), et alors il me donne 1/(a-b) = inf.
Quel est le problème et son remède.
Merci.

[Matthieu Brucher]

Quelle est la valeur de std::numeric_limits<long double>::epsilon() ?

[zaim]

Alors

Code :

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

std::numeric_limits<long double>::epsilon()

me donne 1.0842e-19.
Alors c'est grave Docteur ?
ah !! j'oubliai de préciser que je travail sous SuSE 9.0 avec g++. sur un pc.

[Matthieu Brucher]

Peut-être dû à la surcharge de l'opérateur << pour les long double ?
Tu as regardé au débuggeur la valeur de c et d ?

[Laurent Gomila]

Et en augmentant le nombre de décimales pour l'affichage ?

http://c.developpez.com/faq/cpp/?pag...stream_iomanip

[Jean-Marc.Bourguet]

Dans les formats les plus courants -- et les autres sont vraiment exceptionnels -- ne sont representables exactement par des flottants (float, double, long double en C++) que des nombres de la forme

(-1)^s m 2^e

ou ^ denote l'exponentiation
s vaut 0 ou 1
m est un entier naturel
e est un entier relatif

Donc la plupart des nombres ecrits sous forme decimale ne sont pas representables exactement. Il ne faut donc pas t'attendre a voir un flottant valoir precisement 0.999999999.

Dans le cas present, tu peux t'arranger pour que ce soit la valeur affichee en jouant sur le nombre de chiffres apres la virgule (voir par exemple le lien de Laurent).