home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #26 / NN_1992_26.iso / spool / comp / sys / apollo / 3953 < prev    next >
Encoding:
Internet Message Format  |  1992-11-13  |  3.8 KB
  1. Path: sparky!uunet!caen!spool.mu.edu!m.mscs.mu.edu!tone
  2. From: tone@m.mscs.mu.edu (Peter Tonellato)
  3. Newsgroups: comp.sys.apollo
  4. Subject: Machine Constants for Numerical Packages. (DN10000)
  5. Date: 13 Nov 1992 15:54:48 GMT
  6. Organization: Marquette University - Department MSCS
  7. Lines: 97
  8. Distribution: world
  9. Message-ID: <1e0j48INN4fs@spool.mu.edu>
  10. Reply-To: tone@m.mscs.mu.edu (Peter Tonellato)
  11. NNTP-Posting-Host: m.mscs.mu.edu
  12. Originator: tone@m.mscs.mu.edu
  13.  
  14. Dear Netters:
  15.  
  16. Typical numerical packages include subroutines that determine the standard
  17. machine constants used for machine-epsilon accurate numerical computations.
  18.  
  19. A common subroutine (e.g. used in MINPACK) is DPMPAR (or SPMPAR for single
  20. precision.) This subroutine has been designed so that appropriate entries
  21. for many architectures are available. To use the routine on your home
  22. machine one need only uncomment the appropriate entry. 
  23.  
  24. My question is can someone please send me the correct machine constants for
  25. the DN1000 PRISM architecture? I would very much appreciate it and will submit
  26. the entry to "netlib" for inclusion in future releases of the public domain
  27. numerical packages they offer.
  28.  
  29. Thanks Very Much,
  30.  
  31. Peter J. Tonellato
  32.  
  33. US Mail: Department of Mathematics, Statistics and Computer Science
  34.          Marquette University
  35.          Milwaukee, Wisconsin 53233
  36. Voice:   414-288-5228 (Office)
  37.          414-288-7573 (Message)
  38. FAX:     414-288-5472 (FAX)
  39. Arpa:    tone@m.mscs.mu.edu   
  40.  
  41. ***************************************************************************
  42. DPMPAR.F:
  43. ***************************************************************************
  44.  
  45.       DOUBLE PRECISION FUNCTION DPMPAR(I)
  46.       INTEGER I
  47. C     **********
  48. C
  49. C     FUNCTION DPMPAR
  50. C
  51. C     THIS FUNCTION PROVIDES DOUBLE PRECISION MACHINE PARAMETERS
  52. C     WHEN THE APPROPRIATE SET OF DATA STATEMENTS IS ACTIVATED (BY
  53. C     REMOVING THE C FROM COLUMN 1) AND ALL OTHER DATA STATEMENTS ARE
  54. C     RENDERED INACTIVE. MOST OF THE PARAMETER VALUES WERE OBTAINED
  55. C     FROM THE CORRESPONDING BELL LABORATORIES PORT LIBRARY FUNCTION.
  56. C
  57. C     THE FUNCTION STATEMENT IS
  58. C
  59. C       DOUBLE PRECISION FUNCTION DPMPAR(I)
  60. C
  61. C     WHERE
  62. C
  63. C       I IS AN INTEGER INPUT VARIABLE SET TO 1, 2, OR 3 WHICH
  64. C         SELECTS THE DESIRED MACHINE PARAMETER. IF THE MACHINE HAS
  65. C         T BASE B DIGITS AND ITS SMALLEST AND LARGEST EXPONENTS ARE
  66. C         EMIN AND EMAX, RESPECTIVELY, THEN THESE PARAMETERS ARE
  67. C
  68. C         DPMPAR(1) = B**(1 - T), THE MACHINE PRECISION,
  69. C
  70. C         DPMPAR(2) = B**(EMIN - 1), THE SMALLEST MAGNITUDE,
  71. C
  72. C         DPMPAR(3) = B**EMAX*(1 - B**(-T)), THE LARGEST MAGNITUDE.
  73. C
  74. C     ARGONNE NATIONAL LABORATORY. MINPACK PROJECT. MARCH 1980.
  75. C     BURTON S. GARBOW, KENNETH E. HILLSTROM, JORGE J. MORE
  76. C
  77. C     **********
  78.       INTEGER MCHEPS(4)
  79.       INTEGER MINMAG(4)
  80.       INTEGER MAXMAG(4)
  81.       DOUBLE PRECISION DMACH(3)
  82.       EQUIVALENCE (DMACH(1),MCHEPS(1))
  83.       EQUIVALENCE (DMACH(2),MINMAG(1))
  84.       EQUIVALENCE (DMACH(3),MAXMAG(1))
  85. C
  86. C     MACHINE CONSTANTS FOR THE VAX-11.
  87. C
  88. C     DATA MCHEPS(1),MCHEPS(2) /   9472,  0 /
  89. C     DATA MINMAG(1),MINMAG(2) /    128,  0 /
  90. C     DATA MAXMAG(1),MAXMAG(2) / -32769, -1 /
  91. C
  92. C     MACHINE CONSTANTS FOR IEEE ARITHMETIC MACHINES, SUCH AS THE AT&T
  93. C     3B SERIES AND MOTOROLA 68000 BASED MACHINES (E.G. SUN 3 AND AT&T
  94. C     PC 7300), IN WHICH THE MOST SIGNIFICANT BYTE IS STORED FIRST.
  95. C     Sun machine constants are provided by Dr.George F. Corliss (8/92)
  96. C     
  97.       DATA MCHEPS(1),MCHEPS(2) / 1018167296,          0 /
  98.       DATA MINMAG(1),MINMAG(2) / 2146435071,         -1 /
  99.       DATA MAXMAG(1),MAXMAG(2) /    1048576,          0 / 
  100. C
  101.       DPMPAR = DMACH(I)
  102.       RETURN
  103. C
  104. C     LAST CARD OF FUNCTION DPMPAR.
  105. C
  106.       END
  107.  
  108. ***************************************************************************
  109. ***************************************************************************
  110.  
  111.