Discussion :

[gil38]

Bonjour,
Je developpe sur visual C++ 2003 et j'ai un probleme de précision de calcul.
Voici un bout de code :
double VarDbl1,VarDbl2;
double VarDblRes;

VarDbl1 = 5.2;
VarDbl2 = 4;

VarDblRes = VarDbl1 + VarDbl2 ;

Le resultat n'est pas 9.2 comme on pourrait si attendre mais 9.1999999999999993 !!!
Je vous laisse imaginer ce que cela peut donner dans une equation complexe.
Je ne vois pas comment resoudre ce pb.
Y a t il une option dans le compilateur ?

Merci d'avance pour votre aide

[Laurent Gomila]

Ces imprécisions sont inhérentes à l'encodage des nombres flottants. Si tu veux plus de précision il faudra te tourner vers, au choix :

- Le type long double qui peut faire 8 octets sur certaines architectures
- Des classes spéciales gérant des flottants avec plus de précision
- Des classes de nombres en virgule fixe
- Le calcul formel

[jobherzt]

tu as aussi l'option d'utiliser des intervalles pour representer les floats.. comme ca tu garantit le resultat dans un intervalle.

mais je te signale quand meme que l'erreur sur le resultat de ton calcul est inferieur a 10^-15, meme si ca fait drole c'est quand meme minuscule. apres, dans ce cas precis je t'accorde que c'est quand meme etrange...

[Ulmo]

Apres, dans ce cas precis je t'accorde que c'est quand meme etrange...

Pas tant que ça si tu remarques que 0.2 ne tombe pas juste (en binaire).

[jobherzt]

bien vu

[gil38]

Merci pour vos réponses.

En fait mon probleme vient du fait que j'ai besoin d'arrondir cette valeur au 0.5 le plus proche.
j'utilise cette petite methode

double AR;

AR = 10*Val- (WORD)(10*Val);
if (AR >= 0.5)
Val = (float)((((float)(WORD)(10*Val)) + 0.5)/10);
else
Val = (float)(((float)(WORD)(10*Val))/10);

return Val;

A cause de cet imprecision l'arrondi ne se fait pas toujours du bon coté.

Pour minimiser ce pb j'ajoute à ma valeur un epsilon (0.0001) et cela me donne des résultats plus correcte.

[Laurent Gomila]

Et avec ce code ?

Code :

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

return static_cast<int>(Val * 2 + 0.5) * 0.5;

[poukill]

Laurent vient de te montrer:
1) que ton code de 7 lignes tient en une seule!
2) la bonne façon de coder en C++ (static_cast, etc...)
3) qui sait, la solution à ton problème...

En tous cas, c'est un code vraiment mieux

[gil38]

Ok merci pour la leçon de prog
je ne connaissais pas cette syntaxe
et en plus cela a fait disparaitre les problemes de precision
Mon pb est résolu !