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
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
|
*DMPREC
DOUBLE PRECISION FUNCTION DMPREC()
C***BEGIN PROLOGUE DPREC
C***REFER TO DODR,DODRC
C***ROUTINES CALLED (NONE)
C***DATE WRITTEN 860529 (YYMMDD)
C***REVISION DATE 920304 (YYMMDD)
C***PURPOSE DETERMINE MACHINE PRECISION FOR TARGET MACHINE AND COMPILER
C ASSUMING FLOATING-POINT NUMBERS ARE REPRESENTED IN THE
C T-DIGIT, BASE-B FORM
C SIGN (B**E)*( (X(1)/B) + ... + (X(T)/B**T) )
C WHERE 0 .LE. X(I) .LT. B FOR I=1,...,T, AND
C 0 .LT. X(1).
C TO ALTER THIS FUNCTION FOR A PARTICULAR TARGET MACHINE,
C EITHER
C ACTIVATE THE DESIRED SET OF DATA STATEMENTS BY
C REMOVING THE C FROM COLUMN 1
C OR
C SET B, TD AND TS USING I1MACH BY ACTIVATING
C THE DECLARATION STATEMENTS FOR I1MACH
C AND THE STATEMENTS PRECEEDING THE FIRST
C EXECUTABLE STATEMENT BELOW.
C***END PROLOGUE DPREC
C...LOCAL SCALARS
DOUBLE PRECISION
+ B
INTEGER
+ TD,TS
C...EXTERNAL FUNCTIONS
C INTEGER
C + I1MACH
C EXTERNAL
C + I1MACH
C...VARIABLE DEFINITIONS (ALPHABETICALLY)
C DOUBLE PRECISION B
C THE BASE OF THE TARGET MACHINE.
C (MAY BE DEFINED USING I1MACH(10).)
C INTEGER TD
C THE NUMBER OF BASE-B DIGITS IN DOUBLE PRECISION.
C (MAY BE DEFINED USING I1MACH(14).)
C INTEGER TS
C THE NUMBER OF BASE-B DIGITS IN SINGLE PRECISION.
C (MAY BE DEFINED USING I1MACH(11).)
C MACHINE CONSTANTS FOR COMPUTERS FOLLOWING IEEE ARITHMETIC STANDARD
C (E.G., MOTOROLA 68000 BASED MACHINES SUCH AS SUN AND SPARC
C WORKSTATIONS, AND AT&T PC 7300; AND 8087 BASED MICROS SUCH AS THE
C IBM PC AND THE AT&T 6300).
DATA B / 2 /
DATA TS / 24 /
DATA TD / 53 /
C MACHINE CONSTANTS FOR THE BURROUGHS 1700 SYSTEM.
C DATA B / 2 /
C DATA TS / 24 /
C DATA TD / 60 /
C MACHINE CONSTANTS FOR THE BURROUGHS 5700 SYSTEM
C THE BURROUGHS 6700/7700 SYSTEMS
C DATA B / 8 /
C DATA TS / 13 /
C DATA TD / 26 /
C MACHINE CONSTANTS FOR THE CDC 6000/7000 (FTN5 COMPILER)
C THE CYBER 170/180 SERIES UNDER NOS
C DATA B / 2 /
C DATA TS / 48 /
C DATA TD / 96 /
C MACHINE CONSTANTS FOR THE CDC 6000/7000 (FTN COMPILER)
C THE CYBER 170/180 SERIES UNDER NOS/VE
C THE CYBER 200 SERIES
C DATA B / 2 /
C DATA TS / 47 /
C DATA TD / 94 /
C MACHINE CONSTANTS FOR THE CRAY
C DATA B / 2 /
C DATA TS / 47 /
C DATA TD / 94 /
C MACHINE CONSTANTS FOR THE DATA GENERAL ECLIPSE S/200
C DATA B / 16 /
C DATA TS / 6 /
C DATA TD / 14 /
C MACHINE CONSTANTS FOR THE HARRIS COMPUTER
C DATA B / 2 /
C DATA TS / 23 /
C DATA TD / 38 /
C MACHINE CONSTANTS FOR THE HONEYWELL DPS 8/70
C THE HONEYWELL 600/6000 SERIES
C DATA B / 2 /
C DATA TS / 27 /
C DATA TD / 63 /
C MACHINE CONSTANTS FOR THE HP 2100
C (3 WORD DOUBLE PRECISION OPTION WITH FTN4)
C DATA B / 2 /
C DATA TS / 23 /
C DATA TD / 39 /
C MACHINE CONSTANTS FOR THE HP 2100
C (4 WORD DOUBLE PRECISION OPTION WITH FTN4)
C DATA B / 2 /
C DATA TS / 23 /
C DATA TD / 55 /
C MACHINE CONSTANTS FOR THE IBM 360/370 SERIES
C DATA B / 16 /
C DATA TS / 6 /
C DATA TD / 14 /
C MACHINE CONSTANTS FOR THE IBM PC
C DATA B / 2 /
C DATA TS / 24 /
C DATA TD / 53 /
C MACHINE CONSTANTS FOR THE INTERDATA (PERKIN ELMER) 7/32
C INTERDATA (PERKIN ELMER) 8/32
C DATA B / 16 /
C DATA TS / 6 /
C DATA TD / 14 /
C MACHINE CONSTANTS FOR THE PDP-10 (KA PROCESSOR).
C DATA B / 2 /
C DATA TS / 27 /
C DATA TD / 54 /
C MACHINE CONSTANTS FOR THE PDP-10 (KI PROCESSOR).
C DATA B / 2 /
C DATA TS / 27 /
C DATA TD / 62 /
C MACHINE CONSTANTS FOR THE PDP-11 SYSTEM
C DATA B / 2 /
C DATA TS / 24 /
C DATA TD / 56 /
C MACHINE CONSTANTS FOR THE PERKIN-ELMER 3230
C DATA B / 16 /
C DATA TS / 6 /
C DATA TD / 14 /
C MACHINE CONSTANTS FOR THE PRIME 850 AND PRIME 4050
C DATA B / 2 /
C DATA TS / 23 /
C DATA TD / 47 /
C MACHINE CONSTANTS FOR THE SEL SYSTEMS 85/86
C DATA B / 16 /
C DATA TS / 6 /
C DATA TD / 14 /
C MACHINE CONSTANTS FOR SUN AND SPARC WORKSTATIONS
C DATA B / 2 /
C DATA TS / 24 /
C DATA TD / 53 /
C MACHINE CONSTANTS FOR THE UNIVAC 1100 SERIES
C DATA B / 2 /
C DATA TS / 27 /
C DATA TD / 60 /
C MACHINE CONSTANTS FOR THE VAX-11 WITH FORTRAN IV-PLUS COMPILER
C DATA B / 2 /
C DATA TS / 24 /
C DATA TD / 56 /
C MACHINE CONSTANTS FOR THE VAX/VMS SYSTEM WITHOUT G_FLOATING
C DATA B / 2 /
C DATA TS / 24 /
C DATA TD / 56 /
C MACHINE CONSTANTS FOR THE VAX/VMS SYSTEM WITH G_FLOATING
C DATA B / 2 /
C DATA TS / 24 /
C DATA TD / 53 /
C MACHINE CONSTANTS FOR THE XEROX SIGMA 5/7/9
C DATA B / 16 /
C DATA TS / 6 /
C DATA TD / 14 /
C***FIRST EXECUTABLE STATEMENT DMPREC
C B = I1MACH(10)
C TS = I1MACH(11)
C TD = I1MACH(14)
DMPREC = B ** (1-TD)
RETURN
END
|
