Discussion :

[Ladgalen]

Bonjour, voici un code tout simple :

Code :

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

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 :

Code :

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

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 ?

[FR119492]

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

[Ladgalen]

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.

[Sylvain Bergeron]

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.

[FR119492]

Salut!

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

[genteur slayer]

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)