R devel - Mar 1999 - bessel_?.c constants (fwd)

This doesn't seem to have made it to the list, so here it is again.  Sorry
if this results in duplicates.

---------- Forwarded message ----------
Date: Fri, 5 Mar 1999 21:36:44 -0800 (PST)
From: Gregory R. Warnes <warnes@biostat.washington.edu>
To: R devel <r-devel@stat.math.ethz.ch>
Subject: bessel_?.c constants


Hi All,

I've been digging around in src/nmath and have discovered that there are
two sets of "machine constants" in bessel_l.c and bessel_j.c which
have
different values:

bessel_i.c:    static double ensig = 1e16;
bessel_i.c:    static double rtnsig = 1e-4;
bessel_i.c:    static double enmten = 8.9e-308;
bessel_i.c:    static double enten = 1e308;

bessel_j.c:    static double ensig = 1e17;
bessel_j.c:    static double rtnsig = 1e-4;
bessel_j.c:    static double enmten = 1.2e-37;
bessel_j.c:    static double enten = 1e38;

According to the comments heading the files, the first set is correct for
IEEE Double-Precision hardware (the comments give the same values for
these parameters in each file):


         Approximate values for some important machines are:
    
    			    it	  NSIG	  ENTEN	      ENSIG
      [...]    
      IEEE (IBM/XT,
      SUN, etc.) (D.P.)	    53	   16	 1.0D+308    1.0D+16
      [...]

    			   RTNSIG      ENMTEN	   XLARGE
      [...]   
      IEEE (IBM/XT,
      SUN, etc.) (D.P.)	   1.0E-4    8.90D-308	  1.0D+4
      [...]

Is there some reason for using different values in bessel_j.c than in
bessel_l.c?  Or is this just a bug? ;^)

-Greg

PS:  Where is DBL_EPSILON defined, and what is its value?




-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To:
r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

"Gregory R. Warnes" <warnes@biostat.washington.edu> writes:
> This doesn't seem to have made it to the list, so here it is again. 
Sorry
> if this results in duplicates.
It did...
> PS:  Where is DBL_EPSILON defined, and what is its value?
In the compiler libs, here:

/usr/lib/gcc-lib/i386-redhat-linux/2.7.2.3/include/float.h:
#define DBL_EPSILON 2.2204460492503131e-16


-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To:
r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._