Constantes en double précision

Version imprimable

16/02/2011, 22h34
DV666

Constantes en double précision

Bonjour,

J'ai des calculs très volumineux à réaliser avec une extrême précision (11 ième chiffre après la virgule). Ces développements sont importables (sans resaisie) sous la forme :

Q = Q + 0.40989414976 * Cos(1.48302034194 + 26087.9031415742 * T)
Q = Q + 0.05046294199 * Cos(4.47785489540 + 52175.8062831484 * T)
Q = Q + 0.00855346843 * Cos(1.16520322351 + 78263.7094247226 * T)
Q = Q + 0.00165590362 * Cos(4.11969163181 + 104351.612566297 * T) etc..

Malheureusement je viens de me rendre compte que les constantes ne sont pas gérées en double précision !!!

Double précision :: a

a = 4.40250710144

print *, a

a qui est affiché : 4.40250730514

si cette fois je déclare a = 440250.710144d-5

a qui est affiché devient 4.40250710144000 ce qui est bon.

J'ESPERE QUE JE NE VAIS PAS DEVOIR RE-ECRIRE EN PUISSANCE DE 10 LES 10800 CONSTANTES DE MES CALCULS !!!!

Y-a-t-il une solution simple pour obtenir des calculs justes avec mes développements importés ?

Merci d'avance.
17/02/2011, 01h49
FR119492

Salut!

Citation:

print *, a

Pourquoi n'utilises-tu pas des formats pour la lecture et l'affichage de tes valeurs?
Jean-Marc Blanc
17/02/2011, 08h59
DV666

Citation:

Envoyé par FR119492

Salut!

Pourquoi n'utilises-tu pas des formats pour la lecture et l'affichage de tes valeurs?
Jean-Marc Blanc

Bonjour,

Le formatage de sortie n'a pas de lien avec mon problème, car seuls les 6 permiers chiffres après la virgule sont pris en compte sur les constantes écrites sous le format standard. Il semble qu'il faille ré-écrire toutes les constantes sous la la forme xxxxx.yyyyyyd-z pour obtenir tous les chiffres significatifs !!! Même si je formate le format de sortie, le résultat reste le même et les calculs sont faux !
17/02/2011, 09h50
bobbyboy

Bonjour,

Quel compilateur utilisez-vous? Il faut forcer le compilateur à utiliser par défaut des réels double précision via ces options :

sous ifort "-r8"
sous gfortran "-fdefault-real-8"
17/02/2011, 12h07
DV666

Citation:

Envoyé par bobbyboy

Bonjour,

Quel compilateur utilisez-vous? Il faut forcer le compilateur à utiliser par défaut des réels double précision via ces options :

sous ifort "-r8"
sous gfortran "-fdefault-real-8"

j'utilise INTEL VISUAL FORTRAN sous environnement VISUAL STUDIO 2008.
je vais regarder si il y a des options de compilation dans ce sens ... Merci.
17/02/2011, 21h34
DV666

Une constante double précision doit finalemnt obligatoirement être écrite en virgule flottante, le e étant remplacé par un d.

Exemples :

0d0
0.d0
1.d0
1d0
3.1415d0
31415d-4
1.6d-19
1d12
-36.d0

eh oui il va falloir ajouter un d0 à toutes les constantes .... !!!
17/02/2011, 21h55
Sylvain Bergeron
Si tu dois modifier tout ton code pour beaucoup de constantes, tu devrais regarder du côté des kinds de Fortran 90.

Au lieu d'ajouter d0 à la fin de chaque constante, tu dois plutôt ajouter _"kind". Par exemple :
Code:

1 2 3 4 5 6 integer, parameter :: rk = 8 ! ou = selected_real_kind(15) ... real (rk) :: Variable ... Variable = 1234.56789_rk * ...
A priori, ça peut sembler très lourd, mais tu auras l'avantage de pouvoir changer très facilement la précision de ton code en changeant uniquement la valeur de rk...