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
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
|
*
*
PROGRAM PZGSEPDRIVER
*
* -- ScaLAPACK routine (version 1.7) --
* University of Tennessee, Knoxville, Oak Ridge National Laboratory,
* and University of California, Berkeley.
* May 1, 1997
*
* Parallel COMPLEX*16 Hermitian eigenproblem test driver
*
* The user should modify TOTMEM to indicate the maximum amount of
* memory in bytes her system has. Remember to leave room in memory
* for operating system, the BLACS buffer, etc. INTSIZ and DBLSIZ
* indicate the length in bytes on the given platform for an integer
* and a double precision real.
* For example, on our system with 8 MB of memory, TOTMEM=6500000
* (leaves 1.5 MB for OS, code, BLACS buffer, etc), the length of a
* DOUBLE is 8, and an integer takes up 4 bytes. Some playing around
* to discover what the maximum value you can set MEMSIZ to may be
* required.
* All arrays used by factorization and solve are allocated out of
* big array called MEM.
*
* The full tester requires approximately (5 n + 5 n^2/p + slop)
* COMPLEX*16 words and 6*n integer words.
* So, TOTMEM should be set to at least 1.1 * 8 * (5n + 5n^2/p)
*
* WHAT WE TEST
* ============
*
* This routine tests PZHEGVX, the expert driver for the parallel
* Hermitian eigenvalue problem. We would like to cover all
* possible combinations of: matrix size, process configuration
* (nprow and npcol), block size (nb), matrix type (??), range
* of eigenvalue (all, by value, by position), sorting options,
* and upper vs. lower storage.
*
* We intend to provide two types of test input files, an
* installation test and a thorough test.
*
* We also intend that the reports be meaningful. Our input file
* will allow multiple requests where each request is a cross product
* of the following sets:
* matrix sizes: n
* process configuration triples: nprow, npcol, nb
* matrix types:
* eigenvalue requests: all, by value, by position
* storage (upper vs. lower): uplo
*
* TERMS:
* Request - means a set of tests, which is the cross product of
* a set of specifications from the input file.
* Test - one element in the cross product, i.e. a specific input
* size and type, process configuration, etc.
*
* .. Parameters ..
*
INTEGER TOTMEM, ZPLXSZ, NIN
PARAMETER ( TOTMEM = 2000000, ZPLXSZ = 16, NIN = 11 )
INTEGER MEMSIZ
PARAMETER ( MEMSIZ = TOTMEM / ZPLXSZ )
* ..
* .. Local Scalars ..
CHARACTER HETERO
CHARACTER*80 SUMMRY, USRINFO
INTEGER CONTEXT, IAM, INFO, ISIEEE, MAXNODES, NNOCHECK,
$ NOUT, NPASSED, NPROCS, NSKIPPED, NTESTS
* ..
* .. Local Arrays ..
*
INTEGER ISEED( 4 )
COMPLEX*16 MEM( MEMSIZ )
* ..
* .. External Functions ..
DOUBLE PRECISION DLAMCH
EXTERNAL DLAMCH
* ..
* .. External Subroutines ..
*
EXTERNAL BLACS_EXIT, BLACS_GET, BLACS_GRIDEXIT,
$ BLACS_GRIDINIT, BLACS_PINFO, BLACS_SETUP,
$ IGAMN2D, PDLACHKIEEE, PDLASNBT, PZGSEPREQ
* ..
* .. Executable Statements ..
*
* Get starting information
*
CALL BLACS_PINFO( IAM, NPROCS )
*
*
IF( IAM.EQ.0 ) THEN
*
* Open file and skip data file header
*
OPEN( UNIT = NIN, FILE = 'SEP.dat', STATUS = 'OLD' )
READ( NIN, FMT = * )SUMMRY
SUMMRY = ' '
*
* Read in user-supplied info about machine type, compiler, etc.
*
READ( NIN, FMT = 9999 )USRINFO
*
* Read name and unit number for summary output file
*
READ( NIN, FMT = * )SUMMRY
READ( NIN, FMT = * )NOUT
IF( NOUT.NE.0 .AND. NOUT.NE.6 )
$ OPEN( NOUT, FILE = SUMMRY, STATUS = 'UNKNOWN' )
READ( NIN, FMT = * )MAXNODES
READ( NIN, FMT = * )HETERO
END IF
*
IF( NPROCS.LT.1 ) THEN
CALL BLACS_SETUP( IAM, MAXNODES )
NPROCS = MAXNODES
END IF
*
CALL BLACS_GET( -1, 0, CONTEXT )
CALL BLACS_GRIDINIT( CONTEXT, 'R', 1, NPROCS )
*
CALL PDLASNBT( ISIEEE )
*
CALL IGAMN2D( CONTEXT, 'a', ' ', 1, 1, ISIEEE, 1, 1, 1, -1, -1,
$ 0 )
*
IF( ( ISIEEE.NE.0 ) ) THEN
IF( IAM.EQ.0 ) THEN
WRITE( NOUT, FMT = 9998 )
WRITE( NOUT, FMT = 9997 )
WRITE( NOUT, FMT = 9996 )
WRITE( NOUT, FMT = 9995 )
WRITE( NOUT, FMT = 9994 )
WRITE( NOUT, FMT = 9993 )
WRITE( NOUT, FMT = 9992 )
WRITE( NOUT, FMT = 9991 )
WRITE( NOUT, FMT = 9990 )
END IF
*
CALL PDLACHKIEEE( ISIEEE, DLAMCH( 'O' ), DLAMCH( 'U' ) )
*
CALL IGAMN2D( CONTEXT, 'a', ' ', 1, 1, ISIEEE, 1, 1, 1, -1, -1,
$ 0 )
*
IF( ISIEEE.EQ.0 ) THEN
IF( IAM.EQ.0 ) THEN
WRITE( NOUT, FMT = 9989 )
WRITE( NOUT, FMT = 9988 )
WRITE( NOUT, FMT = 9987 )
END IF
GO TO 20
END IF
*
IF( IAM.EQ.0 ) THEN
WRITE( NOUT, FMT = 9986 )
END IF
*
END IF
IF( IAM.EQ.0 ) THEN
WRITE( NOUT, FMT = 9999 )
$ 'SCALAPACK Hermitian Eigendecomposition routines.'
WRITE( NOUT, FMT = 9999 )USRINFO
WRITE( NOUT, FMT = 9999 )' '
WRITE( NOUT, FMT = 9999 )'Running tests of the parallel ' //
$ 'generalized ' // 'Hermitian eigenvalue routine: PZHEGVX.'
WRITE( NOUT, FMT = 9999 )'A scaled residual check, ' //
$ 'will be computed'
WRITE( NOUT, FMT = 9999 )
WRITE( NOUT, FMT = 9999 )'An explanation of the ' //
$ 'input/output parameters follows:'
WRITE( NOUT, FMT = 9999 )'RESULT : passed; or ' //
$ 'an indication of which eigen request test failed'
WRITE( NOUT, FMT = 9999 )
$ 'N : The number of rows and columns ' //
$ 'of the matrix A.'
WRITE( NOUT, FMT = 9999 )
$ 'P : The number of process rows.'
WRITE( NOUT, FMT = 9999 )
$ 'Q : The number of process columns.'
WRITE( NOUT, FMT = 9999 )
$ 'NB : The size of the square blocks' //
$ ' the matrix A is split into.'
WRITE( NOUT, FMT = 9999 )
$ 'THRESH : If a residual value is less ' //
$ 'than THRESH, RESULT is flagged as PASSED.'
WRITE( NOUT, FMT = 9999 )
$ ' : the QTQ norm is allowed to exceed THRESH' //
$ ' for those eigenvectors'
WRITE( NOUT, FMT = 9999 )' : which could not be ' //
$ 'reorthogonalized for lack of workspace.'
WRITE( NOUT, FMT = 9999 )
$ 'TYP : matrix type (see pZGSEPtst.f).'
WRITE( NOUT, FMT = 9999 )
$ 'IBTYPE : Generalized eigenproblem type' //
$ ' (see pZHEGVx.f)'
WRITE( NOUT, FMT = 9999 )'SUB : Subtests ' //
$ '(see pZGSEPtst).f'
WRITE( NOUT, FMT = 9999 )'CHK : The scaled residual'
WRITE( NOUT, FMT = 9999 )' '
END IF
*
NTESTS = 0
NPASSED = 0
NSKIPPED = 0
NNOCHECK = 0
*
IF( IAM.EQ.0 ) THEN
WRITE( NOUT, FMT = 9979 )
WRITE( NOUT, FMT = 9978 )
END IF
*
10 CONTINUE
*
ISEED( 1 ) = 139
ISEED( 2 ) = 1139
ISEED( 3 ) = 2139
ISEED( 4 ) = 3139
*
CALL PZGSEPREQ( NIN, MEM, MEMSIZ, NOUT, ISEED, NTESTS, NSKIPPED,
$ NNOCHECK, NPASSED, INFO )
IF( INFO.EQ.0 )
$ GO TO 10
*
IF( IAM.EQ.0 ) THEN
WRITE( NOUT, FMT = 9985 )NTESTS
WRITE( NOUT, FMT = 9984 )NPASSED
WRITE( NOUT, FMT = 9983 )NNOCHECK
WRITE( NOUT, FMT = 9982 )NSKIPPED
WRITE( NOUT, FMT = 9981 )NTESTS - NPASSED - NSKIPPED -
$ NNOCHECK
WRITE( NOUT, FMT = * )
WRITE( NOUT, FMT = * )
WRITE( NOUT, FMT = 9980 )
END IF
*
* Uncomment this line on SUN systems to avoid the useless print out
*
* CALL IEEE_FLAGS( 'clear', 'exception', 'underflow', '')
*
*
*
20 CONTINUE
IF( IAM.EQ.0 ) THEN
CLOSE ( NIN )
IF( NOUT.NE.6 .AND. NOUT.NE.0 )
$ CLOSE ( NOUT )
END IF
*
CALL BLACS_GRIDEXIT( CONTEXT )
*
CALL BLACS_EXIT( 0 )
STOP
*
*
9999 FORMAT( A )
9998 FORMAT( ' I am about to check to make sure that overflow' )
9997 FORMAT( ' is handled in the ieee default manner. If this' )
9996 FORMAT( ' is the last output you see, you should assume' )
9995 FORMAT( ' that overflow caused a floating point exception.' )
9994 FORMAT( ' In that case, we recommend that you add -DNO_IEEE' )
9993 FORMAT( ' to the CDEFS line in SLmake.inc.' )
9992 FORMAT( ' Alternatively, you could set CDEFS in SLmake.inc ' )
9991 FORMAT( ' to enable the default ieee behaviour, However, this' )
9990 FORMAT( ' may result in good or very bad performance.' )
9989 FORMAT( ' Either signed zeroes or signed infinities ' )
9988 FORMAT( ' work incorrectly or your system. Change your' )
9987 FORMAT( ' SLmake.inc as suggested above.' )
*
9986 FORMAT( ' Your system appears to handle ieee overflow.' )
*
9985 FORMAT( 'Finished ', I6, ' tests, with the following results:' )
9984 FORMAT( I5, ' tests completed and passed residual checks.' )
9983 FORMAT( I5, ' tests completed without checking.' )
9982 FORMAT( I5, ' tests skipped for lack of memory.' )
9981 FORMAT( I5, ' tests completed and failed.' )
9980 FORMAT( 'END OF TESTS.' )
9979 FORMAT( ' N NB P Q TYP IBTYPE SUB WALL CPU ',
$ ' CHK CHECK' )
9978 FORMAT( ' ----- --- --- --- --- ------ --- -------- --------',
$ ' --------- -----' )
*
* End of PZGSEPDRIVER
*
END
|
