home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #3 / NN_1993_3.iso / spool / sci / math / numanal / 3977 < prev    next >
Encoding:
Text File  |  1993-01-28  |  1.9 KB  |  49 lines
  1. Newsgroups: sci.math.num-analysis
  2. Path: sparky!uunet!ftpbox!news.acns.nwu.edu!zaphod.mps.ohio-state.edu!caen!sol.ctr.columbia.edu!howland.reston.ans.net!spool.mu.edu!sgiblab!munnari.oz.au!cs.mu.OZ.AU!mundoe!kbriggs
  3. From: kbriggs@mundoe.maths.mu.OZ.AU (Keith Briggs rba8 7088)
  4. Subject: LAPACK machine constants
  5. Message-ID: <9302717.22929@mulga.cs.mu.OZ.AU>
  6. Keywords: LAPACK, machine constants
  7. Sender: news@cs.mu.OZ.AU
  8. Organization: Computer Science, University of Melbourne, Australia
  9. Date: Wed, 27 Jan 1993 06:31:02 GMT
  10. Lines: 37
  11.  
  12. The routine DLAMCH in LAPACK is supposed to return machine constants.
  13. Here is the output when compiled and run with xlf on a IBM RISC/6000:
  14.  
  15.  eps      = 0.111022302462515654E-15  (machine eps)
  16.  sfmin    = 0.22250738585072015E-307  (safe minimum)
  17.  base     =  2.00000000000000000     
  18.  eps*base = 0.222044604925031308E-15 
  19.  t        =  53.0000000000000000      (digits)
  20.  rnd      =  1.00000000000000000     
  21.  emin     = -1021.00000000000000     
  22.  rmin     = 0.22250738585072015E-307 
  23.  emax     =  1024.00000000000000     
  24.  rmax     = 0.17976931348623157E+309 
  25.  
  26. No problems!   However, after converting the type declarations to 
  27. quadruple precision we get
  28.  
  29.  eps      =-0.1112536929253600691545116358666202E-307 
  30.  sfmin    =-0.46566128730773925781250000000000000E-09 
  31.  base     = -2147483648.0000000000000000000000000     
  32.  eps*base = 0.2389154863368240330214629234492972E-298 
  33.  t        =  34.000000000000000000000000000000000     
  34.  rnd      = 0.00000000000000000000000000000000000E+00 
  35.  emin     = 0.00000000000000000000000000000000000E+00 
  36.  rmin     =-0.46566128730773925781250000000000000E-09 
  37.  emax     =  1.0000000000000000000000000000000000     
  38.  rmax     = -2147483648.0000000000000000000000000     
  39.  
  40.  
  41. The weird (and unusable) results are caused by the IBM representation of
  42. quadruple as two doubles.    In other words, DLAMCH has not been written
  43. in a sufficiently portable way.
  44.  
  45. DOES ANYBODY KNOW THE CORRECT VALUES?
  46.  
  47. Keith.
  48.  
  49.