Discussion :

[vinzzzz]

Hello
Lorsque j'execute ceci:

Code :

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

1
2
3
 
double ts = atof("0.6");
printf("ts: %.10f\nts: %.40f\n", ddd, r0, ts, ts);

J'obtient cela:

ts: 0.6000000000
ts: 0.5999999999999999777955395074968691915274

Bref, je constate qu'il y a une perte énorme (tout est relatif ms ds mon cas c'est énorme) de précision avec la fonction atof!!

Donc, dans la mesure où je dois récupérer des arguments de la ligne de commande, et les convertir en double, quelqu'un aurait-il une idée de la manière dont je dois m'y prendre afin de garder une présision optimale?

Merci d'avance

[Jean-Marc.Bourguet]

Le format double est incapable de representer exactement 0.6. Voir les archives, on a explique cela en long et en large. Voir aussi:
http://docs.sun.com/source/806-3568/ncg_goldberg.html

[vinzzzz]

Je sais bien. Le truc c'est que lorsqu'on utilise par exemple le code suivant:

Code :

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

1
2
3
 
double ts = 0.6 ;
printf("ts: %.10f\nts: %.40f\n", ddd, r0, ts, ts)

On obtient bien une valeur satisfaisante pour ts. Je voulais simplement savoir si il y avais un moyens d'obtenir une meilleur précision que celle obtenue avec atof pour ce que je souhaite faire...

[Jean-Marc.Bourguet]

Je sais bien. Le truc c'est que lorsqu'on utilise par exemple le code suivant:

Code :

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

1
2
3
 
double ts = 0.6 ;
printf("ts: %.10f\nts: %.40f\n", ddd, r0, ts, ts)

On obtient bien une valeur satisfaisante pour ts. Je voulais simplement savoir si il y avais un moyens d'obtenir une meilleur précision que celle obtenue avec atof pour ce que je souhaite faire...

J'obtiens:

ts: 0.6000000000
ts: 0.5999999999999999777955395074968691915274

Je ne vois pas pourquoi c'est plus acceptable que ce qu'il y a ci-dessus.

[vinzzzz]

Ah tiens moi j'obtient une longue série de 0...

ts: 0.6000000000
ts: 0.6000000000000000000000000000000000000000

[droggo]

Jai,

De toute manière, pourquoi veux-tu afficher 40 décimales pour un double, alors que ce format donne au mieux 16 (ou 17, je ne sais plus ) chiffres significatifs ?
(sans oublier que dans la gamme possible, les valeurs représentées exactement sont rares. )

Et quand tu fais

Code :

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

printf("ts: %.10f\nts: %.40f\n", ddd, r0, ts, ts)

ton compilateur ne te signale pas que le nombre d'arguments passés ne correspond pas à ceux demandés dans la chaîne de format ?

[Jean-Marc.Bourguet]

Ah tiens moi j'obtient une longue série de 0...

Intéressant. Est-ce que tu peux afficher les deux avec

Code :

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

1
2
3
4
5
6
7
8
9
10
void printhex(double d)
{
   static char tohex[]="0123456789ABCDEF";
   unsigned char* ptr = (unsigned char*)&d;
   int i;
   for (i = 0; i < sizeof d; ++i) {
      putchar(tohex[(ptr[i]>>4) & 0xF);
      putchar(tohex[(ptr[i]>>0) & 0xF);
   }
}