Esope source files

   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 :