Discussion :

[arcanis]

Plop,

Autre probleme etrange :

J'ai une fonction myvec_norme() qui me renvoie la norme d'un vecteur, enregistree dans un double.
Le soucis, c'est que si je fais :
myvec_norme(1,1,1) * myvec_norme(1,1,1)
(la norme d'un vecteur (1,1,1) est de racine de trois)
J'obtient 2,9999999996.
Et pourtant, myvec_norme(1,1,1)==sqrt(3) m'affiche vrai (via gdb).

Une idee ?

[diogene]

Et pourtant, myvec_norme(1,1,1)==sqrt(3) m'affiche vrai (via gdb).

Oui, mais il n'y a aucune chance que sqrt(3) donne la racine carrée de 3 de façon exacte. sqrt(3) ne peut en donner qu'une valeur approchée (certes de très près).

J'obtient 2,9999999996.

C'est vrai que pour un calcul en double, on est assez loin de la bonne valeur (il manque des 9).

[Bobsleigh]

Ta fonction myvec_norme() doit utiliser un calcul du genre sqrt(a^2+b^2+c^2), ici sqrt(1+1+1). Il est donc normal d'avoir au final myvec_norme(1,1,1)==sqrt(3) car ça reviens à faire sqrt(1+1+1)==sqrt(3), ce qui est naturellement vrai. Si tu fais sqrt(3)*sqrt(3), tu as de bonne chances d'obtenir également 2,9999999996, si tu utilise la même précision.