Discussion :

[lambertmular]

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 :

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

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 :

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

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

[__dardanos__]

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 :

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

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 :

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

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

[lambertmular]

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

[__dardanos__]

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