Double precision et print

Bonjour,

Je suis en train d’écrire un programme ou j'ai besoin d'une grande précision (double precision). Cependant je constate une chose bizarre lorsque je fais un print de mes variables. Par exemple j'ai un vecteur complexe (dimension 3) :

Code:

---

beta=(/(1.246663751424521,0.249332750284904),(0,0),(0,0)/)*1d4

---

Si je fais un print de ce nombre a partir de la 8eme décimale le nombre est différent de celui que j'ai rentré. pourtant j'ai bien déclaré

Code:

---

complex(kind=8),dimension(3) :: beta

---

Je pense bien que c'est un problème d'arrondis et de double precision mais je ne vois pas ou j'ai commis l'erreur.
Merci pour votre aide

Salut,
C'est parce que les composantes des nombres complexes sont déclarées en simple précision, puis converties pour initialiser le nombre (qui ensuite est multiplié par 1.d4). Il me semble que la syntaxe est :

Code:

---

1 2 3 4 5 6 7 8 9 10

program main implicit none complex(kind=8) :: c1, c2   c1 = (1.246663751424521, 0.249332750284904) * 1.d4 print '(2e24.16)', c1   c2 = (0.1246663751424521D5, 0.249332750284904D4) print '(2e24.16)', c2 end program main

---

Ce qui donne :

Code:

---

1 2	0.1246663808822632E+05 0.2493327558040619E+04 ! -> pb 0.1246663751424521E+05 0.2493327502849040E+04 ! -> ok

---

Merci bien! C'est vrai seulement pour les complexes ou aussi pour les arrays?

C'est valable pour déclarer tout nombre (réel ou complexe), qu'il soit scalaire ou dans un tableau.