Erreur de précision

Bonjour à tous, voici grosso modo mon code :

Code:

---

1 2 3 4 5 6 7 8 9

double a = 720; double b = 152; double c = a/b; if (c==5){ cout << "Youhou" << std::endl; } if (abs(c-5) < 1e-16){ cout << "Bouhouhou" << std::endl; }

---

et j'obtiens sur ma machine : "Youhou" et "Bouhouhou" mais lorsque j'envoie le code sur mon serveur, je le compile (machine linux), je n'obtiens que "Bouhouhou"

Existe-t-il des règles à respecter pour bien écrire un calcul et ne pas avoir des erreurs d'approximation ?

Bonsoir,

Faits touts tes calculs à "epsilon" près et tu n'auras pas de problème. C'est ainsi qu'on doit faire lorsque l'on travail avec des réels.

En fait, ce qui m'ennuie, c'est que je fais après un ceil(c) par après. Sachant que mon c vaut 5.000...0001 le résultat vaudra 6 (sur ma machine, 5). Là je ne vois pas comment faire à epsilon près. Si quelqu'un a une idée de bonne pratique.

Je suis obligé de rester en "double". Là je suis tombé sur ce cas particulier, pour toutes autres valeurs, ça marche bien, sauf lorsqu'on tombe sur un entier. Enfin j'ai pas testé toutes les possibilités.

Bonjour,
Désolé de mettre mon grain de sel, mais 720 / 152 = 4,74 à 0,01 près.
C'est étrange que tu obtiennes un c si près de 5.
Après un ceil(), ok, mais pas avant.

Tu es sûr de tes valeurs ?

j'accepte ce grain de sel, je me suis trompé en tapant. C'est 760 à la place de 720. Autant pour moi.

J'ai fixé le problème avec des nombres à virgule fixe.

Ok. :ccool:

Reste que tant que tu es en virgule flottante le test c==5 ne sera jamais fiable car ces nombres sont approximatifs par construction. Comme le faisais remarquer Flob90, tu ne peut te fier qu'à des inégalités comme abs(c-5)<epsilon.

Concernant ceil:

Citation:

---

Returns the smallest integral value that is not less than x.

---

. Si l'erreur t'amène à 5-epsilon, ça donne 5, sinon 6. Par ailleurs, comme ça retourne toujours un double, tu ne peux toujours pas tester l'égalité car rien ne te dit que la représentation IEE du 5 ou 6 tombe juste sur 5 ou 6...l'imprécision des nombre à virgule flottante encore==>abs(c-5)<epsilon encore et toujours.
EDIT: il faut bien sûr prendre la valeur absolue de c-5. Tu demandais la bonne pratique: tu l'avais déjà implémentée quelque part dans ton code!