Discussion :

[michubuntu]

Bonjour à tous !
Voila j'ai deux questions,
La première :

J'ai un code en FORTRAN 77 qui écrit (WRITE) la valeur de reel (4 Octets)
avec le FORMAT E16.9, Or il me semble que la précision sur quatres Octets est au maximum de 7 chiffres significatifs donc un E14.7 devrai etre suffisant non ?
Si c'est bien le cas les valeur écritents au dela du 7 ieme chiffre sont aléatoire ?

La deuxième :
Sans parler de l'algorithme (qui est douteux je sais), le code lit un fichier ASCII contenant des valeurs flottantes et les stoquent dans une variable de type CHARATER disons string
ensuite ses valeurs sont converties en reel via l'instruction READ de type :
read (string,*) flottant

la forme est stupide (je sais) mais ce qui m'etonne c'est que la valeur flottante n'est pas tout à fait la même que celle affichée par le string.
par exemple :

string=' 1.234000E+02 '
read (string,*) flottant

flottant va par exemple prendre la valeur :
1.234012E+02

Quelqu'un à une explication ?

Merci

[michubuntu]

Je suis un peu déçu de ne pas avoir de réponse.

Peut être ai-je pal posé mes questions ou je suis simplement trop impatient ....

Dans tous les cas merci de m'aider à comprendre

[Fortran90]

salut

Faut être patient dans la vie hein...

Alors un réel simple précision est stocké sur 32 bits (4 octets ou bytes)
8 sont réservés à l'exposant, et 23 à la mantisse (là où on met les chiffres qu'on voit après la virgule en gros)
23+8=31...bien sur parceque le denier bit est réservé au signe.
Alors avec 23 bits on monte jusque 2^22 en maximum représentatif d'où la limite des sept chiffres significatifs.
si tu en écris plus c'est au bon vouloir du compilo. Exemple, avec en input X=0.1111111111111
avec un print en E50.25 (soyons fous) on a
pgi donne 0.1111111100000000000000000E+00
intel donne 0.1111111119389533996582031E+00, tout comme xlf et f77 (gnu). Peut être utilsent-ils le même algorithme de gestion de précision. En tout cas ce qui est sur c'est que ces chiffres (au delà du 7e) ne sont pas dans la mantisse de X.

Pour ta deuxième question, c'est plus étrange en deçà de la limite des 7 chiffres.

Si tu as CHAR='0.123456789E+02'
tu fais READ(CHAR,*)NUM
alors NUM=0.1234567XXE+02
et XX <> 89 ... of course, mais pour sur on a au moins nos 1234567 ...

As tu un exemple concret pour ton soucis ( à moins de 7 chiffres significatifs) ?

[michubuntu]

Ha merci pour ta réponse.
Sinon pour l'histoire du peuso transtypage, c'était du meme ordre d'idée

j'affichais plus de chiffre que la précision maximale d'ou l'ajout...

Voila encore merci et bonne soirée
--