Epsilon machine et son calcul

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

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

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

Bonjour,

Citation:

---

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:

---

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

---