Discussion :

[feynman]

bonjour,
j aime savoir comment calculer l'epsilon machine pour double precision en fortran 77
merci pour l aide

[FR119492]

Salut!
Tu trouves tout ce qu'il te faut dans le BLAS.
Jean-Marc Blanc

[FR119492]

Salut!
L'epsilon n'est pas une caractéristique du Fortran, mais du processeur. Comme il est défini par la norme IEEE 754, il n'est pas nécessaire de le calculer: il suffit de le lire dans la norme.
Jean-Marc Blanc

[Ehouarn]

Bonjour,

j aime savoir comment calculer l'epsilon machine pour double precision en fortran 77

Et bien en recherchant la plus petite valeur dont la somme avec 1 diffère de 1 ... par exemple par dichotomie.

Une alternative plus simple et efficace est d'accepter de se projeter dans le 21ème siècle et d'utiliser un peu de Fortran 90; la fonction intrinsèque "epsilon()" est justement là pour ça:

Code :

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

1
2
3
4
5
6
7
8
9
10
 
program show_epsilon
implicit none
double precision :: d
real :: r
 
write(*,*) "epsilon(double precision)=",epsilon(d)
write(*,*) "epsilon(simple precision)=",epsilon(r)
 
end program show_epsilon