Stupéfaction dans un format !

Bonjour, voici un code tout simple :

Code:

---

1 2 3 4 5 6 7 8 9

program test implicit none   double precision :: a = 10.d0, b = 45.3543d12, c = 74.d0   write(*,"(F8.3,1pD16.4,F8.3)") a, b, c write(*,"(F8.3,D16.4,F8.3)") a, b, c   end program test

---

Dans le premier write je met dans le format le 1p pour avoir un chiffre avant la virgule. Jusque là rien d'anormal. Regardez la sortie fait un peu peur 8O :

Code:

---

1 2 3 4

> a.out 10.000 4.5354D+13 740.000 10.000 0.4535D+14 74.000

---

J'ai la même chose avec ifort et gfortran. C'est un truc écrit dans la norme ou c'est un bug de compilateur 8O ?

Salut!
Je ne vois pas ce qui te fait peur. Dans le manuel de référence de mon Compaq Fortran, je lis:

• For F editing, the external value equals the internal value of the I/O list item multiplied by 10^k. This changes the magnitude of the data.
• For E and D editing, the external decimal field of the I/O list item is multiplied by 10^k, and k is substracted from the exponent. this changes the form of the data.

Jean-Marc Blanc

Ce qui me surprend beaucoup c'est la 3eme colonne. Sur la deuxième colonne (où j'applique le format 1p) aucun souci. Le format est correct.

Par contre le chiffre de la 3eme colonne est faux ce qui est plus embétant.

C'est que tu ne comprends pas la règle pour les scaling factor P. Le format P n'est pas un préfixe applicable à un masque D ou E, c'est un masque qui change la représentation à partir de là. Le nouveau comportement s'applique donc jusqu'à la fin du format, ou jusqu'à ce qu'un nouveau P soit défini.

Salut!

Citation:

---

Par contre le chiffre de la 3eme colonne est faux ce qui est plus embétant

---

Ce n'est pas un chiffre, mais un nombre: un chiffre est seulement un caractère typographique.
Jean-Marc Blanc

si j'ai bien compris, le "1p s'applique jusqu'à la fin de la définition du format donc
"(F8.3,1pD16.4,F8.3)" = "(F8.3,1pD16.4,1pF8.3)"
et 1pF8.3 nous donne en fait 10^1*<nombre à afficher>
et comme on demande un format F, il n'affiche pas l'exposant... et donc le résultat semble faux... c'est bien ça?

au fait, au lieu de "1pD16.4" tu ne peux pas utiliser "ES16.4"? (par contre l'exposant est affiché avec un E et non un D)

Merci pour les explications. :aie:

Par rapport à l'écriture, je ne suis pas obligé de mettre un chiffre avant la virgule. Dans ce que je faisais ça m'arrangeait bien, je trouvais ça plus lisible. Mais vu les conséquences je préfère me passer de ce genre de détail.