1 PROGRAM ZBLAT2
2 *
3 * Test program for the COMPLEX*16 Level 2 Blas.
4 *
5 * The program must be driven by a short data file. The first 18 records
6 * of the file are read using list-directed input, the last 17 records
7 * are read using the format ( A6, L2 ). An annotated example of a data
8 * file can be obtained by deleting the first 3 characters from the
9 * following 35 lines:
10 * 'ZBLAT2.SUMM' NAME OF SUMMARY OUTPUT FILE
11 * 6 UNIT NUMBER OF SUMMARY FILE
12 * 'CBLA2T.SNAP' NAME OF SNAPSHOT OUTPUT FILE
13 * -1 UNIT NUMBER OF SNAPSHOT FILE (NOT USED IF .LT. 0)
14 * F LOGICAL FLAG, T TO REWIND SNAPSHOT FILE AFTER EACH RECORD.
15 * F LOGICAL FLAG, T TO STOP ON FAILURES.
16 * T LOGICAL FLAG, T TO TEST ERROR EXITS.
17 * 16.0 THRESHOLD VALUE OF TEST RATIO
18 * 6 NUMBER OF VALUES OF N
19 * 0 1 2 3 5 9 VALUES OF N
20 * 4 NUMBER OF VALUES OF K
21 * 0 1 2 4 VALUES OF K
22 * 4 NUMBER OF VALUES OF INCX AND INCY
23 * 1 2 -1 -2 VALUES OF INCX AND INCY
24 * 3 NUMBER OF VALUES OF ALPHA
25 * (0.0,0.0) (1.0,0.0) (0.7,-0.9) VALUES OF ALPHA
26 * 3 NUMBER OF VALUES OF BETA
27 * (0.0,0.0) (1.0,0.0) (1.3,-1.1) VALUES OF BETA
28 * ZGEMV T PUT F FOR NO TEST. SAME COLUMNS.
29 * ZGBMV T PUT F FOR NO TEST. SAME COLUMNS.
30 * ZHEMV T PUT F FOR NO TEST. SAME COLUMNS.
31 * ZHBMV T PUT F FOR NO TEST. SAME COLUMNS.
32 * ZHPMV T PUT F FOR NO TEST. SAME COLUMNS.
33 * ZTRMV T PUT F FOR NO TEST. SAME COLUMNS.
34 * ZTBMV T PUT F FOR NO TEST. SAME COLUMNS.
35 * ZTPMV T PUT F FOR NO TEST. SAME COLUMNS.
36 * ZTRSV T PUT F FOR NO TEST. SAME COLUMNS.
37 * ZTBSV T PUT F FOR NO TEST. SAME COLUMNS.
38 * ZTPSV T PUT F FOR NO TEST. SAME COLUMNS.
39 * ZGERC T PUT F FOR NO TEST. SAME COLUMNS.
40 * ZGERU T PUT F FOR NO TEST. SAME COLUMNS.
41 * ZHER T PUT F FOR NO TEST. SAME COLUMNS.
42 * ZHPR T PUT F FOR NO TEST. SAME COLUMNS.
43 * ZHER2 T PUT F FOR NO TEST. SAME COLUMNS.
44 * ZHPR2 T PUT F FOR NO TEST. SAME COLUMNS.
45 *
46 * See:
47 *
48 * Dongarra J. J., Du Croz J. J., Hammarling S. and Hanson R. J..
49 * An extended set of Fortran Basic Linear Algebra Subprograms.
50 *
51 * Technical Memoranda Nos. 41 (revision 3) and 81, Mathematics
52 * and Computer Science Division, Argonne National Laboratory,
53 * 9700 South Cass Avenue, Argonne, Illinois 60439, US.
54 *
55 * Or
56 *
57 * NAG Technical Reports TR3/87 and TR4/87, Numerical Algorithms
58 * Group Ltd., NAG Central Office, 256 Banbury Road, Oxford
59 * OX2 7DE, UK, and Numerical Algorithms Group Inc., 1101 31st
60 * Street, Suite 100, Downers Grove, Illinois 60515-1263, USA.
61 *
62 *
63 * -- Written on 10-August-1987.
64 * Richard Hanson, Sandia National Labs.
65 * Jeremy Du Croz, NAG Central Office.
66 *
67 * .. Parameters ..
68 INTEGER NIN
69 PARAMETER ( NIN = 5 )
70 INTEGER NSUBS
71 PARAMETER ( NSUBS = 17 )
72 COMPLEX*16 ZERO, ONE
73 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
74 $ ONE = ( 1.0D0, 0.0D0 ) )
75 DOUBLE PRECISION RZERO, RHALF, RONE
76 PARAMETER ( RZERO = 0.0D0, RHALF = 0.5D0, RONE = 1.0D0 )
77 INTEGER NMAX, INCMAX
78 PARAMETER ( NMAX = 65, INCMAX = 2 )
79 INTEGER NINMAX, NIDMAX, NKBMAX, NALMAX, NBEMAX
80 PARAMETER ( NINMAX = 7, NIDMAX = 9, NKBMAX = 7,
81 $ NALMAX = 7, NBEMAX = 7 )
82 * .. Local Scalars ..
83 DOUBLE PRECISION EPS, ERR, THRESH
84 INTEGER I, ISNUM, J, N, NALF, NBET, NIDIM, NINC, NKB,
85 $ NOUT, NTRA
86 LOGICAL FATAL, LTESTT, REWI, SAME, SFATAL, TRACE,
87 $ TSTERR
88 CHARACTER*1 TRANS
89 CHARACTER*6 SNAMET
90 CHARACTER*32 SNAPS, SUMMRY
91 * .. Local Arrays ..
92 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ),
93 $ ALF( NALMAX ), AS( NMAX*NMAX ), BET( NBEMAX ),
94 $ X( NMAX ), XS( NMAX*INCMAX ),
95 $ XX( NMAX*INCMAX ), Y( NMAX ),
96 $ YS( NMAX*INCMAX ), YT( NMAX ),
97 $ YY( NMAX*INCMAX ), Z( 2*NMAX )
98 DOUBLE PRECISION G( NMAX )
99 INTEGER IDIM( NIDMAX ), INC( NINMAX ), KB( NKBMAX )
100 LOGICAL LTEST( NSUBS )
101 CHARACTER*6 SNAMES( NSUBS )
102 * .. External Functions ..
103 DOUBLE PRECISION DDIFF
104 LOGICAL LZE
105 EXTERNAL DDIFF, LZE
106 * .. External Subroutines ..
107 EXTERNAL ZCHK1, ZCHK2, ZCHK3, ZCHK4, ZCHK5, ZCHK6,
108 $ ZCHKE, ZMVCH
109 * .. Intrinsic Functions ..
110 INTRINSIC ABS, MAX, MIN
111 * .. Scalars in Common ..
112 INTEGER INFOT, NOUTC
113 LOGICAL LERR, OK
114 CHARACTER*6 SRNAMT
115 * .. Common blocks ..
116 COMMON /INFOC/INFOT, NOUTC, OK, LERR
117 COMMON /SRNAMC/SRNAMT
118 * .. Data statements ..
119 DATA SNAMES/'ZGEMV ', 'ZGBMV ', 'ZHEMV ', 'ZHBMV ',
120 $ 'ZHPMV ', 'ZTRMV ', 'ZTBMV ', 'ZTPMV ',
121 $ 'ZTRSV ', 'ZTBSV ', 'ZTPSV ', 'ZGERC ',
122 $ 'ZGERU ', 'ZHER ', 'ZHPR ', 'ZHER2 ',
123 $ 'ZHPR2 '/
124 * .. Executable Statements ..
125 *
126 * Read name and unit number for summary output file and open file.
127 *
128 READ( NIN, FMT = * )SUMMRY
129 READ( NIN, FMT = * )NOUT
130 OPEN( NOUT, FILE = SUMMRY, STATUS = 'NEW' )
131 NOUTC = NOUT
132 *
133 * Read name and unit number for snapshot output file and open file.
134 *
135 READ( NIN, FMT = * )SNAPS
136 READ( NIN, FMT = * )NTRA
137 TRACE = NTRA.GE.0
138 IF( TRACE )THEN
139 OPEN( NTRA, FILE = SNAPS, STATUS = 'NEW' )
140 END IF
141 * Read the flag that directs rewinding of the snapshot file.
142 READ( NIN, FMT = * )REWI
143 REWI = REWI.AND.TRACE
144 * Read the flag that directs stopping on any failure.
145 READ( NIN, FMT = * )SFATAL
146 * Read the flag that indicates whether error exits are to be tested.
147 READ( NIN, FMT = * )TSTERR
148 * Read the threshold value of the test ratio
149 READ( NIN, FMT = * )THRESH
150 *
151 * Read and check the parameter values for the tests.
152 *
153 * Values of N
154 READ( NIN, FMT = * )NIDIM
155 IF( NIDIM.LT.1.OR.NIDIM.GT.NIDMAX )THEN
156 WRITE( NOUT, FMT = 9997 )'N', NIDMAX
157 GO TO 230
158 END IF
159 READ( NIN, FMT = * )( IDIM( I ), I = 1, NIDIM )
160 DO 10 I = 1, NIDIM
161 IF( IDIM( I ).LT.0.OR.IDIM( I ).GT.NMAX )THEN
162 WRITE( NOUT, FMT = 9996 )NMAX
163 GO TO 230
164 END IF
165 10 CONTINUE
166 * Values of K
167 READ( NIN, FMT = * )NKB
168 IF( NKB.LT.1.OR.NKB.GT.NKBMAX )THEN
169 WRITE( NOUT, FMT = 9997 )'K', NKBMAX
170 GO TO 230
171 END IF
172 READ( NIN, FMT = * )( KB( I ), I = 1, NKB )
173 DO 20 I = 1, NKB
174 IF( KB( I ).LT.0 )THEN
175 WRITE( NOUT, FMT = 9995 )
176 GO TO 230
177 END IF
178 20 CONTINUE
179 * Values of INCX and INCY
180 READ( NIN, FMT = * )NINC
181 IF( NINC.LT.1.OR.NINC.GT.NINMAX )THEN
182 WRITE( NOUT, FMT = 9997 )'INCX AND INCY', NINMAX
183 GO TO 230
184 END IF
185 READ( NIN, FMT = * )( INC( I ), I = 1, NINC )
186 DO 30 I = 1, NINC
187 IF( INC( I ).EQ.0.OR.ABS( INC( I ) ).GT.INCMAX )THEN
188 WRITE( NOUT, FMT = 9994 )INCMAX
189 GO TO 230
190 END IF
191 30 CONTINUE
192 * Values of ALPHA
193 READ( NIN, FMT = * )NALF
194 IF( NALF.LT.1.OR.NALF.GT.NALMAX )THEN
195 WRITE( NOUT, FMT = 9997 )'ALPHA', NALMAX
196 GO TO 230
197 END IF
198 READ( NIN, FMT = * )( ALF( I ), I = 1, NALF )
199 * Values of BETA
200 READ( NIN, FMT = * )NBET
201 IF( NBET.LT.1.OR.NBET.GT.NBEMAX )THEN
202 WRITE( NOUT, FMT = 9997 )'BETA', NBEMAX
203 GO TO 230
204 END IF
205 READ( NIN, FMT = * )( BET( I ), I = 1, NBET )
206 *
207 * Report values of parameters.
208 *
209 WRITE( NOUT, FMT = 9993 )
210 WRITE( NOUT, FMT = 9992 )( IDIM( I ), I = 1, NIDIM )
211 WRITE( NOUT, FMT = 9991 )( KB( I ), I = 1, NKB )
212 WRITE( NOUT, FMT = 9990 )( INC( I ), I = 1, NINC )
213 WRITE( NOUT, FMT = 9989 )( ALF( I ), I = 1, NALF )
214 WRITE( NOUT, FMT = 9988 )( BET( I ), I = 1, NBET )
215 IF( .NOT.TSTERR )THEN
216 WRITE( NOUT, FMT = * )
217 WRITE( NOUT, FMT = 9980 )
218 END IF
219 WRITE( NOUT, FMT = * )
220 WRITE( NOUT, FMT = 9999 )THRESH
221 WRITE( NOUT, FMT = * )
222 *
223 * Read names of subroutines and flags which indicate
224 * whether they are to be tested.
225 *
226 DO 40 I = 1, NSUBS
227 LTEST( I ) = .FALSE.
228 40 CONTINUE
229 50 READ( NIN, FMT = 9984, END = 80 )SNAMET, LTESTT
230 DO 60 I = 1, NSUBS
231 IF( SNAMET.EQ.SNAMES( I ) )
232 $ GO TO 70
233 60 CONTINUE
234 WRITE( NOUT, FMT = 9986 )SNAMET
235 STOP
236 70 LTEST( I ) = LTESTT
237 GO TO 50
238 *
239 80 CONTINUE
240 CLOSE ( NIN )
241 *
242 * Compute EPS (the machine precision).
243 *
244 EPS = RONE
245 90 CONTINUE
246 IF( DDIFF( RONE + EPS, RONE ).EQ.RZERO )
247 $ GO TO 100
248 EPS = RHALF*EPS
249 GO TO 90
250 100 CONTINUE
251 EPS = EPS + EPS
252 WRITE( NOUT, FMT = 9998 )EPS
253 *
254 * Check the reliability of ZMVCH using exact data.
255 *
256 N = MIN( 32, NMAX )
257 DO 120 J = 1, N
258 DO 110 I = 1, N
259 A( I, J ) = MAX( I - J + 1, 0 )
260 110 CONTINUE
261 X( J ) = J
262 Y( J ) = ZERO
263 120 CONTINUE
264 DO 130 J = 1, N
265 YY( J ) = J*( ( J + 1 )*J )/2 - ( ( J + 1 )*J*( J - 1 ) )/3
266 130 CONTINUE
267 * YY holds the exact result. On exit from ZMVCH YT holds
268 * the result computed by ZMVCH.
269 TRANS = 'N'
270 CALL ZMVCH( TRANS, N, N, ONE, A, NMAX, X, 1, ZERO, Y, 1, YT, G,
271 $ YY, EPS, ERR, FATAL, NOUT, .TRUE. )
272 SAME = LZE( YY, YT, N )
273 IF( .NOT.SAME.OR.ERR.NE.RZERO )THEN
274 WRITE( NOUT, FMT = 9985 )TRANS, SAME, ERR
275 STOP
276 END IF
277 TRANS = 'T'
278 CALL ZMVCH( TRANS, N, N, ONE, A, NMAX, X, -1, ZERO, Y, -1, YT, G,
279 $ YY, EPS, ERR, FATAL, NOUT, .TRUE. )
280 SAME = LZE( YY, YT, N )
281 IF( .NOT.SAME.OR.ERR.NE.RZERO )THEN
282 WRITE( NOUT, FMT = 9985 )TRANS, SAME, ERR
283 STOP
284 END IF
285 *
286 * Test each subroutine in turn.
287 *
288 DO 210 ISNUM = 1, NSUBS
289 WRITE( NOUT, FMT = * )
290 IF( .NOT.LTEST( ISNUM ) )THEN
291 * Subprogram is not to be tested.
292 WRITE( NOUT, FMT = 9983 )SNAMES( ISNUM )
293 ELSE
294 SRNAMT = SNAMES( ISNUM )
295 * Test error exits.
296 IF( TSTERR )THEN
297 CALL ZCHKE( ISNUM, SNAMES( ISNUM ), NOUT )
298 WRITE( NOUT, FMT = * )
299 END IF
300 * Test computations.
301 INFOT = 0
302 OK = .TRUE.
303 FATAL = .FALSE.
304 GO TO ( 140, 140, 150, 150, 150, 160, 160,
305 $ 160, 160, 160, 160, 170, 170, 180,
306 $ 180, 190, 190 )ISNUM
307 * Test ZGEMV, 01, and ZGBMV, 02.
308 140 CALL ZCHK1( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
309 $ REWI, FATAL, NIDIM, IDIM, NKB, KB, NALF, ALF,
310 $ NBET, BET, NINC, INC, NMAX, INCMAX, A, AA, AS,
311 $ X, XX, XS, Y, YY, YS, YT, G )
312 GO TO 200
313 * Test ZHEMV, 03, ZHBMV, 04, and ZHPMV, 05.
314 150 CALL ZCHK2( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
315 $ REWI, FATAL, NIDIM, IDIM, NKB, KB, NALF, ALF,
316 $ NBET, BET, NINC, INC, NMAX, INCMAX, A, AA, AS,
317 $ X, XX, XS, Y, YY, YS, YT, G )
318 GO TO 200
319 * Test ZTRMV, 06, ZTBMV, 07, ZTPMV, 08,
320 * ZTRSV, 09, ZTBSV, 10, and ZTPSV, 11.
321 160 CALL ZCHK3( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
322 $ REWI, FATAL, NIDIM, IDIM, NKB, KB, NINC, INC,
323 $ NMAX, INCMAX, A, AA, AS, Y, YY, YS, YT, G, Z )
324 GO TO 200
325 * Test ZGERC, 12, ZGERU, 13.
326 170 CALL ZCHK4( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
327 $ REWI, FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC,
328 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
329 $ YT, G, Z )
330 GO TO 200
331 * Test ZHER, 14, and ZHPR, 15.
332 180 CALL ZCHK5( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
333 $ REWI, FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC,
334 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
335 $ YT, G, Z )
336 GO TO 200
337 * Test ZHER2, 16, and ZHPR2, 17.
338 190 CALL ZCHK6( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
339 $ REWI, FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC,
340 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
341 $ YT, G, Z )
342 *
343 200 IF( FATAL.AND.SFATAL )
344 $ GO TO 220
345 END IF
346 210 CONTINUE
347 WRITE( NOUT, FMT = 9982 )
348 GO TO 240
349 *
350 220 CONTINUE
351 WRITE( NOUT, FMT = 9981 )
352 GO TO 240
353 *
354 230 CONTINUE
355 WRITE( NOUT, FMT = 9987 )
356 *
357 240 CONTINUE
358 IF( TRACE )
359 $ CLOSE ( NTRA )
360 CLOSE ( NOUT )
361 STOP
362 *
363 9999 FORMAT( ' ROUTINES PASS COMPUTATIONAL TESTS IF TEST RATIO IS LES',
364 $ 'S THAN', F8.2 )
365 9998 FORMAT( ' RELATIVE MACHINE PRECISION IS TAKEN TO BE', 1P, D9.1 )
366 9997 FORMAT( ' NUMBER OF VALUES OF ', A, ' IS LESS THAN 1 OR GREATER ',
367 $ 'THAN ', I2 )
368 9996 FORMAT( ' VALUE OF N IS LESS THAN 0 OR GREATER THAN ', I2 )
369 9995 FORMAT( ' VALUE OF K IS LESS THAN 0' )
370 9994 FORMAT( ' ABSOLUTE VALUE OF INCX OR INCY IS 0 OR GREATER THAN ',
371 $ I2 )
372 9993 FORMAT( ' TESTS OF THE COMPLEX*16 LEVEL 2 BLAS', //' THE F',
373 $ 'OLLOWING PARAMETER VALUES WILL BE USED:' )
374 9992 FORMAT( ' FOR N ', 9I6 )
375 9991 FORMAT( ' FOR K ', 7I6 )
376 9990 FORMAT( ' FOR INCX AND INCY ', 7I6 )
377 9989 FORMAT( ' FOR ALPHA ',
378 $ 7( '(', F4.1, ',', F4.1, ') ', : ) )
379 9988 FORMAT( ' FOR BETA ',
380 $ 7( '(', F4.1, ',', F4.1, ') ', : ) )
381 9987 FORMAT( ' AMEND DATA FILE OR INCREASE ARRAY SIZES IN PROGRAM',
382 $ /' ******* TESTS ABANDONED *******' )
383 9986 FORMAT( ' SUBPROGRAM NAME ', A6, ' NOT RECOGNIZED', /' ******* T',
384 $ 'ESTS ABANDONED *******' )
385 9985 FORMAT( ' ERROR IN ZMVCH - IN-LINE DOT PRODUCTS ARE BEING EVALU',
386 $ 'ATED WRONGLY.', /' ZMVCH WAS CALLED WITH TRANS = ', A1,
387 $ ' AND RETURNED SAME = ', L1, ' AND ERR = ', F12.3, '.', /
388 $ ' THIS MAY BE DUE TO FAULTS IN THE ARITHMETIC OR THE COMPILER.'
389 $ , /' ******* TESTS ABANDONED *******' )
390 9984 FORMAT( A6, L2 )
391 9983 FORMAT( 1X, A6, ' WAS NOT TESTED' )
392 9982 FORMAT( /' END OF TESTS' )
393 9981 FORMAT( /' ******* FATAL ERROR - TESTS ABANDONED *******' )
394 9980 FORMAT( ' ERROR-EXITS WILL NOT BE TESTED' )
395 *
396 * End of ZBLAT2.
397 *
398 END
399 SUBROUTINE ZCHK1( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
400 $ FATAL, NIDIM, IDIM, NKB, KB, NALF, ALF, NBET,
401 $ BET, NINC, INC, NMAX, INCMAX, A, AA, AS, X, XX,
402 $ XS, Y, YY, YS, YT, G )
403 *
404 * Tests ZGEMV and ZGBMV.
405 *
406 * Auxiliary routine for test program for Level 2 Blas.
407 *
408 * -- Written on 10-August-1987.
409 * Richard Hanson, Sandia National Labs.
410 * Jeremy Du Croz, NAG Central Office.
411 *
412 * .. Parameters ..
413 COMPLEX*16 ZERO, HALF
414 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
415 $ HALF = ( 0.5D0, 0.0D0 ) )
416 DOUBLE PRECISION RZERO
417 PARAMETER ( RZERO = 0.0D0 )
418 * .. Scalar Arguments ..
419 DOUBLE PRECISION EPS, THRESH
420 INTEGER INCMAX, NALF, NBET, NIDIM, NINC, NKB, NMAX,
421 $ NOUT, NTRA
422 LOGICAL FATAL, REWI, TRACE
423 CHARACTER*6 SNAME
424 * .. Array Arguments ..
425 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ), ALF( NALF ),
426 $ AS( NMAX*NMAX ), BET( NBET ), X( NMAX ),
427 $ XS( NMAX*INCMAX ), XX( NMAX*INCMAX ),
428 $ Y( NMAX ), YS( NMAX*INCMAX ), YT( NMAX ),
429 $ YY( NMAX*INCMAX )
430 DOUBLE PRECISION G( NMAX )
431 INTEGER IDIM( NIDIM ), INC( NINC ), KB( NKB )
432 * .. Local Scalars ..
433 COMPLEX*16 ALPHA, ALS, BETA, BLS, TRANSL
434 DOUBLE PRECISION ERR, ERRMAX
435 INTEGER I, IA, IB, IC, IKU, IM, IN, INCX, INCXS, INCY,
436 $ INCYS, IX, IY, KL, KLS, KU, KUS, LAA, LDA,
437 $ LDAS, LX, LY, M, ML, MS, N, NARGS, NC, ND, NK,
438 $ NL, NS
439 LOGICAL BANDED, FULL, NULL, RESET, SAME, TRAN
440 CHARACTER*1 TRANS, TRANSS
441 CHARACTER*3 ICH
442 * .. Local Arrays ..
443 LOGICAL ISAME( 13 )
444 * .. External Functions ..
445 LOGICAL LZE, LZERES
446 EXTERNAL LZE, LZERES
447 * .. External Subroutines ..
448 EXTERNAL ZGBMV, ZGEMV, ZMAKE, ZMVCH
449 * .. Intrinsic Functions ..
450 INTRINSIC ABS, MAX, MIN
451 * .. Scalars in Common ..
452 INTEGER INFOT, NOUTC
453 LOGICAL LERR, OK
454 * .. Common blocks ..
455 COMMON /INFOC/INFOT, NOUTC, OK, LERR
456 * .. Data statements ..
457 DATA ICH/'NTC'/
458 * .. Executable Statements ..
459 FULL = SNAME( 3: 3 ).EQ.'E'
460 BANDED = SNAME( 3: 3 ).EQ.'B'
461 * Define the number of arguments.
462 IF( FULL )THEN
463 NARGS = 11
464 ELSE IF( BANDED )THEN
465 NARGS = 13
466 END IF
467 *
468 NC = 0
469 RESET = .TRUE.
470 ERRMAX = RZERO
471 *
472 DO 120 IN = 1, NIDIM
473 N = IDIM( IN )
474 ND = N/2 + 1
475 *
476 DO 110 IM = 1, 2
477 IF( IM.EQ.1 )
478 $ M = MAX( N - ND, 0 )
479 IF( IM.EQ.2 )
480 $ M = MIN( N + ND, NMAX )
481 *
482 IF( BANDED )THEN
483 NK = NKB
484 ELSE
485 NK = 1
486 END IF
487 DO 100 IKU = 1, NK
488 IF( BANDED )THEN
489 KU = KB( IKU )
490 KL = MAX( KU - 1, 0 )
491 ELSE
492 KU = N - 1
493 KL = M - 1
494 END IF
495 * Set LDA to 1 more than minimum value if room.
496 IF( BANDED )THEN
497 LDA = KL + KU + 1
498 ELSE
499 LDA = M
500 END IF
501 IF( LDA.LT.NMAX )
502 $ LDA = LDA + 1
503 * Skip tests if not enough room.
504 IF( LDA.GT.NMAX )
505 $ GO TO 100
506 LAA = LDA*N
507 NULL = N.LE.0.OR.M.LE.0
508 *
509 * Generate the matrix A.
510 *
511 TRANSL = ZERO
512 CALL ZMAKE( SNAME( 2: 3 ), ' ', ' ', M, N, A, NMAX, AA,
513 $ LDA, KL, KU, RESET, TRANSL )
514 *
515 DO 90 IC = 1, 3
516 TRANS = ICH( IC: IC )
517 TRAN = TRANS.EQ.'T'.OR.TRANS.EQ.'C'
518 *
519 IF( TRAN )THEN
520 ML = N
521 NL = M
522 ELSE
523 ML = M
524 NL = N
525 END IF
526 *
527 DO 80 IX = 1, NINC
528 INCX = INC( IX )
529 LX = ABS( INCX )*NL
530 *
531 * Generate the vector X.
532 *
533 TRANSL = HALF
534 CALL ZMAKE( 'GE', ' ', ' ', 1, NL, X, 1, XX,
535 $ ABS( INCX ), 0, NL - 1, RESET, TRANSL )
536 IF( NL.GT.1 )THEN
537 X( NL/2 ) = ZERO
538 XX( 1 + ABS( INCX )*( NL/2 - 1 ) ) = ZERO
539 END IF
540 *
541 DO 70 IY = 1, NINC
542 INCY = INC( IY )
543 LY = ABS( INCY )*ML
544 *
545 DO 60 IA = 1, NALF
546 ALPHA = ALF( IA )
547 *
548 DO 50 IB = 1, NBET
549 BETA = BET( IB )
550 *
551 * Generate the vector Y.
552 *
553 TRANSL = ZERO
554 CALL ZMAKE( 'GE', ' ', ' ', 1, ML, Y, 1,
555 $ YY, ABS( INCY ), 0, ML - 1,
556 $ RESET, TRANSL )
557 *
558 NC = NC + 1
559 *
560 * Save every datum before calling the
561 * subroutine.
562 *
563 TRANSS = TRANS
564 MS = M
565 NS = N
566 KLS = KL
567 KUS = KU
568 ALS = ALPHA
569 DO 10 I = 1, LAA
570 AS( I ) = AA( I )
571 10 CONTINUE
572 LDAS = LDA
573 DO 20 I = 1, LX
574 XS( I ) = XX( I )
575 20 CONTINUE
576 INCXS = INCX
577 BLS = BETA
578 DO 30 I = 1, LY
579 YS( I ) = YY( I )
580 30 CONTINUE
581 INCYS = INCY
582 *
583 * Call the subroutine.
584 *
585 IF( FULL )THEN
586 IF( TRACE )
587 $ WRITE( NTRA, FMT = 9994 )NC, SNAME,
588 $ TRANS, M, N, ALPHA, LDA, INCX, BETA,
589 $ INCY
590 IF( REWI )
591 $ REWIND NTRA
592 CALL ZGEMV( TRANS, M, N, ALPHA, AA,
593 $ LDA, XX, INCX, BETA, YY,
594 $ INCY )
595 ELSE IF( BANDED )THEN
596 IF( TRACE )
597 $ WRITE( NTRA, FMT = 9995 )NC, SNAME,
598 $ TRANS, M, N, KL, KU, ALPHA, LDA,
599 $ INCX, BETA, INCY
600 IF( REWI )
601 $ REWIND NTRA
602 CALL ZGBMV( TRANS, M, N, KL, KU, ALPHA,
603 $ AA, LDA, XX, INCX, BETA,
604 $ YY, INCY )
605 END IF
606 *
607 * Check if error-exit was taken incorrectly.
608 *
609 IF( .NOT.OK )THEN
610 WRITE( NOUT, FMT = 9993 )
611 FATAL = .TRUE.
612 GO TO 130
613 END IF
614 *
615 * See what data changed inside subroutines.
616 *
617 ISAME( 1 ) = TRANS.EQ.TRANSS
618 ISAME( 2 ) = MS.EQ.M
619 ISAME( 3 ) = NS.EQ.N
620 IF( FULL )THEN
621 ISAME( 4 ) = ALS.EQ.ALPHA
622 ISAME( 5 ) = LZE( AS, AA, LAA )
623 ISAME( 6 ) = LDAS.EQ.LDA
624 ISAME( 7 ) = LZE( XS, XX, LX )
625 ISAME( 8 ) = INCXS.EQ.INCX
626 ISAME( 9 ) = BLS.EQ.BETA
627 IF( NULL )THEN
628 ISAME( 10 ) = LZE( YS, YY, LY )
629 ELSE
630 ISAME( 10 ) = LZERES( 'GE', ' ', 1,
631 $ ML, YS, YY,
632 $ ABS( INCY ) )
633 END IF
634 ISAME( 11 ) = INCYS.EQ.INCY
635 ELSE IF( BANDED )THEN
636 ISAME( 4 ) = KLS.EQ.KL
637 ISAME( 5 ) = KUS.EQ.KU
638 ISAME( 6 ) = ALS.EQ.ALPHA
639 ISAME( 7 ) = LZE( AS, AA, LAA )
640 ISAME( 8 ) = LDAS.EQ.LDA
641 ISAME( 9 ) = LZE( XS, XX, LX )
642 ISAME( 10 ) = INCXS.EQ.INCX
643 ISAME( 11 ) = BLS.EQ.BETA
644 IF( NULL )THEN
645 ISAME( 12 ) = LZE( YS, YY, LY )
646 ELSE
647 ISAME( 12 ) = LZERES( 'GE', ' ', 1,
648 $ ML, YS, YY,
649 $ ABS( INCY ) )
650 END IF
651 ISAME( 13 ) = INCYS.EQ.INCY
652 END IF
653 *
654 * If data was incorrectly changed, report
655 * and return.
656 *
657 SAME = .TRUE.
658 DO 40 I = 1, NARGS
659 SAME = SAME.AND.ISAME( I )
660 IF( .NOT.ISAME( I ) )
661 $ WRITE( NOUT, FMT = 9998 )I
662 40 CONTINUE
663 IF( .NOT.SAME )THEN
664 FATAL = .TRUE.
665 GO TO 130
666 END IF
667 *
668 IF( .NOT.NULL )THEN
669 *
670 * Check the result.
671 *
672 CALL ZMVCH( TRANS, M, N, ALPHA, A,
673 $ NMAX, X, INCX, BETA, Y,
674 $ INCY, YT, G, YY, EPS, ERR,
675 $ FATAL, NOUT, .TRUE. )
676 ERRMAX = MAX( ERRMAX, ERR )
677 * If got really bad answer, report and
678 * return.
679 IF( FATAL )
680 $ GO TO 130
681 ELSE
682 * Avoid repeating tests with M.le.0 or
683 * N.le.0.
684 GO TO 110
685 END IF
686 *
687 50 CONTINUE
688 *
689 60 CONTINUE
690 *
691 70 CONTINUE
692 *
693 80 CONTINUE
694 *
695 90 CONTINUE
696 *
697 100 CONTINUE
698 *
699 110 CONTINUE
700 *
701 120 CONTINUE
702 *
703 * Report result.
704 *
705 IF( ERRMAX.LT.THRESH )THEN
706 WRITE( NOUT, FMT = 9999 )SNAME, NC
707 ELSE
708 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
709 END IF
710 GO TO 140
711 *
712 130 CONTINUE
713 WRITE( NOUT, FMT = 9996 )SNAME
714 IF( FULL )THEN
715 WRITE( NOUT, FMT = 9994 )NC, SNAME, TRANS, M, N, ALPHA, LDA,
716 $ INCX, BETA, INCY
717 ELSE IF( BANDED )THEN
718 WRITE( NOUT, FMT = 9995 )NC, SNAME, TRANS, M, N, KL, KU,
719 $ ALPHA, LDA, INCX, BETA, INCY
720 END IF
721 *
722 140 CONTINUE
723 RETURN
724 *
725 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
726 $ 'S)' )
727 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
728 $ 'ANGED INCORRECTLY *******' )
729 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
730 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
731 $ ' - SUSPECT *******' )
732 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
733 9995 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 4( I3, ',' ), '(',
734 $ F4.1, ',', F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',',
735 $ F4.1, '), Y,', I2, ') .' )
736 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 2( I3, ',' ), '(',
737 $ F4.1, ',', F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',',
738 $ F4.1, '), Y,', I2, ') .' )
739 9993 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
740 $ '******' )
741 *
742 * End of ZCHK1.
743 *
744 END
745 SUBROUTINE ZCHK2( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
746 $ FATAL, NIDIM, IDIM, NKB, KB, NALF, ALF, NBET,
747 $ BET, NINC, INC, NMAX, INCMAX, A, AA, AS, X, XX,
748 $ XS, Y, YY, YS, YT, G )
749 *
750 * Tests ZHEMV, ZHBMV and ZHPMV.
751 *
752 * Auxiliary routine for test program for Level 2 Blas.
753 *
754 * -- Written on 10-August-1987.
755 * Richard Hanson, Sandia National Labs.
756 * Jeremy Du Croz, NAG Central Office.
757 *
758 * .. Parameters ..
759 COMPLEX*16 ZERO, HALF
760 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
761 $ HALF = ( 0.5D0, 0.0D0 ) )
762 DOUBLE PRECISION RZERO
763 PARAMETER ( RZERO = 0.0D0 )
764 * .. Scalar Arguments ..
765 DOUBLE PRECISION EPS, THRESH
766 INTEGER INCMAX, NALF, NBET, NIDIM, NINC, NKB, NMAX,
767 $ NOUT, NTRA
768 LOGICAL FATAL, REWI, TRACE
769 CHARACTER*6 SNAME
770 * .. Array Arguments ..
771 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ), ALF( NALF ),
772 $ AS( NMAX*NMAX ), BET( NBET ), X( NMAX ),
773 $ XS( NMAX*INCMAX ), XX( NMAX*INCMAX ),
774 $ Y( NMAX ), YS( NMAX*INCMAX ), YT( NMAX ),
775 $ YY( NMAX*INCMAX )
776 DOUBLE PRECISION G( NMAX )
777 INTEGER IDIM( NIDIM ), INC( NINC ), KB( NKB )
778 * .. Local Scalars ..
779 COMPLEX*16 ALPHA, ALS, BETA, BLS, TRANSL
780 DOUBLE PRECISION ERR, ERRMAX
781 INTEGER I, IA, IB, IC, IK, IN, INCX, INCXS, INCY,
782 $ INCYS, IX, IY, K, KS, LAA, LDA, LDAS, LX, LY,
783 $ N, NARGS, NC, NK, NS
784 LOGICAL BANDED, FULL, NULL, PACKED, RESET, SAME
785 CHARACTER*1 UPLO, UPLOS
786 CHARACTER*2 ICH
787 * .. Local Arrays ..
788 LOGICAL ISAME( 13 )
789 * .. External Functions ..
790 LOGICAL LZE, LZERES
791 EXTERNAL LZE, LZERES
792 * .. External Subroutines ..
793 EXTERNAL ZHBMV, ZHEMV, ZHPMV, ZMAKE, ZMVCH
794 * .. Intrinsic Functions ..
795 INTRINSIC ABS, MAX
796 * .. Scalars in Common ..
797 INTEGER INFOT, NOUTC
798 LOGICAL LERR, OK
799 * .. Common blocks ..
800 COMMON /INFOC/INFOT, NOUTC, OK, LERR
801 * .. Data statements ..
802 DATA ICH/'UL'/
803 * .. Executable Statements ..
804 FULL = SNAME( 3: 3 ).EQ.'E'
805 BANDED = SNAME( 3: 3 ).EQ.'B'
806 PACKED = SNAME( 3: 3 ).EQ.'P'
807 * Define the number of arguments.
808 IF( FULL )THEN
809 NARGS = 10
810 ELSE IF( BANDED )THEN
811 NARGS = 11
812 ELSE IF( PACKED )THEN
813 NARGS = 9
814 END IF
815 *
816 NC = 0
817 RESET = .TRUE.
818 ERRMAX = RZERO
819 *
820 DO 110 IN = 1, NIDIM
821 N = IDIM( IN )
822 *
823 IF( BANDED )THEN
824 NK = NKB
825 ELSE
826 NK = 1
827 END IF
828 DO 100 IK = 1, NK
829 IF( BANDED )THEN
830 K = KB( IK )
831 ELSE
832 K = N - 1
833 END IF
834 * Set LDA to 1 more than minimum value if room.
835 IF( BANDED )THEN
836 LDA = K + 1
837 ELSE
838 LDA = N
839 END IF
840 IF( LDA.LT.NMAX )
841 $ LDA = LDA + 1
842 * Skip tests if not enough room.
843 IF( LDA.GT.NMAX )
844 $ GO TO 100
845 IF( PACKED )THEN
846 LAA = ( N*( N + 1 ) )/2
847 ELSE
848 LAA = LDA*N
849 END IF
850 NULL = N.LE.0
851 *
852 DO 90 IC = 1, 2
853 UPLO = ICH( IC: IC )
854 *
855 * Generate the matrix A.
856 *
857 TRANSL = ZERO
858 CALL ZMAKE( SNAME( 2: 3 ), UPLO, ' ', N, N, A, NMAX, AA,
859 $ LDA, K, K, RESET, TRANSL )
860 *
861 DO 80 IX = 1, NINC
862 INCX = INC( IX )
863 LX = ABS( INCX )*N
864 *
865 * Generate the vector X.
866 *
867 TRANSL = HALF
868 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX,
869 $ ABS( INCX ), 0, N - 1, RESET, TRANSL )
870 IF( N.GT.1 )THEN
871 X( N/2 ) = ZERO
872 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
873 END IF
874 *
875 DO 70 IY = 1, NINC
876 INCY = INC( IY )
877 LY = ABS( INCY )*N
878 *
879 DO 60 IA = 1, NALF
880 ALPHA = ALF( IA )
881 *
882 DO 50 IB = 1, NBET
883 BETA = BET( IB )
884 *
885 * Generate the vector Y.
886 *
887 TRANSL = ZERO
888 CALL ZMAKE( 'GE', ' ', ' ', 1, N, Y, 1, YY,
889 $ ABS( INCY ), 0, N - 1, RESET,
890 $ TRANSL )
891 *
892 NC = NC + 1
893 *
894 * Save every datum before calling the
895 * subroutine.
896 *
897 UPLOS = UPLO
898 NS = N
899 KS = K
900 ALS = ALPHA
901 DO 10 I = 1, LAA
902 AS( I ) = AA( I )
903 10 CONTINUE
904 LDAS = LDA
905 DO 20 I = 1, LX
906 XS( I ) = XX( I )
907 20 CONTINUE
908 INCXS = INCX
909 BLS = BETA
910 DO 30 I = 1, LY
911 YS( I ) = YY( I )
912 30 CONTINUE
913 INCYS = INCY
914 *
915 * Call the subroutine.
916 *
917 IF( FULL )THEN
918 IF( TRACE )
919 $ WRITE( NTRA, FMT = 9993 )NC, SNAME,
920 $ UPLO, N, ALPHA, LDA, INCX, BETA, INCY
921 IF( REWI )
922 $ REWIND NTRA
923 CALL ZHEMV( UPLO, N, ALPHA, AA, LDA, XX,
924 $ INCX, BETA, YY, INCY )
925 ELSE IF( BANDED )THEN
926 IF( TRACE )
927 $ WRITE( NTRA, FMT = 9994 )NC, SNAME,
928 $ UPLO, N, K, ALPHA, LDA, INCX, BETA,
929 $ INCY
930 IF( REWI )
931 $ REWIND NTRA
932 CALL ZHBMV( UPLO, N, K, ALPHA, AA, LDA,
933 $ XX, INCX, BETA, YY, INCY )
934 ELSE IF( PACKED )THEN
935 IF( TRACE )
936 $ WRITE( NTRA, FMT = 9995 )NC, SNAME,
937 $ UPLO, N, ALPHA, INCX, BETA, INCY
938 IF( REWI )
939 $ REWIND NTRA
940 CALL ZHPMV( UPLO, N, ALPHA, AA, XX, INCX,
941 $ BETA, YY, INCY )
942 END IF
943 *
944 * Check if error-exit was taken incorrectly.
945 *
946 IF( .NOT.OK )THEN
947 WRITE( NOUT, FMT = 9992 )
948 FATAL = .TRUE.
949 GO TO 120
950 END IF
951 *
952 * See what data changed inside subroutines.
953 *
954 ISAME( 1 ) = UPLO.EQ.UPLOS
955 ISAME( 2 ) = NS.EQ.N
956 IF( FULL )THEN
957 ISAME( 3 ) = ALS.EQ.ALPHA
958 ISAME( 4 ) = LZE( AS, AA, LAA )
959 ISAME( 5 ) = LDAS.EQ.LDA
960 ISAME( 6 ) = LZE( XS, XX, LX )
961 ISAME( 7 ) = INCXS.EQ.INCX
962 ISAME( 8 ) = BLS.EQ.BETA
963 IF( NULL )THEN
964 ISAME( 9 ) = LZE( YS, YY, LY )
965 ELSE
966 ISAME( 9 ) = LZERES( 'GE', ' ', 1, N,
967 $ YS, YY, ABS( INCY ) )
968 END IF
969 ISAME( 10 ) = INCYS.EQ.INCY
970 ELSE IF( BANDED )THEN
971 ISAME( 3 ) = KS.EQ.K
972 ISAME( 4 ) = ALS.EQ.ALPHA
973 ISAME( 5 ) = LZE( AS, AA, LAA )
974 ISAME( 6 ) = LDAS.EQ.LDA
975 ISAME( 7 ) = LZE( XS, XX, LX )
976 ISAME( 8 ) = INCXS.EQ.INCX
977 ISAME( 9 ) = BLS.EQ.BETA
978 IF( NULL )THEN
979 ISAME( 10 ) = LZE( YS, YY, LY )
980 ELSE
981 ISAME( 10 ) = LZERES( 'GE', ' ', 1, N,
982 $ YS, YY, ABS( INCY ) )
983 END IF
984 ISAME( 11 ) = INCYS.EQ.INCY
985 ELSE IF( PACKED )THEN
986 ISAME( 3 ) = ALS.EQ.ALPHA
987 ISAME( 4 ) = LZE( AS, AA, LAA )
988 ISAME( 5 ) = LZE( XS, XX, LX )
989 ISAME( 6 ) = INCXS.EQ.INCX
990 ISAME( 7 ) = BLS.EQ.BETA
991 IF( NULL )THEN
992 ISAME( 8 ) = LZE( YS, YY, LY )
993 ELSE
994 ISAME( 8 ) = LZERES( 'GE', ' ', 1, N,
995 $ YS, YY, ABS( INCY ) )
996 END IF
997 ISAME( 9 ) = INCYS.EQ.INCY
998 END IF
999 *
1000 * If data was incorrectly changed, report and
1001 * return.
1002 *
1003 SAME = .TRUE.
1004 DO 40 I = 1, NARGS
1005 SAME = SAME.AND.ISAME( I )
1006 IF( .NOT.ISAME( I ) )
1007 $ WRITE( NOUT, FMT = 9998 )I
1008 40 CONTINUE
1009 IF( .NOT.SAME )THEN
1010 FATAL = .TRUE.
1011 GO TO 120
1012 END IF
1013 *
1014 IF( .NOT.NULL )THEN
1015 *
1016 * Check the result.
1017 *
1018 CALL ZMVCH( 'N', N, N, ALPHA, A, NMAX, X,
1019 $ INCX, BETA, Y, INCY, YT, G,
1020 $ YY, EPS, ERR, FATAL, NOUT,
1021 $ .TRUE. )
1022 ERRMAX = MAX( ERRMAX, ERR )
1023 * If got really bad answer, report and
1024 * return.
1025 IF( FATAL )
1026 $ GO TO 120
1027 ELSE
1028 * Avoid repeating tests with N.le.0
1029 GO TO 110
1030 END IF
1031 *
1032 50 CONTINUE
1033 *
1034 60 CONTINUE
1035 *
1036 70 CONTINUE
1037 *
1038 80 CONTINUE
1039 *
1040 90 CONTINUE
1041 *
1042 100 CONTINUE
1043 *
1044 110 CONTINUE
1045 *
1046 * Report result.
1047 *
1048 IF( ERRMAX.LT.THRESH )THEN
1049 WRITE( NOUT, FMT = 9999 )SNAME, NC
1050 ELSE
1051 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
1052 END IF
1053 GO TO 130
1054 *
1055 120 CONTINUE
1056 WRITE( NOUT, FMT = 9996 )SNAME
1057 IF( FULL )THEN
1058 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, N, ALPHA, LDA, INCX,
1059 $ BETA, INCY
1060 ELSE IF( BANDED )THEN
1061 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, N, K, ALPHA, LDA,
1062 $ INCX, BETA, INCY
1063 ELSE IF( PACKED )THEN
1064 WRITE( NOUT, FMT = 9995 )NC, SNAME, UPLO, N, ALPHA, INCX,
1065 $ BETA, INCY
1066 END IF
1067 *
1068 130 CONTINUE
1069 RETURN
1070 *
1071 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
1072 $ 'S)' )
1073 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1074 $ 'ANGED INCORRECTLY *******' )
1075 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
1076 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
1077 $ ' - SUSPECT *******' )
1078 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
1079 9995 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
1080 $ F4.1, '), AP, X,', I2, ',(', F4.1, ',', F4.1, '), Y,', I2,
1081 $ ') .' )
1082 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 2( I3, ',' ), '(',
1083 $ F4.1, ',', F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',',
1084 $ F4.1, '), Y,', I2, ') .' )
1085 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
1086 $ F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',', F4.1, '), ',
1087 $ 'Y,', I2, ') .' )
1088 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1089 $ '******' )
1090 *
1091 * End of ZCHK2.
1092 *
1093 END
1094 SUBROUTINE ZCHK3( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
1095 $ FATAL, NIDIM, IDIM, NKB, KB, NINC, INC, NMAX,
1096 $ INCMAX, A, AA, AS, X, XX, XS, XT, G, Z )
1097 *
1098 * Tests ZTRMV, ZTBMV, ZTPMV, ZTRSV, ZTBSV and ZTPSV.
1099 *
1100 * Auxiliary routine for test program for Level 2 Blas.
1101 *
1102 * -- Written on 10-August-1987.
1103 * Richard Hanson, Sandia National Labs.
1104 * Jeremy Du Croz, NAG Central Office.
1105 *
1106 * .. Parameters ..
1107 COMPLEX*16 ZERO, HALF, ONE
1108 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
1109 $ HALF = ( 0.5D0, 0.0D0 ),
1110 $ ONE = ( 1.0D0, 0.0D0 ) )
1111 DOUBLE PRECISION RZERO
1112 PARAMETER ( RZERO = 0.0D0 )
1113 * .. Scalar Arguments ..
1114 DOUBLE PRECISION EPS, THRESH
1115 INTEGER INCMAX, NIDIM, NINC, NKB, NMAX, NOUT, NTRA
1116 LOGICAL FATAL, REWI, TRACE
1117 CHARACTER*6 SNAME
1118 * .. Array Arguments ..
1119 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ),
1120 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ),
1121 $ XT( NMAX ), XX( NMAX*INCMAX ), Z( NMAX )
1122 DOUBLE PRECISION G( NMAX )
1123 INTEGER IDIM( NIDIM ), INC( NINC ), KB( NKB )
1124 * .. Local Scalars ..
1125 COMPLEX*16 TRANSL
1126 DOUBLE PRECISION ERR, ERRMAX
1127 INTEGER I, ICD, ICT, ICU, IK, IN, INCX, INCXS, IX, K,
1128 $ KS, LAA, LDA, LDAS, LX, N, NARGS, NC, NK, NS
1129 LOGICAL BANDED, FULL, NULL, PACKED, RESET, SAME
1130 CHARACTER*1 DIAG, DIAGS, TRANS, TRANSS, UPLO, UPLOS
1131 CHARACTER*2 ICHD, ICHU
1132 CHARACTER*3 ICHT
1133 * .. Local Arrays ..
1134 LOGICAL ISAME( 13 )
1135 * .. External Functions ..
1136 LOGICAL LZE, LZERES
1137 EXTERNAL LZE, LZERES
1138 * .. External Subroutines ..
1139 EXTERNAL ZMAKE, ZMVCH, ZTBMV, ZTBSV, ZTPMV, ZTPSV,
1140 $ ZTRMV, ZTRSV
1141 * .. Intrinsic Functions ..
1142 INTRINSIC ABS, MAX
1143 * .. Scalars in Common ..
1144 INTEGER INFOT, NOUTC
1145 LOGICAL LERR, OK
1146 * .. Common blocks ..
1147 COMMON /INFOC/INFOT, NOUTC, OK, LERR
1148 * .. Data statements ..
1149 DATA ICHU/'UL'/, ICHT/'NTC'/, ICHD/'UN'/
1150 * .. Executable Statements ..
1151 FULL = SNAME( 3: 3 ).EQ.'R'
1152 BANDED = SNAME( 3: 3 ).EQ.'B'
1153 PACKED = SNAME( 3: 3 ).EQ.'P'
1154 * Define the number of arguments.
1155 IF( FULL )THEN
1156 NARGS = 8
1157 ELSE IF( BANDED )THEN
1158 NARGS = 9
1159 ELSE IF( PACKED )THEN
1160 NARGS = 7
1161 END IF
1162 *
1163 NC = 0
1164 RESET = .TRUE.
1165 ERRMAX = RZERO
1166 * Set up zero vector for ZMVCH.
1167 DO 10 I = 1, NMAX
1168 Z( I ) = ZERO
1169 10 CONTINUE
1170 *
1171 DO 110 IN = 1, NIDIM
1172 N = IDIM( IN )
1173 *
1174 IF( BANDED )THEN
1175 NK = NKB
1176 ELSE
1177 NK = 1
1178 END IF
1179 DO 100 IK = 1, NK
1180 IF( BANDED )THEN
1181 K = KB( IK )
1182 ELSE
1183 K = N - 1
1184 END IF
1185 * Set LDA to 1 more than minimum value if room.
1186 IF( BANDED )THEN
1187 LDA = K + 1
1188 ELSE
1189 LDA = N
1190 END IF
1191 IF( LDA.LT.NMAX )
1192 $ LDA = LDA + 1
1193 * Skip tests if not enough room.
1194 IF( LDA.GT.NMAX )
1195 $ GO TO 100
1196 IF( PACKED )THEN
1197 LAA = ( N*( N + 1 ) )/2
1198 ELSE
1199 LAA = LDA*N
1200 END IF
1201 NULL = N.LE.0
1202 *
1203 DO 90 ICU = 1, 2
1204 UPLO = ICHU( ICU: ICU )
1205 *
1206 DO 80 ICT = 1, 3
1207 TRANS = ICHT( ICT: ICT )
1208 *
1209 DO 70 ICD = 1, 2
1210 DIAG = ICHD( ICD: ICD )
1211 *
1212 * Generate the matrix A.
1213 *
1214 TRANSL = ZERO
1215 CALL ZMAKE( SNAME( 2: 3 ), UPLO, DIAG, N, N, A,
1216 $ NMAX, AA, LDA, K, K, RESET, TRANSL )
1217 *
1218 DO 60 IX = 1, NINC
1219 INCX = INC( IX )
1220 LX = ABS( INCX )*N
1221 *
1222 * Generate the vector X.
1223 *
1224 TRANSL = HALF
1225 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX,
1226 $ ABS( INCX ), 0, N - 1, RESET,
1227 $ TRANSL )
1228 IF( N.GT.1 )THEN
1229 X( N/2 ) = ZERO
1230 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
1231 END IF
1232 *
1233 NC = NC + 1
1234 *
1235 * Save every datum before calling the subroutine.
1236 *
1237 UPLOS = UPLO
1238 TRANSS = TRANS
1239 DIAGS = DIAG
1240 NS = N
1241 KS = K
1242 DO 20 I = 1, LAA
1243 AS( I ) = AA( I )
1244 20 CONTINUE
1245 LDAS = LDA
1246 DO 30 I = 1, LX
1247 XS( I ) = XX( I )
1248 30 CONTINUE
1249 INCXS = INCX
1250 *
1251 * Call the subroutine.
1252 *
1253 IF( SNAME( 4: 5 ).EQ.'MV' )THEN
1254 IF( FULL )THEN
1255 IF( TRACE )
1256 $ WRITE( NTRA, FMT = 9993 )NC, SNAME,
1257 $ UPLO, TRANS, DIAG, N, LDA, INCX
1258 IF( REWI )
1259 $ REWIND NTRA
1260 CALL ZTRMV( UPLO, TRANS, DIAG, N, AA, LDA,
1261 $ XX, INCX )
1262 ELSE IF( BANDED )THEN
1263 IF( TRACE )
1264 $ WRITE( NTRA, FMT = 9994 )NC, SNAME,
1265 $ UPLO, TRANS, DIAG, N, K, LDA, INCX
1266 IF( REWI )
1267 $ REWIND NTRA
1268 CALL ZTBMV( UPLO, TRANS, DIAG, N, K, AA,
1269 $ LDA, XX, INCX )
1270 ELSE IF( PACKED )THEN
1271 IF( TRACE )
1272 $ WRITE( NTRA, FMT = 9995 )NC, SNAME,
1273 $ UPLO, TRANS, DIAG, N, INCX
1274 IF( REWI )
1275 $ REWIND NTRA
1276 CALL ZTPMV( UPLO, TRANS, DIAG, N, AA, XX,
1277 $ INCX )
1278 END IF
1279 ELSE IF( SNAME( 4: 5 ).EQ.'SV' )THEN
1280 IF( FULL )THEN
1281 IF( TRACE )
1282 $ WRITE( NTRA, FMT = 9993 )NC, SNAME,
1283 $ UPLO, TRANS, DIAG, N, LDA, INCX
1284 IF( REWI )
1285 $ REWIND NTRA
1286 CALL ZTRSV( UPLO, TRANS, DIAG, N, AA, LDA,
1287 $ XX, INCX )
1288 ELSE IF( BANDED )THEN
1289 IF( TRACE )
1290 $ WRITE( NTRA, FMT = 9994 )NC, SNAME,
1291 $ UPLO, TRANS, DIAG, N, K, LDA, INCX
1292 IF( REWI )
1293 $ REWIND NTRA
1294 CALL ZTBSV( UPLO, TRANS, DIAG, N, K, AA,
1295 $ LDA, XX, INCX )
1296 ELSE IF( PACKED )THEN
1297 IF( TRACE )
1298 $ WRITE( NTRA, FMT = 9995 )NC, SNAME,
1299 $ UPLO, TRANS, DIAG, N, INCX
1300 IF( REWI )
1301 $ REWIND NTRA
1302 CALL ZTPSV( UPLO, TRANS, DIAG, N, AA, XX,
1303 $ INCX )
1304 END IF
1305 END IF
1306 *
1307 * Check if error-exit was taken incorrectly.
1308 *
1309 IF( .NOT.OK )THEN
1310 WRITE( NOUT, FMT = 9992 )
1311 FATAL = .TRUE.
1312 GO TO 120
1313 END IF
1314 *
1315 * See what data changed inside subroutines.
1316 *
1317 ISAME( 1 ) = UPLO.EQ.UPLOS
1318 ISAME( 2 ) = TRANS.EQ.TRANSS
1319 ISAME( 3 ) = DIAG.EQ.DIAGS
1320 ISAME( 4 ) = NS.EQ.N
1321 IF( FULL )THEN
1322 ISAME( 5 ) = LZE( AS, AA, LAA )
1323 ISAME( 6 ) = LDAS.EQ.LDA
1324 IF( NULL )THEN
1325 ISAME( 7 ) = LZE( XS, XX, LX )
1326 ELSE
1327 ISAME( 7 ) = LZERES( 'GE', ' ', 1, N, XS,
1328 $ XX, ABS( INCX ) )
1329 END IF
1330 ISAME( 8 ) = INCXS.EQ.INCX
1331 ELSE IF( BANDED )THEN
1332 ISAME( 5 ) = KS.EQ.K
1333 ISAME( 6 ) = LZE( AS, AA, LAA )
1334 ISAME( 7 ) = LDAS.EQ.LDA
1335 IF( NULL )THEN
1336 ISAME( 8 ) = LZE( XS, XX, LX )
1337 ELSE
1338 ISAME( 8 ) = LZERES( 'GE', ' ', 1, N, XS,
1339 $ XX, ABS( INCX ) )
1340 END IF
1341 ISAME( 9 ) = INCXS.EQ.INCX
1342 ELSE IF( PACKED )THEN
1343 ISAME( 5 ) = LZE( AS, AA, LAA )
1344 IF( NULL )THEN
1345 ISAME( 6 ) = LZE( XS, XX, LX )
1346 ELSE
1347 ISAME( 6 ) = LZERES( 'GE', ' ', 1, N, XS,
1348 $ XX, ABS( INCX ) )
1349 END IF
1350 ISAME( 7 ) = INCXS.EQ.INCX
1351 END IF
1352 *
1353 * If data was incorrectly changed, report and
1354 * return.
1355 *
1356 SAME = .TRUE.
1357 DO 40 I = 1, NARGS
1358 SAME = SAME.AND.ISAME( I )
1359 IF( .NOT.ISAME( I ) )
1360 $ WRITE( NOUT, FMT = 9998 )I
1361 40 CONTINUE
1362 IF( .NOT.SAME )THEN
1363 FATAL = .TRUE.
1364 GO TO 120
1365 END IF
1366 *
1367 IF( .NOT.NULL )THEN
1368 IF( SNAME( 4: 5 ).EQ.'MV' )THEN
1369 *
1370 * Check the result.
1371 *
1372 CALL ZMVCH( TRANS, N, N, ONE, A, NMAX, X,
1373 $ INCX, ZERO, Z, INCX, XT, G,
1374 $ XX, EPS, ERR, FATAL, NOUT,
1375 $ .TRUE. )
1376 ELSE IF( SNAME( 4: 5 ).EQ.'SV' )THEN
1377 *
1378 * Compute approximation to original vector.
1379 *
1380 DO 50 I = 1, N
1381 Z( I ) = XX( 1 + ( I - 1 )*
1382 $ ABS( INCX ) )
1383 XX( 1 + ( I - 1 )*ABS( INCX ) )
1384 $ = X( I )
1385 50 CONTINUE
1386 CALL ZMVCH( TRANS, N, N, ONE, A, NMAX, Z,
1387 $ INCX, ZERO, X, INCX, XT, G,
1388 $ XX, EPS, ERR, FATAL, NOUT,
1389 $ .FALSE. )
1390 END IF
1391 ERRMAX = MAX( ERRMAX, ERR )
1392 * If got really bad answer, report and return.
1393 IF( FATAL )
1394 $ GO TO 120
1395 ELSE
1396 * Avoid repeating tests with N.le.0.
1397 GO TO 110
1398 END IF
1399 *
1400 60 CONTINUE
1401 *
1402 70 CONTINUE
1403 *
1404 80 CONTINUE
1405 *
1406 90 CONTINUE
1407 *
1408 100 CONTINUE
1409 *
1410 110 CONTINUE
1411 *
1412 * Report result.
1413 *
1414 IF( ERRMAX.LT.THRESH )THEN
1415 WRITE( NOUT, FMT = 9999 )SNAME, NC
1416 ELSE
1417 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
1418 END IF
1419 GO TO 130
1420 *
1421 120 CONTINUE
1422 WRITE( NOUT, FMT = 9996 )SNAME
1423 IF( FULL )THEN
1424 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, TRANS, DIAG, N, LDA,
1425 $ INCX
1426 ELSE IF( BANDED )THEN
1427 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, TRANS, DIAG, N, K,
1428 $ LDA, INCX
1429 ELSE IF( PACKED )THEN
1430 WRITE( NOUT, FMT = 9995 )NC, SNAME, UPLO, TRANS, DIAG, N, INCX
1431 END IF
1432 *
1433 130 CONTINUE
1434 RETURN
1435 *
1436 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
1437 $ 'S)' )
1438 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1439 $ 'ANGED INCORRECTLY *******' )
1440 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
1441 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
1442 $ ' - SUSPECT *******' )
1443 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
1444 9995 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), I3, ', AP, ',
1445 $ 'X,', I2, ') .' )
1446 9994 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), 2( I3, ',' ),
1447 $ ' A,', I3, ', X,', I2, ') .' )
1448 9993 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), I3, ', A,',
1449 $ I3, ', X,', I2, ') .' )
1450 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1451 $ '******' )
1452 *
1453 * End of ZCHK3.
1454 *
1455 END
1456 SUBROUTINE ZCHK4( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
1457 $ FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC, NMAX,
1458 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G,
1459 $ Z )
1460 *
1461 * Tests ZGERC and ZGERU.
1462 *
1463 * Auxiliary routine for test program for Level 2 Blas.
1464 *
1465 * -- Written on 10-August-1987.
1466 * Richard Hanson, Sandia National Labs.
1467 * Jeremy Du Croz, NAG Central Office.
1468 *
1469 * .. Parameters ..
1470 COMPLEX*16 ZERO, HALF, ONE
1471 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
1472 $ HALF = ( 0.5D0, 0.0D0 ),
1473 $ ONE = ( 1.0D0, 0.0D0 ) )
1474 DOUBLE PRECISION RZERO
1475 PARAMETER ( RZERO = 0.0D0 )
1476 * .. Scalar Arguments ..
1477 DOUBLE PRECISION EPS, THRESH
1478 INTEGER INCMAX, NALF, NIDIM, NINC, NMAX, NOUT, NTRA
1479 LOGICAL FATAL, REWI, TRACE
1480 CHARACTER*6 SNAME
1481 * .. Array Arguments ..
1482 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ), ALF( NALF ),
1483 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ),
1484 $ XX( NMAX*INCMAX ), Y( NMAX ),
1485 $ YS( NMAX*INCMAX ), YT( NMAX ),
1486 $ YY( NMAX*INCMAX ), Z( NMAX )
1487 DOUBLE PRECISION G( NMAX )
1488 INTEGER IDIM( NIDIM ), INC( NINC )
1489 * .. Local Scalars ..
1490 COMPLEX*16 ALPHA, ALS, TRANSL
1491 DOUBLE PRECISION ERR, ERRMAX
1492 INTEGER I, IA, IM, IN, INCX, INCXS, INCY, INCYS, IX,
1493 $ IY, J, LAA, LDA, LDAS, LX, LY, M, MS, N, NARGS,
1494 $ NC, ND, NS
1495 LOGICAL CONJ, NULL, RESET, SAME
1496 * .. Local Arrays ..
1497 COMPLEX*16 W( 1 )
1498 LOGICAL ISAME( 13 )
1499 * .. External Functions ..
1500 LOGICAL LZE, LZERES
1501 EXTERNAL LZE, LZERES
1502 * .. External Subroutines ..
1503 EXTERNAL ZGERC, ZGERU, ZMAKE, ZMVCH
1504 * .. Intrinsic Functions ..
1505 INTRINSIC ABS, DCONJG, MAX, MIN
1506 * .. Scalars in Common ..
1507 INTEGER INFOT, NOUTC
1508 LOGICAL LERR, OK
1509 * .. Common blocks ..
1510 COMMON /INFOC/INFOT, NOUTC, OK, LERR
1511 * .. Executable Statements ..
1512 CONJ = SNAME( 5: 5 ).EQ.'C'
1513 * Define the number of arguments.
1514 NARGS = 9
1515 *
1516 NC = 0
1517 RESET = .TRUE.
1518 ERRMAX = RZERO
1519 *
1520 DO 120 IN = 1, NIDIM
1521 N = IDIM( IN )
1522 ND = N/2 + 1
1523 *
1524 DO 110 IM = 1, 2
1525 IF( IM.EQ.1 )
1526 $ M = MAX( N - ND, 0 )
1527 IF( IM.EQ.2 )
1528 $ M = MIN( N + ND, NMAX )
1529 *
1530 * Set LDA to 1 more than minimum value if room.
1531 LDA = M
1532 IF( LDA.LT.NMAX )
1533 $ LDA = LDA + 1
1534 * Skip tests if not enough room.
1535 IF( LDA.GT.NMAX )
1536 $ GO TO 110
1537 LAA = LDA*N
1538 NULL = N.LE.0.OR.M.LE.0
1539 *
1540 DO 100 IX = 1, NINC
1541 INCX = INC( IX )
1542 LX = ABS( INCX )*M
1543 *
1544 * Generate the vector X.
1545 *
1546 TRANSL = HALF
1547 CALL ZMAKE( 'GE', ' ', ' ', 1, M, X, 1, XX, ABS( INCX ),
1548 $ 0, M - 1, RESET, TRANSL )
1549 IF( M.GT.1 )THEN
1550 X( M/2 ) = ZERO
1551 XX( 1 + ABS( INCX )*( M/2 - 1 ) ) = ZERO
1552 END IF
1553 *
1554 DO 90 IY = 1, NINC
1555 INCY = INC( IY )
1556 LY = ABS( INCY )*N
1557 *
1558 * Generate the vector Y.
1559 *
1560 TRANSL = ZERO
1561 CALL ZMAKE( 'GE', ' ', ' ', 1, N, Y, 1, YY,
1562 $ ABS( INCY ), 0, N - 1, RESET, TRANSL )
1563 IF( N.GT.1 )THEN
1564 Y( N/2 ) = ZERO
1565 YY( 1 + ABS( INCY )*( N/2 - 1 ) ) = ZERO
1566 END IF
1567 *
1568 DO 80 IA = 1, NALF
1569 ALPHA = ALF( IA )
1570 *
1571 * Generate the matrix A.
1572 *
1573 TRANSL = ZERO
1574 CALL ZMAKE( SNAME( 2: 3 ), ' ', ' ', M, N, A, NMAX,
1575 $ AA, LDA, M - 1, N - 1, RESET, TRANSL )
1576 *
1577 NC = NC + 1
1578 *
1579 * Save every datum before calling the subroutine.
1580 *
1581 MS = M
1582 NS = N
1583 ALS = ALPHA
1584 DO 10 I = 1, LAA
1585 AS( I ) = AA( I )
1586 10 CONTINUE
1587 LDAS = LDA
1588 DO 20 I = 1, LX
1589 XS( I ) = XX( I )
1590 20 CONTINUE
1591 INCXS = INCX
1592 DO 30 I = 1, LY
1593 YS( I ) = YY( I )
1594 30 CONTINUE
1595 INCYS = INCY
1596 *
1597 * Call the subroutine.
1598 *
1599 IF( TRACE )
1600 $ WRITE( NTRA, FMT = 9994 )NC, SNAME, M, N,
1601 $ ALPHA, INCX, INCY, LDA
1602 IF( CONJ )THEN
1603 IF( REWI )
1604 $ REWIND NTRA
1605 CALL ZGERC( M, N, ALPHA, XX, INCX, YY, INCY, AA,
1606 $ LDA )
1607 ELSE
1608 IF( REWI )
1609 $ REWIND NTRA
1610 CALL ZGERU( M, N, ALPHA, XX, INCX, YY, INCY, AA,
1611 $ LDA )
1612 END IF
1613 *
1614 * Check if error-exit was taken incorrectly.
1615 *
1616 IF( .NOT.OK )THEN
1617 WRITE( NOUT, FMT = 9993 )
1618 FATAL = .TRUE.
1619 GO TO 140
1620 END IF
1621 *
1622 * See what data changed inside subroutine.
1623 *
1624 ISAME( 1 ) = MS.EQ.M
1625 ISAME( 2 ) = NS.EQ.N
1626 ISAME( 3 ) = ALS.EQ.ALPHA
1627 ISAME( 4 ) = LZE( XS, XX, LX )
1628 ISAME( 5 ) = INCXS.EQ.INCX
1629 ISAME( 6 ) = LZE( YS, YY, LY )
1630 ISAME( 7 ) = INCYS.EQ.INCY
1631 IF( NULL )THEN
1632 ISAME( 8 ) = LZE( AS, AA, LAA )
1633 ELSE
1634 ISAME( 8 ) = LZERES( 'GE', ' ', M, N, AS, AA,
1635 $ LDA )
1636 END IF
1637 ISAME( 9 ) = LDAS.EQ.LDA
1638 *
1639 * If data was incorrectly changed, report and return.
1640 *
1641 SAME = .TRUE.
1642 DO 40 I = 1, NARGS
1643 SAME = SAME.AND.ISAME( I )
1644 IF( .NOT.ISAME( I ) )
1645 $ WRITE( NOUT, FMT = 9998 )I
1646 40 CONTINUE
1647 IF( .NOT.SAME )THEN
1648 FATAL = .TRUE.
1649 GO TO 140
1650 END IF
1651 *
1652 IF( .NOT.NULL )THEN
1653 *
1654 * Check the result column by column.
1655 *
1656 IF( INCX.GT.0 )THEN
1657 DO 50 I = 1, M
1658 Z( I ) = X( I )
1659 50 CONTINUE
1660 ELSE
1661 DO 60 I = 1, M
1662 Z( I ) = X( M - I + 1 )
1663 60 CONTINUE
1664 END IF
1665 DO 70 J = 1, N
1666 IF( INCY.GT.0 )THEN
1667 W( 1 ) = Y( J )
1668 ELSE
1669 W( 1 ) = Y( N - J + 1 )
1670 END IF
1671 IF( CONJ )
1672 $ W( 1 ) = DCONJG( W( 1 ) )
1673 CALL ZMVCH( 'N', M, 1, ALPHA, Z, NMAX, W, 1,
1674 $ ONE, A( 1, J ), 1, YT, G,
1675 $ AA( 1 + ( J - 1 )*LDA ), EPS,
1676 $ ERR, FATAL, NOUT, .TRUE. )
1677 ERRMAX = MAX( ERRMAX, ERR )
1678 * If got really bad answer, report and return.
1679 IF( FATAL )
1680 $ GO TO 130
1681 70 CONTINUE
1682 ELSE
1683 * Avoid repeating tests with M.le.0 or N.le.0.
1684 GO TO 110
1685 END IF
1686 *
1687 80 CONTINUE
1688 *
1689 90 CONTINUE
1690 *
1691 100 CONTINUE
1692 *
1693 110 CONTINUE
1694 *
1695 120 CONTINUE
1696 *
1697 * Report result.
1698 *
1699 IF( ERRMAX.LT.THRESH )THEN
1700 WRITE( NOUT, FMT = 9999 )SNAME, NC
1701 ELSE
1702 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
1703 END IF
1704 GO TO 150
1705 *
1706 130 CONTINUE
1707 WRITE( NOUT, FMT = 9995 )J
1708 *
1709 140 CONTINUE
1710 WRITE( NOUT, FMT = 9996 )SNAME
1711 WRITE( NOUT, FMT = 9994 )NC, SNAME, M, N, ALPHA, INCX, INCY, LDA
1712 *
1713 150 CONTINUE
1714 RETURN
1715 *
1716 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
1717 $ 'S)' )
1718 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1719 $ 'ANGED INCORRECTLY *******' )
1720 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
1721 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
1722 $ ' - SUSPECT *******' )
1723 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
1724 9995 FORMAT( ' THESE ARE THE RESULTS FOR COLUMN ', I3 )
1725 9994 FORMAT( 1X, I6, ': ', A6, '(', 2( I3, ',' ), '(', F4.1, ',', F4.1,
1726 $ '), X,', I2, ', Y,', I2, ', A,', I3, ') ',
1727 $ ' .' )
1728 9993 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1729 $ '******' )
1730 *
1731 * End of ZCHK4.
1732 *
1733 END
1734 SUBROUTINE ZCHK5( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
1735 $ FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC, NMAX,
1736 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G,
1737 $ Z )
1738 *
1739 * Tests ZHER and ZHPR.
1740 *
1741 * Auxiliary routine for test program for Level 2 Blas.
1742 *
1743 * -- Written on 10-August-1987.
1744 * Richard Hanson, Sandia National Labs.
1745 * Jeremy Du Croz, NAG Central Office.
1746 *
1747 * .. Parameters ..
1748 COMPLEX*16 ZERO, HALF, ONE
1749 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
1750 $ HALF = ( 0.5D0, 0.0D0 ),
1751 $ ONE = ( 1.0D0, 0.0D0 ) )
1752 DOUBLE PRECISION RZERO
1753 PARAMETER ( RZERO = 0.0D0 )
1754 * .. Scalar Arguments ..
1755 DOUBLE PRECISION EPS, THRESH
1756 INTEGER INCMAX, NALF, NIDIM, NINC, NMAX, NOUT, NTRA
1757 LOGICAL FATAL, REWI, TRACE
1758 CHARACTER*6 SNAME
1759 * .. Array Arguments ..
1760 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ), ALF( NALF ),
1761 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ),
1762 $ XX( NMAX*INCMAX ), Y( NMAX ),
1763 $ YS( NMAX*INCMAX ), YT( NMAX ),
1764 $ YY( NMAX*INCMAX ), Z( NMAX )
1765 DOUBLE PRECISION G( NMAX )
1766 INTEGER IDIM( NIDIM ), INC( NINC )
1767 * .. Local Scalars ..
1768 COMPLEX*16 ALPHA, TRANSL
1769 DOUBLE PRECISION ERR, ERRMAX, RALPHA, RALS
1770 INTEGER I, IA, IC, IN, INCX, INCXS, IX, J, JA, JJ, LAA,
1771 $ LDA, LDAS, LJ, LX, N, NARGS, NC, NS
1772 LOGICAL FULL, NULL, PACKED, RESET, SAME, UPPER
1773 CHARACTER*1 UPLO, UPLOS
1774 CHARACTER*2 ICH
1775 * .. Local Arrays ..
1776 COMPLEX*16 W( 1 )
1777 LOGICAL ISAME( 13 )
1778 * .. External Functions ..
1779 LOGICAL LZE, LZERES
1780 EXTERNAL LZE, LZERES
1781 * .. External Subroutines ..
1782 EXTERNAL ZHER, ZHPR, ZMAKE, ZMVCH
1783 * .. Intrinsic Functions ..
1784 INTRINSIC ABS, DBLE, DCMPLX, DCONJG, MAX
1785 * .. Scalars in Common ..
1786 INTEGER INFOT, NOUTC
1787 LOGICAL LERR, OK
1788 * .. Common blocks ..
1789 COMMON /INFOC/INFOT, NOUTC, OK, LERR
1790 * .. Data statements ..
1791 DATA ICH/'UL'/
1792 * .. Executable Statements ..
1793 FULL = SNAME( 3: 3 ).EQ.'E'
1794 PACKED = SNAME( 3: 3 ).EQ.'P'
1795 * Define the number of arguments.
1796 IF( FULL )THEN
1797 NARGS = 7
1798 ELSE IF( PACKED )THEN
1799 NARGS = 6
1800 END IF
1801 *
1802 NC = 0
1803 RESET = .TRUE.
1804 ERRMAX = RZERO
1805 *
1806 DO 100 IN = 1, NIDIM
1807 N = IDIM( IN )
1808 * Set LDA to 1 more than minimum value if room.
1809 LDA = N
1810 IF( LDA.LT.NMAX )
1811 $ LDA = LDA + 1
1812 * Skip tests if not enough room.
1813 IF( LDA.GT.NMAX )
1814 $ GO TO 100
1815 IF( PACKED )THEN
1816 LAA = ( N*( N + 1 ) )/2
1817 ELSE
1818 LAA = LDA*N
1819 END IF
1820 *
1821 DO 90 IC = 1, 2
1822 UPLO = ICH( IC: IC )
1823 UPPER = UPLO.EQ.'U'
1824 *
1825 DO 80 IX = 1, NINC
1826 INCX = INC( IX )
1827 LX = ABS( INCX )*N
1828 *
1829 * Generate the vector X.
1830 *
1831 TRANSL = HALF
1832 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX, ABS( INCX ),
1833 $ 0, N - 1, RESET, TRANSL )
1834 IF( N.GT.1 )THEN
1835 X( N/2 ) = ZERO
1836 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
1837 END IF
1838 *
1839 DO 70 IA = 1, NALF
1840 RALPHA = DBLE( ALF( IA ) )
1841 ALPHA = DCMPLX( RALPHA, RZERO )
1842 NULL = N.LE.0.OR.RALPHA.EQ.RZERO
1843 *
1844 * Generate the matrix A.
1845 *
1846 TRANSL = ZERO
1847 CALL ZMAKE( SNAME( 2: 3 ), UPLO, ' ', N, N, A, NMAX,
1848 $ AA, LDA, N - 1, N - 1, RESET, TRANSL )
1849 *
1850 NC = NC + 1
1851 *
1852 * Save every datum before calling the subroutine.
1853 *
1854 UPLOS = UPLO
1855 NS = N
1856 RALS = RALPHA
1857 DO 10 I = 1, LAA
1858 AS( I ) = AA( I )
1859 10 CONTINUE
1860 LDAS = LDA
1861 DO 20 I = 1, LX
1862 XS( I ) = XX( I )
1863 20 CONTINUE
1864 INCXS = INCX
1865 *
1866 * Call the subroutine.
1867 *
1868 IF( FULL )THEN
1869 IF( TRACE )
1870 $ WRITE( NTRA, FMT = 9993 )NC, SNAME, UPLO, N,
1871 $ RALPHA, INCX, LDA
1872 IF( REWI )
1873 $ REWIND NTRA
1874 CALL ZHER( UPLO, N, RALPHA, XX, INCX, AA, LDA )
1875 ELSE IF( PACKED )THEN
1876 IF( TRACE )
1877 $ WRITE( NTRA, FMT = 9994 )NC, SNAME, UPLO, N,
1878 $ RALPHA, INCX
1879 IF( REWI )
1880 $ REWIND NTRA
1881 CALL ZHPR( UPLO, N, RALPHA, XX, INCX, AA )
1882 END IF
1883 *
1884 * Check if error-exit was taken incorrectly.
1885 *
1886 IF( .NOT.OK )THEN
1887 WRITE( NOUT, FMT = 9992 )
1888 FATAL = .TRUE.
1889 GO TO 120
1890 END IF
1891 *
1892 * See what data changed inside subroutines.
1893 *
1894 ISAME( 1 ) = UPLO.EQ.UPLOS
1895 ISAME( 2 ) = NS.EQ.N
1896 ISAME( 3 ) = RALS.EQ.RALPHA
1897 ISAME( 4 ) = LZE( XS, XX, LX )
1898 ISAME( 5 ) = INCXS.EQ.INCX
1899 IF( NULL )THEN
1900 ISAME( 6 ) = LZE( AS, AA, LAA )
1901 ELSE
1902 ISAME( 6 ) = LZERES( SNAME( 2: 3 ), UPLO, N, N, AS,
1903 $ AA, LDA )
1904 END IF
1905 IF( .NOT.PACKED )THEN
1906 ISAME( 7 ) = LDAS.EQ.LDA
1907 END IF
1908 *
1909 * If data was incorrectly changed, report and return.
1910 *
1911 SAME = .TRUE.
1912 DO 30 I = 1, NARGS
1913 SAME = SAME.AND.ISAME( I )
1914 IF( .NOT.ISAME( I ) )
1915 $ WRITE( NOUT, FMT = 9998 )I
1916 30 CONTINUE
1917 IF( .NOT.SAME )THEN
1918 FATAL = .TRUE.
1919 GO TO 120
1920 END IF
1921 *
1922 IF( .NOT.NULL )THEN
1923 *
1924 * Check the result column by column.
1925 *
1926 IF( INCX.GT.0 )THEN
1927 DO 40 I = 1, N
1928 Z( I ) = X( I )
1929 40 CONTINUE
1930 ELSE
1931 DO 50 I = 1, N
1932 Z( I ) = X( N - I + 1 )
1933 50 CONTINUE
1934 END IF
1935 JA = 1
1936 DO 60 J = 1, N
1937 W( 1 ) = DCONJG( Z( J ) )
1938 IF( UPPER )THEN
1939 JJ = 1
1940 LJ = J
1941 ELSE
1942 JJ = J
1943 LJ = N - J + 1
1944 END IF
1945 CALL ZMVCH( 'N', LJ, 1, ALPHA, Z( JJ ), LJ, W,
1946 $ 1, ONE, A( JJ, J ), 1, YT, G,
1947 $ AA( JA ), EPS, ERR, FATAL, NOUT,
1948 $ .TRUE. )
1949 IF( FULL )THEN
1950 IF( UPPER )THEN
1951 JA = JA + LDA
1952 ELSE
1953 JA = JA + LDA + 1
1954 END IF
1955 ELSE
1956 JA = JA + LJ
1957 END IF
1958 ERRMAX = MAX( ERRMAX, ERR )
1959 * If got really bad answer, report and return.
1960 IF( FATAL )
1961 $ GO TO 110
1962 60 CONTINUE
1963 ELSE
1964 * Avoid repeating tests if N.le.0.
1965 IF( N.LE.0 )
1966 $ GO TO 100
1967 END IF
1968 *
1969 70 CONTINUE
1970 *
1971 80 CONTINUE
1972 *
1973 90 CONTINUE
1974 *
1975 100 CONTINUE
1976 *
1977 * Report result.
1978 *
1979 IF( ERRMAX.LT.THRESH )THEN
1980 WRITE( NOUT, FMT = 9999 )SNAME, NC
1981 ELSE
1982 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
1983 END IF
1984 GO TO 130
1985 *
1986 110 CONTINUE
1987 WRITE( NOUT, FMT = 9995 )J
1988 *
1989 120 CONTINUE
1990 WRITE( NOUT, FMT = 9996 )SNAME
1991 IF( FULL )THEN
1992 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, N, RALPHA, INCX, LDA
1993 ELSE IF( PACKED )THEN
1994 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, N, RALPHA, INCX
1995 END IF
1996 *
1997 130 CONTINUE
1998 RETURN
1999 *
2000 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
2001 $ 'S)' )
2002 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
2003 $ 'ANGED INCORRECTLY *******' )
2004 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
2005 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
2006 $ ' - SUSPECT *******' )
2007 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
2008 9995 FORMAT( ' THESE ARE THE RESULTS FOR COLUMN ', I3 )
2009 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', X,',
2010 $ I2, ', AP) .' )
2011 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', X,',
2012 $ I2, ', A,', I3, ') .' )
2013 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
2014 $ '******' )
2015 *
2016 * End of ZCHK5.
2017 *
2018 END
2019 SUBROUTINE ZCHK6( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
2020 $ FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC, NMAX,
2021 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G,
2022 $ Z )
2023 *
2024 * Tests ZHER2 and ZHPR2.
2025 *
2026 * Auxiliary routine for test program for Level 2 Blas.
2027 *
2028 * -- Written on 10-August-1987.
2029 * Richard Hanson, Sandia National Labs.
2030 * Jeremy Du Croz, NAG Central Office.
2031 *
2032 * .. Parameters ..
2033 COMPLEX*16 ZERO, HALF, ONE
2034 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
2035 $ HALF = ( 0.5D0, 0.0D0 ),
2036 $ ONE = ( 1.0D0, 0.0D0 ) )
2037 DOUBLE PRECISION RZERO
2038 PARAMETER ( RZERO = 0.0D0 )
2039 * .. Scalar Arguments ..
2040 DOUBLE PRECISION EPS, THRESH
2041 INTEGER INCMAX, NALF, NIDIM, NINC, NMAX, NOUT, NTRA
2042 LOGICAL FATAL, REWI, TRACE
2043 CHARACTER*6 SNAME
2044 * .. Array Arguments ..
2045 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ), ALF( NALF ),
2046 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ),
2047 $ XX( NMAX*INCMAX ), Y( NMAX ),
2048 $ YS( NMAX*INCMAX ), YT( NMAX ),
2049 $ YY( NMAX*INCMAX ), Z( NMAX, 2 )
2050 DOUBLE PRECISION G( NMAX )
2051 INTEGER IDIM( NIDIM ), INC( NINC )
2052 * .. Local Scalars ..
2053 COMPLEX*16 ALPHA, ALS, TRANSL
2054 DOUBLE PRECISION ERR, ERRMAX
2055 INTEGER I, IA, IC, IN, INCX, INCXS, INCY, INCYS, IX,
2056 $ IY, J, JA, JJ, LAA, LDA, LDAS, LJ, LX, LY, N,
2057 $ NARGS, NC, NS
2058 LOGICAL FULL, NULL, PACKED, RESET, SAME, UPPER
2059 CHARACTER*1 UPLO, UPLOS
2060 CHARACTER*2 ICH
2061 * .. Local Arrays ..
2062 COMPLEX*16 W( 2 )
2063 LOGICAL ISAME( 13 )
2064 * .. External Functions ..
2065 LOGICAL LZE, LZERES
2066 EXTERNAL LZE, LZERES
2067 * .. External Subroutines ..
2068 EXTERNAL ZHER2, ZHPR2, ZMAKE, ZMVCH
2069 * .. Intrinsic Functions ..
2070 INTRINSIC ABS, DCONJG, MAX
2071 * .. Scalars in Common ..
2072 INTEGER INFOT, NOUTC
2073 LOGICAL LERR, OK
2074 * .. Common blocks ..
2075 COMMON /INFOC/INFOT, NOUTC, OK, LERR
2076 * .. Data statements ..
2077 DATA ICH/'UL'/
2078 * .. Executable Statements ..
2079 FULL = SNAME( 3: 3 ).EQ.'E'
2080 PACKED = SNAME( 3: 3 ).EQ.'P'
2081 * Define the number of arguments.
2082 IF( FULL )THEN
2083 NARGS = 9
2084 ELSE IF( PACKED )THEN
2085 NARGS = 8
2086 END IF
2087 *
2088 NC = 0
2089 RESET = .TRUE.
2090 ERRMAX = RZERO
2091 *
2092 DO 140 IN = 1, NIDIM
2093 N = IDIM( IN )
2094 * Set LDA to 1 more than minimum value if room.
2095 LDA = N
2096 IF( LDA.LT.NMAX )
2097 $ LDA = LDA + 1
2098 * Skip tests if not enough room.
2099 IF( LDA.GT.NMAX )
2100 $ GO TO 140
2101 IF( PACKED )THEN
2102 LAA = ( N*( N + 1 ) )/2
2103 ELSE
2104 LAA = LDA*N
2105 END IF
2106 *
2107 DO 130 IC = 1, 2
2108 UPLO = ICH( IC: IC )
2109 UPPER = UPLO.EQ.'U'
2110 *
2111 DO 120 IX = 1, NINC
2112 INCX = INC( IX )
2113 LX = ABS( INCX )*N
2114 *
2115 * Generate the vector X.
2116 *
2117 TRANSL = HALF
2118 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX, ABS( INCX ),
2119 $ 0, N - 1, RESET, TRANSL )
2120 IF( N.GT.1 )THEN
2121 X( N/2 ) = ZERO
2122 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
2123 END IF
2124 *
2125 DO 110 IY = 1, NINC
2126 INCY = INC( IY )
2127 LY = ABS( INCY )*N
2128 *
2129 * Generate the vector Y.
2130 *
2131 TRANSL = ZERO
2132 CALL ZMAKE( 'GE', ' ', ' ', 1, N, Y, 1, YY,
2133 $ ABS( INCY ), 0, N - 1, RESET, TRANSL )
2134 IF( N.GT.1 )THEN
2135 Y( N/2 ) = ZERO
2136 YY( 1 + ABS( INCY )*( N/2 - 1 ) ) = ZERO
2137 END IF
2138 *
2139 DO 100 IA = 1, NALF
2140 ALPHA = ALF( IA )
2141 NULL = N.LE.0.OR.ALPHA.EQ.ZERO
2142 *
2143 * Generate the matrix A.
2144 *
2145 TRANSL = ZERO
2146 CALL ZMAKE( SNAME( 2: 3 ), UPLO, ' ', N, N, A,
2147 $ NMAX, AA, LDA, N - 1, N - 1, RESET,
2148 $ TRANSL )
2149 *
2150 NC = NC + 1
2151 *
2152 * Save every datum before calling the subroutine.
2153 *
2154 UPLOS = UPLO
2155 NS = N
2156 ALS = ALPHA
2157 DO 10 I = 1, LAA
2158 AS( I ) = AA( I )
2159 10 CONTINUE
2160 LDAS = LDA
2161 DO 20 I = 1, LX
2162 XS( I ) = XX( I )
2163 20 CONTINUE
2164 INCXS = INCX
2165 DO 30 I = 1, LY
2166 YS( I ) = YY( I )
2167 30 CONTINUE
2168 INCYS = INCY
2169 *
2170 * Call the subroutine.
2171 *
2172 IF( FULL )THEN
2173 IF( TRACE )
2174 $ WRITE( NTRA, FMT = 9993 )NC, SNAME, UPLO, N,
2175 $ ALPHA, INCX, INCY, LDA
2176 IF( REWI )
2177 $ REWIND NTRA
2178 CALL ZHER2( UPLO, N, ALPHA, XX, INCX, YY, INCY,
2179 $ AA, LDA )
2180 ELSE IF( PACKED )THEN
2181 IF( TRACE )
2182 $ WRITE( NTRA, FMT = 9994 )NC, SNAME, UPLO, N,
2183 $ ALPHA, INCX, INCY
2184 IF( REWI )
2185 $ REWIND NTRA
2186 CALL ZHPR2( UPLO, N, ALPHA, XX, INCX, YY, INCY,
2187 $ AA )
2188 END IF
2189 *
2190 * Check if error-exit was taken incorrectly.
2191 *
2192 IF( .NOT.OK )THEN
2193 WRITE( NOUT, FMT = 9992 )
2194 FATAL = .TRUE.
2195 GO TO 160
2196 END IF
2197 *
2198 * See what data changed inside subroutines.
2199 *
2200 ISAME( 1 ) = UPLO.EQ.UPLOS
2201 ISAME( 2 ) = NS.EQ.N
2202 ISAME( 3 ) = ALS.EQ.ALPHA
2203 ISAME( 4 ) = LZE( XS, XX, LX )
2204 ISAME( 5 ) = INCXS.EQ.INCX
2205 ISAME( 6 ) = LZE( YS, YY, LY )
2206 ISAME( 7 ) = INCYS.EQ.INCY
2207 IF( NULL )THEN
2208 ISAME( 8 ) = LZE( AS, AA, LAA )
2209 ELSE
2210 ISAME( 8 ) = LZERES( SNAME( 2: 3 ), UPLO, N, N,
2211 $ AS, AA, LDA )
2212 END IF
2213 IF( .NOT.PACKED )THEN
2214 ISAME( 9 ) = LDAS.EQ.LDA
2215 END IF
2216 *
2217 * If data was incorrectly changed, report and return.
2218 *
2219 SAME = .TRUE.
2220 DO 40 I = 1, NARGS
2221 SAME = SAME.AND.ISAME( I )
2222 IF( .NOT.ISAME( I ) )
2223 $ WRITE( NOUT, FMT = 9998 )I
2224 40 CONTINUE
2225 IF( .NOT.SAME )THEN
2226 FATAL = .TRUE.
2227 GO TO 160
2228 END IF
2229 *
2230 IF( .NOT.NULL )THEN
2231 *
2232 * Check the result column by column.
2233 *
2234 IF( INCX.GT.0 )THEN
2235 DO 50 I = 1, N
2236 Z( I, 1 ) = X( I )
2237 50 CONTINUE
2238 ELSE
2239 DO 60 I = 1, N
2240 Z( I, 1 ) = X( N - I + 1 )
2241 60 CONTINUE
2242 END IF
2243 IF( INCY.GT.0 )THEN
2244 DO 70 I = 1, N
2245 Z( I, 2 ) = Y( I )
2246 70 CONTINUE
2247 ELSE
2248 DO 80 I = 1, N
2249 Z( I, 2 ) = Y( N - I + 1 )
2250 80 CONTINUE
2251 END IF
2252 JA = 1
2253 DO 90 J = 1, N
2254 W( 1 ) = ALPHA*DCONJG( Z( J, 2 ) )
2255 W( 2 ) = DCONJG( ALPHA )*DCONJG( Z( J, 1 ) )
2256 IF( UPPER )THEN
2257 JJ = 1
2258 LJ = J
2259 ELSE
2260 JJ = J
2261 LJ = N - J + 1
2262 END IF
2263 CALL ZMVCH( 'N', LJ, 2, ONE, Z( JJ, 1 ),
2264 $ NMAX, W, 1, ONE, A( JJ, J ), 1,
2265 $ YT, G, AA( JA ), EPS, ERR, FATAL,
2266 $ NOUT, .TRUE. )
2267 IF( FULL )THEN
2268 IF( UPPER )THEN
2269 JA = JA + LDA
2270 ELSE
2271 JA = JA + LDA + 1
2272 END IF
2273 ELSE
2274 JA = JA + LJ
2275 END IF
2276 ERRMAX = MAX( ERRMAX, ERR )
2277 * If got really bad answer, report and return.
2278 IF( FATAL )
2279 $ GO TO 150
2280 90 CONTINUE
2281 ELSE
2282 * Avoid repeating tests with N.le.0.
2283 IF( N.LE.0 )
2284 $ GO TO 140
2285 END IF
2286 *
2287 100 CONTINUE
2288 *
2289 110 CONTINUE
2290 *
2291 120 CONTINUE
2292 *
2293 130 CONTINUE
2294 *
2295 140 CONTINUE
2296 *
2297 * Report result.
2298 *
2299 IF( ERRMAX.LT.THRESH )THEN
2300 WRITE( NOUT, FMT = 9999 )SNAME, NC
2301 ELSE
2302 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
2303 END IF
2304 GO TO 170
2305 *
2306 150 CONTINUE
2307 WRITE( NOUT, FMT = 9995 )J
2308 *
2309 160 CONTINUE
2310 WRITE( NOUT, FMT = 9996 )SNAME
2311 IF( FULL )THEN
2312 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, N, ALPHA, INCX,
2313 $ INCY, LDA
2314 ELSE IF( PACKED )THEN
2315 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, N, ALPHA, INCX, INCY
2316 END IF
2317 *
2318 170 CONTINUE
2319 RETURN
2320 *
2321 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
2322 $ 'S)' )
2323 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
2324 $ 'ANGED INCORRECTLY *******' )
2325 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
2326 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
2327 $ ' - SUSPECT *******' )
2328 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
2329 9995 FORMAT( ' THESE ARE THE RESULTS FOR COLUMN ', I3 )
2330 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
2331 $ F4.1, '), X,', I2, ', Y,', I2, ', AP) ',
2332 $ ' .' )
2333 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
2334 $ F4.1, '), X,', I2, ', Y,', I2, ', A,', I3, ') ',
2335 $ ' .' )
2336 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
2337 $ '******' )
2338 *
2339 * End of ZCHK6.
2340 *
2341 END
2342 SUBROUTINE ZCHKE( ISNUM, SRNAMT, NOUT )
2343 *
2344 * Tests the error exits from the Level 2 Blas.
2345 * Requires a special version of the error-handling routine XERBLA.
2346 * ALPHA, RALPHA, BETA, A, X and Y should not need to be defined.
2347 *
2348 * Auxiliary routine for test program for Level 2 Blas.
2349 *
2350 * -- Written on 10-August-1987.
2351 * Richard Hanson, Sandia National Labs.
2352 * Jeremy Du Croz, NAG Central Office.
2353 *
2354 * .. Scalar Arguments ..
2355 INTEGER ISNUM, NOUT
2356 CHARACTER*6 SRNAMT
2357 * .. Scalars in Common ..
2358 INTEGER INFOT, NOUTC
2359 LOGICAL LERR, OK
2360 * .. Local Scalars ..
2361 COMPLEX*16 ALPHA, BETA
2362 DOUBLE PRECISION RALPHA
2363 * .. Local Arrays ..
2364 COMPLEX*16 A( 1, 1 ), X( 1 ), Y( 1 )
2365 * .. External Subroutines ..
2366 EXTERNAL CHKXER, ZGBMV, ZGEMV, ZGERC, ZGERU, ZHBMV,
2367 $ ZHEMV, ZHER, ZHER2, ZHPMV, ZHPR, ZHPR2, ZTBMV,
2368 $ ZTBSV, ZTPMV, ZTPSV, ZTRMV, ZTRSV
2369 * .. Common blocks ..
2370 COMMON /INFOC/INFOT, NOUTC, OK, LERR
2371 * .. Executable Statements ..
2372 * OK is set to .FALSE. by the special version of XERBLA or by CHKXER
2373 * if anything is wrong.
2374 OK = .TRUE.
2375 * LERR is set to .TRUE. by the special version of XERBLA each time
2376 * it is called, and is then tested and re-set by CHKXER.
2377 LERR = .FALSE.
2378 GO TO ( 10, 20, 30, 40, 50, 60, 70, 80,
2379 $ 90, 100, 110, 120, 130, 140, 150, 160,
2380 $ 170 )ISNUM
2381 10 INFOT = 1
2382 CALL ZGEMV( '/', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2383 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2384 INFOT = 2
2385 CALL ZGEMV( 'N', -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2386 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2387 INFOT = 3
2388 CALL ZGEMV( 'N', 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2389 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2390 INFOT = 6
2391 CALL ZGEMV( 'N', 2, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2392 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2393 INFOT = 8
2394 CALL ZGEMV( 'N', 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2395 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2396 INFOT = 11
2397 CALL ZGEMV( 'N', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2398 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2399 GO TO 180
2400 20 INFOT = 1
2401 CALL ZGBMV( '/', 0, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2402 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2403 INFOT = 2
2404 CALL ZGBMV( 'N', -1, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2405 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2406 INFOT = 3
2407 CALL ZGBMV( 'N', 0, -1, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2408 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2409 INFOT = 4
2410 CALL ZGBMV( 'N', 0, 0, -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2411 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2412 INFOT = 5
2413 CALL ZGBMV( 'N', 2, 0, 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2414 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2415 INFOT = 8
2416 CALL ZGBMV( 'N', 0, 0, 1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2417 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2418 INFOT = 10
2419 CALL ZGBMV( 'N', 0, 0, 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2420 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2421 INFOT = 13
2422 CALL ZGBMV( 'N', 0, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2423 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2424 GO TO 180
2425 30 INFOT = 1
2426 CALL ZHEMV( '/', 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2427 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2428 INFOT = 2
2429 CALL ZHEMV( 'U', -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2430 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2431 INFOT = 5
2432 CALL ZHEMV( 'U', 2, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2433 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2434 INFOT = 7
2435 CALL ZHEMV( 'U', 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2436 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2437 INFOT = 10
2438 CALL ZHEMV( 'U', 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2439 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2440 GO TO 180
2441 40 INFOT = 1
2442 CALL ZHBMV( '/', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2443 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2444 INFOT = 2
2445 CALL ZHBMV( 'U', -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2446 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2447 INFOT = 3
2448 CALL ZHBMV( 'U', 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2449 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2450 INFOT = 6
2451 CALL ZHBMV( 'U', 0, 1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2452 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2453 INFOT = 8
2454 CALL ZHBMV( 'U', 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2455 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2456 INFOT = 11
2457 CALL ZHBMV( 'U', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2458 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2459 GO TO 180
2460 50 INFOT = 1
2461 CALL ZHPMV( '/', 0, ALPHA, A, X, 1, BETA, Y, 1 )
2462 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2463 INFOT = 2
2464 CALL ZHPMV( 'U', -1, ALPHA, A, X, 1, BETA, Y, 1 )
2465 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2466 INFOT = 6
2467 CALL ZHPMV( 'U', 0, ALPHA, A, X, 0, BETA, Y, 1 )
2468 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2469 INFOT = 9
2470 CALL ZHPMV( 'U', 0, ALPHA, A, X, 1, BETA, Y, 0 )
2471 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2472 GO TO 180
2473 60 INFOT = 1
2474 CALL ZTRMV( '/', 'N', 'N', 0, A, 1, X, 1 )
2475 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2476 INFOT = 2
2477 CALL ZTRMV( 'U', '/', 'N', 0, A, 1, X, 1 )
2478 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2479 INFOT = 3
2480 CALL ZTRMV( 'U', 'N', '/', 0, A, 1, X, 1 )
2481 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2482 INFOT = 4
2483 CALL ZTRMV( 'U', 'N', 'N', -1, A, 1, X, 1 )
2484 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2485 INFOT = 6
2486 CALL ZTRMV( 'U', 'N', 'N', 2, A, 1, X, 1 )
2487 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2488 INFOT = 8
2489 CALL ZTRMV( 'U', 'N', 'N', 0, A, 1, X, 0 )
2490 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2491 GO TO 180
2492 70 INFOT = 1
2493 CALL ZTBMV( '/', 'N', 'N', 0, 0, A, 1, X, 1 )
2494 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2495 INFOT = 2
2496 CALL ZTBMV( 'U', '/', 'N', 0, 0, A, 1, X, 1 )
2497 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2498 INFOT = 3
2499 CALL ZTBMV( 'U', 'N', '/', 0, 0, A, 1, X, 1 )
2500 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2501 INFOT = 4
2502 CALL ZTBMV( 'U', 'N', 'N', -1, 0, A, 1, X, 1 )
2503 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2504 INFOT = 5
2505 CALL ZTBMV( 'U', 'N', 'N', 0, -1, A, 1, X, 1 )
2506 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2507 INFOT = 7
2508 CALL ZTBMV( 'U', 'N', 'N', 0, 1, A, 1, X, 1 )
2509 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2510 INFOT = 9
2511 CALL ZTBMV( 'U', 'N', 'N', 0, 0, A, 1, X, 0 )
2512 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2513 GO TO 180
2514 80 INFOT = 1
2515 CALL ZTPMV( '/', 'N', 'N', 0, A, X, 1 )
2516 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2517 INFOT = 2
2518 CALL ZTPMV( 'U', '/', 'N', 0, A, X, 1 )
2519 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2520 INFOT = 3
2521 CALL ZTPMV( 'U', 'N', '/', 0, A, X, 1 )
2522 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2523 INFOT = 4
2524 CALL ZTPMV( 'U', 'N', 'N', -1, A, X, 1 )
2525 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2526 INFOT = 7
2527 CALL ZTPMV( 'U', 'N', 'N', 0, A, X, 0 )
2528 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2529 GO TO 180
2530 90 INFOT = 1
2531 CALL ZTRSV( '/', 'N', 'N', 0, A, 1, X, 1 )
2532 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2533 INFOT = 2
2534 CALL ZTRSV( 'U', '/', 'N', 0, A, 1, X, 1 )
2535 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2536 INFOT = 3
2537 CALL ZTRSV( 'U', 'N', '/', 0, A, 1, X, 1 )
2538 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2539 INFOT = 4
2540 CALL ZTRSV( 'U', 'N', 'N', -1, A, 1, X, 1 )
2541 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2542 INFOT = 6
2543 CALL ZTRSV( 'U', 'N', 'N', 2, A, 1, X, 1 )
2544 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2545 INFOT = 8
2546 CALL ZTRSV( 'U', 'N', 'N', 0, A, 1, X, 0 )
2547 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2548 GO TO 180
2549 100 INFOT = 1
2550 CALL ZTBSV( '/', 'N', 'N', 0, 0, A, 1, X, 1 )
2551 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2552 INFOT = 2
2553 CALL ZTBSV( 'U', '/', 'N', 0, 0, A, 1, X, 1 )
2554 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2555 INFOT = 3
2556 CALL ZTBSV( 'U', 'N', '/', 0, 0, A, 1, X, 1 )
2557 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2558 INFOT = 4
2559 CALL ZTBSV( 'U', 'N', 'N', -1, 0, A, 1, X, 1 )
2560 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2561 INFOT = 5
2562 CALL ZTBSV( 'U', 'N', 'N', 0, -1, A, 1, X, 1 )
2563 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2564 INFOT = 7
2565 CALL ZTBSV( 'U', 'N', 'N', 0, 1, A, 1, X, 1 )
2566 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2567 INFOT = 9
2568 CALL ZTBSV( 'U', 'N', 'N', 0, 0, A, 1, X, 0 )
2569 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2570 GO TO 180
2571 110 INFOT = 1
2572 CALL ZTPSV( '/', 'N', 'N', 0, A, X, 1 )
2573 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2574 INFOT = 2
2575 CALL ZTPSV( 'U', '/', 'N', 0, A, X, 1 )
2576 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2577 INFOT = 3
2578 CALL ZTPSV( 'U', 'N', '/', 0, A, X, 1 )
2579 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2580 INFOT = 4
2581 CALL ZTPSV( 'U', 'N', 'N', -1, A, X, 1 )
2582 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2583 INFOT = 7
2584 CALL ZTPSV( 'U', 'N', 'N', 0, A, X, 0 )
2585 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2586 GO TO 180
2587 120 INFOT = 1
2588 CALL ZGERC( -1, 0, ALPHA, X, 1, Y, 1, A, 1 )
2589 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2590 INFOT = 2
2591 CALL ZGERC( 0, -1, ALPHA, X, 1, Y, 1, A, 1 )
2592 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2593 INFOT = 5
2594 CALL ZGERC( 0, 0, ALPHA, X, 0, Y, 1, A, 1 )
2595 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2596 INFOT = 7
2597 CALL ZGERC( 0, 0, ALPHA, X, 1, Y, 0, A, 1 )
2598 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2599 INFOT = 9
2600 CALL ZGERC( 2, 0, ALPHA, X, 1, Y, 1, A, 1 )
2601 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2602 GO TO 180
2603 130 INFOT = 1
2604 CALL ZGERU( -1, 0, ALPHA, X, 1, Y, 1, A, 1 )
2605 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2606 INFOT = 2
2607 CALL ZGERU( 0, -1, ALPHA, X, 1, Y, 1, A, 1 )
2608 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2609 INFOT = 5
2610 CALL ZGERU( 0, 0, ALPHA, X, 0, Y, 1, A, 1 )
2611 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2612 INFOT = 7
2613 CALL ZGERU( 0, 0, ALPHA, X, 1, Y, 0, A, 1 )
2614 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2615 INFOT = 9
2616 CALL ZGERU( 2, 0, ALPHA, X, 1, Y, 1, A, 1 )
2617 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2618 GO TO 180
2619 140 INFOT = 1
2620 CALL ZHER( '/', 0, RALPHA, X, 1, A, 1 )
2621 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2622 INFOT = 2
2623 CALL ZHER( 'U', -1, RALPHA, X, 1, A, 1 )
2624 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2625 INFOT = 5
2626 CALL ZHER( 'U', 0, RALPHA, X, 0, A, 1 )
2627 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2628 INFOT = 7
2629 CALL ZHER( 'U', 2, RALPHA, X, 1, A, 1 )
2630 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2631 GO TO 180
2632 150 INFOT = 1
2633 CALL ZHPR( '/', 0, RALPHA, X, 1, A )
2634 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2635 INFOT = 2
2636 CALL ZHPR( 'U', -1, RALPHA, X, 1, A )
2637 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2638 INFOT = 5
2639 CALL ZHPR( 'U', 0, RALPHA, X, 0, A )
2640 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2641 GO TO 180
2642 160 INFOT = 1
2643 CALL ZHER2( '/', 0, ALPHA, X, 1, Y, 1, A, 1 )
2644 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2645 INFOT = 2
2646 CALL ZHER2( 'U', -1, ALPHA, X, 1, Y, 1, A, 1 )
2647 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2648 INFOT = 5
2649 CALL ZHER2( 'U', 0, ALPHA, X, 0, Y, 1, A, 1 )
2650 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2651 INFOT = 7
2652 CALL ZHER2( 'U', 0, ALPHA, X, 1, Y, 0, A, 1 )
2653 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2654 INFOT = 9
2655 CALL ZHER2( 'U', 2, ALPHA, X, 1, Y, 1, A, 1 )
2656 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2657 GO TO 180
2658 170 INFOT = 1
2659 CALL ZHPR2( '/', 0, ALPHA, X, 1, Y, 1, A )
2660 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2661 INFOT = 2
2662 CALL ZHPR2( 'U', -1, ALPHA, X, 1, Y, 1, A )
2663 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2664 INFOT = 5
2665 CALL ZHPR2( 'U', 0, ALPHA, X, 0, Y, 1, A )
2666 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2667 INFOT = 7
2668 CALL ZHPR2( 'U', 0, ALPHA, X, 1, Y, 0, A )
2669 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2670 *
2671 180 IF( OK )THEN
2672 WRITE( NOUT, FMT = 9999 )SRNAMT
2673 ELSE
2674 WRITE( NOUT, FMT = 9998 )SRNAMT
2675 END IF
2676 RETURN
2677 *
2678 9999 FORMAT( ' ', A6, ' PASSED THE TESTS OF ERROR-EXITS' )
2679 9998 FORMAT( ' ******* ', A6, ' FAILED THE TESTS OF ERROR-EXITS *****',
2680 $ '**' )
2681 *
2682 * End of ZCHKE.
2683 *
2684 END
2685 SUBROUTINE ZMAKE( TYPE, UPLO, DIAG, M, N, A, NMAX, AA, LDA, KL,
2686 $ KU, RESET, TRANSL )
2687 *
2688 * Generates values for an M by N matrix A within the bandwidth
2689 * defined by KL and KU.
2690 * Stores the values in the array AA in the data structure required
2691 * by the routine, with unwanted elements set to rogue value.
2692 *
2693 * TYPE is 'GE', 'GB', 'HE', 'HB', 'HP', 'TR', 'TB' OR 'TP'.
2694 *
2695 * Auxiliary routine for test program for Level 2 Blas.
2696 *
2697 * -- Written on 10-August-1987.
2698 * Richard Hanson, Sandia National Labs.
2699 * Jeremy Du Croz, NAG Central Office.
2700 *
2701 * .. Parameters ..
2702 COMPLEX*16 ZERO, ONE
2703 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
2704 $ ONE = ( 1.0D0, 0.0D0 ) )
2705 COMPLEX*16 ROGUE
2706 PARAMETER ( ROGUE = ( -1.0D10, 1.0D10 ) )
2707 DOUBLE PRECISION RZERO
2708 PARAMETER ( RZERO = 0.0D0 )
2709 DOUBLE PRECISION RROGUE
2710 PARAMETER ( RROGUE = -1.0D10 )
2711 * .. Scalar Arguments ..
2712 COMPLEX*16 TRANSL
2713 INTEGER KL, KU, LDA, M, N, NMAX
2714 LOGICAL RESET
2715 CHARACTER*1 DIAG, UPLO
2716 CHARACTER*2 TYPE
2717 * .. Array Arguments ..
2718 COMPLEX*16 A( NMAX, * ), AA( * )
2719 * .. Local Scalars ..
2720 INTEGER I, I1, I2, I3, IBEG, IEND, IOFF, J, JJ, KK
2721 LOGICAL GEN, LOWER, SYM, TRI, UNIT, UPPER
2722 * .. External Functions ..
2723 COMPLEX*16 ZBEG
2724 EXTERNAL ZBEG
2725 * .. Intrinsic Functions ..
2726 INTRINSIC DBLE, DCMPLX, DCONJG, MAX, MIN
2727 * .. Executable Statements ..
2728 GEN = TYPE( 1: 1 ).EQ.'G'
2729 SYM = TYPE( 1: 1 ).EQ.'H'
2730 TRI = TYPE( 1: 1 ).EQ.'T'
2731 UPPER = ( SYM.OR.TRI ).AND.UPLO.EQ.'U'
2732 LOWER = ( SYM.OR.TRI ).AND.UPLO.EQ.'L'
2733 UNIT = TRI.AND.DIAG.EQ.'U'
2734 *
2735 * Generate data in array A.
2736 *
2737 DO 20 J = 1, N
2738 DO 10 I = 1, M
2739 IF( GEN.OR.( UPPER.AND.I.LE.J ).OR.( LOWER.AND.I.GE.J ) )
2740 $ THEN
2741 IF( ( I.LE.J.AND.J - I.LE.KU ).OR.
2742 $ ( I.GE.J.AND.I - J.LE.KL ) )THEN
2743 A( I, J ) = ZBEG( RESET ) + TRANSL
2744 ELSE
2745 A( I, J ) = ZERO
2746 END IF
2747 IF( I.NE.J )THEN
2748 IF( SYM )THEN
2749 A( J, I ) = DCONJG( A( I, J ) )
2750 ELSE IF( TRI )THEN
2751 A( J, I ) = ZERO
2752 END IF
2753 END IF
2754 END IF
2755 10 CONTINUE
2756 IF( SYM )
2757 $ A( J, J ) = DCMPLX( DBLE( A( J, J ) ), RZERO )
2758 IF( TRI )
2759 $ A( J, J ) = A( J, J ) + ONE
2760 IF( UNIT )
2761 $ A( J, J ) = ONE
2762 20 CONTINUE
2763 *
2764 * Store elements in array AS in data structure required by routine.
2765 *
2766 IF( TYPE.EQ.'GE' )THEN
2767 DO 50 J = 1, N
2768 DO 30 I = 1, M
2769 AA( I + ( J - 1 )*LDA ) = A( I, J )
2770 30 CONTINUE
2771 DO 40 I = M + 1, LDA
2772 AA( I + ( J - 1 )*LDA ) = ROGUE
2773 40 CONTINUE
2774 50 CONTINUE
2775 ELSE IF( TYPE.EQ.'GB' )THEN
2776 DO 90 J = 1, N
2777 DO 60 I1 = 1, KU + 1 - J
2778 AA( I1 + ( J - 1 )*LDA ) = ROGUE
2779 60 CONTINUE
2780 DO 70 I2 = I1, MIN( KL + KU + 1, KU + 1 + M - J )
2781 AA( I2 + ( J - 1 )*LDA ) = A( I2 + J - KU - 1, J )
2782 70 CONTINUE
2783 DO 80 I3 = I2, LDA
2784 AA( I3 + ( J - 1 )*LDA ) = ROGUE
2785 80 CONTINUE
2786 90 CONTINUE
2787 ELSE IF( TYPE.EQ.'HE'.OR.TYPE.EQ.'TR' )THEN
2788 DO 130 J = 1, N
2789 IF( UPPER )THEN
2790 IBEG = 1
2791 IF( UNIT )THEN
2792 IEND = J - 1
2793 ELSE
2794 IEND = J
2795 END IF
2796 ELSE
2797 IF( UNIT )THEN
2798 IBEG = J + 1
2799 ELSE
2800 IBEG = J
2801 END IF
2802 IEND = N
2803 END IF
2804 DO 100 I = 1, IBEG - 1
2805 AA( I + ( J - 1 )*LDA ) = ROGUE
2806 100 CONTINUE
2807 DO 110 I = IBEG, IEND
2808 AA( I + ( J - 1 )*LDA ) = A( I, J )
2809 110 CONTINUE
2810 DO 120 I = IEND + 1, LDA
2811 AA( I + ( J - 1 )*LDA ) = ROGUE
2812 120 CONTINUE
2813 IF( SYM )THEN
2814 JJ = J + ( J - 1 )*LDA
2815 AA( JJ ) = DCMPLX( DBLE( AA( JJ ) ), RROGUE )
2816 END IF
2817 130 CONTINUE
2818 ELSE IF( TYPE.EQ.'HB'.OR.TYPE.EQ.'TB' )THEN
2819 DO 170 J = 1, N
2820 IF( UPPER )THEN
2821 KK = KL + 1
2822 IBEG = MAX( 1, KL + 2 - J )
2823 IF( UNIT )THEN
2824 IEND = KL
2825 ELSE
2826 IEND = KL + 1
2827 END IF
2828 ELSE
2829 KK = 1
2830 IF( UNIT )THEN
2831 IBEG = 2
2832 ELSE
2833 IBEG = 1
2834 END IF
2835 IEND = MIN( KL + 1, 1 + M - J )
2836 END IF
2837 DO 140 I = 1, IBEG - 1
2838 AA( I + ( J - 1 )*LDA ) = ROGUE
2839 140 CONTINUE
2840 DO 150 I = IBEG, IEND
2841 AA( I + ( J - 1 )*LDA ) = A( I + J - KK, J )
2842 150 CONTINUE
2843 DO 160 I = IEND + 1, LDA
2844 AA( I + ( J - 1 )*LDA ) = ROGUE
2845 160 CONTINUE
2846 IF( SYM )THEN
2847 JJ = KK + ( J - 1 )*LDA
2848 AA( JJ ) = DCMPLX( DBLE( AA( JJ ) ), RROGUE )
2849 END IF
2850 170 CONTINUE
2851 ELSE IF( TYPE.EQ.'HP'.OR.TYPE.EQ.'TP' )THEN
2852 IOFF = 0
2853 DO 190 J = 1, N
2854 IF( UPPER )THEN
2855 IBEG = 1
2856 IEND = J
2857 ELSE
2858 IBEG = J
2859 IEND = N
2860 END IF
2861 DO 180 I = IBEG, IEND
2862 IOFF = IOFF + 1
2863 AA( IOFF ) = A( I, J )
2864 IF( I.EQ.J )THEN
2865 IF( UNIT )
2866 $ AA( IOFF ) = ROGUE
2867 IF( SYM )
2868 $ AA( IOFF ) = DCMPLX( DBLE( AA( IOFF ) ), RROGUE )
2869 END IF
2870 180 CONTINUE
2871 190 CONTINUE
2872 END IF
2873 RETURN
2874 *
2875 * End of ZMAKE.
2876 *
2877 END
2878 SUBROUTINE ZMVCH( TRANS, M, N, ALPHA, A, NMAX, X, INCX, BETA, Y,
2879 $ INCY, YT, G, YY, EPS, ERR, FATAL, NOUT, MV )
2880 *
2881 * Checks the results of the computational tests.
2882 *
2883 * Auxiliary routine for test program for Level 2 Blas.
2884 *
2885 * -- Written on 10-August-1987.
2886 * Richard Hanson, Sandia National Labs.
2887 * Jeremy Du Croz, NAG Central Office.
2888 *
2889 * .. Parameters ..
2890 COMPLEX*16 ZERO
2891 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ) )
2892 DOUBLE PRECISION RZERO, RONE
2893 PARAMETER ( RZERO = 0.0D0, RONE = 1.0D0 )
2894 * .. Scalar Arguments ..
2895 COMPLEX*16 ALPHA, BETA
2896 DOUBLE PRECISION EPS, ERR
2897 INTEGER INCX, INCY, M, N, NMAX, NOUT
2898 LOGICAL FATAL, MV
2899 CHARACTER*1 TRANS
2900 * .. Array Arguments ..
2901 COMPLEX*16 A( NMAX, * ), X( * ), Y( * ), YT( * ), YY( * )
2902 DOUBLE PRECISION G( * )
2903 * .. Local Scalars ..
2904 COMPLEX*16 C
2905 DOUBLE PRECISION ERRI
2906 INTEGER I, INCXL, INCYL, IY, J, JX, KX, KY, ML, NL
2907 LOGICAL CTRAN, TRAN
2908 * .. Intrinsic Functions ..
2909 INTRINSIC ABS, DBLE, DCONJG, DIMAG, MAX, SQRT
2910 * .. Statement Functions ..
2911 DOUBLE PRECISION ABS1
2912 * .. Statement Function definitions ..
2913 ABS1( C ) = ABS( DBLE( C ) ) + ABS( DIMAG( C ) )
2914 * .. Executable Statements ..
2915 TRAN = TRANS.EQ.'T'
2916 CTRAN = TRANS.EQ.'C'
2917 IF( TRAN.OR.CTRAN )THEN
2918 ML = N
2919 NL = M
2920 ELSE
2921 ML = M
2922 NL = N
2923 END IF
2924 IF( INCX.LT.0 )THEN
2925 KX = NL
2926 INCXL = -1
2927 ELSE
2928 KX = 1
2929 INCXL = 1
2930 END IF
2931 IF( INCY.LT.0 )THEN
2932 KY = ML
2933 INCYL = -1
2934 ELSE
2935 KY = 1
2936 INCYL = 1
2937 END IF
2938 *
2939 * Compute expected result in YT using data in A, X and Y.
2940 * Compute gauges in G.
2941 *
2942 IY = KY
2943 DO 40 I = 1, ML
2944 YT( IY ) = ZERO
2945 G( IY ) = RZERO
2946 JX = KX
2947 IF( TRAN )THEN
2948 DO 10 J = 1, NL
2949 YT( IY ) = YT( IY ) + A( J, I )*X( JX )
2950 G( IY ) = G( IY ) + ABS1( A( J, I ) )*ABS1( X( JX ) )
2951 JX = JX + INCXL
2952 10 CONTINUE
2953 ELSE IF( CTRAN )THEN
2954 DO 20 J = 1, NL
2955 YT( IY ) = YT( IY ) + DCONJG( A( J, I ) )*X( JX )
2956 G( IY ) = G( IY ) + ABS1( A( J, I ) )*ABS1( X( JX ) )
2957 JX = JX + INCXL
2958 20 CONTINUE
2959 ELSE
2960 DO 30 J = 1, NL
2961 YT( IY ) = YT( IY ) + A( I, J )*X( JX )
2962 G( IY ) = G( IY ) + ABS1( A( I, J ) )*ABS1( X( JX ) )
2963 JX = JX + INCXL
2964 30 CONTINUE
2965 END IF
2966 YT( IY ) = ALPHA*YT( IY ) + BETA*Y( IY )
2967 G( IY ) = ABS1( ALPHA )*G( IY ) + ABS1( BETA )*ABS1( Y( IY ) )
2968 IY = IY + INCYL
2969 40 CONTINUE
2970 *
2971 * Compute the error ratio for this result.
2972 *
2973 ERR = ZERO
2974 DO 50 I = 1, ML
2975 ERRI = ABS( YT( I ) - YY( 1 + ( I - 1 )*ABS( INCY ) ) )/EPS
2976 IF( G( I ).NE.RZERO )
2977 $ ERRI = ERRI/G( I )
2978 ERR = MAX( ERR, ERRI )
2979 IF( ERR*SQRT( EPS ).GE.RONE )
2980 $ GO TO 60
2981 50 CONTINUE
2982 * If the loop completes, all results are at least half accurate.
2983 GO TO 80
2984 *
2985 * Report fatal error.
2986 *
2987 60 FATAL = .TRUE.
2988 WRITE( NOUT, FMT = 9999 )
2989 DO 70 I = 1, ML
2990 IF( MV )THEN
2991 WRITE( NOUT, FMT = 9998 )I, YT( I ),
2992 $ YY( 1 + ( I - 1 )*ABS( INCY ) )
2993 ELSE
2994 WRITE( NOUT, FMT = 9998 )I,
2995 $ YY( 1 + ( I - 1 )*ABS( INCY ) ), YT( I )
2996 END IF
2997 70 CONTINUE
2998 *
2999 80 CONTINUE
3000 RETURN
3001 *
3002 9999 FORMAT( ' ******* FATAL ERROR - COMPUTED RESULT IS LESS THAN HAL',
3003 $ 'F ACCURATE *******', /' EXPECTED RE',
3004 $ 'SULT COMPUTED RESULT' )
3005 9998 FORMAT( 1X, I7, 2( ' (', G15.6, ',', G15.6, ')' ) )
3006 *
3007 * End of ZMVCH.
3008 *
3009 END
3010 LOGICAL FUNCTION LZE( RI, RJ, LR )
3011 *
3012 * Tests if two arrays are identical.
3013 *
3014 * Auxiliary routine for test program for Level 2 Blas.
3015 *
3016 * -- Written on 10-August-1987.
3017 * Richard Hanson, Sandia National Labs.
3018 * Jeremy Du Croz, NAG Central Office.
3019 *
3020 * .. Scalar Arguments ..
3021 INTEGER LR
3022 * .. Array Arguments ..
3023 COMPLEX*16 RI( * ), RJ( * )
3024 * .. Local Scalars ..
3025 INTEGER I
3026 * .. Executable Statements ..
3027 DO 10 I = 1, LR
3028 IF( RI( I ).NE.RJ( I ) )
3029 $ GO TO 20
3030 10 CONTINUE
3031 LZE = .TRUE.
3032 GO TO 30
3033 20 CONTINUE
3034 LZE = .FALSE.
3035 30 RETURN
3036 *
3037 * End of LZE.
3038 *
3039 END
3040 LOGICAL FUNCTION LZERES( TYPE, UPLO, M, N, AA, AS, LDA )
3041 *
3042 * Tests if selected elements in two arrays are equal.
3043 *
3044 * TYPE is 'GE', 'HE' or 'HP'.
3045 *
3046 * Auxiliary routine for test program for Level 2 Blas.
3047 *
3048 * -- Written on 10-August-1987.
3049 * Richard Hanson, Sandia National Labs.
3050 * Jeremy Du Croz, NAG Central Office.
3051 *
3052 * .. Scalar Arguments ..
3053 INTEGER LDA, M, N
3054 CHARACTER*1 UPLO
3055 CHARACTER*2 TYPE
3056 * .. Array Arguments ..
3057 COMPLEX*16 AA( LDA, * ), AS( LDA, * )
3058 * .. Local Scalars ..
3059 INTEGER I, IBEG, IEND, J
3060 LOGICAL UPPER
3061 * .. Executable Statements ..
3062 UPPER = UPLO.EQ.'U'
3063 IF( TYPE.EQ.'GE' )THEN
3064 DO 20 J = 1, N
3065 DO 10 I = M + 1, LDA
3066 IF( AA( I, J ).NE.AS( I, J ) )
3067 $ GO TO 70
3068 10 CONTINUE
3069 20 CONTINUE
3070 ELSE IF( TYPE.EQ.'HE' )THEN
3071 DO 50 J = 1, N
3072 IF( UPPER )THEN
3073 IBEG = 1
3074 IEND = J
3075 ELSE
3076 IBEG = J
3077 IEND = N
3078 END IF
3079 DO 30 I = 1, IBEG - 1
3080 IF( AA( I, J ).NE.AS( I, J ) )
3081 $ GO TO 70
3082 30 CONTINUE
3083 DO 40 I = IEND + 1, LDA
3084 IF( AA( I, J ).NE.AS( I, J ) )
3085 $ GO TO 70
3086 40 CONTINUE
3087 50 CONTINUE
3088 END IF
3089 *
3090 60 CONTINUE
3091 LZERES = .TRUE.
3092 GO TO 80
3093 70 CONTINUE
3094 LZERES = .FALSE.
3095 80 RETURN
3096 *
3097 * End of LZERES.
3098 *
3099 END
3100 COMPLEX*16 FUNCTION ZBEG( RESET )
3101 *
3102 * Generates complex numbers as pairs of random numbers uniformly
3103 * distributed between -0.5 and 0.5.
3104 *
3105 * Auxiliary routine for test program for Level 2 Blas.
3106 *
3107 * -- Written on 10-August-1987.
3108 * Richard Hanson, Sandia National Labs.
3109 * Jeremy Du Croz, NAG Central Office.
3110 *
3111 * .. Scalar Arguments ..
3112 LOGICAL RESET
3113 * .. Local Scalars ..
3114 INTEGER I, IC, J, MI, MJ
3115 * .. Save statement ..
3116 SAVE I, IC, J, MI, MJ
3117 * .. Intrinsic Functions ..
3118 INTRINSIC DCMPLX
3119 * .. Executable Statements ..
3120 IF( RESET )THEN
3121 * Initialize local variables.
3122 MI = 891
3123 MJ = 457
3124 I = 7
3125 J = 7
3126 IC = 0
3127 RESET = .FALSE.
3128 END IF
3129 *
3130 * The sequence of values of I or J is bounded between 1 and 999.
3131 * If initial I or J = 1,2,3,6,7 or 9, the period will be 50.
3132 * If initial I or J = 4 or 8, the period will be 25.
3133 * If initial I or J = 5, the period will be 10.
3134 * IC is used to break up the period by skipping 1 value of I or J
3135 * in 6.
3136 *
3137 IC = IC + 1
3138 10 I = I*MI
3139 J = J*MJ
3140 I = I - 1000*( I/1000 )
3141 J = J - 1000*( J/1000 )
3142 IF( IC.GE.5 )THEN
3143 IC = 0
3144 GO TO 10
3145 END IF
3146 ZBEG = DCMPLX( ( I - 500 )/1001.0D0, ( J - 500 )/1001.0D0 )
3147 RETURN
3148 *
3149 * End of ZBEG.
3150 *
3151 END
3152 DOUBLE PRECISION FUNCTION DDIFF( X, Y )
3153 *
3154 * Auxiliary routine for test program for Level 2 Blas.
3155 *
3156 * -- Written on 10-August-1987.
3157 * Richard Hanson, Sandia National Labs.
3158 *
3159 * .. Scalar Arguments ..
3160 DOUBLE PRECISION X, Y
3161 * .. Executable Statements ..
3162 DDIFF = X - Y
3163 RETURN
3164 *
3165 * End of DDIFF.
3166 *
3167 END
3168 SUBROUTINE CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
3169 *
3170 * Tests whether XERBLA has detected an error when it should.
3171 *
3172 * Auxiliary routine for test program for Level 2 Blas.
3173 *
3174 * -- Written on 10-August-1987.
3175 * Richard Hanson, Sandia National Labs.
3176 * Jeremy Du Croz, NAG Central Office.
3177 *
3178 * .. Scalar Arguments ..
3179 INTEGER INFOT, NOUT
3180 LOGICAL LERR, OK
3181 CHARACTER*6 SRNAMT
3182 * .. Executable Statements ..
3183 IF( .NOT.LERR )THEN
3184 WRITE( NOUT, FMT = 9999 )INFOT, SRNAMT
3185 OK = .FALSE.
3186 END IF
3187 LERR = .FALSE.
3188 RETURN
3189 *
3190 9999 FORMAT( ' ***** ILLEGAL VALUE OF PARAMETER NUMBER ', I2, ' NOT D',
3191 $ 'ETECTED BY ', A6, ' *****' )
3192 *
3193 * End of CHKXER.
3194 *
3195 END
3196 SUBROUTINE XERBLA( SRNAME, INFO )
3197 *
3198 * This is a special version of XERBLA to be used only as part of
3199 * the test program for testing error exits from the Level 2 BLAS
3200 * routines.
3201 *
3202 * XERBLA is an error handler for the Level 2 BLAS routines.
3203 *
3204 * It is called by the Level 2 BLAS routines if an input parameter is
3205 * invalid.
3206 *
3207 * Auxiliary routine for test program for Level 2 Blas.
3208 *
3209 * -- Written on 10-August-1987.
3210 * Richard Hanson, Sandia National Labs.
3211 * Jeremy Du Croz, NAG Central Office.
3212 *
3213 * .. Scalar Arguments ..
3214 INTEGER INFO
3215 CHARACTER*6 SRNAME
3216 * .. Scalars in Common ..
3217 INTEGER INFOT, NOUT
3218 LOGICAL LERR, OK
3219 CHARACTER*6 SRNAMT
3220 * .. Common blocks ..
3221 COMMON /INFOC/INFOT, NOUT, OK, LERR
3222 COMMON /SRNAMC/SRNAMT
3223 * .. Executable Statements ..
3224 LERR = .TRUE.
3225 IF( INFO.NE.INFOT )THEN
3226 IF( INFOT.NE.0 )THEN
3227 WRITE( NOUT, FMT = 9999 )INFO, INFOT
3228 ELSE
3229 WRITE( NOUT, FMT = 9997 )INFO
3230 END IF
3231 OK = .FALSE.
3232 END IF
3233 IF( SRNAME.NE.SRNAMT )THEN
3234 WRITE( NOUT, FMT = 9998 )SRNAME, SRNAMT
3235 OK = .FALSE.
3236 END IF
3237 RETURN
3238 *
3239 9999 FORMAT( ' ******* XERBLA WAS CALLED WITH INFO = ', I6, ' INSTEAD',
3240 $ ' OF ', I2, ' *******' )
3241 9998 FORMAT( ' ******* XERBLA WAS CALLED WITH SRNAME = ', A6, ' INSTE',
3242 $ 'AD OF ', A6, ' *******' )
3243 9997 FORMAT( ' ******* XERBLA WAS CALLED WITH INFO = ', I6,
3244 $ ' *******' )
3245 *
3246 * End of XERBLA
3247 *
3248 END
3249
2 *
3 * Test program for the COMPLEX*16 Level 2 Blas.
4 *
5 * The program must be driven by a short data file. The first 18 records
6 * of the file are read using list-directed input, the last 17 records
7 * are read using the format ( A6, L2 ). An annotated example of a data
8 * file can be obtained by deleting the first 3 characters from the
9 * following 35 lines:
10 * 'ZBLAT2.SUMM' NAME OF SUMMARY OUTPUT FILE
11 * 6 UNIT NUMBER OF SUMMARY FILE
12 * 'CBLA2T.SNAP' NAME OF SNAPSHOT OUTPUT FILE
13 * -1 UNIT NUMBER OF SNAPSHOT FILE (NOT USED IF .LT. 0)
14 * F LOGICAL FLAG, T TO REWIND SNAPSHOT FILE AFTER EACH RECORD.
15 * F LOGICAL FLAG, T TO STOP ON FAILURES.
16 * T LOGICAL FLAG, T TO TEST ERROR EXITS.
17 * 16.0 THRESHOLD VALUE OF TEST RATIO
18 * 6 NUMBER OF VALUES OF N
19 * 0 1 2 3 5 9 VALUES OF N
20 * 4 NUMBER OF VALUES OF K
21 * 0 1 2 4 VALUES OF K
22 * 4 NUMBER OF VALUES OF INCX AND INCY
23 * 1 2 -1 -2 VALUES OF INCX AND INCY
24 * 3 NUMBER OF VALUES OF ALPHA
25 * (0.0,0.0) (1.0,0.0) (0.7,-0.9) VALUES OF ALPHA
26 * 3 NUMBER OF VALUES OF BETA
27 * (0.0,0.0) (1.0,0.0) (1.3,-1.1) VALUES OF BETA
28 * ZGEMV T PUT F FOR NO TEST. SAME COLUMNS.
29 * ZGBMV T PUT F FOR NO TEST. SAME COLUMNS.
30 * ZHEMV T PUT F FOR NO TEST. SAME COLUMNS.
31 * ZHBMV T PUT F FOR NO TEST. SAME COLUMNS.
32 * ZHPMV T PUT F FOR NO TEST. SAME COLUMNS.
33 * ZTRMV T PUT F FOR NO TEST. SAME COLUMNS.
34 * ZTBMV T PUT F FOR NO TEST. SAME COLUMNS.
35 * ZTPMV T PUT F FOR NO TEST. SAME COLUMNS.
36 * ZTRSV T PUT F FOR NO TEST. SAME COLUMNS.
37 * ZTBSV T PUT F FOR NO TEST. SAME COLUMNS.
38 * ZTPSV T PUT F FOR NO TEST. SAME COLUMNS.
39 * ZGERC T PUT F FOR NO TEST. SAME COLUMNS.
40 * ZGERU T PUT F FOR NO TEST. SAME COLUMNS.
41 * ZHER T PUT F FOR NO TEST. SAME COLUMNS.
42 * ZHPR T PUT F FOR NO TEST. SAME COLUMNS.
43 * ZHER2 T PUT F FOR NO TEST. SAME COLUMNS.
44 * ZHPR2 T PUT F FOR NO TEST. SAME COLUMNS.
45 *
46 * See:
47 *
48 * Dongarra J. J., Du Croz J. J., Hammarling S. and Hanson R. J..
49 * An extended set of Fortran Basic Linear Algebra Subprograms.
50 *
51 * Technical Memoranda Nos. 41 (revision 3) and 81, Mathematics
52 * and Computer Science Division, Argonne National Laboratory,
53 * 9700 South Cass Avenue, Argonne, Illinois 60439, US.
54 *
55 * Or
56 *
57 * NAG Technical Reports TR3/87 and TR4/87, Numerical Algorithms
58 * Group Ltd., NAG Central Office, 256 Banbury Road, Oxford
59 * OX2 7DE, UK, and Numerical Algorithms Group Inc., 1101 31st
60 * Street, Suite 100, Downers Grove, Illinois 60515-1263, USA.
61 *
62 *
63 * -- Written on 10-August-1987.
64 * Richard Hanson, Sandia National Labs.
65 * Jeremy Du Croz, NAG Central Office.
66 *
67 * .. Parameters ..
68 INTEGER NIN
69 PARAMETER ( NIN = 5 )
70 INTEGER NSUBS
71 PARAMETER ( NSUBS = 17 )
72 COMPLEX*16 ZERO, ONE
73 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
74 $ ONE = ( 1.0D0, 0.0D0 ) )
75 DOUBLE PRECISION RZERO, RHALF, RONE
76 PARAMETER ( RZERO = 0.0D0, RHALF = 0.5D0, RONE = 1.0D0 )
77 INTEGER NMAX, INCMAX
78 PARAMETER ( NMAX = 65, INCMAX = 2 )
79 INTEGER NINMAX, NIDMAX, NKBMAX, NALMAX, NBEMAX
80 PARAMETER ( NINMAX = 7, NIDMAX = 9, NKBMAX = 7,
81 $ NALMAX = 7, NBEMAX = 7 )
82 * .. Local Scalars ..
83 DOUBLE PRECISION EPS, ERR, THRESH
84 INTEGER I, ISNUM, J, N, NALF, NBET, NIDIM, NINC, NKB,
85 $ NOUT, NTRA
86 LOGICAL FATAL, LTESTT, REWI, SAME, SFATAL, TRACE,
87 $ TSTERR
88 CHARACTER*1 TRANS
89 CHARACTER*6 SNAMET
90 CHARACTER*32 SNAPS, SUMMRY
91 * .. Local Arrays ..
92 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ),
93 $ ALF( NALMAX ), AS( NMAX*NMAX ), BET( NBEMAX ),
94 $ X( NMAX ), XS( NMAX*INCMAX ),
95 $ XX( NMAX*INCMAX ), Y( NMAX ),
96 $ YS( NMAX*INCMAX ), YT( NMAX ),
97 $ YY( NMAX*INCMAX ), Z( 2*NMAX )
98 DOUBLE PRECISION G( NMAX )
99 INTEGER IDIM( NIDMAX ), INC( NINMAX ), KB( NKBMAX )
100 LOGICAL LTEST( NSUBS )
101 CHARACTER*6 SNAMES( NSUBS )
102 * .. External Functions ..
103 DOUBLE PRECISION DDIFF
104 LOGICAL LZE
105 EXTERNAL DDIFF, LZE
106 * .. External Subroutines ..
107 EXTERNAL ZCHK1, ZCHK2, ZCHK3, ZCHK4, ZCHK5, ZCHK6,
108 $ ZCHKE, ZMVCH
109 * .. Intrinsic Functions ..
110 INTRINSIC ABS, MAX, MIN
111 * .. Scalars in Common ..
112 INTEGER INFOT, NOUTC
113 LOGICAL LERR, OK
114 CHARACTER*6 SRNAMT
115 * .. Common blocks ..
116 COMMON /INFOC/INFOT, NOUTC, OK, LERR
117 COMMON /SRNAMC/SRNAMT
118 * .. Data statements ..
119 DATA SNAMES/'ZGEMV ', 'ZGBMV ', 'ZHEMV ', 'ZHBMV ',
120 $ 'ZHPMV ', 'ZTRMV ', 'ZTBMV ', 'ZTPMV ',
121 $ 'ZTRSV ', 'ZTBSV ', 'ZTPSV ', 'ZGERC ',
122 $ 'ZGERU ', 'ZHER ', 'ZHPR ', 'ZHER2 ',
123 $ 'ZHPR2 '/
124 * .. Executable Statements ..
125 *
126 * Read name and unit number for summary output file and open file.
127 *
128 READ( NIN, FMT = * )SUMMRY
129 READ( NIN, FMT = * )NOUT
130 OPEN( NOUT, FILE = SUMMRY, STATUS = 'NEW' )
131 NOUTC = NOUT
132 *
133 * Read name and unit number for snapshot output file and open file.
134 *
135 READ( NIN, FMT = * )SNAPS
136 READ( NIN, FMT = * )NTRA
137 TRACE = NTRA.GE.0
138 IF( TRACE )THEN
139 OPEN( NTRA, FILE = SNAPS, STATUS = 'NEW' )
140 END IF
141 * Read the flag that directs rewinding of the snapshot file.
142 READ( NIN, FMT = * )REWI
143 REWI = REWI.AND.TRACE
144 * Read the flag that directs stopping on any failure.
145 READ( NIN, FMT = * )SFATAL
146 * Read the flag that indicates whether error exits are to be tested.
147 READ( NIN, FMT = * )TSTERR
148 * Read the threshold value of the test ratio
149 READ( NIN, FMT = * )THRESH
150 *
151 * Read and check the parameter values for the tests.
152 *
153 * Values of N
154 READ( NIN, FMT = * )NIDIM
155 IF( NIDIM.LT.1.OR.NIDIM.GT.NIDMAX )THEN
156 WRITE( NOUT, FMT = 9997 )'N', NIDMAX
157 GO TO 230
158 END IF
159 READ( NIN, FMT = * )( IDIM( I ), I = 1, NIDIM )
160 DO 10 I = 1, NIDIM
161 IF( IDIM( I ).LT.0.OR.IDIM( I ).GT.NMAX )THEN
162 WRITE( NOUT, FMT = 9996 )NMAX
163 GO TO 230
164 END IF
165 10 CONTINUE
166 * Values of K
167 READ( NIN, FMT = * )NKB
168 IF( NKB.LT.1.OR.NKB.GT.NKBMAX )THEN
169 WRITE( NOUT, FMT = 9997 )'K', NKBMAX
170 GO TO 230
171 END IF
172 READ( NIN, FMT = * )( KB( I ), I = 1, NKB )
173 DO 20 I = 1, NKB
174 IF( KB( I ).LT.0 )THEN
175 WRITE( NOUT, FMT = 9995 )
176 GO TO 230
177 END IF
178 20 CONTINUE
179 * Values of INCX and INCY
180 READ( NIN, FMT = * )NINC
181 IF( NINC.LT.1.OR.NINC.GT.NINMAX )THEN
182 WRITE( NOUT, FMT = 9997 )'INCX AND INCY', NINMAX
183 GO TO 230
184 END IF
185 READ( NIN, FMT = * )( INC( I ), I = 1, NINC )
186 DO 30 I = 1, NINC
187 IF( INC( I ).EQ.0.OR.ABS( INC( I ) ).GT.INCMAX )THEN
188 WRITE( NOUT, FMT = 9994 )INCMAX
189 GO TO 230
190 END IF
191 30 CONTINUE
192 * Values of ALPHA
193 READ( NIN, FMT = * )NALF
194 IF( NALF.LT.1.OR.NALF.GT.NALMAX )THEN
195 WRITE( NOUT, FMT = 9997 )'ALPHA', NALMAX
196 GO TO 230
197 END IF
198 READ( NIN, FMT = * )( ALF( I ), I = 1, NALF )
199 * Values of BETA
200 READ( NIN, FMT = * )NBET
201 IF( NBET.LT.1.OR.NBET.GT.NBEMAX )THEN
202 WRITE( NOUT, FMT = 9997 )'BETA', NBEMAX
203 GO TO 230
204 END IF
205 READ( NIN, FMT = * )( BET( I ), I = 1, NBET )
206 *
207 * Report values of parameters.
208 *
209 WRITE( NOUT, FMT = 9993 )
210 WRITE( NOUT, FMT = 9992 )( IDIM( I ), I = 1, NIDIM )
211 WRITE( NOUT, FMT = 9991 )( KB( I ), I = 1, NKB )
212 WRITE( NOUT, FMT = 9990 )( INC( I ), I = 1, NINC )
213 WRITE( NOUT, FMT = 9989 )( ALF( I ), I = 1, NALF )
214 WRITE( NOUT, FMT = 9988 )( BET( I ), I = 1, NBET )
215 IF( .NOT.TSTERR )THEN
216 WRITE( NOUT, FMT = * )
217 WRITE( NOUT, FMT = 9980 )
218 END IF
219 WRITE( NOUT, FMT = * )
220 WRITE( NOUT, FMT = 9999 )THRESH
221 WRITE( NOUT, FMT = * )
222 *
223 * Read names of subroutines and flags which indicate
224 * whether they are to be tested.
225 *
226 DO 40 I = 1, NSUBS
227 LTEST( I ) = .FALSE.
228 40 CONTINUE
229 50 READ( NIN, FMT = 9984, END = 80 )SNAMET, LTESTT
230 DO 60 I = 1, NSUBS
231 IF( SNAMET.EQ.SNAMES( I ) )
232 $ GO TO 70
233 60 CONTINUE
234 WRITE( NOUT, FMT = 9986 )SNAMET
235 STOP
236 70 LTEST( I ) = LTESTT
237 GO TO 50
238 *
239 80 CONTINUE
240 CLOSE ( NIN )
241 *
242 * Compute EPS (the machine precision).
243 *
244 EPS = RONE
245 90 CONTINUE
246 IF( DDIFF( RONE + EPS, RONE ).EQ.RZERO )
247 $ GO TO 100
248 EPS = RHALF*EPS
249 GO TO 90
250 100 CONTINUE
251 EPS = EPS + EPS
252 WRITE( NOUT, FMT = 9998 )EPS
253 *
254 * Check the reliability of ZMVCH using exact data.
255 *
256 N = MIN( 32, NMAX )
257 DO 120 J = 1, N
258 DO 110 I = 1, N
259 A( I, J ) = MAX( I - J + 1, 0 )
260 110 CONTINUE
261 X( J ) = J
262 Y( J ) = ZERO
263 120 CONTINUE
264 DO 130 J = 1, N
265 YY( J ) = J*( ( J + 1 )*J )/2 - ( ( J + 1 )*J*( J - 1 ) )/3
266 130 CONTINUE
267 * YY holds the exact result. On exit from ZMVCH YT holds
268 * the result computed by ZMVCH.
269 TRANS = 'N'
270 CALL ZMVCH( TRANS, N, N, ONE, A, NMAX, X, 1, ZERO, Y, 1, YT, G,
271 $ YY, EPS, ERR, FATAL, NOUT, .TRUE. )
272 SAME = LZE( YY, YT, N )
273 IF( .NOT.SAME.OR.ERR.NE.RZERO )THEN
274 WRITE( NOUT, FMT = 9985 )TRANS, SAME, ERR
275 STOP
276 END IF
277 TRANS = 'T'
278 CALL ZMVCH( TRANS, N, N, ONE, A, NMAX, X, -1, ZERO, Y, -1, YT, G,
279 $ YY, EPS, ERR, FATAL, NOUT, .TRUE. )
280 SAME = LZE( YY, YT, N )
281 IF( .NOT.SAME.OR.ERR.NE.RZERO )THEN
282 WRITE( NOUT, FMT = 9985 )TRANS, SAME, ERR
283 STOP
284 END IF
285 *
286 * Test each subroutine in turn.
287 *
288 DO 210 ISNUM = 1, NSUBS
289 WRITE( NOUT, FMT = * )
290 IF( .NOT.LTEST( ISNUM ) )THEN
291 * Subprogram is not to be tested.
292 WRITE( NOUT, FMT = 9983 )SNAMES( ISNUM )
293 ELSE
294 SRNAMT = SNAMES( ISNUM )
295 * Test error exits.
296 IF( TSTERR )THEN
297 CALL ZCHKE( ISNUM, SNAMES( ISNUM ), NOUT )
298 WRITE( NOUT, FMT = * )
299 END IF
300 * Test computations.
301 INFOT = 0
302 OK = .TRUE.
303 FATAL = .FALSE.
304 GO TO ( 140, 140, 150, 150, 150, 160, 160,
305 $ 160, 160, 160, 160, 170, 170, 180,
306 $ 180, 190, 190 )ISNUM
307 * Test ZGEMV, 01, and ZGBMV, 02.
308 140 CALL ZCHK1( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
309 $ REWI, FATAL, NIDIM, IDIM, NKB, KB, NALF, ALF,
310 $ NBET, BET, NINC, INC, NMAX, INCMAX, A, AA, AS,
311 $ X, XX, XS, Y, YY, YS, YT, G )
312 GO TO 200
313 * Test ZHEMV, 03, ZHBMV, 04, and ZHPMV, 05.
314 150 CALL ZCHK2( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
315 $ REWI, FATAL, NIDIM, IDIM, NKB, KB, NALF, ALF,
316 $ NBET, BET, NINC, INC, NMAX, INCMAX, A, AA, AS,
317 $ X, XX, XS, Y, YY, YS, YT, G )
318 GO TO 200
319 * Test ZTRMV, 06, ZTBMV, 07, ZTPMV, 08,
320 * ZTRSV, 09, ZTBSV, 10, and ZTPSV, 11.
321 160 CALL ZCHK3( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
322 $ REWI, FATAL, NIDIM, IDIM, NKB, KB, NINC, INC,
323 $ NMAX, INCMAX, A, AA, AS, Y, YY, YS, YT, G, Z )
324 GO TO 200
325 * Test ZGERC, 12, ZGERU, 13.
326 170 CALL ZCHK4( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
327 $ REWI, FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC,
328 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
329 $ YT, G, Z )
330 GO TO 200
331 * Test ZHER, 14, and ZHPR, 15.
332 180 CALL ZCHK5( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
333 $ REWI, FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC,
334 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
335 $ YT, G, Z )
336 GO TO 200
337 * Test ZHER2, 16, and ZHPR2, 17.
338 190 CALL ZCHK6( SNAMES( ISNUM ), EPS, THRESH, NOUT, NTRA, TRACE,
339 $ REWI, FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC,
340 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
341 $ YT, G, Z )
342 *
343 200 IF( FATAL.AND.SFATAL )
344 $ GO TO 220
345 END IF
346 210 CONTINUE
347 WRITE( NOUT, FMT = 9982 )
348 GO TO 240
349 *
350 220 CONTINUE
351 WRITE( NOUT, FMT = 9981 )
352 GO TO 240
353 *
354 230 CONTINUE
355 WRITE( NOUT, FMT = 9987 )
356 *
357 240 CONTINUE
358 IF( TRACE )
359 $ CLOSE ( NTRA )
360 CLOSE ( NOUT )
361 STOP
362 *
363 9999 FORMAT( ' ROUTINES PASS COMPUTATIONAL TESTS IF TEST RATIO IS LES',
364 $ 'S THAN', F8.2 )
365 9998 FORMAT( ' RELATIVE MACHINE PRECISION IS TAKEN TO BE', 1P, D9.1 )
366 9997 FORMAT( ' NUMBER OF VALUES OF ', A, ' IS LESS THAN 1 OR GREATER ',
367 $ 'THAN ', I2 )
368 9996 FORMAT( ' VALUE OF N IS LESS THAN 0 OR GREATER THAN ', I2 )
369 9995 FORMAT( ' VALUE OF K IS LESS THAN 0' )
370 9994 FORMAT( ' ABSOLUTE VALUE OF INCX OR INCY IS 0 OR GREATER THAN ',
371 $ I2 )
372 9993 FORMAT( ' TESTS OF THE COMPLEX*16 LEVEL 2 BLAS', //' THE F',
373 $ 'OLLOWING PARAMETER VALUES WILL BE USED:' )
374 9992 FORMAT( ' FOR N ', 9I6 )
375 9991 FORMAT( ' FOR K ', 7I6 )
376 9990 FORMAT( ' FOR INCX AND INCY ', 7I6 )
377 9989 FORMAT( ' FOR ALPHA ',
378 $ 7( '(', F4.1, ',', F4.1, ') ', : ) )
379 9988 FORMAT( ' FOR BETA ',
380 $ 7( '(', F4.1, ',', F4.1, ') ', : ) )
381 9987 FORMAT( ' AMEND DATA FILE OR INCREASE ARRAY SIZES IN PROGRAM',
382 $ /' ******* TESTS ABANDONED *******' )
383 9986 FORMAT( ' SUBPROGRAM NAME ', A6, ' NOT RECOGNIZED', /' ******* T',
384 $ 'ESTS ABANDONED *******' )
385 9985 FORMAT( ' ERROR IN ZMVCH - IN-LINE DOT PRODUCTS ARE BEING EVALU',
386 $ 'ATED WRONGLY.', /' ZMVCH WAS CALLED WITH TRANS = ', A1,
387 $ ' AND RETURNED SAME = ', L1, ' AND ERR = ', F12.3, '.', /
388 $ ' THIS MAY BE DUE TO FAULTS IN THE ARITHMETIC OR THE COMPILER.'
389 $ , /' ******* TESTS ABANDONED *******' )
390 9984 FORMAT( A6, L2 )
391 9983 FORMAT( 1X, A6, ' WAS NOT TESTED' )
392 9982 FORMAT( /' END OF TESTS' )
393 9981 FORMAT( /' ******* FATAL ERROR - TESTS ABANDONED *******' )
394 9980 FORMAT( ' ERROR-EXITS WILL NOT BE TESTED' )
395 *
396 * End of ZBLAT2.
397 *
398 END
399 SUBROUTINE ZCHK1( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
400 $ FATAL, NIDIM, IDIM, NKB, KB, NALF, ALF, NBET,
401 $ BET, NINC, INC, NMAX, INCMAX, A, AA, AS, X, XX,
402 $ XS, Y, YY, YS, YT, G )
403 *
404 * Tests ZGEMV and ZGBMV.
405 *
406 * Auxiliary routine for test program for Level 2 Blas.
407 *
408 * -- Written on 10-August-1987.
409 * Richard Hanson, Sandia National Labs.
410 * Jeremy Du Croz, NAG Central Office.
411 *
412 * .. Parameters ..
413 COMPLEX*16 ZERO, HALF
414 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
415 $ HALF = ( 0.5D0, 0.0D0 ) )
416 DOUBLE PRECISION RZERO
417 PARAMETER ( RZERO = 0.0D0 )
418 * .. Scalar Arguments ..
419 DOUBLE PRECISION EPS, THRESH
420 INTEGER INCMAX, NALF, NBET, NIDIM, NINC, NKB, NMAX,
421 $ NOUT, NTRA
422 LOGICAL FATAL, REWI, TRACE
423 CHARACTER*6 SNAME
424 * .. Array Arguments ..
425 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ), ALF( NALF ),
426 $ AS( NMAX*NMAX ), BET( NBET ), X( NMAX ),
427 $ XS( NMAX*INCMAX ), XX( NMAX*INCMAX ),
428 $ Y( NMAX ), YS( NMAX*INCMAX ), YT( NMAX ),
429 $ YY( NMAX*INCMAX )
430 DOUBLE PRECISION G( NMAX )
431 INTEGER IDIM( NIDIM ), INC( NINC ), KB( NKB )
432 * .. Local Scalars ..
433 COMPLEX*16 ALPHA, ALS, BETA, BLS, TRANSL
434 DOUBLE PRECISION ERR, ERRMAX
435 INTEGER I, IA, IB, IC, IKU, IM, IN, INCX, INCXS, INCY,
436 $ INCYS, IX, IY, KL, KLS, KU, KUS, LAA, LDA,
437 $ LDAS, LX, LY, M, ML, MS, N, NARGS, NC, ND, NK,
438 $ NL, NS
439 LOGICAL BANDED, FULL, NULL, RESET, SAME, TRAN
440 CHARACTER*1 TRANS, TRANSS
441 CHARACTER*3 ICH
442 * .. Local Arrays ..
443 LOGICAL ISAME( 13 )
444 * .. External Functions ..
445 LOGICAL LZE, LZERES
446 EXTERNAL LZE, LZERES
447 * .. External Subroutines ..
448 EXTERNAL ZGBMV, ZGEMV, ZMAKE, ZMVCH
449 * .. Intrinsic Functions ..
450 INTRINSIC ABS, MAX, MIN
451 * .. Scalars in Common ..
452 INTEGER INFOT, NOUTC
453 LOGICAL LERR, OK
454 * .. Common blocks ..
455 COMMON /INFOC/INFOT, NOUTC, OK, LERR
456 * .. Data statements ..
457 DATA ICH/'NTC'/
458 * .. Executable Statements ..
459 FULL = SNAME( 3: 3 ).EQ.'E'
460 BANDED = SNAME( 3: 3 ).EQ.'B'
461 * Define the number of arguments.
462 IF( FULL )THEN
463 NARGS = 11
464 ELSE IF( BANDED )THEN
465 NARGS = 13
466 END IF
467 *
468 NC = 0
469 RESET = .TRUE.
470 ERRMAX = RZERO
471 *
472 DO 120 IN = 1, NIDIM
473 N = IDIM( IN )
474 ND = N/2 + 1
475 *
476 DO 110 IM = 1, 2
477 IF( IM.EQ.1 )
478 $ M = MAX( N - ND, 0 )
479 IF( IM.EQ.2 )
480 $ M = MIN( N + ND, NMAX )
481 *
482 IF( BANDED )THEN
483 NK = NKB
484 ELSE
485 NK = 1
486 END IF
487 DO 100 IKU = 1, NK
488 IF( BANDED )THEN
489 KU = KB( IKU )
490 KL = MAX( KU - 1, 0 )
491 ELSE
492 KU = N - 1
493 KL = M - 1
494 END IF
495 * Set LDA to 1 more than minimum value if room.
496 IF( BANDED )THEN
497 LDA = KL + KU + 1
498 ELSE
499 LDA = M
500 END IF
501 IF( LDA.LT.NMAX )
502 $ LDA = LDA + 1
503 * Skip tests if not enough room.
504 IF( LDA.GT.NMAX )
505 $ GO TO 100
506 LAA = LDA*N
507 NULL = N.LE.0.OR.M.LE.0
508 *
509 * Generate the matrix A.
510 *
511 TRANSL = ZERO
512 CALL ZMAKE( SNAME( 2: 3 ), ' ', ' ', M, N, A, NMAX, AA,
513 $ LDA, KL, KU, RESET, TRANSL )
514 *
515 DO 90 IC = 1, 3
516 TRANS = ICH( IC: IC )
517 TRAN = TRANS.EQ.'T'.OR.TRANS.EQ.'C'
518 *
519 IF( TRAN )THEN
520 ML = N
521 NL = M
522 ELSE
523 ML = M
524 NL = N
525 END IF
526 *
527 DO 80 IX = 1, NINC
528 INCX = INC( IX )
529 LX = ABS( INCX )*NL
530 *
531 * Generate the vector X.
532 *
533 TRANSL = HALF
534 CALL ZMAKE( 'GE', ' ', ' ', 1, NL, X, 1, XX,
535 $ ABS( INCX ), 0, NL - 1, RESET, TRANSL )
536 IF( NL.GT.1 )THEN
537 X( NL/2 ) = ZERO
538 XX( 1 + ABS( INCX )*( NL/2 - 1 ) ) = ZERO
539 END IF
540 *
541 DO 70 IY = 1, NINC
542 INCY = INC( IY )
543 LY = ABS( INCY )*ML
544 *
545 DO 60 IA = 1, NALF
546 ALPHA = ALF( IA )
547 *
548 DO 50 IB = 1, NBET
549 BETA = BET( IB )
550 *
551 * Generate the vector Y.
552 *
553 TRANSL = ZERO
554 CALL ZMAKE( 'GE', ' ', ' ', 1, ML, Y, 1,
555 $ YY, ABS( INCY ), 0, ML - 1,
556 $ RESET, TRANSL )
557 *
558 NC = NC + 1
559 *
560 * Save every datum before calling the
561 * subroutine.
562 *
563 TRANSS = TRANS
564 MS = M
565 NS = N
566 KLS = KL
567 KUS = KU
568 ALS = ALPHA
569 DO 10 I = 1, LAA
570 AS( I ) = AA( I )
571 10 CONTINUE
572 LDAS = LDA
573 DO 20 I = 1, LX
574 XS( I ) = XX( I )
575 20 CONTINUE
576 INCXS = INCX
577 BLS = BETA
578 DO 30 I = 1, LY
579 YS( I ) = YY( I )
580 30 CONTINUE
581 INCYS = INCY
582 *
583 * Call the subroutine.
584 *
585 IF( FULL )THEN
586 IF( TRACE )
587 $ WRITE( NTRA, FMT = 9994 )NC, SNAME,
588 $ TRANS, M, N, ALPHA, LDA, INCX, BETA,
589 $ INCY
590 IF( REWI )
591 $ REWIND NTRA
592 CALL ZGEMV( TRANS, M, N, ALPHA, AA,
593 $ LDA, XX, INCX, BETA, YY,
594 $ INCY )
595 ELSE IF( BANDED )THEN
596 IF( TRACE )
597 $ WRITE( NTRA, FMT = 9995 )NC, SNAME,
598 $ TRANS, M, N, KL, KU, ALPHA, LDA,
599 $ INCX, BETA, INCY
600 IF( REWI )
601 $ REWIND NTRA
602 CALL ZGBMV( TRANS, M, N, KL, KU, ALPHA,
603 $ AA, LDA, XX, INCX, BETA,
604 $ YY, INCY )
605 END IF
606 *
607 * Check if error-exit was taken incorrectly.
608 *
609 IF( .NOT.OK )THEN
610 WRITE( NOUT, FMT = 9993 )
611 FATAL = .TRUE.
612 GO TO 130
613 END IF
614 *
615 * See what data changed inside subroutines.
616 *
617 ISAME( 1 ) = TRANS.EQ.TRANSS
618 ISAME( 2 ) = MS.EQ.M
619 ISAME( 3 ) = NS.EQ.N
620 IF( FULL )THEN
621 ISAME( 4 ) = ALS.EQ.ALPHA
622 ISAME( 5 ) = LZE( AS, AA, LAA )
623 ISAME( 6 ) = LDAS.EQ.LDA
624 ISAME( 7 ) = LZE( XS, XX, LX )
625 ISAME( 8 ) = INCXS.EQ.INCX
626 ISAME( 9 ) = BLS.EQ.BETA
627 IF( NULL )THEN
628 ISAME( 10 ) = LZE( YS, YY, LY )
629 ELSE
630 ISAME( 10 ) = LZERES( 'GE', ' ', 1,
631 $ ML, YS, YY,
632 $ ABS( INCY ) )
633 END IF
634 ISAME( 11 ) = INCYS.EQ.INCY
635 ELSE IF( BANDED )THEN
636 ISAME( 4 ) = KLS.EQ.KL
637 ISAME( 5 ) = KUS.EQ.KU
638 ISAME( 6 ) = ALS.EQ.ALPHA
639 ISAME( 7 ) = LZE( AS, AA, LAA )
640 ISAME( 8 ) = LDAS.EQ.LDA
641 ISAME( 9 ) = LZE( XS, XX, LX )
642 ISAME( 10 ) = INCXS.EQ.INCX
643 ISAME( 11 ) = BLS.EQ.BETA
644 IF( NULL )THEN
645 ISAME( 12 ) = LZE( YS, YY, LY )
646 ELSE
647 ISAME( 12 ) = LZERES( 'GE', ' ', 1,
648 $ ML, YS, YY,
649 $ ABS( INCY ) )
650 END IF
651 ISAME( 13 ) = INCYS.EQ.INCY
652 END IF
653 *
654 * If data was incorrectly changed, report
655 * and return.
656 *
657 SAME = .TRUE.
658 DO 40 I = 1, NARGS
659 SAME = SAME.AND.ISAME( I )
660 IF( .NOT.ISAME( I ) )
661 $ WRITE( NOUT, FMT = 9998 )I
662 40 CONTINUE
663 IF( .NOT.SAME )THEN
664 FATAL = .TRUE.
665 GO TO 130
666 END IF
667 *
668 IF( .NOT.NULL )THEN
669 *
670 * Check the result.
671 *
672 CALL ZMVCH( TRANS, M, N, ALPHA, A,
673 $ NMAX, X, INCX, BETA, Y,
674 $ INCY, YT, G, YY, EPS, ERR,
675 $ FATAL, NOUT, .TRUE. )
676 ERRMAX = MAX( ERRMAX, ERR )
677 * If got really bad answer, report and
678 * return.
679 IF( FATAL )
680 $ GO TO 130
681 ELSE
682 * Avoid repeating tests with M.le.0 or
683 * N.le.0.
684 GO TO 110
685 END IF
686 *
687 50 CONTINUE
688 *
689 60 CONTINUE
690 *
691 70 CONTINUE
692 *
693 80 CONTINUE
694 *
695 90 CONTINUE
696 *
697 100 CONTINUE
698 *
699 110 CONTINUE
700 *
701 120 CONTINUE
702 *
703 * Report result.
704 *
705 IF( ERRMAX.LT.THRESH )THEN
706 WRITE( NOUT, FMT = 9999 )SNAME, NC
707 ELSE
708 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
709 END IF
710 GO TO 140
711 *
712 130 CONTINUE
713 WRITE( NOUT, FMT = 9996 )SNAME
714 IF( FULL )THEN
715 WRITE( NOUT, FMT = 9994 )NC, SNAME, TRANS, M, N, ALPHA, LDA,
716 $ INCX, BETA, INCY
717 ELSE IF( BANDED )THEN
718 WRITE( NOUT, FMT = 9995 )NC, SNAME, TRANS, M, N, KL, KU,
719 $ ALPHA, LDA, INCX, BETA, INCY
720 END IF
721 *
722 140 CONTINUE
723 RETURN
724 *
725 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
726 $ 'S)' )
727 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
728 $ 'ANGED INCORRECTLY *******' )
729 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
730 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
731 $ ' - SUSPECT *******' )
732 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
733 9995 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 4( I3, ',' ), '(',
734 $ F4.1, ',', F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',',
735 $ F4.1, '), Y,', I2, ') .' )
736 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 2( I3, ',' ), '(',
737 $ F4.1, ',', F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',',
738 $ F4.1, '), Y,', I2, ') .' )
739 9993 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
740 $ '******' )
741 *
742 * End of ZCHK1.
743 *
744 END
745 SUBROUTINE ZCHK2( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
746 $ FATAL, NIDIM, IDIM, NKB, KB, NALF, ALF, NBET,
747 $ BET, NINC, INC, NMAX, INCMAX, A, AA, AS, X, XX,
748 $ XS, Y, YY, YS, YT, G )
749 *
750 * Tests ZHEMV, ZHBMV and ZHPMV.
751 *
752 * Auxiliary routine for test program for Level 2 Blas.
753 *
754 * -- Written on 10-August-1987.
755 * Richard Hanson, Sandia National Labs.
756 * Jeremy Du Croz, NAG Central Office.
757 *
758 * .. Parameters ..
759 COMPLEX*16 ZERO, HALF
760 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
761 $ HALF = ( 0.5D0, 0.0D0 ) )
762 DOUBLE PRECISION RZERO
763 PARAMETER ( RZERO = 0.0D0 )
764 * .. Scalar Arguments ..
765 DOUBLE PRECISION EPS, THRESH
766 INTEGER INCMAX, NALF, NBET, NIDIM, NINC, NKB, NMAX,
767 $ NOUT, NTRA
768 LOGICAL FATAL, REWI, TRACE
769 CHARACTER*6 SNAME
770 * .. Array Arguments ..
771 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ), ALF( NALF ),
772 $ AS( NMAX*NMAX ), BET( NBET ), X( NMAX ),
773 $ XS( NMAX*INCMAX ), XX( NMAX*INCMAX ),
774 $ Y( NMAX ), YS( NMAX*INCMAX ), YT( NMAX ),
775 $ YY( NMAX*INCMAX )
776 DOUBLE PRECISION G( NMAX )
777 INTEGER IDIM( NIDIM ), INC( NINC ), KB( NKB )
778 * .. Local Scalars ..
779 COMPLEX*16 ALPHA, ALS, BETA, BLS, TRANSL
780 DOUBLE PRECISION ERR, ERRMAX
781 INTEGER I, IA, IB, IC, IK, IN, INCX, INCXS, INCY,
782 $ INCYS, IX, IY, K, KS, LAA, LDA, LDAS, LX, LY,
783 $ N, NARGS, NC, NK, NS
784 LOGICAL BANDED, FULL, NULL, PACKED, RESET, SAME
785 CHARACTER*1 UPLO, UPLOS
786 CHARACTER*2 ICH
787 * .. Local Arrays ..
788 LOGICAL ISAME( 13 )
789 * .. External Functions ..
790 LOGICAL LZE, LZERES
791 EXTERNAL LZE, LZERES
792 * .. External Subroutines ..
793 EXTERNAL ZHBMV, ZHEMV, ZHPMV, ZMAKE, ZMVCH
794 * .. Intrinsic Functions ..
795 INTRINSIC ABS, MAX
796 * .. Scalars in Common ..
797 INTEGER INFOT, NOUTC
798 LOGICAL LERR, OK
799 * .. Common blocks ..
800 COMMON /INFOC/INFOT, NOUTC, OK, LERR
801 * .. Data statements ..
802 DATA ICH/'UL'/
803 * .. Executable Statements ..
804 FULL = SNAME( 3: 3 ).EQ.'E'
805 BANDED = SNAME( 3: 3 ).EQ.'B'
806 PACKED = SNAME( 3: 3 ).EQ.'P'
807 * Define the number of arguments.
808 IF( FULL )THEN
809 NARGS = 10
810 ELSE IF( BANDED )THEN
811 NARGS = 11
812 ELSE IF( PACKED )THEN
813 NARGS = 9
814 END IF
815 *
816 NC = 0
817 RESET = .TRUE.
818 ERRMAX = RZERO
819 *
820 DO 110 IN = 1, NIDIM
821 N = IDIM( IN )
822 *
823 IF( BANDED )THEN
824 NK = NKB
825 ELSE
826 NK = 1
827 END IF
828 DO 100 IK = 1, NK
829 IF( BANDED )THEN
830 K = KB( IK )
831 ELSE
832 K = N - 1
833 END IF
834 * Set LDA to 1 more than minimum value if room.
835 IF( BANDED )THEN
836 LDA = K + 1
837 ELSE
838 LDA = N
839 END IF
840 IF( LDA.LT.NMAX )
841 $ LDA = LDA + 1
842 * Skip tests if not enough room.
843 IF( LDA.GT.NMAX )
844 $ GO TO 100
845 IF( PACKED )THEN
846 LAA = ( N*( N + 1 ) )/2
847 ELSE
848 LAA = LDA*N
849 END IF
850 NULL = N.LE.0
851 *
852 DO 90 IC = 1, 2
853 UPLO = ICH( IC: IC )
854 *
855 * Generate the matrix A.
856 *
857 TRANSL = ZERO
858 CALL ZMAKE( SNAME( 2: 3 ), UPLO, ' ', N, N, A, NMAX, AA,
859 $ LDA, K, K, RESET, TRANSL )
860 *
861 DO 80 IX = 1, NINC
862 INCX = INC( IX )
863 LX = ABS( INCX )*N
864 *
865 * Generate the vector X.
866 *
867 TRANSL = HALF
868 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX,
869 $ ABS( INCX ), 0, N - 1, RESET, TRANSL )
870 IF( N.GT.1 )THEN
871 X( N/2 ) = ZERO
872 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
873 END IF
874 *
875 DO 70 IY = 1, NINC
876 INCY = INC( IY )
877 LY = ABS( INCY )*N
878 *
879 DO 60 IA = 1, NALF
880 ALPHA = ALF( IA )
881 *
882 DO 50 IB = 1, NBET
883 BETA = BET( IB )
884 *
885 * Generate the vector Y.
886 *
887 TRANSL = ZERO
888 CALL ZMAKE( 'GE', ' ', ' ', 1, N, Y, 1, YY,
889 $ ABS( INCY ), 0, N - 1, RESET,
890 $ TRANSL )
891 *
892 NC = NC + 1
893 *
894 * Save every datum before calling the
895 * subroutine.
896 *
897 UPLOS = UPLO
898 NS = N
899 KS = K
900 ALS = ALPHA
901 DO 10 I = 1, LAA
902 AS( I ) = AA( I )
903 10 CONTINUE
904 LDAS = LDA
905 DO 20 I = 1, LX
906 XS( I ) = XX( I )
907 20 CONTINUE
908 INCXS = INCX
909 BLS = BETA
910 DO 30 I = 1, LY
911 YS( I ) = YY( I )
912 30 CONTINUE
913 INCYS = INCY
914 *
915 * Call the subroutine.
916 *
917 IF( FULL )THEN
918 IF( TRACE )
919 $ WRITE( NTRA, FMT = 9993 )NC, SNAME,
920 $ UPLO, N, ALPHA, LDA, INCX, BETA, INCY
921 IF( REWI )
922 $ REWIND NTRA
923 CALL ZHEMV( UPLO, N, ALPHA, AA, LDA, XX,
924 $ INCX, BETA, YY, INCY )
925 ELSE IF( BANDED )THEN
926 IF( TRACE )
927 $ WRITE( NTRA, FMT = 9994 )NC, SNAME,
928 $ UPLO, N, K, ALPHA, LDA, INCX, BETA,
929 $ INCY
930 IF( REWI )
931 $ REWIND NTRA
932 CALL ZHBMV( UPLO, N, K, ALPHA, AA, LDA,
933 $ XX, INCX, BETA, YY, INCY )
934 ELSE IF( PACKED )THEN
935 IF( TRACE )
936 $ WRITE( NTRA, FMT = 9995 )NC, SNAME,
937 $ UPLO, N, ALPHA, INCX, BETA, INCY
938 IF( REWI )
939 $ REWIND NTRA
940 CALL ZHPMV( UPLO, N, ALPHA, AA, XX, INCX,
941 $ BETA, YY, INCY )
942 END IF
943 *
944 * Check if error-exit was taken incorrectly.
945 *
946 IF( .NOT.OK )THEN
947 WRITE( NOUT, FMT = 9992 )
948 FATAL = .TRUE.
949 GO TO 120
950 END IF
951 *
952 * See what data changed inside subroutines.
953 *
954 ISAME( 1 ) = UPLO.EQ.UPLOS
955 ISAME( 2 ) = NS.EQ.N
956 IF( FULL )THEN
957 ISAME( 3 ) = ALS.EQ.ALPHA
958 ISAME( 4 ) = LZE( AS, AA, LAA )
959 ISAME( 5 ) = LDAS.EQ.LDA
960 ISAME( 6 ) = LZE( XS, XX, LX )
961 ISAME( 7 ) = INCXS.EQ.INCX
962 ISAME( 8 ) = BLS.EQ.BETA
963 IF( NULL )THEN
964 ISAME( 9 ) = LZE( YS, YY, LY )
965 ELSE
966 ISAME( 9 ) = LZERES( 'GE', ' ', 1, N,
967 $ YS, YY, ABS( INCY ) )
968 END IF
969 ISAME( 10 ) = INCYS.EQ.INCY
970 ELSE IF( BANDED )THEN
971 ISAME( 3 ) = KS.EQ.K
972 ISAME( 4 ) = ALS.EQ.ALPHA
973 ISAME( 5 ) = LZE( AS, AA, LAA )
974 ISAME( 6 ) = LDAS.EQ.LDA
975 ISAME( 7 ) = LZE( XS, XX, LX )
976 ISAME( 8 ) = INCXS.EQ.INCX
977 ISAME( 9 ) = BLS.EQ.BETA
978 IF( NULL )THEN
979 ISAME( 10 ) = LZE( YS, YY, LY )
980 ELSE
981 ISAME( 10 ) = LZERES( 'GE', ' ', 1, N,
982 $ YS, YY, ABS( INCY ) )
983 END IF
984 ISAME( 11 ) = INCYS.EQ.INCY
985 ELSE IF( PACKED )THEN
986 ISAME( 3 ) = ALS.EQ.ALPHA
987 ISAME( 4 ) = LZE( AS, AA, LAA )
988 ISAME( 5 ) = LZE( XS, XX, LX )
989 ISAME( 6 ) = INCXS.EQ.INCX
990 ISAME( 7 ) = BLS.EQ.BETA
991 IF( NULL )THEN
992 ISAME( 8 ) = LZE( YS, YY, LY )
993 ELSE
994 ISAME( 8 ) = LZERES( 'GE', ' ', 1, N,
995 $ YS, YY, ABS( INCY ) )
996 END IF
997 ISAME( 9 ) = INCYS.EQ.INCY
998 END IF
999 *
1000 * If data was incorrectly changed, report and
1001 * return.
1002 *
1003 SAME = .TRUE.
1004 DO 40 I = 1, NARGS
1005 SAME = SAME.AND.ISAME( I )
1006 IF( .NOT.ISAME( I ) )
1007 $ WRITE( NOUT, FMT = 9998 )I
1008 40 CONTINUE
1009 IF( .NOT.SAME )THEN
1010 FATAL = .TRUE.
1011 GO TO 120
1012 END IF
1013 *
1014 IF( .NOT.NULL )THEN
1015 *
1016 * Check the result.
1017 *
1018 CALL ZMVCH( 'N', N, N, ALPHA, A, NMAX, X,
1019 $ INCX, BETA, Y, INCY, YT, G,
1020 $ YY, EPS, ERR, FATAL, NOUT,
1021 $ .TRUE. )
1022 ERRMAX = MAX( ERRMAX, ERR )
1023 * If got really bad answer, report and
1024 * return.
1025 IF( FATAL )
1026 $ GO TO 120
1027 ELSE
1028 * Avoid repeating tests with N.le.0
1029 GO TO 110
1030 END IF
1031 *
1032 50 CONTINUE
1033 *
1034 60 CONTINUE
1035 *
1036 70 CONTINUE
1037 *
1038 80 CONTINUE
1039 *
1040 90 CONTINUE
1041 *
1042 100 CONTINUE
1043 *
1044 110 CONTINUE
1045 *
1046 * Report result.
1047 *
1048 IF( ERRMAX.LT.THRESH )THEN
1049 WRITE( NOUT, FMT = 9999 )SNAME, NC
1050 ELSE
1051 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
1052 END IF
1053 GO TO 130
1054 *
1055 120 CONTINUE
1056 WRITE( NOUT, FMT = 9996 )SNAME
1057 IF( FULL )THEN
1058 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, N, ALPHA, LDA, INCX,
1059 $ BETA, INCY
1060 ELSE IF( BANDED )THEN
1061 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, N, K, ALPHA, LDA,
1062 $ INCX, BETA, INCY
1063 ELSE IF( PACKED )THEN
1064 WRITE( NOUT, FMT = 9995 )NC, SNAME, UPLO, N, ALPHA, INCX,
1065 $ BETA, INCY
1066 END IF
1067 *
1068 130 CONTINUE
1069 RETURN
1070 *
1071 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
1072 $ 'S)' )
1073 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1074 $ 'ANGED INCORRECTLY *******' )
1075 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
1076 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
1077 $ ' - SUSPECT *******' )
1078 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
1079 9995 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
1080 $ F4.1, '), AP, X,', I2, ',(', F4.1, ',', F4.1, '), Y,', I2,
1081 $ ') .' )
1082 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 2( I3, ',' ), '(',
1083 $ F4.1, ',', F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',',
1084 $ F4.1, '), Y,', I2, ') .' )
1085 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
1086 $ F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',', F4.1, '), ',
1087 $ 'Y,', I2, ') .' )
1088 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1089 $ '******' )
1090 *
1091 * End of ZCHK2.
1092 *
1093 END
1094 SUBROUTINE ZCHK3( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
1095 $ FATAL, NIDIM, IDIM, NKB, KB, NINC, INC, NMAX,
1096 $ INCMAX, A, AA, AS, X, XX, XS, XT, G, Z )
1097 *
1098 * Tests ZTRMV, ZTBMV, ZTPMV, ZTRSV, ZTBSV and ZTPSV.
1099 *
1100 * Auxiliary routine for test program for Level 2 Blas.
1101 *
1102 * -- Written on 10-August-1987.
1103 * Richard Hanson, Sandia National Labs.
1104 * Jeremy Du Croz, NAG Central Office.
1105 *
1106 * .. Parameters ..
1107 COMPLEX*16 ZERO, HALF, ONE
1108 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
1109 $ HALF = ( 0.5D0, 0.0D0 ),
1110 $ ONE = ( 1.0D0, 0.0D0 ) )
1111 DOUBLE PRECISION RZERO
1112 PARAMETER ( RZERO = 0.0D0 )
1113 * .. Scalar Arguments ..
1114 DOUBLE PRECISION EPS, THRESH
1115 INTEGER INCMAX, NIDIM, NINC, NKB, NMAX, NOUT, NTRA
1116 LOGICAL FATAL, REWI, TRACE
1117 CHARACTER*6 SNAME
1118 * .. Array Arguments ..
1119 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ),
1120 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ),
1121 $ XT( NMAX ), XX( NMAX*INCMAX ), Z( NMAX )
1122 DOUBLE PRECISION G( NMAX )
1123 INTEGER IDIM( NIDIM ), INC( NINC ), KB( NKB )
1124 * .. Local Scalars ..
1125 COMPLEX*16 TRANSL
1126 DOUBLE PRECISION ERR, ERRMAX
1127 INTEGER I, ICD, ICT, ICU, IK, IN, INCX, INCXS, IX, K,
1128 $ KS, LAA, LDA, LDAS, LX, N, NARGS, NC, NK, NS
1129 LOGICAL BANDED, FULL, NULL, PACKED, RESET, SAME
1130 CHARACTER*1 DIAG, DIAGS, TRANS, TRANSS, UPLO, UPLOS
1131 CHARACTER*2 ICHD, ICHU
1132 CHARACTER*3 ICHT
1133 * .. Local Arrays ..
1134 LOGICAL ISAME( 13 )
1135 * .. External Functions ..
1136 LOGICAL LZE, LZERES
1137 EXTERNAL LZE, LZERES
1138 * .. External Subroutines ..
1139 EXTERNAL ZMAKE, ZMVCH, ZTBMV, ZTBSV, ZTPMV, ZTPSV,
1140 $ ZTRMV, ZTRSV
1141 * .. Intrinsic Functions ..
1142 INTRINSIC ABS, MAX
1143 * .. Scalars in Common ..
1144 INTEGER INFOT, NOUTC
1145 LOGICAL LERR, OK
1146 * .. Common blocks ..
1147 COMMON /INFOC/INFOT, NOUTC, OK, LERR
1148 * .. Data statements ..
1149 DATA ICHU/'UL'/, ICHT/'NTC'/, ICHD/'UN'/
1150 * .. Executable Statements ..
1151 FULL = SNAME( 3: 3 ).EQ.'R'
1152 BANDED = SNAME( 3: 3 ).EQ.'B'
1153 PACKED = SNAME( 3: 3 ).EQ.'P'
1154 * Define the number of arguments.
1155 IF( FULL )THEN
1156 NARGS = 8
1157 ELSE IF( BANDED )THEN
1158 NARGS = 9
1159 ELSE IF( PACKED )THEN
1160 NARGS = 7
1161 END IF
1162 *
1163 NC = 0
1164 RESET = .TRUE.
1165 ERRMAX = RZERO
1166 * Set up zero vector for ZMVCH.
1167 DO 10 I = 1, NMAX
1168 Z( I ) = ZERO
1169 10 CONTINUE
1170 *
1171 DO 110 IN = 1, NIDIM
1172 N = IDIM( IN )
1173 *
1174 IF( BANDED )THEN
1175 NK = NKB
1176 ELSE
1177 NK = 1
1178 END IF
1179 DO 100 IK = 1, NK
1180 IF( BANDED )THEN
1181 K = KB( IK )
1182 ELSE
1183 K = N - 1
1184 END IF
1185 * Set LDA to 1 more than minimum value if room.
1186 IF( BANDED )THEN
1187 LDA = K + 1
1188 ELSE
1189 LDA = N
1190 END IF
1191 IF( LDA.LT.NMAX )
1192 $ LDA = LDA + 1
1193 * Skip tests if not enough room.
1194 IF( LDA.GT.NMAX )
1195 $ GO TO 100
1196 IF( PACKED )THEN
1197 LAA = ( N*( N + 1 ) )/2
1198 ELSE
1199 LAA = LDA*N
1200 END IF
1201 NULL = N.LE.0
1202 *
1203 DO 90 ICU = 1, 2
1204 UPLO = ICHU( ICU: ICU )
1205 *
1206 DO 80 ICT = 1, 3
1207 TRANS = ICHT( ICT: ICT )
1208 *
1209 DO 70 ICD = 1, 2
1210 DIAG = ICHD( ICD: ICD )
1211 *
1212 * Generate the matrix A.
1213 *
1214 TRANSL = ZERO
1215 CALL ZMAKE( SNAME( 2: 3 ), UPLO, DIAG, N, N, A,
1216 $ NMAX, AA, LDA, K, K, RESET, TRANSL )
1217 *
1218 DO 60 IX = 1, NINC
1219 INCX = INC( IX )
1220 LX = ABS( INCX )*N
1221 *
1222 * Generate the vector X.
1223 *
1224 TRANSL = HALF
1225 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX,
1226 $ ABS( INCX ), 0, N - 1, RESET,
1227 $ TRANSL )
1228 IF( N.GT.1 )THEN
1229 X( N/2 ) = ZERO
1230 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
1231 END IF
1232 *
1233 NC = NC + 1
1234 *
1235 * Save every datum before calling the subroutine.
1236 *
1237 UPLOS = UPLO
1238 TRANSS = TRANS
1239 DIAGS = DIAG
1240 NS = N
1241 KS = K
1242 DO 20 I = 1, LAA
1243 AS( I ) = AA( I )
1244 20 CONTINUE
1245 LDAS = LDA
1246 DO 30 I = 1, LX
1247 XS( I ) = XX( I )
1248 30 CONTINUE
1249 INCXS = INCX
1250 *
1251 * Call the subroutine.
1252 *
1253 IF( SNAME( 4: 5 ).EQ.'MV' )THEN
1254 IF( FULL )THEN
1255 IF( TRACE )
1256 $ WRITE( NTRA, FMT = 9993 )NC, SNAME,
1257 $ UPLO, TRANS, DIAG, N, LDA, INCX
1258 IF( REWI )
1259 $ REWIND NTRA
1260 CALL ZTRMV( UPLO, TRANS, DIAG, N, AA, LDA,
1261 $ XX, INCX )
1262 ELSE IF( BANDED )THEN
1263 IF( TRACE )
1264 $ WRITE( NTRA, FMT = 9994 )NC, SNAME,
1265 $ UPLO, TRANS, DIAG, N, K, LDA, INCX
1266 IF( REWI )
1267 $ REWIND NTRA
1268 CALL ZTBMV( UPLO, TRANS, DIAG, N, K, AA,
1269 $ LDA, XX, INCX )
1270 ELSE IF( PACKED )THEN
1271 IF( TRACE )
1272 $ WRITE( NTRA, FMT = 9995 )NC, SNAME,
1273 $ UPLO, TRANS, DIAG, N, INCX
1274 IF( REWI )
1275 $ REWIND NTRA
1276 CALL ZTPMV( UPLO, TRANS, DIAG, N, AA, XX,
1277 $ INCX )
1278 END IF
1279 ELSE IF( SNAME( 4: 5 ).EQ.'SV' )THEN
1280 IF( FULL )THEN
1281 IF( TRACE )
1282 $ WRITE( NTRA, FMT = 9993 )NC, SNAME,
1283 $ UPLO, TRANS, DIAG, N, LDA, INCX
1284 IF( REWI )
1285 $ REWIND NTRA
1286 CALL ZTRSV( UPLO, TRANS, DIAG, N, AA, LDA,
1287 $ XX, INCX )
1288 ELSE IF( BANDED )THEN
1289 IF( TRACE )
1290 $ WRITE( NTRA, FMT = 9994 )NC, SNAME,
1291 $ UPLO, TRANS, DIAG, N, K, LDA, INCX
1292 IF( REWI )
1293 $ REWIND NTRA
1294 CALL ZTBSV( UPLO, TRANS, DIAG, N, K, AA,
1295 $ LDA, XX, INCX )
1296 ELSE IF( PACKED )THEN
1297 IF( TRACE )
1298 $ WRITE( NTRA, FMT = 9995 )NC, SNAME,
1299 $ UPLO, TRANS, DIAG, N, INCX
1300 IF( REWI )
1301 $ REWIND NTRA
1302 CALL ZTPSV( UPLO, TRANS, DIAG, N, AA, XX,
1303 $ INCX )
1304 END IF
1305 END IF
1306 *
1307 * Check if error-exit was taken incorrectly.
1308 *
1309 IF( .NOT.OK )THEN
1310 WRITE( NOUT, FMT = 9992 )
1311 FATAL = .TRUE.
1312 GO TO 120
1313 END IF
1314 *
1315 * See what data changed inside subroutines.
1316 *
1317 ISAME( 1 ) = UPLO.EQ.UPLOS
1318 ISAME( 2 ) = TRANS.EQ.TRANSS
1319 ISAME( 3 ) = DIAG.EQ.DIAGS
1320 ISAME( 4 ) = NS.EQ.N
1321 IF( FULL )THEN
1322 ISAME( 5 ) = LZE( AS, AA, LAA )
1323 ISAME( 6 ) = LDAS.EQ.LDA
1324 IF( NULL )THEN
1325 ISAME( 7 ) = LZE( XS, XX, LX )
1326 ELSE
1327 ISAME( 7 ) = LZERES( 'GE', ' ', 1, N, XS,
1328 $ XX, ABS( INCX ) )
1329 END IF
1330 ISAME( 8 ) = INCXS.EQ.INCX
1331 ELSE IF( BANDED )THEN
1332 ISAME( 5 ) = KS.EQ.K
1333 ISAME( 6 ) = LZE( AS, AA, LAA )
1334 ISAME( 7 ) = LDAS.EQ.LDA
1335 IF( NULL )THEN
1336 ISAME( 8 ) = LZE( XS, XX, LX )
1337 ELSE
1338 ISAME( 8 ) = LZERES( 'GE', ' ', 1, N, XS,
1339 $ XX, ABS( INCX ) )
1340 END IF
1341 ISAME( 9 ) = INCXS.EQ.INCX
1342 ELSE IF( PACKED )THEN
1343 ISAME( 5 ) = LZE( AS, AA, LAA )
1344 IF( NULL )THEN
1345 ISAME( 6 ) = LZE( XS, XX, LX )
1346 ELSE
1347 ISAME( 6 ) = LZERES( 'GE', ' ', 1, N, XS,
1348 $ XX, ABS( INCX ) )
1349 END IF
1350 ISAME( 7 ) = INCXS.EQ.INCX
1351 END IF
1352 *
1353 * If data was incorrectly changed, report and
1354 * return.
1355 *
1356 SAME = .TRUE.
1357 DO 40 I = 1, NARGS
1358 SAME = SAME.AND.ISAME( I )
1359 IF( .NOT.ISAME( I ) )
1360 $ WRITE( NOUT, FMT = 9998 )I
1361 40 CONTINUE
1362 IF( .NOT.SAME )THEN
1363 FATAL = .TRUE.
1364 GO TO 120
1365 END IF
1366 *
1367 IF( .NOT.NULL )THEN
1368 IF( SNAME( 4: 5 ).EQ.'MV' )THEN
1369 *
1370 * Check the result.
1371 *
1372 CALL ZMVCH( TRANS, N, N, ONE, A, NMAX, X,
1373 $ INCX, ZERO, Z, INCX, XT, G,
1374 $ XX, EPS, ERR, FATAL, NOUT,
1375 $ .TRUE. )
1376 ELSE IF( SNAME( 4: 5 ).EQ.'SV' )THEN
1377 *
1378 * Compute approximation to original vector.
1379 *
1380 DO 50 I = 1, N
1381 Z( I ) = XX( 1 + ( I - 1 )*
1382 $ ABS( INCX ) )
1383 XX( 1 + ( I - 1 )*ABS( INCX ) )
1384 $ = X( I )
1385 50 CONTINUE
1386 CALL ZMVCH( TRANS, N, N, ONE, A, NMAX, Z,
1387 $ INCX, ZERO, X, INCX, XT, G,
1388 $ XX, EPS, ERR, FATAL, NOUT,
1389 $ .FALSE. )
1390 END IF
1391 ERRMAX = MAX( ERRMAX, ERR )
1392 * If got really bad answer, report and return.
1393 IF( FATAL )
1394 $ GO TO 120
1395 ELSE
1396 * Avoid repeating tests with N.le.0.
1397 GO TO 110
1398 END IF
1399 *
1400 60 CONTINUE
1401 *
1402 70 CONTINUE
1403 *
1404 80 CONTINUE
1405 *
1406 90 CONTINUE
1407 *
1408 100 CONTINUE
1409 *
1410 110 CONTINUE
1411 *
1412 * Report result.
1413 *
1414 IF( ERRMAX.LT.THRESH )THEN
1415 WRITE( NOUT, FMT = 9999 )SNAME, NC
1416 ELSE
1417 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
1418 END IF
1419 GO TO 130
1420 *
1421 120 CONTINUE
1422 WRITE( NOUT, FMT = 9996 )SNAME
1423 IF( FULL )THEN
1424 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, TRANS, DIAG, N, LDA,
1425 $ INCX
1426 ELSE IF( BANDED )THEN
1427 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, TRANS, DIAG, N, K,
1428 $ LDA, INCX
1429 ELSE IF( PACKED )THEN
1430 WRITE( NOUT, FMT = 9995 )NC, SNAME, UPLO, TRANS, DIAG, N, INCX
1431 END IF
1432 *
1433 130 CONTINUE
1434 RETURN
1435 *
1436 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
1437 $ 'S)' )
1438 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1439 $ 'ANGED INCORRECTLY *******' )
1440 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
1441 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
1442 $ ' - SUSPECT *******' )
1443 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
1444 9995 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), I3, ', AP, ',
1445 $ 'X,', I2, ') .' )
1446 9994 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), 2( I3, ',' ),
1447 $ ' A,', I3, ', X,', I2, ') .' )
1448 9993 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), I3, ', A,',
1449 $ I3, ', X,', I2, ') .' )
1450 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1451 $ '******' )
1452 *
1453 * End of ZCHK3.
1454 *
1455 END
1456 SUBROUTINE ZCHK4( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
1457 $ FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC, NMAX,
1458 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G,
1459 $ Z )
1460 *
1461 * Tests ZGERC and ZGERU.
1462 *
1463 * Auxiliary routine for test program for Level 2 Blas.
1464 *
1465 * -- Written on 10-August-1987.
1466 * Richard Hanson, Sandia National Labs.
1467 * Jeremy Du Croz, NAG Central Office.
1468 *
1469 * .. Parameters ..
1470 COMPLEX*16 ZERO, HALF, ONE
1471 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
1472 $ HALF = ( 0.5D0, 0.0D0 ),
1473 $ ONE = ( 1.0D0, 0.0D0 ) )
1474 DOUBLE PRECISION RZERO
1475 PARAMETER ( RZERO = 0.0D0 )
1476 * .. Scalar Arguments ..
1477 DOUBLE PRECISION EPS, THRESH
1478 INTEGER INCMAX, NALF, NIDIM, NINC, NMAX, NOUT, NTRA
1479 LOGICAL FATAL, REWI, TRACE
1480 CHARACTER*6 SNAME
1481 * .. Array Arguments ..
1482 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ), ALF( NALF ),
1483 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ),
1484 $ XX( NMAX*INCMAX ), Y( NMAX ),
1485 $ YS( NMAX*INCMAX ), YT( NMAX ),
1486 $ YY( NMAX*INCMAX ), Z( NMAX )
1487 DOUBLE PRECISION G( NMAX )
1488 INTEGER IDIM( NIDIM ), INC( NINC )
1489 * .. Local Scalars ..
1490 COMPLEX*16 ALPHA, ALS, TRANSL
1491 DOUBLE PRECISION ERR, ERRMAX
1492 INTEGER I, IA, IM, IN, INCX, INCXS, INCY, INCYS, IX,
1493 $ IY, J, LAA, LDA, LDAS, LX, LY, M, MS, N, NARGS,
1494 $ NC, ND, NS
1495 LOGICAL CONJ, NULL, RESET, SAME
1496 * .. Local Arrays ..
1497 COMPLEX*16 W( 1 )
1498 LOGICAL ISAME( 13 )
1499 * .. External Functions ..
1500 LOGICAL LZE, LZERES
1501 EXTERNAL LZE, LZERES
1502 * .. External Subroutines ..
1503 EXTERNAL ZGERC, ZGERU, ZMAKE, ZMVCH
1504 * .. Intrinsic Functions ..
1505 INTRINSIC ABS, DCONJG, MAX, MIN
1506 * .. Scalars in Common ..
1507 INTEGER INFOT, NOUTC
1508 LOGICAL LERR, OK
1509 * .. Common blocks ..
1510 COMMON /INFOC/INFOT, NOUTC, OK, LERR
1511 * .. Executable Statements ..
1512 CONJ = SNAME( 5: 5 ).EQ.'C'
1513 * Define the number of arguments.
1514 NARGS = 9
1515 *
1516 NC = 0
1517 RESET = .TRUE.
1518 ERRMAX = RZERO
1519 *
1520 DO 120 IN = 1, NIDIM
1521 N = IDIM( IN )
1522 ND = N/2 + 1
1523 *
1524 DO 110 IM = 1, 2
1525 IF( IM.EQ.1 )
1526 $ M = MAX( N - ND, 0 )
1527 IF( IM.EQ.2 )
1528 $ M = MIN( N + ND, NMAX )
1529 *
1530 * Set LDA to 1 more than minimum value if room.
1531 LDA = M
1532 IF( LDA.LT.NMAX )
1533 $ LDA = LDA + 1
1534 * Skip tests if not enough room.
1535 IF( LDA.GT.NMAX )
1536 $ GO TO 110
1537 LAA = LDA*N
1538 NULL = N.LE.0.OR.M.LE.0
1539 *
1540 DO 100 IX = 1, NINC
1541 INCX = INC( IX )
1542 LX = ABS( INCX )*M
1543 *
1544 * Generate the vector X.
1545 *
1546 TRANSL = HALF
1547 CALL ZMAKE( 'GE', ' ', ' ', 1, M, X, 1, XX, ABS( INCX ),
1548 $ 0, M - 1, RESET, TRANSL )
1549 IF( M.GT.1 )THEN
1550 X( M/2 ) = ZERO
1551 XX( 1 + ABS( INCX )*( M/2 - 1 ) ) = ZERO
1552 END IF
1553 *
1554 DO 90 IY = 1, NINC
1555 INCY = INC( IY )
1556 LY = ABS( INCY )*N
1557 *
1558 * Generate the vector Y.
1559 *
1560 TRANSL = ZERO
1561 CALL ZMAKE( 'GE', ' ', ' ', 1, N, Y, 1, YY,
1562 $ ABS( INCY ), 0, N - 1, RESET, TRANSL )
1563 IF( N.GT.1 )THEN
1564 Y( N/2 ) = ZERO
1565 YY( 1 + ABS( INCY )*( N/2 - 1 ) ) = ZERO
1566 END IF
1567 *
1568 DO 80 IA = 1, NALF
1569 ALPHA = ALF( IA )
1570 *
1571 * Generate the matrix A.
1572 *
1573 TRANSL = ZERO
1574 CALL ZMAKE( SNAME( 2: 3 ), ' ', ' ', M, N, A, NMAX,
1575 $ AA, LDA, M - 1, N - 1, RESET, TRANSL )
1576 *
1577 NC = NC + 1
1578 *
1579 * Save every datum before calling the subroutine.
1580 *
1581 MS = M
1582 NS = N
1583 ALS = ALPHA
1584 DO 10 I = 1, LAA
1585 AS( I ) = AA( I )
1586 10 CONTINUE
1587 LDAS = LDA
1588 DO 20 I = 1, LX
1589 XS( I ) = XX( I )
1590 20 CONTINUE
1591 INCXS = INCX
1592 DO 30 I = 1, LY
1593 YS( I ) = YY( I )
1594 30 CONTINUE
1595 INCYS = INCY
1596 *
1597 * Call the subroutine.
1598 *
1599 IF( TRACE )
1600 $ WRITE( NTRA, FMT = 9994 )NC, SNAME, M, N,
1601 $ ALPHA, INCX, INCY, LDA
1602 IF( CONJ )THEN
1603 IF( REWI )
1604 $ REWIND NTRA
1605 CALL ZGERC( M, N, ALPHA, XX, INCX, YY, INCY, AA,
1606 $ LDA )
1607 ELSE
1608 IF( REWI )
1609 $ REWIND NTRA
1610 CALL ZGERU( M, N, ALPHA, XX, INCX, YY, INCY, AA,
1611 $ LDA )
1612 END IF
1613 *
1614 * Check if error-exit was taken incorrectly.
1615 *
1616 IF( .NOT.OK )THEN
1617 WRITE( NOUT, FMT = 9993 )
1618 FATAL = .TRUE.
1619 GO TO 140
1620 END IF
1621 *
1622 * See what data changed inside subroutine.
1623 *
1624 ISAME( 1 ) = MS.EQ.M
1625 ISAME( 2 ) = NS.EQ.N
1626 ISAME( 3 ) = ALS.EQ.ALPHA
1627 ISAME( 4 ) = LZE( XS, XX, LX )
1628 ISAME( 5 ) = INCXS.EQ.INCX
1629 ISAME( 6 ) = LZE( YS, YY, LY )
1630 ISAME( 7 ) = INCYS.EQ.INCY
1631 IF( NULL )THEN
1632 ISAME( 8 ) = LZE( AS, AA, LAA )
1633 ELSE
1634 ISAME( 8 ) = LZERES( 'GE', ' ', M, N, AS, AA,
1635 $ LDA )
1636 END IF
1637 ISAME( 9 ) = LDAS.EQ.LDA
1638 *
1639 * If data was incorrectly changed, report and return.
1640 *
1641 SAME = .TRUE.
1642 DO 40 I = 1, NARGS
1643 SAME = SAME.AND.ISAME( I )
1644 IF( .NOT.ISAME( I ) )
1645 $ WRITE( NOUT, FMT = 9998 )I
1646 40 CONTINUE
1647 IF( .NOT.SAME )THEN
1648 FATAL = .TRUE.
1649 GO TO 140
1650 END IF
1651 *
1652 IF( .NOT.NULL )THEN
1653 *
1654 * Check the result column by column.
1655 *
1656 IF( INCX.GT.0 )THEN
1657 DO 50 I = 1, M
1658 Z( I ) = X( I )
1659 50 CONTINUE
1660 ELSE
1661 DO 60 I = 1, M
1662 Z( I ) = X( M - I + 1 )
1663 60 CONTINUE
1664 END IF
1665 DO 70 J = 1, N
1666 IF( INCY.GT.0 )THEN
1667 W( 1 ) = Y( J )
1668 ELSE
1669 W( 1 ) = Y( N - J + 1 )
1670 END IF
1671 IF( CONJ )
1672 $ W( 1 ) = DCONJG( W( 1 ) )
1673 CALL ZMVCH( 'N', M, 1, ALPHA, Z, NMAX, W, 1,
1674 $ ONE, A( 1, J ), 1, YT, G,
1675 $ AA( 1 + ( J - 1 )*LDA ), EPS,
1676 $ ERR, FATAL, NOUT, .TRUE. )
1677 ERRMAX = MAX( ERRMAX, ERR )
1678 * If got really bad answer, report and return.
1679 IF( FATAL )
1680 $ GO TO 130
1681 70 CONTINUE
1682 ELSE
1683 * Avoid repeating tests with M.le.0 or N.le.0.
1684 GO TO 110
1685 END IF
1686 *
1687 80 CONTINUE
1688 *
1689 90 CONTINUE
1690 *
1691 100 CONTINUE
1692 *
1693 110 CONTINUE
1694 *
1695 120 CONTINUE
1696 *
1697 * Report result.
1698 *
1699 IF( ERRMAX.LT.THRESH )THEN
1700 WRITE( NOUT, FMT = 9999 )SNAME, NC
1701 ELSE
1702 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
1703 END IF
1704 GO TO 150
1705 *
1706 130 CONTINUE
1707 WRITE( NOUT, FMT = 9995 )J
1708 *
1709 140 CONTINUE
1710 WRITE( NOUT, FMT = 9996 )SNAME
1711 WRITE( NOUT, FMT = 9994 )NC, SNAME, M, N, ALPHA, INCX, INCY, LDA
1712 *
1713 150 CONTINUE
1714 RETURN
1715 *
1716 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
1717 $ 'S)' )
1718 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1719 $ 'ANGED INCORRECTLY *******' )
1720 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
1721 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
1722 $ ' - SUSPECT *******' )
1723 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
1724 9995 FORMAT( ' THESE ARE THE RESULTS FOR COLUMN ', I3 )
1725 9994 FORMAT( 1X, I6, ': ', A6, '(', 2( I3, ',' ), '(', F4.1, ',', F4.1,
1726 $ '), X,', I2, ', Y,', I2, ', A,', I3, ') ',
1727 $ ' .' )
1728 9993 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1729 $ '******' )
1730 *
1731 * End of ZCHK4.
1732 *
1733 END
1734 SUBROUTINE ZCHK5( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
1735 $ FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC, NMAX,
1736 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G,
1737 $ Z )
1738 *
1739 * Tests ZHER and ZHPR.
1740 *
1741 * Auxiliary routine for test program for Level 2 Blas.
1742 *
1743 * -- Written on 10-August-1987.
1744 * Richard Hanson, Sandia National Labs.
1745 * Jeremy Du Croz, NAG Central Office.
1746 *
1747 * .. Parameters ..
1748 COMPLEX*16 ZERO, HALF, ONE
1749 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
1750 $ HALF = ( 0.5D0, 0.0D0 ),
1751 $ ONE = ( 1.0D0, 0.0D0 ) )
1752 DOUBLE PRECISION RZERO
1753 PARAMETER ( RZERO = 0.0D0 )
1754 * .. Scalar Arguments ..
1755 DOUBLE PRECISION EPS, THRESH
1756 INTEGER INCMAX, NALF, NIDIM, NINC, NMAX, NOUT, NTRA
1757 LOGICAL FATAL, REWI, TRACE
1758 CHARACTER*6 SNAME
1759 * .. Array Arguments ..
1760 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ), ALF( NALF ),
1761 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ),
1762 $ XX( NMAX*INCMAX ), Y( NMAX ),
1763 $ YS( NMAX*INCMAX ), YT( NMAX ),
1764 $ YY( NMAX*INCMAX ), Z( NMAX )
1765 DOUBLE PRECISION G( NMAX )
1766 INTEGER IDIM( NIDIM ), INC( NINC )
1767 * .. Local Scalars ..
1768 COMPLEX*16 ALPHA, TRANSL
1769 DOUBLE PRECISION ERR, ERRMAX, RALPHA, RALS
1770 INTEGER I, IA, IC, IN, INCX, INCXS, IX, J, JA, JJ, LAA,
1771 $ LDA, LDAS, LJ, LX, N, NARGS, NC, NS
1772 LOGICAL FULL, NULL, PACKED, RESET, SAME, UPPER
1773 CHARACTER*1 UPLO, UPLOS
1774 CHARACTER*2 ICH
1775 * .. Local Arrays ..
1776 COMPLEX*16 W( 1 )
1777 LOGICAL ISAME( 13 )
1778 * .. External Functions ..
1779 LOGICAL LZE, LZERES
1780 EXTERNAL LZE, LZERES
1781 * .. External Subroutines ..
1782 EXTERNAL ZHER, ZHPR, ZMAKE, ZMVCH
1783 * .. Intrinsic Functions ..
1784 INTRINSIC ABS, DBLE, DCMPLX, DCONJG, MAX
1785 * .. Scalars in Common ..
1786 INTEGER INFOT, NOUTC
1787 LOGICAL LERR, OK
1788 * .. Common blocks ..
1789 COMMON /INFOC/INFOT, NOUTC, OK, LERR
1790 * .. Data statements ..
1791 DATA ICH/'UL'/
1792 * .. Executable Statements ..
1793 FULL = SNAME( 3: 3 ).EQ.'E'
1794 PACKED = SNAME( 3: 3 ).EQ.'P'
1795 * Define the number of arguments.
1796 IF( FULL )THEN
1797 NARGS = 7
1798 ELSE IF( PACKED )THEN
1799 NARGS = 6
1800 END IF
1801 *
1802 NC = 0
1803 RESET = .TRUE.
1804 ERRMAX = RZERO
1805 *
1806 DO 100 IN = 1, NIDIM
1807 N = IDIM( IN )
1808 * Set LDA to 1 more than minimum value if room.
1809 LDA = N
1810 IF( LDA.LT.NMAX )
1811 $ LDA = LDA + 1
1812 * Skip tests if not enough room.
1813 IF( LDA.GT.NMAX )
1814 $ GO TO 100
1815 IF( PACKED )THEN
1816 LAA = ( N*( N + 1 ) )/2
1817 ELSE
1818 LAA = LDA*N
1819 END IF
1820 *
1821 DO 90 IC = 1, 2
1822 UPLO = ICH( IC: IC )
1823 UPPER = UPLO.EQ.'U'
1824 *
1825 DO 80 IX = 1, NINC
1826 INCX = INC( IX )
1827 LX = ABS( INCX )*N
1828 *
1829 * Generate the vector X.
1830 *
1831 TRANSL = HALF
1832 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX, ABS( INCX ),
1833 $ 0, N - 1, RESET, TRANSL )
1834 IF( N.GT.1 )THEN
1835 X( N/2 ) = ZERO
1836 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
1837 END IF
1838 *
1839 DO 70 IA = 1, NALF
1840 RALPHA = DBLE( ALF( IA ) )
1841 ALPHA = DCMPLX( RALPHA, RZERO )
1842 NULL = N.LE.0.OR.RALPHA.EQ.RZERO
1843 *
1844 * Generate the matrix A.
1845 *
1846 TRANSL = ZERO
1847 CALL ZMAKE( SNAME( 2: 3 ), UPLO, ' ', N, N, A, NMAX,
1848 $ AA, LDA, N - 1, N - 1, RESET, TRANSL )
1849 *
1850 NC = NC + 1
1851 *
1852 * Save every datum before calling the subroutine.
1853 *
1854 UPLOS = UPLO
1855 NS = N
1856 RALS = RALPHA
1857 DO 10 I = 1, LAA
1858 AS( I ) = AA( I )
1859 10 CONTINUE
1860 LDAS = LDA
1861 DO 20 I = 1, LX
1862 XS( I ) = XX( I )
1863 20 CONTINUE
1864 INCXS = INCX
1865 *
1866 * Call the subroutine.
1867 *
1868 IF( FULL )THEN
1869 IF( TRACE )
1870 $ WRITE( NTRA, FMT = 9993 )NC, SNAME, UPLO, N,
1871 $ RALPHA, INCX, LDA
1872 IF( REWI )
1873 $ REWIND NTRA
1874 CALL ZHER( UPLO, N, RALPHA, XX, INCX, AA, LDA )
1875 ELSE IF( PACKED )THEN
1876 IF( TRACE )
1877 $ WRITE( NTRA, FMT = 9994 )NC, SNAME, UPLO, N,
1878 $ RALPHA, INCX
1879 IF( REWI )
1880 $ REWIND NTRA
1881 CALL ZHPR( UPLO, N, RALPHA, XX, INCX, AA )
1882 END IF
1883 *
1884 * Check if error-exit was taken incorrectly.
1885 *
1886 IF( .NOT.OK )THEN
1887 WRITE( NOUT, FMT = 9992 )
1888 FATAL = .TRUE.
1889 GO TO 120
1890 END IF
1891 *
1892 * See what data changed inside subroutines.
1893 *
1894 ISAME( 1 ) = UPLO.EQ.UPLOS
1895 ISAME( 2 ) = NS.EQ.N
1896 ISAME( 3 ) = RALS.EQ.RALPHA
1897 ISAME( 4 ) = LZE( XS, XX, LX )
1898 ISAME( 5 ) = INCXS.EQ.INCX
1899 IF( NULL )THEN
1900 ISAME( 6 ) = LZE( AS, AA, LAA )
1901 ELSE
1902 ISAME( 6 ) = LZERES( SNAME( 2: 3 ), UPLO, N, N, AS,
1903 $ AA, LDA )
1904 END IF
1905 IF( .NOT.PACKED )THEN
1906 ISAME( 7 ) = LDAS.EQ.LDA
1907 END IF
1908 *
1909 * If data was incorrectly changed, report and return.
1910 *
1911 SAME = .TRUE.
1912 DO 30 I = 1, NARGS
1913 SAME = SAME.AND.ISAME( I )
1914 IF( .NOT.ISAME( I ) )
1915 $ WRITE( NOUT, FMT = 9998 )I
1916 30 CONTINUE
1917 IF( .NOT.SAME )THEN
1918 FATAL = .TRUE.
1919 GO TO 120
1920 END IF
1921 *
1922 IF( .NOT.NULL )THEN
1923 *
1924 * Check the result column by column.
1925 *
1926 IF( INCX.GT.0 )THEN
1927 DO 40 I = 1, N
1928 Z( I ) = X( I )
1929 40 CONTINUE
1930 ELSE
1931 DO 50 I = 1, N
1932 Z( I ) = X( N - I + 1 )
1933 50 CONTINUE
1934 END IF
1935 JA = 1
1936 DO 60 J = 1, N
1937 W( 1 ) = DCONJG( Z( J ) )
1938 IF( UPPER )THEN
1939 JJ = 1
1940 LJ = J
1941 ELSE
1942 JJ = J
1943 LJ = N - J + 1
1944 END IF
1945 CALL ZMVCH( 'N', LJ, 1, ALPHA, Z( JJ ), LJ, W,
1946 $ 1, ONE, A( JJ, J ), 1, YT, G,
1947 $ AA( JA ), EPS, ERR, FATAL, NOUT,
1948 $ .TRUE. )
1949 IF( FULL )THEN
1950 IF( UPPER )THEN
1951 JA = JA + LDA
1952 ELSE
1953 JA = JA + LDA + 1
1954 END IF
1955 ELSE
1956 JA = JA + LJ
1957 END IF
1958 ERRMAX = MAX( ERRMAX, ERR )
1959 * If got really bad answer, report and return.
1960 IF( FATAL )
1961 $ GO TO 110
1962 60 CONTINUE
1963 ELSE
1964 * Avoid repeating tests if N.le.0.
1965 IF( N.LE.0 )
1966 $ GO TO 100
1967 END IF
1968 *
1969 70 CONTINUE
1970 *
1971 80 CONTINUE
1972 *
1973 90 CONTINUE
1974 *
1975 100 CONTINUE
1976 *
1977 * Report result.
1978 *
1979 IF( ERRMAX.LT.THRESH )THEN
1980 WRITE( NOUT, FMT = 9999 )SNAME, NC
1981 ELSE
1982 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
1983 END IF
1984 GO TO 130
1985 *
1986 110 CONTINUE
1987 WRITE( NOUT, FMT = 9995 )J
1988 *
1989 120 CONTINUE
1990 WRITE( NOUT, FMT = 9996 )SNAME
1991 IF( FULL )THEN
1992 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, N, RALPHA, INCX, LDA
1993 ELSE IF( PACKED )THEN
1994 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, N, RALPHA, INCX
1995 END IF
1996 *
1997 130 CONTINUE
1998 RETURN
1999 *
2000 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
2001 $ 'S)' )
2002 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
2003 $ 'ANGED INCORRECTLY *******' )
2004 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
2005 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
2006 $ ' - SUSPECT *******' )
2007 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
2008 9995 FORMAT( ' THESE ARE THE RESULTS FOR COLUMN ', I3 )
2009 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', X,',
2010 $ I2, ', AP) .' )
2011 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', X,',
2012 $ I2, ', A,', I3, ') .' )
2013 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
2014 $ '******' )
2015 *
2016 * End of ZCHK5.
2017 *
2018 END
2019 SUBROUTINE ZCHK6( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
2020 $ FATAL, NIDIM, IDIM, NALF, ALF, NINC, INC, NMAX,
2021 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G,
2022 $ Z )
2023 *
2024 * Tests ZHER2 and ZHPR2.
2025 *
2026 * Auxiliary routine for test program for Level 2 Blas.
2027 *
2028 * -- Written on 10-August-1987.
2029 * Richard Hanson, Sandia National Labs.
2030 * Jeremy Du Croz, NAG Central Office.
2031 *
2032 * .. Parameters ..
2033 COMPLEX*16 ZERO, HALF, ONE
2034 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
2035 $ HALF = ( 0.5D0, 0.0D0 ),
2036 $ ONE = ( 1.0D0, 0.0D0 ) )
2037 DOUBLE PRECISION RZERO
2038 PARAMETER ( RZERO = 0.0D0 )
2039 * .. Scalar Arguments ..
2040 DOUBLE PRECISION EPS, THRESH
2041 INTEGER INCMAX, NALF, NIDIM, NINC, NMAX, NOUT, NTRA
2042 LOGICAL FATAL, REWI, TRACE
2043 CHARACTER*6 SNAME
2044 * .. Array Arguments ..
2045 COMPLEX*16 A( NMAX, NMAX ), AA( NMAX*NMAX ), ALF( NALF ),
2046 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ),
2047 $ XX( NMAX*INCMAX ), Y( NMAX ),
2048 $ YS( NMAX*INCMAX ), YT( NMAX ),
2049 $ YY( NMAX*INCMAX ), Z( NMAX, 2 )
2050 DOUBLE PRECISION G( NMAX )
2051 INTEGER IDIM( NIDIM ), INC( NINC )
2052 * .. Local Scalars ..
2053 COMPLEX*16 ALPHA, ALS, TRANSL
2054 DOUBLE PRECISION ERR, ERRMAX
2055 INTEGER I, IA, IC, IN, INCX, INCXS, INCY, INCYS, IX,
2056 $ IY, J, JA, JJ, LAA, LDA, LDAS, LJ, LX, LY, N,
2057 $ NARGS, NC, NS
2058 LOGICAL FULL, NULL, PACKED, RESET, SAME, UPPER
2059 CHARACTER*1 UPLO, UPLOS
2060 CHARACTER*2 ICH
2061 * .. Local Arrays ..
2062 COMPLEX*16 W( 2 )
2063 LOGICAL ISAME( 13 )
2064 * .. External Functions ..
2065 LOGICAL LZE, LZERES
2066 EXTERNAL LZE, LZERES
2067 * .. External Subroutines ..
2068 EXTERNAL ZHER2, ZHPR2, ZMAKE, ZMVCH
2069 * .. Intrinsic Functions ..
2070 INTRINSIC ABS, DCONJG, MAX
2071 * .. Scalars in Common ..
2072 INTEGER INFOT, NOUTC
2073 LOGICAL LERR, OK
2074 * .. Common blocks ..
2075 COMMON /INFOC/INFOT, NOUTC, OK, LERR
2076 * .. Data statements ..
2077 DATA ICH/'UL'/
2078 * .. Executable Statements ..
2079 FULL = SNAME( 3: 3 ).EQ.'E'
2080 PACKED = SNAME( 3: 3 ).EQ.'P'
2081 * Define the number of arguments.
2082 IF( FULL )THEN
2083 NARGS = 9
2084 ELSE IF( PACKED )THEN
2085 NARGS = 8
2086 END IF
2087 *
2088 NC = 0
2089 RESET = .TRUE.
2090 ERRMAX = RZERO
2091 *
2092 DO 140 IN = 1, NIDIM
2093 N = IDIM( IN )
2094 * Set LDA to 1 more than minimum value if room.
2095 LDA = N
2096 IF( LDA.LT.NMAX )
2097 $ LDA = LDA + 1
2098 * Skip tests if not enough room.
2099 IF( LDA.GT.NMAX )
2100 $ GO TO 140
2101 IF( PACKED )THEN
2102 LAA = ( N*( N + 1 ) )/2
2103 ELSE
2104 LAA = LDA*N
2105 END IF
2106 *
2107 DO 130 IC = 1, 2
2108 UPLO = ICH( IC: IC )
2109 UPPER = UPLO.EQ.'U'
2110 *
2111 DO 120 IX = 1, NINC
2112 INCX = INC( IX )
2113 LX = ABS( INCX )*N
2114 *
2115 * Generate the vector X.
2116 *
2117 TRANSL = HALF
2118 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX, ABS( INCX ),
2119 $ 0, N - 1, RESET, TRANSL )
2120 IF( N.GT.1 )THEN
2121 X( N/2 ) = ZERO
2122 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
2123 END IF
2124 *
2125 DO 110 IY = 1, NINC
2126 INCY = INC( IY )
2127 LY = ABS( INCY )*N
2128 *
2129 * Generate the vector Y.
2130 *
2131 TRANSL = ZERO
2132 CALL ZMAKE( 'GE', ' ', ' ', 1, N, Y, 1, YY,
2133 $ ABS( INCY ), 0, N - 1, RESET, TRANSL )
2134 IF( N.GT.1 )THEN
2135 Y( N/2 ) = ZERO
2136 YY( 1 + ABS( INCY )*( N/2 - 1 ) ) = ZERO
2137 END IF
2138 *
2139 DO 100 IA = 1, NALF
2140 ALPHA = ALF( IA )
2141 NULL = N.LE.0.OR.ALPHA.EQ.ZERO
2142 *
2143 * Generate the matrix A.
2144 *
2145 TRANSL = ZERO
2146 CALL ZMAKE( SNAME( 2: 3 ), UPLO, ' ', N, N, A,
2147 $ NMAX, AA, LDA, N - 1, N - 1, RESET,
2148 $ TRANSL )
2149 *
2150 NC = NC + 1
2151 *
2152 * Save every datum before calling the subroutine.
2153 *
2154 UPLOS = UPLO
2155 NS = N
2156 ALS = ALPHA
2157 DO 10 I = 1, LAA
2158 AS( I ) = AA( I )
2159 10 CONTINUE
2160 LDAS = LDA
2161 DO 20 I = 1, LX
2162 XS( I ) = XX( I )
2163 20 CONTINUE
2164 INCXS = INCX
2165 DO 30 I = 1, LY
2166 YS( I ) = YY( I )
2167 30 CONTINUE
2168 INCYS = INCY
2169 *
2170 * Call the subroutine.
2171 *
2172 IF( FULL )THEN
2173 IF( TRACE )
2174 $ WRITE( NTRA, FMT = 9993 )NC, SNAME, UPLO, N,
2175 $ ALPHA, INCX, INCY, LDA
2176 IF( REWI )
2177 $ REWIND NTRA
2178 CALL ZHER2( UPLO, N, ALPHA, XX, INCX, YY, INCY,
2179 $ AA, LDA )
2180 ELSE IF( PACKED )THEN
2181 IF( TRACE )
2182 $ WRITE( NTRA, FMT = 9994 )NC, SNAME, UPLO, N,
2183 $ ALPHA, INCX, INCY
2184 IF( REWI )
2185 $ REWIND NTRA
2186 CALL ZHPR2( UPLO, N, ALPHA, XX, INCX, YY, INCY,
2187 $ AA )
2188 END IF
2189 *
2190 * Check if error-exit was taken incorrectly.
2191 *
2192 IF( .NOT.OK )THEN
2193 WRITE( NOUT, FMT = 9992 )
2194 FATAL = .TRUE.
2195 GO TO 160
2196 END IF
2197 *
2198 * See what data changed inside subroutines.
2199 *
2200 ISAME( 1 ) = UPLO.EQ.UPLOS
2201 ISAME( 2 ) = NS.EQ.N
2202 ISAME( 3 ) = ALS.EQ.ALPHA
2203 ISAME( 4 ) = LZE( XS, XX, LX )
2204 ISAME( 5 ) = INCXS.EQ.INCX
2205 ISAME( 6 ) = LZE( YS, YY, LY )
2206 ISAME( 7 ) = INCYS.EQ.INCY
2207 IF( NULL )THEN
2208 ISAME( 8 ) = LZE( AS, AA, LAA )
2209 ELSE
2210 ISAME( 8 ) = LZERES( SNAME( 2: 3 ), UPLO, N, N,
2211 $ AS, AA, LDA )
2212 END IF
2213 IF( .NOT.PACKED )THEN
2214 ISAME( 9 ) = LDAS.EQ.LDA
2215 END IF
2216 *
2217 * If data was incorrectly changed, report and return.
2218 *
2219 SAME = .TRUE.
2220 DO 40 I = 1, NARGS
2221 SAME = SAME.AND.ISAME( I )
2222 IF( .NOT.ISAME( I ) )
2223 $ WRITE( NOUT, FMT = 9998 )I
2224 40 CONTINUE
2225 IF( .NOT.SAME )THEN
2226 FATAL = .TRUE.
2227 GO TO 160
2228 END IF
2229 *
2230 IF( .NOT.NULL )THEN
2231 *
2232 * Check the result column by column.
2233 *
2234 IF( INCX.GT.0 )THEN
2235 DO 50 I = 1, N
2236 Z( I, 1 ) = X( I )
2237 50 CONTINUE
2238 ELSE
2239 DO 60 I = 1, N
2240 Z( I, 1 ) = X( N - I + 1 )
2241 60 CONTINUE
2242 END IF
2243 IF( INCY.GT.0 )THEN
2244 DO 70 I = 1, N
2245 Z( I, 2 ) = Y( I )
2246 70 CONTINUE
2247 ELSE
2248 DO 80 I = 1, N
2249 Z( I, 2 ) = Y( N - I + 1 )
2250 80 CONTINUE
2251 END IF
2252 JA = 1
2253 DO 90 J = 1, N
2254 W( 1 ) = ALPHA*DCONJG( Z( J, 2 ) )
2255 W( 2 ) = DCONJG( ALPHA )*DCONJG( Z( J, 1 ) )
2256 IF( UPPER )THEN
2257 JJ = 1
2258 LJ = J
2259 ELSE
2260 JJ = J
2261 LJ = N - J + 1
2262 END IF
2263 CALL ZMVCH( 'N', LJ, 2, ONE, Z( JJ, 1 ),
2264 $ NMAX, W, 1, ONE, A( JJ, J ), 1,
2265 $ YT, G, AA( JA ), EPS, ERR, FATAL,
2266 $ NOUT, .TRUE. )
2267 IF( FULL )THEN
2268 IF( UPPER )THEN
2269 JA = JA + LDA
2270 ELSE
2271 JA = JA + LDA + 1
2272 END IF
2273 ELSE
2274 JA = JA + LJ
2275 END IF
2276 ERRMAX = MAX( ERRMAX, ERR )
2277 * If got really bad answer, report and return.
2278 IF( FATAL )
2279 $ GO TO 150
2280 90 CONTINUE
2281 ELSE
2282 * Avoid repeating tests with N.le.0.
2283 IF( N.LE.0 )
2284 $ GO TO 140
2285 END IF
2286 *
2287 100 CONTINUE
2288 *
2289 110 CONTINUE
2290 *
2291 120 CONTINUE
2292 *
2293 130 CONTINUE
2294 *
2295 140 CONTINUE
2296 *
2297 * Report result.
2298 *
2299 IF( ERRMAX.LT.THRESH )THEN
2300 WRITE( NOUT, FMT = 9999 )SNAME, NC
2301 ELSE
2302 WRITE( NOUT, FMT = 9997 )SNAME, NC, ERRMAX
2303 END IF
2304 GO TO 170
2305 *
2306 150 CONTINUE
2307 WRITE( NOUT, FMT = 9995 )J
2308 *
2309 160 CONTINUE
2310 WRITE( NOUT, FMT = 9996 )SNAME
2311 IF( FULL )THEN
2312 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, N, ALPHA, INCX,
2313 $ INCY, LDA
2314 ELSE IF( PACKED )THEN
2315 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, N, ALPHA, INCX, INCY
2316 END IF
2317 *
2318 170 CONTINUE
2319 RETURN
2320 *
2321 9999 FORMAT( ' ', A6, ' PASSED THE COMPUTATIONAL TESTS (', I6, ' CALL',
2322 $ 'S)' )
2323 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
2324 $ 'ANGED INCORRECTLY *******' )
2325 9997 FORMAT( ' ', A6, ' COMPLETED THE COMPUTATIONAL TESTS (', I6, ' C',
2326 $ 'ALLS)', /' ******* BUT WITH MAXIMUM TEST RATIO', F8.2,
2327 $ ' - SUSPECT *******' )
2328 9996 FORMAT( ' ******* ', A6, ' FAILED ON CALL NUMBER:' )
2329 9995 FORMAT( ' THESE ARE THE RESULTS FOR COLUMN ', I3 )
2330 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
2331 $ F4.1, '), X,', I2, ', Y,', I2, ', AP) ',
2332 $ ' .' )
2333 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
2334 $ F4.1, '), X,', I2, ', Y,', I2, ', A,', I3, ') ',
2335 $ ' .' )
2336 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
2337 $ '******' )
2338 *
2339 * End of ZCHK6.
2340 *
2341 END
2342 SUBROUTINE ZCHKE( ISNUM, SRNAMT, NOUT )
2343 *
2344 * Tests the error exits from the Level 2 Blas.
2345 * Requires a special version of the error-handling routine XERBLA.
2346 * ALPHA, RALPHA, BETA, A, X and Y should not need to be defined.
2347 *
2348 * Auxiliary routine for test program for Level 2 Blas.
2349 *
2350 * -- Written on 10-August-1987.
2351 * Richard Hanson, Sandia National Labs.
2352 * Jeremy Du Croz, NAG Central Office.
2353 *
2354 * .. Scalar Arguments ..
2355 INTEGER ISNUM, NOUT
2356 CHARACTER*6 SRNAMT
2357 * .. Scalars in Common ..
2358 INTEGER INFOT, NOUTC
2359 LOGICAL LERR, OK
2360 * .. Local Scalars ..
2361 COMPLEX*16 ALPHA, BETA
2362 DOUBLE PRECISION RALPHA
2363 * .. Local Arrays ..
2364 COMPLEX*16 A( 1, 1 ), X( 1 ), Y( 1 )
2365 * .. External Subroutines ..
2366 EXTERNAL CHKXER, ZGBMV, ZGEMV, ZGERC, ZGERU, ZHBMV,
2367 $ ZHEMV, ZHER, ZHER2, ZHPMV, ZHPR, ZHPR2, ZTBMV,
2368 $ ZTBSV, ZTPMV, ZTPSV, ZTRMV, ZTRSV
2369 * .. Common blocks ..
2370 COMMON /INFOC/INFOT, NOUTC, OK, LERR
2371 * .. Executable Statements ..
2372 * OK is set to .FALSE. by the special version of XERBLA or by CHKXER
2373 * if anything is wrong.
2374 OK = .TRUE.
2375 * LERR is set to .TRUE. by the special version of XERBLA each time
2376 * it is called, and is then tested and re-set by CHKXER.
2377 LERR = .FALSE.
2378 GO TO ( 10, 20, 30, 40, 50, 60, 70, 80,
2379 $ 90, 100, 110, 120, 130, 140, 150, 160,
2380 $ 170 )ISNUM
2381 10 INFOT = 1
2382 CALL ZGEMV( '/', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2383 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2384 INFOT = 2
2385 CALL ZGEMV( 'N', -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2386 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2387 INFOT = 3
2388 CALL ZGEMV( 'N', 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2389 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2390 INFOT = 6
2391 CALL ZGEMV( 'N', 2, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2392 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2393 INFOT = 8
2394 CALL ZGEMV( 'N', 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2395 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2396 INFOT = 11
2397 CALL ZGEMV( 'N', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2398 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2399 GO TO 180
2400 20 INFOT = 1
2401 CALL ZGBMV( '/', 0, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2402 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2403 INFOT = 2
2404 CALL ZGBMV( 'N', -1, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2405 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2406 INFOT = 3
2407 CALL ZGBMV( 'N', 0, -1, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2408 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2409 INFOT = 4
2410 CALL ZGBMV( 'N', 0, 0, -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2411 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2412 INFOT = 5
2413 CALL ZGBMV( 'N', 2, 0, 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2414 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2415 INFOT = 8
2416 CALL ZGBMV( 'N', 0, 0, 1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2417 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2418 INFOT = 10
2419 CALL ZGBMV( 'N', 0, 0, 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2420 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2421 INFOT = 13
2422 CALL ZGBMV( 'N', 0, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2423 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2424 GO TO 180
2425 30 INFOT = 1
2426 CALL ZHEMV( '/', 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2427 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2428 INFOT = 2
2429 CALL ZHEMV( 'U', -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2430 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2431 INFOT = 5
2432 CALL ZHEMV( 'U', 2, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2433 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2434 INFOT = 7
2435 CALL ZHEMV( 'U', 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2436 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2437 INFOT = 10
2438 CALL ZHEMV( 'U', 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2439 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2440 GO TO 180
2441 40 INFOT = 1
2442 CALL ZHBMV( '/', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2443 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2444 INFOT = 2
2445 CALL ZHBMV( 'U', -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2446 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2447 INFOT = 3
2448 CALL ZHBMV( 'U', 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2449 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2450 INFOT = 6
2451 CALL ZHBMV( 'U', 0, 1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2452 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2453 INFOT = 8
2454 CALL ZHBMV( 'U', 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2455 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2456 INFOT = 11
2457 CALL ZHBMV( 'U', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2458 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2459 GO TO 180
2460 50 INFOT = 1
2461 CALL ZHPMV( '/', 0, ALPHA, A, X, 1, BETA, Y, 1 )
2462 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2463 INFOT = 2
2464 CALL ZHPMV( 'U', -1, ALPHA, A, X, 1, BETA, Y, 1 )
2465 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2466 INFOT = 6
2467 CALL ZHPMV( 'U', 0, ALPHA, A, X, 0, BETA, Y, 1 )
2468 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2469 INFOT = 9
2470 CALL ZHPMV( 'U', 0, ALPHA, A, X, 1, BETA, Y, 0 )
2471 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2472 GO TO 180
2473 60 INFOT = 1
2474 CALL ZTRMV( '/', 'N', 'N', 0, A, 1, X, 1 )
2475 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2476 INFOT = 2
2477 CALL ZTRMV( 'U', '/', 'N', 0, A, 1, X, 1 )
2478 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2479 INFOT = 3
2480 CALL ZTRMV( 'U', 'N', '/', 0, A, 1, X, 1 )
2481 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2482 INFOT = 4
2483 CALL ZTRMV( 'U', 'N', 'N', -1, A, 1, X, 1 )
2484 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2485 INFOT = 6
2486 CALL ZTRMV( 'U', 'N', 'N', 2, A, 1, X, 1 )
2487 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2488 INFOT = 8
2489 CALL ZTRMV( 'U', 'N', 'N', 0, A, 1, X, 0 )
2490 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2491 GO TO 180
2492 70 INFOT = 1
2493 CALL ZTBMV( '/', 'N', 'N', 0, 0, A, 1, X, 1 )
2494 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2495 INFOT = 2
2496 CALL ZTBMV( 'U', '/', 'N', 0, 0, A, 1, X, 1 )
2497 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2498 INFOT = 3
2499 CALL ZTBMV( 'U', 'N', '/', 0, 0, A, 1, X, 1 )
2500 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2501 INFOT = 4
2502 CALL ZTBMV( 'U', 'N', 'N', -1, 0, A, 1, X, 1 )
2503 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2504 INFOT = 5
2505 CALL ZTBMV( 'U', 'N', 'N', 0, -1, A, 1, X, 1 )
2506 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2507 INFOT = 7
2508 CALL ZTBMV( 'U', 'N', 'N', 0, 1, A, 1, X, 1 )
2509 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2510 INFOT = 9
2511 CALL ZTBMV( 'U', 'N', 'N', 0, 0, A, 1, X, 0 )
2512 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2513 GO TO 180
2514 80 INFOT = 1
2515 CALL ZTPMV( '/', 'N', 'N', 0, A, X, 1 )
2516 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2517 INFOT = 2
2518 CALL ZTPMV( 'U', '/', 'N', 0, A, X, 1 )
2519 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2520 INFOT = 3
2521 CALL ZTPMV( 'U', 'N', '/', 0, A, X, 1 )
2522 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2523 INFOT = 4
2524 CALL ZTPMV( 'U', 'N', 'N', -1, A, X, 1 )
2525 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2526 INFOT = 7
2527 CALL ZTPMV( 'U', 'N', 'N', 0, A, X, 0 )
2528 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2529 GO TO 180
2530 90 INFOT = 1
2531 CALL ZTRSV( '/', 'N', 'N', 0, A, 1, X, 1 )
2532 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2533 INFOT = 2
2534 CALL ZTRSV( 'U', '/', 'N', 0, A, 1, X, 1 )
2535 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2536 INFOT = 3
2537 CALL ZTRSV( 'U', 'N', '/', 0, A, 1, X, 1 )
2538 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2539 INFOT = 4
2540 CALL ZTRSV( 'U', 'N', 'N', -1, A, 1, X, 1 )
2541 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2542 INFOT = 6
2543 CALL ZTRSV( 'U', 'N', 'N', 2, A, 1, X, 1 )
2544 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2545 INFOT = 8
2546 CALL ZTRSV( 'U', 'N', 'N', 0, A, 1, X, 0 )
2547 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2548 GO TO 180
2549 100 INFOT = 1
2550 CALL ZTBSV( '/', 'N', 'N', 0, 0, A, 1, X, 1 )
2551 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2552 INFOT = 2
2553 CALL ZTBSV( 'U', '/', 'N', 0, 0, A, 1, X, 1 )
2554 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2555 INFOT = 3
2556 CALL ZTBSV( 'U', 'N', '/', 0, 0, A, 1, X, 1 )
2557 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2558 INFOT = 4
2559 CALL ZTBSV( 'U', 'N', 'N', -1, 0, A, 1, X, 1 )
2560 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2561 INFOT = 5
2562 CALL ZTBSV( 'U', 'N', 'N', 0, -1, A, 1, X, 1 )
2563 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2564 INFOT = 7
2565 CALL ZTBSV( 'U', 'N', 'N', 0, 1, A, 1, X, 1 )
2566 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2567 INFOT = 9
2568 CALL ZTBSV( 'U', 'N', 'N', 0, 0, A, 1, X, 0 )
2569 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2570 GO TO 180
2571 110 INFOT = 1
2572 CALL ZTPSV( '/', 'N', 'N', 0, A, X, 1 )
2573 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2574 INFOT = 2
2575 CALL ZTPSV( 'U', '/', 'N', 0, A, X, 1 )
2576 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2577 INFOT = 3
2578 CALL ZTPSV( 'U', 'N', '/', 0, A, X, 1 )
2579 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2580 INFOT = 4
2581 CALL ZTPSV( 'U', 'N', 'N', -1, A, X, 1 )
2582 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2583 INFOT = 7
2584 CALL ZTPSV( 'U', 'N', 'N', 0, A, X, 0 )
2585 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2586 GO TO 180
2587 120 INFOT = 1
2588 CALL ZGERC( -1, 0, ALPHA, X, 1, Y, 1, A, 1 )
2589 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2590 INFOT = 2
2591 CALL ZGERC( 0, -1, ALPHA, X, 1, Y, 1, A, 1 )
2592 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2593 INFOT = 5
2594 CALL ZGERC( 0, 0, ALPHA, X, 0, Y, 1, A, 1 )
2595 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2596 INFOT = 7
2597 CALL ZGERC( 0, 0, ALPHA, X, 1, Y, 0, A, 1 )
2598 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2599 INFOT = 9
2600 CALL ZGERC( 2, 0, ALPHA, X, 1, Y, 1, A, 1 )
2601 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2602 GO TO 180
2603 130 INFOT = 1
2604 CALL ZGERU( -1, 0, ALPHA, X, 1, Y, 1, A, 1 )
2605 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2606 INFOT = 2
2607 CALL ZGERU( 0, -1, ALPHA, X, 1, Y, 1, A, 1 )
2608 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2609 INFOT = 5
2610 CALL ZGERU( 0, 0, ALPHA, X, 0, Y, 1, A, 1 )
2611 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2612 INFOT = 7
2613 CALL ZGERU( 0, 0, ALPHA, X, 1, Y, 0, A, 1 )
2614 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2615 INFOT = 9
2616 CALL ZGERU( 2, 0, ALPHA, X, 1, Y, 1, A, 1 )
2617 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2618 GO TO 180
2619 140 INFOT = 1
2620 CALL ZHER( '/', 0, RALPHA, X, 1, A, 1 )
2621 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2622 INFOT = 2
2623 CALL ZHER( 'U', -1, RALPHA, X, 1, A, 1 )
2624 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2625 INFOT = 5
2626 CALL ZHER( 'U', 0, RALPHA, X, 0, A, 1 )
2627 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2628 INFOT = 7
2629 CALL ZHER( 'U', 2, RALPHA, X, 1, A, 1 )
2630 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2631 GO TO 180
2632 150 INFOT = 1
2633 CALL ZHPR( '/', 0, RALPHA, X, 1, A )
2634 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2635 INFOT = 2
2636 CALL ZHPR( 'U', -1, RALPHA, X, 1, A )
2637 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2638 INFOT = 5
2639 CALL ZHPR( 'U', 0, RALPHA, X, 0, A )
2640 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2641 GO TO 180
2642 160 INFOT = 1
2643 CALL ZHER2( '/', 0, ALPHA, X, 1, Y, 1, A, 1 )
2644 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2645 INFOT = 2
2646 CALL ZHER2( 'U', -1, ALPHA, X, 1, Y, 1, A, 1 )
2647 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2648 INFOT = 5
2649 CALL ZHER2( 'U', 0, ALPHA, X, 0, Y, 1, A, 1 )
2650 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2651 INFOT = 7
2652 CALL ZHER2( 'U', 0, ALPHA, X, 1, Y, 0, A, 1 )
2653 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2654 INFOT = 9
2655 CALL ZHER2( 'U', 2, ALPHA, X, 1, Y, 1, A, 1 )
2656 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2657 GO TO 180
2658 170 INFOT = 1
2659 CALL ZHPR2( '/', 0, ALPHA, X, 1, Y, 1, A )
2660 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2661 INFOT = 2
2662 CALL ZHPR2( 'U', -1, ALPHA, X, 1, Y, 1, A )
2663 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2664 INFOT = 5
2665 CALL ZHPR2( 'U', 0, ALPHA, X, 0, Y, 1, A )
2666 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2667 INFOT = 7
2668 CALL ZHPR2( 'U', 0, ALPHA, X, 1, Y, 0, A )
2669 CALL CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
2670 *
2671 180 IF( OK )THEN
2672 WRITE( NOUT, FMT = 9999 )SRNAMT
2673 ELSE
2674 WRITE( NOUT, FMT = 9998 )SRNAMT
2675 END IF
2676 RETURN
2677 *
2678 9999 FORMAT( ' ', A6, ' PASSED THE TESTS OF ERROR-EXITS' )
2679 9998 FORMAT( ' ******* ', A6, ' FAILED THE TESTS OF ERROR-EXITS *****',
2680 $ '**' )
2681 *
2682 * End of ZCHKE.
2683 *
2684 END
2685 SUBROUTINE ZMAKE( TYPE, UPLO, DIAG, M, N, A, NMAX, AA, LDA, KL,
2686 $ KU, RESET, TRANSL )
2687 *
2688 * Generates values for an M by N matrix A within the bandwidth
2689 * defined by KL and KU.
2690 * Stores the values in the array AA in the data structure required
2691 * by the routine, with unwanted elements set to rogue value.
2692 *
2693 * TYPE is 'GE', 'GB', 'HE', 'HB', 'HP', 'TR', 'TB' OR 'TP'.
2694 *
2695 * Auxiliary routine for test program for Level 2 Blas.
2696 *
2697 * -- Written on 10-August-1987.
2698 * Richard Hanson, Sandia National Labs.
2699 * Jeremy Du Croz, NAG Central Office.
2700 *
2701 * .. Parameters ..
2702 COMPLEX*16 ZERO, ONE
2703 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ),
2704 $ ONE = ( 1.0D0, 0.0D0 ) )
2705 COMPLEX*16 ROGUE
2706 PARAMETER ( ROGUE = ( -1.0D10, 1.0D10 ) )
2707 DOUBLE PRECISION RZERO
2708 PARAMETER ( RZERO = 0.0D0 )
2709 DOUBLE PRECISION RROGUE
2710 PARAMETER ( RROGUE = -1.0D10 )
2711 * .. Scalar Arguments ..
2712 COMPLEX*16 TRANSL
2713 INTEGER KL, KU, LDA, M, N, NMAX
2714 LOGICAL RESET
2715 CHARACTER*1 DIAG, UPLO
2716 CHARACTER*2 TYPE
2717 * .. Array Arguments ..
2718 COMPLEX*16 A( NMAX, * ), AA( * )
2719 * .. Local Scalars ..
2720 INTEGER I, I1, I2, I3, IBEG, IEND, IOFF, J, JJ, KK
2721 LOGICAL GEN, LOWER, SYM, TRI, UNIT, UPPER
2722 * .. External Functions ..
2723 COMPLEX*16 ZBEG
2724 EXTERNAL ZBEG
2725 * .. Intrinsic Functions ..
2726 INTRINSIC DBLE, DCMPLX, DCONJG, MAX, MIN
2727 * .. Executable Statements ..
2728 GEN = TYPE( 1: 1 ).EQ.'G'
2729 SYM = TYPE( 1: 1 ).EQ.'H'
2730 TRI = TYPE( 1: 1 ).EQ.'T'
2731 UPPER = ( SYM.OR.TRI ).AND.UPLO.EQ.'U'
2732 LOWER = ( SYM.OR.TRI ).AND.UPLO.EQ.'L'
2733 UNIT = TRI.AND.DIAG.EQ.'U'
2734 *
2735 * Generate data in array A.
2736 *
2737 DO 20 J = 1, N
2738 DO 10 I = 1, M
2739 IF( GEN.OR.( UPPER.AND.I.LE.J ).OR.( LOWER.AND.I.GE.J ) )
2740 $ THEN
2741 IF( ( I.LE.J.AND.J - I.LE.KU ).OR.
2742 $ ( I.GE.J.AND.I - J.LE.KL ) )THEN
2743 A( I, J ) = ZBEG( RESET ) + TRANSL
2744 ELSE
2745 A( I, J ) = ZERO
2746 END IF
2747 IF( I.NE.J )THEN
2748 IF( SYM )THEN
2749 A( J, I ) = DCONJG( A( I, J ) )
2750 ELSE IF( TRI )THEN
2751 A( J, I ) = ZERO
2752 END IF
2753 END IF
2754 END IF
2755 10 CONTINUE
2756 IF( SYM )
2757 $ A( J, J ) = DCMPLX( DBLE( A( J, J ) ), RZERO )
2758 IF( TRI )
2759 $ A( J, J ) = A( J, J ) + ONE
2760 IF( UNIT )
2761 $ A( J, J ) = ONE
2762 20 CONTINUE
2763 *
2764 * Store elements in array AS in data structure required by routine.
2765 *
2766 IF( TYPE.EQ.'GE' )THEN
2767 DO 50 J = 1, N
2768 DO 30 I = 1, M
2769 AA( I + ( J - 1 )*LDA ) = A( I, J )
2770 30 CONTINUE
2771 DO 40 I = M + 1, LDA
2772 AA( I + ( J - 1 )*LDA ) = ROGUE
2773 40 CONTINUE
2774 50 CONTINUE
2775 ELSE IF( TYPE.EQ.'GB' )THEN
2776 DO 90 J = 1, N
2777 DO 60 I1 = 1, KU + 1 - J
2778 AA( I1 + ( J - 1 )*LDA ) = ROGUE
2779 60 CONTINUE
2780 DO 70 I2 = I1, MIN( KL + KU + 1, KU + 1 + M - J )
2781 AA( I2 + ( J - 1 )*LDA ) = A( I2 + J - KU - 1, J )
2782 70 CONTINUE
2783 DO 80 I3 = I2, LDA
2784 AA( I3 + ( J - 1 )*LDA ) = ROGUE
2785 80 CONTINUE
2786 90 CONTINUE
2787 ELSE IF( TYPE.EQ.'HE'.OR.TYPE.EQ.'TR' )THEN
2788 DO 130 J = 1, N
2789 IF( UPPER )THEN
2790 IBEG = 1
2791 IF( UNIT )THEN
2792 IEND = J - 1
2793 ELSE
2794 IEND = J
2795 END IF
2796 ELSE
2797 IF( UNIT )THEN
2798 IBEG = J + 1
2799 ELSE
2800 IBEG = J
2801 END IF
2802 IEND = N
2803 END IF
2804 DO 100 I = 1, IBEG - 1
2805 AA( I + ( J - 1 )*LDA ) = ROGUE
2806 100 CONTINUE
2807 DO 110 I = IBEG, IEND
2808 AA( I + ( J - 1 )*LDA ) = A( I, J )
2809 110 CONTINUE
2810 DO 120 I = IEND + 1, LDA
2811 AA( I + ( J - 1 )*LDA ) = ROGUE
2812 120 CONTINUE
2813 IF( SYM )THEN
2814 JJ = J + ( J - 1 )*LDA
2815 AA( JJ ) = DCMPLX( DBLE( AA( JJ ) ), RROGUE )
2816 END IF
2817 130 CONTINUE
2818 ELSE IF( TYPE.EQ.'HB'.OR.TYPE.EQ.'TB' )THEN
2819 DO 170 J = 1, N
2820 IF( UPPER )THEN
2821 KK = KL + 1
2822 IBEG = MAX( 1, KL + 2 - J )
2823 IF( UNIT )THEN
2824 IEND = KL
2825 ELSE
2826 IEND = KL + 1
2827 END IF
2828 ELSE
2829 KK = 1
2830 IF( UNIT )THEN
2831 IBEG = 2
2832 ELSE
2833 IBEG = 1
2834 END IF
2835 IEND = MIN( KL + 1, 1 + M - J )
2836 END IF
2837 DO 140 I = 1, IBEG - 1
2838 AA( I + ( J - 1 )*LDA ) = ROGUE
2839 140 CONTINUE
2840 DO 150 I = IBEG, IEND
2841 AA( I + ( J - 1 )*LDA ) = A( I + J - KK, J )
2842 150 CONTINUE
2843 DO 160 I = IEND + 1, LDA
2844 AA( I + ( J - 1 )*LDA ) = ROGUE
2845 160 CONTINUE
2846 IF( SYM )THEN
2847 JJ = KK + ( J - 1 )*LDA
2848 AA( JJ ) = DCMPLX( DBLE( AA( JJ ) ), RROGUE )
2849 END IF
2850 170 CONTINUE
2851 ELSE IF( TYPE.EQ.'HP'.OR.TYPE.EQ.'TP' )THEN
2852 IOFF = 0
2853 DO 190 J = 1, N
2854 IF( UPPER )THEN
2855 IBEG = 1
2856 IEND = J
2857 ELSE
2858 IBEG = J
2859 IEND = N
2860 END IF
2861 DO 180 I = IBEG, IEND
2862 IOFF = IOFF + 1
2863 AA( IOFF ) = A( I, J )
2864 IF( I.EQ.J )THEN
2865 IF( UNIT )
2866 $ AA( IOFF ) = ROGUE
2867 IF( SYM )
2868 $ AA( IOFF ) = DCMPLX( DBLE( AA( IOFF ) ), RROGUE )
2869 END IF
2870 180 CONTINUE
2871 190 CONTINUE
2872 END IF
2873 RETURN
2874 *
2875 * End of ZMAKE.
2876 *
2877 END
2878 SUBROUTINE ZMVCH( TRANS, M, N, ALPHA, A, NMAX, X, INCX, BETA, Y,
2879 $ INCY, YT, G, YY, EPS, ERR, FATAL, NOUT, MV )
2880 *
2881 * Checks the results of the computational tests.
2882 *
2883 * Auxiliary routine for test program for Level 2 Blas.
2884 *
2885 * -- Written on 10-August-1987.
2886 * Richard Hanson, Sandia National Labs.
2887 * Jeremy Du Croz, NAG Central Office.
2888 *
2889 * .. Parameters ..
2890 COMPLEX*16 ZERO
2891 PARAMETER ( ZERO = ( 0.0D0, 0.0D0 ) )
2892 DOUBLE PRECISION RZERO, RONE
2893 PARAMETER ( RZERO = 0.0D0, RONE = 1.0D0 )
2894 * .. Scalar Arguments ..
2895 COMPLEX*16 ALPHA, BETA
2896 DOUBLE PRECISION EPS, ERR
2897 INTEGER INCX, INCY, M, N, NMAX, NOUT
2898 LOGICAL FATAL, MV
2899 CHARACTER*1 TRANS
2900 * .. Array Arguments ..
2901 COMPLEX*16 A( NMAX, * ), X( * ), Y( * ), YT( * ), YY( * )
2902 DOUBLE PRECISION G( * )
2903 * .. Local Scalars ..
2904 COMPLEX*16 C
2905 DOUBLE PRECISION ERRI
2906 INTEGER I, INCXL, INCYL, IY, J, JX, KX, KY, ML, NL
2907 LOGICAL CTRAN, TRAN
2908 * .. Intrinsic Functions ..
2909 INTRINSIC ABS, DBLE, DCONJG, DIMAG, MAX, SQRT
2910 * .. Statement Functions ..
2911 DOUBLE PRECISION ABS1
2912 * .. Statement Function definitions ..
2913 ABS1( C ) = ABS( DBLE( C ) ) + ABS( DIMAG( C ) )
2914 * .. Executable Statements ..
2915 TRAN = TRANS.EQ.'T'
2916 CTRAN = TRANS.EQ.'C'
2917 IF( TRAN.OR.CTRAN )THEN
2918 ML = N
2919 NL = M
2920 ELSE
2921 ML = M
2922 NL = N
2923 END IF
2924 IF( INCX.LT.0 )THEN
2925 KX = NL
2926 INCXL = -1
2927 ELSE
2928 KX = 1
2929 INCXL = 1
2930 END IF
2931 IF( INCY.LT.0 )THEN
2932 KY = ML
2933 INCYL = -1
2934 ELSE
2935 KY = 1
2936 INCYL = 1
2937 END IF
2938 *
2939 * Compute expected result in YT using data in A, X and Y.
2940 * Compute gauges in G.
2941 *
2942 IY = KY
2943 DO 40 I = 1, ML
2944 YT( IY ) = ZERO
2945 G( IY ) = RZERO
2946 JX = KX
2947 IF( TRAN )THEN
2948 DO 10 J = 1, NL
2949 YT( IY ) = YT( IY ) + A( J, I )*X( JX )
2950 G( IY ) = G( IY ) + ABS1( A( J, I ) )*ABS1( X( JX ) )
2951 JX = JX + INCXL
2952 10 CONTINUE
2953 ELSE IF( CTRAN )THEN
2954 DO 20 J = 1, NL
2955 YT( IY ) = YT( IY ) + DCONJG( A( J, I ) )*X( JX )
2956 G( IY ) = G( IY ) + ABS1( A( J, I ) )*ABS1( X( JX ) )
2957 JX = JX + INCXL
2958 20 CONTINUE
2959 ELSE
2960 DO 30 J = 1, NL
2961 YT( IY ) = YT( IY ) + A( I, J )*X( JX )
2962 G( IY ) = G( IY ) + ABS1( A( I, J ) )*ABS1( X( JX ) )
2963 JX = JX + INCXL
2964 30 CONTINUE
2965 END IF
2966 YT( IY ) = ALPHA*YT( IY ) + BETA*Y( IY )
2967 G( IY ) = ABS1( ALPHA )*G( IY ) + ABS1( BETA )*ABS1( Y( IY ) )
2968 IY = IY + INCYL
2969 40 CONTINUE
2970 *
2971 * Compute the error ratio for this result.
2972 *
2973 ERR = ZERO
2974 DO 50 I = 1, ML
2975 ERRI = ABS( YT( I ) - YY( 1 + ( I - 1 )*ABS( INCY ) ) )/EPS
2976 IF( G( I ).NE.RZERO )
2977 $ ERRI = ERRI/G( I )
2978 ERR = MAX( ERR, ERRI )
2979 IF( ERR*SQRT( EPS ).GE.RONE )
2980 $ GO TO 60
2981 50 CONTINUE
2982 * If the loop completes, all results are at least half accurate.
2983 GO TO 80
2984 *
2985 * Report fatal error.
2986 *
2987 60 FATAL = .TRUE.
2988 WRITE( NOUT, FMT = 9999 )
2989 DO 70 I = 1, ML
2990 IF( MV )THEN
2991 WRITE( NOUT, FMT = 9998 )I, YT( I ),
2992 $ YY( 1 + ( I - 1 )*ABS( INCY ) )
2993 ELSE
2994 WRITE( NOUT, FMT = 9998 )I,
2995 $ YY( 1 + ( I - 1 )*ABS( INCY ) ), YT( I )
2996 END IF
2997 70 CONTINUE
2998 *
2999 80 CONTINUE
3000 RETURN
3001 *
3002 9999 FORMAT( ' ******* FATAL ERROR - COMPUTED RESULT IS LESS THAN HAL',
3003 $ 'F ACCURATE *******', /' EXPECTED RE',
3004 $ 'SULT COMPUTED RESULT' )
3005 9998 FORMAT( 1X, I7, 2( ' (', G15.6, ',', G15.6, ')' ) )
3006 *
3007 * End of ZMVCH.
3008 *
3009 END
3010 LOGICAL FUNCTION LZE( RI, RJ, LR )
3011 *
3012 * Tests if two arrays are identical.
3013 *
3014 * Auxiliary routine for test program for Level 2 Blas.
3015 *
3016 * -- Written on 10-August-1987.
3017 * Richard Hanson, Sandia National Labs.
3018 * Jeremy Du Croz, NAG Central Office.
3019 *
3020 * .. Scalar Arguments ..
3021 INTEGER LR
3022 * .. Array Arguments ..
3023 COMPLEX*16 RI( * ), RJ( * )
3024 * .. Local Scalars ..
3025 INTEGER I
3026 * .. Executable Statements ..
3027 DO 10 I = 1, LR
3028 IF( RI( I ).NE.RJ( I ) )
3029 $ GO TO 20
3030 10 CONTINUE
3031 LZE = .TRUE.
3032 GO TO 30
3033 20 CONTINUE
3034 LZE = .FALSE.
3035 30 RETURN
3036 *
3037 * End of LZE.
3038 *
3039 END
3040 LOGICAL FUNCTION LZERES( TYPE, UPLO, M, N, AA, AS, LDA )
3041 *
3042 * Tests if selected elements in two arrays are equal.
3043 *
3044 * TYPE is 'GE', 'HE' or 'HP'.
3045 *
3046 * Auxiliary routine for test program for Level 2 Blas.
3047 *
3048 * -- Written on 10-August-1987.
3049 * Richard Hanson, Sandia National Labs.
3050 * Jeremy Du Croz, NAG Central Office.
3051 *
3052 * .. Scalar Arguments ..
3053 INTEGER LDA, M, N
3054 CHARACTER*1 UPLO
3055 CHARACTER*2 TYPE
3056 * .. Array Arguments ..
3057 COMPLEX*16 AA( LDA, * ), AS( LDA, * )
3058 * .. Local Scalars ..
3059 INTEGER I, IBEG, IEND, J
3060 LOGICAL UPPER
3061 * .. Executable Statements ..
3062 UPPER = UPLO.EQ.'U'
3063 IF( TYPE.EQ.'GE' )THEN
3064 DO 20 J = 1, N
3065 DO 10 I = M + 1, LDA
3066 IF( AA( I, J ).NE.AS( I, J ) )
3067 $ GO TO 70
3068 10 CONTINUE
3069 20 CONTINUE
3070 ELSE IF( TYPE.EQ.'HE' )THEN
3071 DO 50 J = 1, N
3072 IF( UPPER )THEN
3073 IBEG = 1
3074 IEND = J
3075 ELSE
3076 IBEG = J
3077 IEND = N
3078 END IF
3079 DO 30 I = 1, IBEG - 1
3080 IF( AA( I, J ).NE.AS( I, J ) )
3081 $ GO TO 70
3082 30 CONTINUE
3083 DO 40 I = IEND + 1, LDA
3084 IF( AA( I, J ).NE.AS( I, J ) )
3085 $ GO TO 70
3086 40 CONTINUE
3087 50 CONTINUE
3088 END IF
3089 *
3090 60 CONTINUE
3091 LZERES = .TRUE.
3092 GO TO 80
3093 70 CONTINUE
3094 LZERES = .FALSE.
3095 80 RETURN
3096 *
3097 * End of LZERES.
3098 *
3099 END
3100 COMPLEX*16 FUNCTION ZBEG( RESET )
3101 *
3102 * Generates complex numbers as pairs of random numbers uniformly
3103 * distributed between -0.5 and 0.5.
3104 *
3105 * Auxiliary routine for test program for Level 2 Blas.
3106 *
3107 * -- Written on 10-August-1987.
3108 * Richard Hanson, Sandia National Labs.
3109 * Jeremy Du Croz, NAG Central Office.
3110 *
3111 * .. Scalar Arguments ..
3112 LOGICAL RESET
3113 * .. Local Scalars ..
3114 INTEGER I, IC, J, MI, MJ
3115 * .. Save statement ..
3116 SAVE I, IC, J, MI, MJ
3117 * .. Intrinsic Functions ..
3118 INTRINSIC DCMPLX
3119 * .. Executable Statements ..
3120 IF( RESET )THEN
3121 * Initialize local variables.
3122 MI = 891
3123 MJ = 457
3124 I = 7
3125 J = 7
3126 IC = 0
3127 RESET = .FALSE.
3128 END IF
3129 *
3130 * The sequence of values of I or J is bounded between 1 and 999.
3131 * If initial I or J = 1,2,3,6,7 or 9, the period will be 50.
3132 * If initial I or J = 4 or 8, the period will be 25.
3133 * If initial I or J = 5, the period will be 10.
3134 * IC is used to break up the period by skipping 1 value of I or J
3135 * in 6.
3136 *
3137 IC = IC + 1
3138 10 I = I*MI
3139 J = J*MJ
3140 I = I - 1000*( I/1000 )
3141 J = J - 1000*( J/1000 )
3142 IF( IC.GE.5 )THEN
3143 IC = 0
3144 GO TO 10
3145 END IF
3146 ZBEG = DCMPLX( ( I - 500 )/1001.0D0, ( J - 500 )/1001.0D0 )
3147 RETURN
3148 *
3149 * End of ZBEG.
3150 *
3151 END
3152 DOUBLE PRECISION FUNCTION DDIFF( X, Y )
3153 *
3154 * Auxiliary routine for test program for Level 2 Blas.
3155 *
3156 * -- Written on 10-August-1987.
3157 * Richard Hanson, Sandia National Labs.
3158 *
3159 * .. Scalar Arguments ..
3160 DOUBLE PRECISION X, Y
3161 * .. Executable Statements ..
3162 DDIFF = X - Y
3163 RETURN
3164 *
3165 * End of DDIFF.
3166 *
3167 END
3168 SUBROUTINE CHKXER( SRNAMT, INFOT, NOUT, LERR, OK )
3169 *
3170 * Tests whether XERBLA has detected an error when it should.
3171 *
3172 * Auxiliary routine for test program for Level 2 Blas.
3173 *
3174 * -- Written on 10-August-1987.
3175 * Richard Hanson, Sandia National Labs.
3176 * Jeremy Du Croz, NAG Central Office.
3177 *
3178 * .. Scalar Arguments ..
3179 INTEGER INFOT, NOUT
3180 LOGICAL LERR, OK
3181 CHARACTER*6 SRNAMT
3182 * .. Executable Statements ..
3183 IF( .NOT.LERR )THEN
3184 WRITE( NOUT, FMT = 9999 )INFOT, SRNAMT
3185 OK = .FALSE.
3186 END IF
3187 LERR = .FALSE.
3188 RETURN
3189 *
3190 9999 FORMAT( ' ***** ILLEGAL VALUE OF PARAMETER NUMBER ', I2, ' NOT D',
3191 $ 'ETECTED BY ', A6, ' *****' )
3192 *
3193 * End of CHKXER.
3194 *
3195 END
3196 SUBROUTINE XERBLA( SRNAME, INFO )
3197 *
3198 * This is a special version of XERBLA to be used only as part of
3199 * the test program for testing error exits from the Level 2 BLAS
3200 * routines.
3201 *
3202 * XERBLA is an error handler for the Level 2 BLAS routines.
3203 *
3204 * It is called by the Level 2 BLAS routines if an input parameter is
3205 * invalid.
3206 *
3207 * Auxiliary routine for test program for Level 2 Blas.
3208 *
3209 * -- Written on 10-August-1987.
3210 * Richard Hanson, Sandia National Labs.
3211 * Jeremy Du Croz, NAG Central Office.
3212 *
3213 * .. Scalar Arguments ..
3214 INTEGER INFO
3215 CHARACTER*6 SRNAME
3216 * .. Scalars in Common ..
3217 INTEGER INFOT, NOUT
3218 LOGICAL LERR, OK
3219 CHARACTER*6 SRNAMT
3220 * .. Common blocks ..
3221 COMMON /INFOC/INFOT, NOUT, OK, LERR
3222 COMMON /SRNAMC/SRNAMT
3223 * .. Executable Statements ..
3224 LERR = .TRUE.
3225 IF( INFO.NE.INFOT )THEN
3226 IF( INFOT.NE.0 )THEN
3227 WRITE( NOUT, FMT = 9999 )INFO, INFOT
3228 ELSE
3229 WRITE( NOUT, FMT = 9997 )INFO
3230 END IF
3231 OK = .FALSE.
3232 END IF
3233 IF( SRNAME.NE.SRNAMT )THEN
3234 WRITE( NOUT, FMT = 9998 )SRNAME, SRNAMT
3235 OK = .FALSE.
3236 END IF
3237 RETURN
3238 *
3239 9999 FORMAT( ' ******* XERBLA WAS CALLED WITH INFO = ', I6, ' INSTEAD',
3240 $ ' OF ', I2, ' *******' )
3241 9998 FORMAT( ' ******* XERBLA WAS CALLED WITH SRNAME = ', A6, ' INSTE',
3242 $ 'AD OF ', A6, ' *******' )
3243 9997 FORMAT( ' ******* XERBLA WAS CALLED WITH INFO = ', I6,
3244 $ ' *******' )
3245 *
3246 * End of XERBLA
3247 *
3248 END
3249