Discussion :

[DV666]

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.

[FR119492]

Salut!

print *, a

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

[DV666]

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 !

[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"

[DV666]

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.

[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 .... !!!