Télécharger epslon.eso

Retour à la liste

Numérotation des lignes :

epslon
  1. C EPSLON SOURCE CHAT 05/01/12 23:42:50 5004
  2. FUNCTION EPSLON ( X )
  3. *
  4. *****************************************************************************
  5. * ESTIMATION DE L ERREUR D ARRONDI DANS CHAQUE MACHINE *
  6. * _________________________________________________________________________ *
  7. * *
  8. * DATE : 19 avril 1995 *
  9. * AUTEUR : Nicolas BENECH *
  10. * _________________________________________________________________________ *
  11. * *
  12. * MODULE(S) APPELANT(S) : QZBASC, QZREDU *
  13. * *
  14. * MODULE(S) APPELE(S) : *
  15. * _________________________________________________________________________ *
  16. * *
  17. * EN ENTREE : *
  18. * - X : unite de base (en general 1) *
  19. * _________________________________________________________________________ *
  20. * *
  21. * EN SORTIE : *
  22. * - EPSLON : erreur d'arrondi *
  23. * _________________________________________________________________________ *
  24. * *
  25. * REMARQUE : Ce sous-programme doit bien fonctionner dans les systemes *
  26. * tels que : *
  27. * - la base de representation des reels n'est pas une *
  28. * puissance de 3; *
  29. * - la quantite A est representee avec la precision *
  30. * utilisee dans le stockage en memoire. *
  31. *****************************************************************************
  32. *
  33. IMPLICIT INTEGER(I-N)
  34. IMPLICIT REAL*8 (A-H,O-Z)
  35.  
  36. *
  37. A = 4.0D0/3.0D0
  38. 10 B = A - 1.0D0
  39. C = B + B + B
  40. EPS = ABS(C-1.0D0)
  41. IF ( EPS . EQ . 0.0D0 ) GOTO 10
  42. EPSLON = EPS*ABS(X)
  43. *
  44. RETURN
  45. END
  46.  
  47.  
  48.  
  49.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales