Discussion :

[lock42]

Bonjour.
Dans le cadre de mon travail, je dois avoir des variables contenant des très petits nombres : 1E-46. Je m'y connais peu en Fortran, mais ayant créé des variables Real*16, il m'aurait semblé que ca marche or le compilateur (compilateur ifort sous Linux 64bit) me répond :

warning #7919: The value was too small when converting to REAL(KIND=4); the result is zero. [1E-46]

Je vous remercie par avance de votre aide

Cordialement

[bobbyboy]

Bonjour,

J'ai exactement la même configuration que toi (Linux 64bits et ifort) et ca passe en real*16.
Comment initialises-tu tes variables ? La bonne syntaxe en double (ou quadruple) precision est :

Code :

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

1
2
3
 
READ(KIND=16) :: x
x = 1.0D-46

Par contre, en travaillant si près de l'erreur machine, attend toi à avoir des comportements bizzares.

Le simple code suivant laisse déja à réfléchir :

Code :

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

1
2
3
4
5
6
 
READ(KIND=16) :: x , y
x = 1.0D-46
y = 1.1D-46
write(*,*) x
write(*,*) y

Avec ifort, tu peux utiliser l'option suivante pour avoir moins de troncatures mais ça ne résoud pas tout :
-fp-model precise

[lock42]

Merci pour la réponse,
Alors j'utilise la chose comme ca :

Code :

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

1
2
 
      REAL(KIND=16) :: WAVE(5000000),INTS(5000000),INTS1(5000000)

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
  DO 3000 BI = 1,NPOINTS,2
   WRITE (2,4000)INTS1(BI)*1E-46,INTS1(BI+1)*1E-46
 
 
 
 
 
3000   CONTINUE
 
 
 
 
4000    FORMAT(2(1PE10.3))

[bobbyboy]

Essaye de mettre des "D" à la place des "E" comme ceci (et partout ailleurs dans le code où tu utilises un format double précision):

Code :

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

1
2
 
WRITE (2,4000)INTS1(BI)*1D-46,INTS1(BI+1)*1D-46

et tiens moi au courant stp

[lock42]

Bien sur que je te tiens au courant
Alors tu viens de devenir mon nouveau meilleur ami
Ca marche super bien !!!
La science et moi te remercions énormément

[Sylvain Bergeron]

Quelques précisions :

• une variable real (16) est une variable "quad" précision.
• La constante 1d-46 est une constante double précision (et non quad).
• La constante 1q-46 est la version quad de la constante, mais cette forme est non standard.
• La constante 1e-46_16 est la forme standard quad.