1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145
|
REAL FUNCTION SPMPAR(I)
INTEGER I
C **********
C
C FUNCTION SPMPAR
C
C THIS FUNCTION PROVIDES SINGLE PRECISION MACHINE PARAMETERS
C WHEN THE APPROPRIATE SET OF DATA STATEMENTS IS ACTIVATED (BY
C REMOVING THE C FROM COLUMN 1) AND ALL OTHER DATA STATEMENTS ARE
C RENDERED INACTIVE. MOST OF THE PARAMETER VALUES WERE OBTAINED
C FROM THE CORRESPONDING BELL LABORATORIES PORT LIBRARY FUNCTION.
C
C THE FUNCTION STATEMENT IS
C
C REAL FUNCTION SPMPAR(I)
C
C WHERE
C
C I IS AN INTEGER INPUT VARIABLE SET TO 1, 2, OR 3 WHICH
C SELECTS THE DESIRED MACHINE PARAMETER. IF THE MACHINE HAS
C T BASE B DIGITS AND ITS SMALLEST AND LARGEST EXPONENTS ARE
C EMIN AND EMAX, RESPECTIVELY, THEN THESE PARAMETERS ARE
C
C SPMPAR(1) = B**(1 - T), THE MACHINE PRECISION,
C
C SPMPAR(2) = B**(EMIN - 1), THE SMALLEST MAGNITUDE,
C
C SPMPAR(3) = B**EMAX*(1 - B**(-T)), THE LARGEST MAGNITUDE.
C
C ARGONNE NATIONAL LABORATORY. MINPACK PROJECT. MARCH 1980.
C BURTON S. GARBOW, KENNETH E. HILLSTROM, JORGE J. MORE
C
C **********
INTEGER MCHEPS(2)
INTEGER MINMAG(2)
INTEGER MAXMAG(2)
REAL RMACH(3)
EQUIVALENCE (RMACH(1),MCHEPS(1))
EQUIVALENCE (RMACH(2),MINMAG(1))
EQUIVALENCE (RMACH(3),MAXMAG(1))
C
C MACHINE CONSTANTS FOR THE IBM 360/370 SERIES,
C THE AMDAHL 470/V6, THE ICL 2900, THE ITEL AS/6,
C THE XEROX SIGMA 5/7/9 AND THE SEL SYSTEMS 85/86.
C
DATA RMACH(1) / Z3C100000 /
DATA RMACH(2) / Z00100000 /
DATA RMACH(3) / Z7FFFFFFF /
C
C MACHINE CONSTANTS FOR THE HONEYWELL 600/6000 SERIES.
C
C DATA RMACH(1) / O716400000000 /
C DATA RMACH(2) / O402400000000 /
C DATA RMACH(3) / O376777777777 /
C
C MACHINE CONSTANTS FOR THE CDC 6000/7000 SERIES.
C
C DATA RMACH(1) / 16414000000000000000B /
C DATA RMACH(2) / 00014000000000000000B /
C DATA RMACH(3) / 37767777777777777777B /
C
C MACHINE CONSTANTS FOR THE PDP-10 (KA OR KI PROCESSOR).
C
C DATA RMACH(1) / "147400000000 /
C DATA RMACH(2) / "000400000000 /
C DATA RMACH(3) / "377777777777 /
C
C MACHINE CONSTANTS FOR THE PDP-11 FORTRAN SUPPORTING
C 32-BIT INTEGERS (EXPRESSED IN INTEGER AND OCTAL).
C
C DATA MCHEPS(1) / 889192448 /
C DATA MINMAG(1) / 8388608 /
C DATA MAXMAG(1) / 2147483647 /
C
C DATA RMACH(1) / O06500000000 /
C DATA RMACH(2) / O00040000000 /
C DATA RMACH(3) / O17777777777 /
C
C MACHINE CONSTANTS FOR THE PDP-11 FORTRAN SUPPORTING
C 16-BIT INTEGERS (EXPRESSED IN INTEGER AND OCTAL).
C
C DATA MCHEPS(1),MCHEPS(2) / 13568, 0 /
C DATA MINMAG(1),MINMAG(2) / 128, 0 /
C DATA MAXMAG(1),MAXMAG(2) / 32767, -1 /
C
C DATA MCHEPS(1),MCHEPS(2) / O032400, O000000 /
C DATA MINMAG(1),MINMAG(2) / O000200, O000000 /
C DATA MAXMAG(1),MAXMAG(2) / O077777, O177777 /
C
C MACHINE CONSTANTS FOR THE BURROUGHS 5700/6700/7700 SYSTEMS.
C
C DATA RMACH(1) / O1301000000000000 /
C DATA RMACH(2) / O1771000000000000 /
C DATA RMACH(3) / O0777777777777777 /
C
C MACHINE CONSTANTS FOR THE BURROUGHS 1700 SYSTEM.
C
C DATA RMACH(1) / Z4EA800000 /
C DATA RMACH(2) / Z400800000 /
C DATA RMACH(3) / Z5FFFFFFFF /
C
C MACHINE CONSTANTS FOR THE UNIVAC 1100 SERIES.
C
C DATA RMACH(1) / O147400000000 /
C DATA RMACH(2) / O000400000000 /
C DATA RMACH(3) / O377777777777 /
C
C MACHINE CONSTANTS FOR THE DATA GENERAL ECLIPSE S/200.
C
C NOTE - IT MAY BE APPROPRIATE TO INCLUDE THE FOLLOWING CARD -
C STATIC RMACH(3)
C
C DATA MINMAG/20K,0/,MAXMAG/77777K,177777K/
C DATA MCHEPS/36020K,0/
C
C MACHINE CONSTANTS FOR THE HARRIS 220.
C
C DATA MCHEPS(1) / '20000000, '00000353 /
C DATA MINMAG(1) / '20000000, '00000201 /
C DATA MAXMAG(1) / '37777777, '00000177 /
C
C MACHINE CONSTANTS FOR THE CRAY-1.
C
C DATA RMACH(1) / 0377224000000000000000B /
C DATA RMACH(2) / 0200034000000000000000B /
C DATA RMACH(3) / 0577777777777777777776B /
C
C MACHINE CONSTANTS FOR THE PRIME 400.
C
C DATA MCHEPS(1) / :10000000153 /
C DATA MINMAG(1) / :10000000000 /
C DATA MAXMAG(1) / :17777777777 /
C
C MACHINE CONSTANTS FOR THE VAX-11.
C
C DATA MCHEPS(1) / 13568 /
C DATA MINMAG(1) / 128 /
C DATA MAXMAG(1) / -32769 /
C
SPMPAR = RMACH(I)
RETURN
C
C LAST CARD OF FUNCTION SPMPAR.
C
END
|