LAPACK  3.9.0
LAPACK: Linear Algebra PACKage
zerrst.f
Go to the documentation of this file.
1 *> \brief \b ZERRST
2 *
3 * @precisions fortran z -> c
4 *
5 * =========== DOCUMENTATION ===========
6 *
7 * Online html documentation available at
8 * http://www.netlib.org/lapack/explore-html/
9 *
10 * Definition:
11 * ===========
12 *
13 * SUBROUTINE ZERRST( PATH, NUNIT )
14 *
15 * .. Scalar Arguments ..
16 * CHARACTER*3 PATH
17 * INTEGER NUNIT
18 * ..
19 *
20 *
21 *> \par Purpose:
22 * =============
23 *>
24 *> \verbatim
25 *>
26 *> ZERRST tests the error exits for ZHETRD, ZUNGTR, CUNMTR, ZHPTRD,
27 *> ZUNGTR, ZUPMTR, ZSTEQR, CSTEIN, ZPTEQR, ZHBTRD,
28 *> ZHEEV, CHEEVX, CHEEVD, ZHBEV, CHBEVX, CHBEVD,
29 *> ZHPEV, CHPEVX, CHPEVD, and ZSTEDC.
30 *> ZHEEVD_2STAGE, ZHEEVR_2STAGE, ZHEEVX_2STAGE,
31 *> ZHEEV_2STAGE, ZHBEV_2STAGE, ZHBEVD_2STAGE,
32 *> ZHBEVX_2STAGE, ZHETRD_2STAGE
33 *> \endverbatim
34 *
35 * Arguments:
36 * ==========
37 *
38 *> \param[in] PATH
39 *> \verbatim
40 *> PATH is CHARACTER*3
41 *> The LAPACK path name for the routines to be tested.
42 *> \endverbatim
43 *>
44 *> \param[in] NUNIT
45 *> \verbatim
46 *> NUNIT is INTEGER
47 *> The unit number for output.
48 *> \endverbatim
49 *
50 * Authors:
51 * ========
52 *
53 *> \author Univ. of Tennessee
54 *> \author Univ. of California Berkeley
55 *> \author Univ. of Colorado Denver
56 *> \author NAG Ltd.
57 *
58 *> \date June 2017
59 *
60 *> \ingroup complex16_eig
61 *
62 * =====================================================================
63  SUBROUTINE zerrst( PATH, NUNIT )
64 *
65 * -- LAPACK test routine (version 3.7.1) --
66 * -- LAPACK is a software package provided by Univ. of Tennessee, --
67 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
68 * June 2017
69 *
70 * .. Scalar Arguments ..
71  CHARACTER*3 PATH
72  INTEGER NUNIT
73 * ..
74 *
75 * =====================================================================
76 *
77 * .. Parameters ..
78  INTEGER NMAX, LIW, LW
79  parameter( nmax = 3, liw = 12*nmax, lw = 20*nmax )
80 * ..
81 * .. Local Scalars ..
82  CHARACTER*2 C2
83  INTEGER I, INFO, J, M, N, NT
84 * ..
85 * .. Local Arrays ..
86  INTEGER I1( NMAX ), I2( NMAX ), I3( NMAX ), IW( LIW )
87  DOUBLE PRECISION D( NMAX ), E( NMAX ), R( LW ), RW( LW ),
88  $ X( NMAX )
89  COMPLEX*16 A( NMAX, NMAX ), C( NMAX, NMAX ),
90  $ Q( NMAX, NMAX ), TAU( NMAX ), W( LW ),
91  $ Z( NMAX, NMAX )
92 * ..
93 * .. External Functions ..
94  LOGICAL LSAMEN
95  EXTERNAL lsamen
96 * ..
97 * .. External Subroutines ..
98  EXTERNAL chkxer, zhbev, zhbevd, zhbevx, zhbtrd, zheev,
99  $ zheevd, zheevr, zheevx, zhetrd, zhpev, zhpevd,
100  $ zhpevx, zhptrd, zpteqr, zstedc, zstein, zsteqr,
101  $ zungtr, zunmtr, zupgtr, zupmtr,
102  $ zheevd_2stage, zheevr_2stage, zheevx_2stage,
103  $ zheev_2stage, zhbev_2stage, zhbevd_2stage,
104  $ zhbevx_2stage, zhetrd_2stage
105 * ..
106 * .. Scalars in Common ..
107  LOGICAL LERR, OK
108  CHARACTER*32 SRNAMT
109  INTEGER INFOT, NOUT
110 * ..
111 * .. Common blocks ..
112  COMMON / infoc / infot, nout, ok, lerr
113  COMMON / srnamc / srnamt
114 * ..
115 * .. Intrinsic Functions ..
116  INTRINSIC dble
117 * ..
118 * .. Executable Statements ..
119 *
120  nout = nunit
121  WRITE( nout, fmt = * )
122  c2 = path( 2: 3 )
123 *
124 * Set the variables to innocuous values.
125 *
126  DO 20 j = 1, nmax
127  DO 10 i = 1, nmax
128  a( i, j ) = 1.d0 / dble( i+j )
129  10 CONTINUE
130  20 CONTINUE
131  DO 30 j = 1, nmax
132  d( j ) = dble( j )
133  e( j ) = 0.0d0
134  i1( j ) = j
135  i2( j ) = j
136  tau( j ) = 1.d0
137  30 CONTINUE
138  ok = .true.
139  nt = 0
140 *
141 * Test error exits for the ST path.
142 *
143  IF( lsamen( 2, c2, 'ST' ) ) THEN
144 *
145 * ZHETRD
146 *
147  srnamt = 'ZHETRD'
148  infot = 1
149  CALL zhetrd( '/', 0, a, 1, d, e, tau, w, 1, info )
150  CALL chkxer( 'ZHETRD', infot, nout, lerr, ok )
151  infot = 2
152  CALL zhetrd( 'U', -1, a, 1, d, e, tau, w, 1, info )
153  CALL chkxer( 'ZHETRD', infot, nout, lerr, ok )
154  infot = 4
155  CALL zhetrd( 'U', 2, a, 1, d, e, tau, w, 1, info )
156  CALL chkxer( 'ZHETRD', infot, nout, lerr, ok )
157  infot = 9
158  CALL zhetrd( 'U', 0, a, 1, d, e, tau, w, 0, info )
159  CALL chkxer( 'ZHETRD', infot, nout, lerr, ok )
160  nt = nt + 4
161 *
162 * ZHETRD_2STAGE
163 *
164  srnamt = 'ZHETRD_2STAGE'
165  infot = 1
166  CALL zhetrd_2stage( '/', 'U', 0, a, 1, d, e, tau,
167  $ c, 1, w, 1, info )
168  CALL chkxer( 'ZHETRD_2STAGE', infot, nout, lerr, ok )
169  infot = 1
170  CALL zhetrd_2stage( 'H', 'U', 0, a, 1, d, e, tau,
171  $ c, 1, w, 1, info )
172  CALL chkxer( 'ZHETRD_2STAGE', infot, nout, lerr, ok )
173  infot = 2
174  CALL zhetrd_2stage( 'N', '/', 0, a, 1, d, e, tau,
175  $ c, 1, w, 1, info )
176  CALL chkxer( 'ZHETRD_2STAGE', infot, nout, lerr, ok )
177  infot = 3
178  CALL zhetrd_2stage( 'N', 'U', -1, a, 1, d, e, tau,
179  $ c, 1, w, 1, info )
180  CALL chkxer( 'ZHETRD_2STAGE', infot, nout, lerr, ok )
181  infot = 5
182  CALL zhetrd_2stage( 'N', 'U', 2, a, 1, d, e, tau,
183  $ c, 1, w, 1, info )
184  CALL chkxer( 'ZHETRD_2STAGE', infot, nout, lerr, ok )
185  infot = 10
186  CALL zhetrd_2stage( 'N', 'U', 0, a, 1, d, e, tau,
187  $ c, 0, w, 1, info )
188  CALL chkxer( 'ZHETRD_2STAGE', infot, nout, lerr, ok )
189  infot = 12
190  CALL zhetrd_2stage( 'N', 'U', 0, a, 1, d, e, tau,
191  $ c, 1, w, 0, info )
192  CALL chkxer( 'ZHETRD_2STAGE', infot, nout, lerr, ok )
193  nt = nt + 7
194 *
195 * ZHETRD_HE2HB
196 *
197  srnamt = 'ZHETRD_HE2HB'
198  infot = 1
199  CALL zhetrd_he2hb( '/', 0, 0, a, 1, c, 1, tau, w, 1, info )
200  CALL chkxer( 'ZHETRD_HE2HB', infot, nout, lerr, ok )
201  infot = 2
202  CALL zhetrd_he2hb( 'U', -1, 0, a, 1, c, 1, tau, w, 1, info )
203  CALL chkxer( 'ZHETRD_HE2HB', infot, nout, lerr, ok )
204  infot = 3
205  CALL zhetrd_he2hb( 'U', 0, -1, a, 1, c, 1, tau, w, 1, info )
206  CALL chkxer( 'ZHETRD_HE2HB', infot, nout, lerr, ok )
207  infot = 5
208  CALL zhetrd_he2hb( 'U', 2, 0, a, 1, c, 1, tau, w, 1, info )
209  CALL chkxer( 'ZHETRD_HE2HB', infot, nout, lerr, ok )
210  infot = 7
211  CALL zhetrd_he2hb( 'U', 0, 2, a, 1, c, 1, tau, w, 1, info )
212  CALL chkxer( 'ZHETRD_HE2HB', infot, nout, lerr, ok )
213  infot = 10
214  CALL zhetrd_he2hb( 'U', 0, 0, a, 1, c, 1, tau, w, 0, info )
215  CALL chkxer( 'ZHETRD_HE2HB', infot, nout, lerr, ok )
216  nt = nt + 6
217 *
218 * ZHETRD_HB2ST
219 *
220  srnamt = 'ZHETRD_HB2ST'
221  infot = 1
222  CALL zhetrd_hb2st( '/', 'N', 'U', 0, 0, a, 1, d, e,
223  $ c, 1, w, 1, info )
224  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
225  infot = 2
226  CALL zhetrd_hb2st( 'Y', '/', 'U', 0, 0, a, 1, d, e,
227  $ c, 1, w, 1, info )
228  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
229  infot = 2
230  CALL zhetrd_hb2st( 'Y', 'H', 'U', 0, 0, a, 1, d, e,
231  $ c, 1, w, 1, info )
232  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
233  infot = 3
234  CALL zhetrd_hb2st( 'Y', 'N', '/', 0, 0, a, 1, d, e,
235  $ c, 1, w, 1, info )
236  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
237  infot = 4
238  CALL zhetrd_hb2st( 'Y', 'N', 'U', -1, 0, a, 1, d, e,
239  $ c, 1, w, 1, info )
240  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
241  infot = 5
242  CALL zhetrd_hb2st( 'Y', 'N', 'U', 0, -1, a, 1, d, e,
243  $ c, 1, w, 1, info )
244  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
245  infot = 7
246  CALL zhetrd_hb2st( 'Y', 'N', 'U', 0, 1, a, 1, d, e,
247  $ c, 1, w, 1, info )
248  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
249  infot = 11
250  CALL zhetrd_hb2st( 'Y', 'N', 'U', 0, 0, a, 1, d, e,
251  $ c, 0, w, 1, info )
252  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
253  infot = 13
254  CALL zhetrd_hb2st( 'Y', 'N', 'U', 0, 0, a, 1, d, e,
255  $ c, 1, w, 0, info )
256  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
257  nt = nt + 9
258 *
259 * ZUNGTR
260 *
261  srnamt = 'ZUNGTR'
262  infot = 1
263  CALL zungtr( '/', 0, a, 1, tau, w, 1, info )
264  CALL chkxer( 'ZUNGTR', infot, nout, lerr, ok )
265  infot = 2
266  CALL zungtr( 'U', -1, a, 1, tau, w, 1, info )
267  CALL chkxer( 'ZUNGTR', infot, nout, lerr, ok )
268  infot = 4
269  CALL zungtr( 'U', 2, a, 1, tau, w, 1, info )
270  CALL chkxer( 'ZUNGTR', infot, nout, lerr, ok )
271  infot = 7
272  CALL zungtr( 'U', 3, a, 3, tau, w, 1, info )
273  CALL chkxer( 'ZUNGTR', infot, nout, lerr, ok )
274  nt = nt + 4
275 *
276 * ZUNMTR
277 *
278  srnamt = 'ZUNMTR'
279  infot = 1
280  CALL zunmtr( '/', 'U', 'N', 0, 0, a, 1, tau, c, 1, w, 1, info )
281  CALL chkxer( 'ZUNMTR', infot, nout, lerr, ok )
282  infot = 2
283  CALL zunmtr( 'L', '/', 'N', 0, 0, a, 1, tau, c, 1, w, 1, info )
284  CALL chkxer( 'ZUNMTR', infot, nout, lerr, ok )
285  infot = 3
286  CALL zunmtr( 'L', 'U', '/', 0, 0, a, 1, tau, c, 1, w, 1, info )
287  CALL chkxer( 'ZUNMTR', infot, nout, lerr, ok )
288  infot = 4
289  CALL zunmtr( 'L', 'U', 'N', -1, 0, a, 1, tau, c, 1, w, 1,
290  $ info )
291  CALL chkxer( 'ZUNMTR', infot, nout, lerr, ok )
292  infot = 5
293  CALL zunmtr( 'L', 'U', 'N', 0, -1, a, 1, tau, c, 1, w, 1,
294  $ info )
295  CALL chkxer( 'ZUNMTR', infot, nout, lerr, ok )
296  infot = 7
297  CALL zunmtr( 'L', 'U', 'N', 2, 0, a, 1, tau, c, 2, w, 1, info )
298  CALL chkxer( 'ZUNMTR', infot, nout, lerr, ok )
299  infot = 7
300  CALL zunmtr( 'R', 'U', 'N', 0, 2, a, 1, tau, c, 1, w, 1, info )
301  CALL chkxer( 'ZUNMTR', infot, nout, lerr, ok )
302  infot = 10
303  CALL zunmtr( 'L', 'U', 'N', 2, 0, a, 2, tau, c, 1, w, 1, info )
304  CALL chkxer( 'ZUNMTR', infot, nout, lerr, ok )
305  infot = 12
306  CALL zunmtr( 'L', 'U', 'N', 0, 2, a, 1, tau, c, 1, w, 1, info )
307  CALL chkxer( 'ZUNMTR', infot, nout, lerr, ok )
308  infot = 12
309  CALL zunmtr( 'R', 'U', 'N', 2, 0, a, 1, tau, c, 2, w, 1, info )
310  CALL chkxer( 'ZUNMTR', infot, nout, lerr, ok )
311  nt = nt + 10
312 *
313 * ZHPTRD
314 *
315  srnamt = 'ZHPTRD'
316  infot = 1
317  CALL zhptrd( '/', 0, a, d, e, tau, info )
318  CALL chkxer( 'ZHPTRD', infot, nout, lerr, ok )
319  infot = 2
320  CALL zhptrd( 'U', -1, a, d, e, tau, info )
321  CALL chkxer( 'ZHPTRD', infot, nout, lerr, ok )
322  nt = nt + 2
323 *
324 * ZUPGTR
325 *
326  srnamt = 'ZUPGTR'
327  infot = 1
328  CALL zupgtr( '/', 0, a, tau, z, 1, w, info )
329  CALL chkxer( 'ZUPGTR', infot, nout, lerr, ok )
330  infot = 2
331  CALL zupgtr( 'U', -1, a, tau, z, 1, w, info )
332  CALL chkxer( 'ZUPGTR', infot, nout, lerr, ok )
333  infot = 6
334  CALL zupgtr( 'U', 2, a, tau, z, 1, w, info )
335  CALL chkxer( 'ZUPGTR', infot, nout, lerr, ok )
336  nt = nt + 3
337 *
338 * ZUPMTR
339 *
340  srnamt = 'ZUPMTR'
341  infot = 1
342  CALL zupmtr( '/', 'U', 'N', 0, 0, a, tau, c, 1, w, info )
343  CALL chkxer( 'ZUPMTR', infot, nout, lerr, ok )
344  infot = 2
345  CALL zupmtr( 'L', '/', 'N', 0, 0, a, tau, c, 1, w, info )
346  CALL chkxer( 'ZUPMTR', infot, nout, lerr, ok )
347  infot = 3
348  CALL zupmtr( 'L', 'U', '/', 0, 0, a, tau, c, 1, w, info )
349  CALL chkxer( 'ZUPMTR', infot, nout, lerr, ok )
350  infot = 4
351  CALL zupmtr( 'L', 'U', 'N', -1, 0, a, tau, c, 1, w, info )
352  CALL chkxer( 'ZUPMTR', infot, nout, lerr, ok )
353  infot = 5
354  CALL zupmtr( 'L', 'U', 'N', 0, -1, a, tau, c, 1, w, info )
355  CALL chkxer( 'ZUPMTR', infot, nout, lerr, ok )
356  infot = 9
357  CALL zupmtr( 'L', 'U', 'N', 2, 0, a, tau, c, 1, w, info )
358  CALL chkxer( 'ZUPMTR', infot, nout, lerr, ok )
359  nt = nt + 6
360 *
361 * ZPTEQR
362 *
363  srnamt = 'ZPTEQR'
364  infot = 1
365  CALL zpteqr( '/', 0, d, e, z, 1, rw, info )
366  CALL chkxer( 'ZPTEQR', infot, nout, lerr, ok )
367  infot = 2
368  CALL zpteqr( 'N', -1, d, e, z, 1, rw, info )
369  CALL chkxer( 'ZPTEQR', infot, nout, lerr, ok )
370  infot = 6
371  CALL zpteqr( 'V', 2, d, e, z, 1, rw, info )
372  CALL chkxer( 'ZPTEQR', infot, nout, lerr, ok )
373  nt = nt + 3
374 *
375 * ZSTEIN
376 *
377  srnamt = 'ZSTEIN'
378  infot = 1
379  CALL zstein( -1, d, e, 0, x, i1, i2, z, 1, rw, iw, i3, info )
380  CALL chkxer( 'ZSTEIN', infot, nout, lerr, ok )
381  infot = 4
382  CALL zstein( 0, d, e, -1, x, i1, i2, z, 1, rw, iw, i3, info )
383  CALL chkxer( 'ZSTEIN', infot, nout, lerr, ok )
384  infot = 4
385  CALL zstein( 0, d, e, 1, x, i1, i2, z, 1, rw, iw, i3, info )
386  CALL chkxer( 'ZSTEIN', infot, nout, lerr, ok )
387  infot = 9
388  CALL zstein( 2, d, e, 0, x, i1, i2, z, 1, rw, iw, i3, info )
389  CALL chkxer( 'ZSTEIN', infot, nout, lerr, ok )
390  nt = nt + 4
391 *
392 * ZSTEQR
393 *
394  srnamt = 'ZSTEQR'
395  infot = 1
396  CALL zsteqr( '/', 0, d, e, z, 1, rw, info )
397  CALL chkxer( 'ZSTEQR', infot, nout, lerr, ok )
398  infot = 2
399  CALL zsteqr( 'N', -1, d, e, z, 1, rw, info )
400  CALL chkxer( 'ZSTEQR', infot, nout, lerr, ok )
401  infot = 6
402  CALL zsteqr( 'V', 2, d, e, z, 1, rw, info )
403  CALL chkxer( 'ZSTEQR', infot, nout, lerr, ok )
404  nt = nt + 3
405 *
406 * ZSTEDC
407 *
408  srnamt = 'ZSTEDC'
409  infot = 1
410  CALL zstedc( '/', 0, d, e, z, 1, w, 1, rw, 1, iw, 1, info )
411  CALL chkxer( 'ZSTEDC', infot, nout, lerr, ok )
412  infot = 2
413  CALL zstedc( 'N', -1, d, e, z, 1, w, 1, rw, 1, iw, 1, info )
414  CALL chkxer( 'ZSTEDC', infot, nout, lerr, ok )
415  infot = 6
416  CALL zstedc( 'V', 2, d, e, z, 1, w, 4, rw, 23, iw, 28, info )
417  CALL chkxer( 'ZSTEDC', infot, nout, lerr, ok )
418  infot = 8
419  CALL zstedc( 'N', 2, d, e, z, 1, w, 0, rw, 1, iw, 1, info )
420  CALL chkxer( 'ZSTEDC', infot, nout, lerr, ok )
421  infot = 8
422  CALL zstedc( 'V', 2, d, e, z, 2, w, 0, rw, 23, iw, 28, info )
423  CALL chkxer( 'ZSTEDC', infot, nout, lerr, ok )
424  infot = 10
425  CALL zstedc( 'N', 2, d, e, z, 1, w, 1, rw, 0, iw, 1, info )
426  CALL chkxer( 'ZSTEDC', infot, nout, lerr, ok )
427  infot = 10
428  CALL zstedc( 'I', 2, d, e, z, 2, w, 1, rw, 1, iw, 12, info )
429  CALL chkxer( 'ZSTEDC', infot, nout, lerr, ok )
430  infot = 10
431  CALL zstedc( 'V', 2, d, e, z, 2, w, 4, rw, 1, iw, 28, info )
432  CALL chkxer( 'ZSTEDC', infot, nout, lerr, ok )
433  infot = 12
434  CALL zstedc( 'N', 2, d, e, z, 1, w, 1, rw, 1, iw, 0, info )
435  CALL chkxer( 'ZSTEDC', infot, nout, lerr, ok )
436  infot = 12
437  CALL zstedc( 'I', 2, d, e, z, 2, w, 1, rw, 23, iw, 0, info )
438  CALL chkxer( 'ZSTEDC', infot, nout, lerr, ok )
439  infot = 12
440  CALL zstedc( 'V', 2, d, e, z, 2, w, 4, rw, 23, iw, 0, info )
441  CALL chkxer( 'ZSTEDC', infot, nout, lerr, ok )
442  nt = nt + 11
443 *
444 * ZHEEVD
445 *
446  srnamt = 'ZHEEVD'
447  infot = 1
448  CALL zheevd( '/', 'U', 0, a, 1, x, w, 1, rw, 1, iw, 1, info )
449  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
450  infot = 2
451  CALL zheevd( 'N', '/', 0, a, 1, x, w, 1, rw, 1, iw, 1, info )
452  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
453  infot = 3
454  CALL zheevd( 'N', 'U', -1, a, 1, x, w, 1, rw, 1, iw, 1, info )
455  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
456  infot = 5
457  CALL zheevd( 'N', 'U', 2, a, 1, x, w, 3, rw, 2, iw, 1, info )
458  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
459  infot = 8
460  CALL zheevd( 'N', 'U', 1, a, 1, x, w, 0, rw, 1, iw, 1, info )
461  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
462  infot = 8
463  CALL zheevd( 'N', 'U', 2, a, 2, x, w, 2, rw, 2, iw, 1, info )
464  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
465  infot = 8
466  CALL zheevd( 'V', 'U', 2, a, 2, x, w, 3, rw, 25, iw, 12, info )
467  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
468  infot = 10
469  CALL zheevd( 'N', 'U', 1, a, 1, x, w, 1, rw, 0, iw, 1, info )
470  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
471  infot = 10
472  CALL zheevd( 'N', 'U', 2, a, 2, x, w, 3, rw, 1, iw, 1, info )
473  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
474  infot = 10
475  CALL zheevd( 'V', 'U', 2, a, 2, x, w, 8, rw, 18, iw, 12, info )
476  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
477  infot = 12
478  CALL zheevd( 'N', 'U', 1, a, 1, x, w, 1, rw, 1, iw, 0, info )
479  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
480  infot = 12
481  CALL zheevd( 'V', 'U', 2, a, 2, x, w, 8, rw, 25, iw, 11, info )
482  CALL chkxer( 'ZHEEVD', infot, nout, lerr, ok )
483  nt = nt + 12
484 *
485 * ZHEEVD_2STAGE
486 *
487  srnamt = 'ZHEEVD_2STAGE'
488  infot = 1
489  CALL zheevd_2stage( '/', 'U', 0, a, 1, x, w, 1,
490  $ rw, 1, iw, 1, info )
491  CALL chkxer( 'ZHEEVD_2STAGE', infot, nout, lerr, ok )
492  infot = 1
493  CALL zheevd_2stage( 'V', 'U', 0, a, 1, x, w, 1,
494  $ rw, 1, iw, 1, info )
495  CALL chkxer( 'ZHEEVD_2STAGE', infot, nout, lerr, ok )
496  infot = 2
497  CALL zheevd_2stage( 'N', '/', 0, a, 1, x, w, 1,
498  $ rw, 1, iw, 1, info )
499  CALL chkxer( 'ZHEEVD_2STAGE', infot, nout, lerr, ok )
500  infot = 3
501  CALL zheevd_2stage( 'N', 'U', -1, a, 1, x, w, 1,
502  $ rw, 1, iw, 1, info )
503  CALL chkxer( 'ZHEEVD_2STAGE', infot, nout, lerr, ok )
504  infot = 5
505  CALL zheevd_2stage( 'N', 'U', 2, a, 1, x, w, 3,
506  $ rw, 2, iw, 1, info )
507  CALL chkxer( 'ZHEEVD_2STAGE', infot, nout, lerr, ok )
508  infot = 8
509  CALL zheevd_2stage( 'N', 'U', 1, a, 1, x, w, 0,
510  $ rw, 1, iw, 1, info )
511  CALL chkxer( 'ZHEEVD_2STAGE', infot, nout, lerr, ok )
512  infot = 8
513  CALL zheevd_2stage( 'N', 'U', 2, a, 2, x, w, 2,
514  $ rw, 2, iw, 1, info )
515  CALL chkxer( 'ZHEEVD_2STAGE', infot, nout, lerr, ok )
516 * INFOT = 8
517 * CALL ZHEEVD_2STAGE( 'V', 'U', 2, A, 2, X, W, 3,
518 * $ RW, 25, IW, 12, INFO )
519 * CALL CHKXER( 'ZHEEVD_2STAGE', INFOT, NOUT, LERR, OK )
520  infot = 10
521  CALL zheevd_2stage( 'N', 'U', 1, a, 1, x, w, 1,
522  $ rw, 0, iw, 1, info )
523  CALL chkxer( 'ZHEEVD_2STAGE', infot, nout, lerr, ok )
524  infot = 10
525  CALL zheevd_2stage( 'N', 'U', 2, a, 2, x, w, 25,
526  $ rw, 1, iw, 1, info )
527  CALL chkxer( 'ZHEEVD_2STAGE', infot, nout, lerr, ok )
528 * INFOT = 10
529 * CALL ZHEEVD_2STAGE( 'V', 'U', 2, A, 2, X, W, 8,
530 * $ RW, 18, IW, 12, INFO )
531 * CALL CHKXER( 'ZHEEVD_2STAGE', INFOT, NOUT, LERR, OK )
532  infot = 12
533  CALL zheevd_2stage( 'N', 'U', 1, a, 1, x, w, 1,
534  $ rw, 1, iw, 0, info )
535  CALL chkxer( 'ZHEEVD_2STAGE', infot, nout, lerr, ok )
536  infot = 12
537 * CALL ZHEEVD_2STAGE( 'V', 'U', 2, A, 2, X, W, 8,
538 * $ RW, 25, IW, 11, INFO )
539 * CALL CHKXER( 'ZHEEVD_2STAGE', INFOT, NOUT, LERR, OK )
540  nt = nt + 10
541 *
542 * ZHEEV
543 *
544  srnamt = 'ZHEEV '
545  infot = 1
546  CALL zheev( '/', 'U', 0, a, 1, x, w, 1, rw, info )
547  CALL chkxer( 'ZHEEV ', infot, nout, lerr, ok )
548  infot = 2
549  CALL zheev( 'N', '/', 0, a, 1, x, w, 1, rw, info )
550  CALL chkxer( 'ZHEEV ', infot, nout, lerr, ok )
551  infot = 3
552  CALL zheev( 'N', 'U', -1, a, 1, x, w, 1, rw, info )
553  CALL chkxer( 'ZHEEV ', infot, nout, lerr, ok )
554  infot = 5
555  CALL zheev( 'N', 'U', 2, a, 1, x, w, 3, rw, info )
556  CALL chkxer( 'ZHEEV ', infot, nout, lerr, ok )
557  infot = 8
558  CALL zheev( 'N', 'U', 2, a, 2, x, w, 2, rw, info )
559  CALL chkxer( 'ZHEEV ', infot, nout, lerr, ok )
560  nt = nt + 5
561 *
562 * ZHEEV_2STAGE
563 *
564  srnamt = 'ZHEEV_2STAGE '
565  infot = 1
566  CALL zheev_2stage( '/', 'U', 0, a, 1, x, w, 1, rw, info )
567  CALL chkxer( 'ZHEEV_2STAGE ', infot, nout, lerr, ok )
568  infot = 1
569  CALL zheev_2stage( 'V', 'U', 0, a, 1, x, w, 1, rw, info )
570  CALL chkxer( 'ZHEEV_2STAGE ', infot, nout, lerr, ok )
571  infot = 2
572  CALL zheev_2stage( 'N', '/', 0, a, 1, x, w, 1, rw, info )
573  CALL chkxer( 'ZHEEV_2STAGE ', infot, nout, lerr, ok )
574  infot = 3
575  CALL zheev_2stage( 'N', 'U', -1, a, 1, x, w, 1, rw, info )
576  CALL chkxer( 'ZHEEV_2STAGE ', infot, nout, lerr, ok )
577  infot = 5
578  CALL zheev_2stage( 'N', 'U', 2, a, 1, x, w, 3, rw, info )
579  CALL chkxer( 'ZHEEV_2STAGE ', infot, nout, lerr, ok )
580  infot = 8
581  CALL zheev_2stage( 'N', 'U', 2, a, 2, x, w, 2, rw, info )
582  CALL chkxer( 'ZHEEV_2STAGE ', infot, nout, lerr, ok )
583  nt = nt + 6
584 *
585 * ZHEEVX
586 *
587  srnamt = 'ZHEEVX'
588  infot = 1
589  CALL zheevx( '/', 'A', 'U', 0, a, 1, 0.0d0, 0.0d0, 0, 0, 0.0d0,
590  $ m, x, z, 1, w, 1, rw, iw, i3, info )
591  CALL chkxer( 'ZHEEVX', infot, nout, lerr, ok )
592  infot = 2
593  CALL zheevx( 'V', '/', 'U', 0, a, 1, 0.0d0, 1.0d0, 1, 0, 0.0d0,
594  $ m, x, z, 1, w, 1, rw, iw, i3, info )
595  CALL chkxer( 'ZHEEVX', infot, nout, lerr, ok )
596  infot = 3
597  CALL zheevx( 'V', 'A', '/', 0, a, 1, 0.0d0, 0.0d0, 0, 0, 0.0d0,
598  $ m, x, z, 1, w, 1, rw, iw, i3, info )
599  infot = 4
600  CALL zheevx( 'V', 'A', 'U', -1, a, 1, 0.0d0, 0.0d0, 0, 0,
601  $ 0.0d0, m, x, z, 1, w, 1, rw, iw, i3, info )
602  CALL chkxer( 'ZHEEVX', infot, nout, lerr, ok )
603  infot = 6
604  CALL zheevx( 'V', 'A', 'U', 2, a, 1, 0.0d0, 0.0d0, 0, 0, 0.0d0,
605  $ m, x, z, 2, w, 3, rw, iw, i3, info )
606  CALL chkxer( 'ZHEEVX', infot, nout, lerr, ok )
607  infot = 8
608  CALL zheevx( 'V', 'V', 'U', 1, a, 1, 0.0d0, 0.0d0, 0, 0, 0.0d0,
609  $ m, x, z, 1, w, 1, rw, iw, i3, info )
610  CALL chkxer( 'ZHEEVX', infot, nout, lerr, ok )
611  infot = 9
612  CALL zheevx( 'V', 'I', 'U', 1, a, 1, 0.0d0, 0.0d0, 0, 0, 0.0d0,
613  $ m, x, z, 1, w, 1, rw, iw, i3, info )
614  CALL chkxer( 'ZHEEVX', infot, nout, lerr, ok )
615  infot = 10
616  CALL zheevx( 'V', 'I', 'U', 2, a, 2, 0.0d0, 0.0d0, 2, 1, 0.0d0,
617  $ m, x, z, 2, w, 3, rw, iw, i3, info )
618  CALL chkxer( 'ZHEEVX', infot, nout, lerr, ok )
619  infot = 15
620  CALL zheevx( 'V', 'A', 'U', 2, a, 2, 0.0d0, 0.0d0, 0, 0, 0.0d0,
621  $ m, x, z, 1, w, 3, rw, iw, i3, info )
622  CALL chkxer( 'ZHEEVX', infot, nout, lerr, ok )
623  infot = 17
624  CALL zheevx( 'V', 'A', 'U', 2, a, 2, 0.0d0, 0.0d0, 0, 0, 0.0d0,
625  $ m, x, z, 2, w, 2, rw, iw, i1, info )
626  CALL chkxer( 'ZHEEVX', infot, nout, lerr, ok )
627  nt = nt + 10
628 *
629 * ZHEEVX_2STAGE
630 *
631  srnamt = 'ZHEEVX_2STAGE'
632  infot = 1
633  CALL zheevx_2stage( '/', 'A', 'U', 0, a, 1,
634  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
635  $ m, x, z, 1, w, 1, rw, iw, i3, info )
636  CALL chkxer( 'ZHEEVX_2STAGE', infot, nout, lerr, ok )
637  infot = 1
638  CALL zheevx_2stage( 'V', 'A', 'U', 0, a, 1,
639  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
640  $ m, x, z, 1, w, 1, rw, iw, i3, info )
641  CALL chkxer( 'ZHEEVX_2STAGE', infot, nout, lerr, ok )
642  infot = 2
643  CALL zheevx_2stage( 'N', '/', 'U', 0, a, 1,
644  $ 0.0d0, 1.0d0, 1, 0, 0.0d0,
645  $ m, x, z, 1, w, 1, rw, iw, i3, info )
646  CALL chkxer( 'ZHEEVX_2STAGE', infot, nout, lerr, ok )
647  infot = 3
648  CALL zheevx_2stage( 'N', 'A', '/', 0, a, 1,
649  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
650  $ m, x, z, 1, w, 1, rw, iw, i3, info )
651  infot = 4
652  CALL zheevx_2stage( 'N', 'A', 'U', -1, a, 1,
653  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
654  $ m, x, z, 1, w, 1, rw, iw, i3, info )
655  CALL chkxer( 'ZHEEVX_2STAGE', infot, nout, lerr, ok )
656  infot = 6
657  CALL zheevx_2stage( 'N', 'A', 'U', 2, a, 1,
658  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
659  $ m, x, z, 2, w, 3, rw, iw, i3, info )
660  CALL chkxer( 'ZHEEVX_2STAGE', infot, nout, lerr, ok )
661  infot = 8
662  CALL zheevx_2stage( 'N', 'V', 'U', 1, a, 1,
663  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
664  $ m, x, z, 1, w, 1, rw, iw, i3, info )
665  CALL chkxer( 'ZHEEVX_2STAGE', infot, nout, lerr, ok )
666  infot = 9
667  CALL zheevx_2stage( 'N', 'I', 'U', 1, a, 1,
668  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
669  $ m, x, z, 1, w, 1, rw, iw, i3, info )
670  CALL chkxer( 'ZHEEVX_2STAGE', infot, nout, lerr, ok )
671  infot = 10
672  CALL zheevx_2stage( 'N', 'I', 'U', 2, a, 2,
673  $ 0.0d0, 0.0d0, 2, 1, 0.0d0,
674  $ m, x, z, 2, w, 3, rw, iw, i3, info )
675  CALL chkxer( 'ZHEEVX_2STAGE', infot, nout, lerr, ok )
676  infot = 15
677  CALL zheevx_2stage( 'N', 'A', 'U', 2, a, 2,
678  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
679  $ m, x, z, 0, w, 3, rw, iw, i3, info )
680  CALL chkxer( 'ZHEEVX_2STAGE', infot, nout, lerr, ok )
681  infot = 17
682  CALL zheevx_2stage( 'N', 'A', 'U', 2, a, 2,
683  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
684  $ m, x, z, 2, w, 0, rw, iw, i1, info )
685  CALL chkxer( 'ZHEEVX_2STAGE', infot, nout, lerr, ok )
686  nt = nt + 11
687 *
688 * ZHEEVR
689 *
690  srnamt = 'ZHEEVR'
691  n = 1
692  infot = 1
693  CALL zheevr( '/', 'A', 'U', 0, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
694  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
695  $ 10*n, info )
696  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
697  infot = 2
698  CALL zheevr( 'V', '/', 'U', 0, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
699  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
700  $ 10*n, info )
701  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
702  infot = 3
703  CALL zheevr( 'V', 'A', '/', -1, a, 1, 0.0d0, 0.0d0, 1, 1,
704  $ 0.0d0, m, r, z, 1, iw, q, 2*n, rw, 24*n,
705  $ iw( 2*n+1 ), 10*n, info )
706  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
707  infot = 4
708  CALL zheevr( 'V', 'A', 'U', -1, a, 1, 0.0d0, 0.0d0, 1, 1,
709  $ 0.0d0, m, r, z, 1, iw, q, 2*n, rw, 24*n,
710  $ iw( 2*n+1 ), 10*n, info )
711  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
712  infot = 6
713  CALL zheevr( 'V', 'A', 'U', 2, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
714  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
715  $ 10*n, info )
716  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
717  infot = 8
718  CALL zheevr( 'V', 'V', 'U', 1, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
719  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
720  $ 10*n, info )
721  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
722  infot = 9
723  CALL zheevr( 'V', 'I', 'U', 1, a, 1, 0.0d0, 0.0d0, 0, 1, 0.0d0,
724  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
725  $ 10*n, info )
726  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
727  infot = 10
728 *
729  CALL zheevr( 'V', 'I', 'U', 2, a, 2, 0.0d0, 0.0d0, 2, 1, 0.0d0,
730  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
731  $ 10*n, info )
732  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
733  infot = 15
734  CALL zheevr( 'V', 'I', 'U', 1, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
735  $ m, r, z, 0, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
736  $ 10*n, info )
737  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
738  infot = 18
739  CALL zheevr( 'V', 'I', 'U', 1, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
740  $ m, r, z, 1, iw, q, 2*n-1, rw, 24*n, iw( 2*n+1 ),
741  $ 10*n, info )
742  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
743  infot = 20
744  CALL zheevr( 'V', 'I', 'U', 1, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
745  $ m, r, z, 1, iw, q, 2*n, rw, 24*n-1, iw( 2*n-1 ),
746  $ 10*n, info )
747  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
748  infot = 22
749  CALL zheevr( 'V', 'I', 'U', 1, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
750  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw, 10*n-1,
751  $ info )
752  CALL chkxer( 'ZHEEVR', infot, nout, lerr, ok )
753  nt = nt + 12
754 *
755 * ZHEEVR_2STAGE
756 *
757  srnamt = 'ZHEEVR_2STAGE'
758  n = 1
759  infot = 1
760  CALL zheevr_2stage( '/', 'A', 'U', 0, a, 1,
761  $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
762  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
763  $ 10*n, info )
764  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
765  infot = 1
766  CALL zheevr_2stage( 'V', 'A', 'U', 0, a, 1,
767  $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
768  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
769  $ 10*n, info )
770  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
771  infot = 2
772  CALL zheevr_2stage( 'N', '/', 'U', 0, a, 1,
773  $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
774  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
775  $ 10*n, info )
776  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
777  infot = 3
778  CALL zheevr_2stage( 'N', 'A', '/', -1, a, 1,
779  $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
780  $ m, r, z, 1, iw, q, 2*n, rw, 24*n,
781  $ iw( 2*n+1 ), 10*n, info )
782  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
783  infot = 4
784  CALL zheevr_2stage( 'N', 'A', 'U', -1, a, 1,
785  $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
786  $ m, r, z, 1, iw, q, 2*n, rw, 24*n,
787  $ iw( 2*n+1 ), 10*n, info )
788  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
789  infot = 6
790  CALL zheevr_2stage( 'N', 'A', 'U', 2, a, 1,
791  $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
792  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
793  $ 10*n, info )
794  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
795  infot = 8
796  CALL zheevr_2stage( 'N', 'V', 'U', 1, a, 1,
797  $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
798  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
799  $ 10*n, info )
800  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
801  infot = 9
802  CALL zheevr_2stage( 'N', 'I', 'U', 1, a, 1,
803  $ 0.0d0, 0.0d0, 0, 1, 0.0d0,
804  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
805  $ 10*n, info )
806  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
807  infot = 10
808  CALL zheevr_2stage( 'N', 'I', 'U', 2, a, 2,
809  $ 0.0d0, 0.0d0, 2, 1, 0.0d0,
810  $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
811  $ 10*n, info )
812  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
813  infot = 15
814  CALL zheevr_2stage( 'N', 'I', 'U', 1, a, 1,
815  $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
816  $ m, r, z, 0, iw, q, 2*n, rw, 24*n, iw( 2*n+1 ),
817  $ 10*n, info )
818  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
819  infot = 18
820  CALL zheevr_2stage( 'N', 'I', 'U', 1, a, 1,
821  $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
822  $ m, r, z, 1, iw, q, 2*n-1, rw, 24*n, iw( 2*n+1 ),
823  $ 10*n, info )
824  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
825  infot = 20
826  CALL zheevr_2stage( 'N', 'I', 'U', 1, a, 1,
827  $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
828  $ m, r, z, 1, iw, q, 26*n, rw, 24*n-1, iw( 2*n-1 ),
829  $ 10*n, info )
830  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
831  infot = 22
832  CALL zheevr_2stage( 'N', 'I', 'U', 1, a, 1,
833  $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
834  $ m, r, z, 1, iw, q, 26*n, rw, 24*n, iw, 10*n-1,
835  $ info )
836  CALL chkxer( 'ZHEEVR_2STAGE', infot, nout, lerr, ok )
837  nt = nt + 13
838 *
839 * ZHPEVD
840 *
841  srnamt = 'ZHPEVD'
842  infot = 1
843  CALL zhpevd( '/', 'U', 0, a, x, z, 1, w, 1, rw, 1, iw, 1,
844  $ info )
845  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
846  infot = 2
847  CALL zhpevd( 'N', '/', 0, a, x, z, 1, w, 1, rw, 1, iw, 1,
848  $ info )
849  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
850  infot = 3
851  CALL zhpevd( 'N', 'U', -1, a, x, z, 1, w, 1, rw, 1, iw, 1,
852  $ info )
853  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
854  infot = 7
855  CALL zhpevd( 'V', 'U', 2, a, x, z, 1, w, 4, rw, 25, iw, 12,
856  $ info )
857  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
858  infot = 9
859  CALL zhpevd( 'N', 'U', 1, a, x, z, 1, w, 0, rw, 1, iw, 1,
860  $ info )
861  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
862  infot = 9
863  CALL zhpevd( 'N', 'U', 2, a, x, z, 2, w, 1, rw, 2, iw, 1,
864  $ info )
865  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
866  infot = 9
867  CALL zhpevd( 'V', 'U', 2, a, x, z, 2, w, 2, rw, 25, iw, 12,
868  $ info )
869  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
870  infot = 11
871  CALL zhpevd( 'N', 'U', 1, a, x, z, 1, w, 1, rw, 0, iw, 1,
872  $ info )
873  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
874  infot = 11
875  CALL zhpevd( 'N', 'U', 2, a, x, z, 2, w, 2, rw, 1, iw, 1,
876  $ info )
877  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
878  infot = 11
879  CALL zhpevd( 'V', 'U', 2, a, x, z, 2, w, 4, rw, 18, iw, 12,
880  $ info )
881  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
882  infot = 13
883  CALL zhpevd( 'N', 'U', 1, a, x, z, 1, w, 1, rw, 1, iw, 0,
884  $ info )
885  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
886  infot = 13
887  CALL zhpevd( 'N', 'U', 2, a, x, z, 2, w, 2, rw, 2, iw, 0,
888  $ info )
889  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
890  infot = 13
891  CALL zhpevd( 'V', 'U', 2, a, x, z, 2, w, 4, rw, 25, iw, 2,
892  $ info )
893  CALL chkxer( 'ZHPEVD', infot, nout, lerr, ok )
894  nt = nt + 13
895 *
896 * ZHPEV
897 *
898  srnamt = 'ZHPEV '
899  infot = 1
900  CALL zhpev( '/', 'U', 0, a, x, z, 1, w, rw, info )
901  CALL chkxer( 'ZHPEV ', infot, nout, lerr, ok )
902  infot = 2
903  CALL zhpev( 'N', '/', 0, a, x, z, 1, w, rw, info )
904  CALL chkxer( 'ZHPEV ', infot, nout, lerr, ok )
905  infot = 3
906  CALL zhpev( 'N', 'U', -1, a, x, z, 1, w, rw, info )
907  CALL chkxer( 'ZHPEV ', infot, nout, lerr, ok )
908  infot = 7
909  CALL zhpev( 'V', 'U', 2, a, x, z, 1, w, rw, info )
910  CALL chkxer( 'ZHPEV ', infot, nout, lerr, ok )
911  nt = nt + 4
912 *
913 * ZHPEVX
914 *
915  srnamt = 'ZHPEVX'
916  infot = 1
917  CALL zhpevx( '/', 'A', 'U', 0, a, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
918  $ x, z, 1, w, rw, iw, i3, info )
919  CALL chkxer( 'ZHPEVX', infot, nout, lerr, ok )
920  infot = 2
921  CALL zhpevx( 'V', '/', 'U', 0, a, 0.0d0, 1.0d0, 1, 0, 0.0d0, m,
922  $ x, z, 1, w, rw, iw, i3, info )
923  CALL chkxer( 'ZHPEVX', infot, nout, lerr, ok )
924  infot = 3
925  CALL zhpevx( 'V', 'A', '/', 0, a, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
926  $ x, z, 1, w, rw, iw, i3, info )
927  CALL chkxer( 'ZHPEVX', infot, nout, lerr, ok )
928  infot = 4
929  CALL zhpevx( 'V', 'A', 'U', -1, a, 0.0d0, 0.0d0, 0, 0, 0.0d0,
930  $ m, x, z, 1, w, rw, iw, i3, info )
931  CALL chkxer( 'ZHPEVX', infot, nout, lerr, ok )
932  infot = 7
933  CALL zhpevx( 'V', 'V', 'U', 1, a, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
934  $ x, z, 1, w, rw, iw, i3, info )
935  CALL chkxer( 'ZHPEVX', infot, nout, lerr, ok )
936  infot = 8
937  CALL zhpevx( 'V', 'I', 'U', 1, a, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
938  $ x, z, 1, w, rw, iw, i3, info )
939  CALL chkxer( 'ZHPEVX', infot, nout, lerr, ok )
940  infot = 9
941  CALL zhpevx( 'V', 'I', 'U', 2, a, 0.0d0, 0.0d0, 2, 1, 0.0d0, m,
942  $ x, z, 2, w, rw, iw, i3, info )
943  CALL chkxer( 'ZHPEVX', infot, nout, lerr, ok )
944  infot = 14
945  CALL zhpevx( 'V', 'A', 'U', 2, a, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
946  $ x, z, 1, w, rw, iw, i3, info )
947  CALL chkxer( 'ZHPEVX', infot, nout, lerr, ok )
948  nt = nt + 8
949 *
950 * Test error exits for the HB path.
951 *
952  ELSE IF( lsamen( 2, c2, 'HB' ) ) THEN
953 *
954 * ZHBTRD
955 *
956  srnamt = 'ZHBTRD'
957  infot = 1
958  CALL zhbtrd( '/', 'U', 0, 0, a, 1, d, e, z, 1, w, info )
959  CALL chkxer( 'ZHBTRD', infot, nout, lerr, ok )
960  infot = 2
961  CALL zhbtrd( 'N', '/', 0, 0, a, 1, d, e, z, 1, w, info )
962  CALL chkxer( 'ZHBTRD', infot, nout, lerr, ok )
963  infot = 3
964  CALL zhbtrd( 'N', 'U', -1, 0, a, 1, d, e, z, 1, w, info )
965  CALL chkxer( 'ZHBTRD', infot, nout, lerr, ok )
966  infot = 4
967  CALL zhbtrd( 'N', 'U', 0, -1, a, 1, d, e, z, 1, w, info )
968  CALL chkxer( 'ZHBTRD', infot, nout, lerr, ok )
969  infot = 6
970  CALL zhbtrd( 'N', 'U', 1, 1, a, 1, d, e, z, 1, w, info )
971  CALL chkxer( 'ZHBTRD', infot, nout, lerr, ok )
972  infot = 10
973  CALL zhbtrd( 'V', 'U', 2, 0, a, 1, d, e, z, 1, w, info )
974  CALL chkxer( 'ZHBTRD', infot, nout, lerr, ok )
975  nt = nt + 6
976 *
977 * ZHETRD_HB2ST
978 *
979  srnamt = 'ZHETRD_HB2ST'
980  infot = 1
981  CALL zhetrd_hb2st( '/', 'N', 'U', 0, 0, a, 1, d, e,
982  $ c, 1, w, 1, info )
983  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
984  infot = 2
985  CALL zhetrd_hb2st( 'N', '/', 'U', 0, 0, a, 1, d, e,
986  $ c, 1, w, 1, info )
987  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
988  infot = 2
989  CALL zhetrd_hb2st( 'N', 'H', 'U', 0, 0, a, 1, d, e,
990  $ c, 1, w, 1, info )
991  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
992  infot = 3
993  CALL zhetrd_hb2st( 'N', 'N', '/', 0, 0, a, 1, d, e,
994  $ c, 1, w, 1, info )
995  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
996  infot = 4
997  CALL zhetrd_hb2st( 'N', 'N', 'U', -1, 0, a, 1, d, e,
998  $ c, 1, w, 1, info )
999  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
1000  infot = 5
1001  CALL zhetrd_hb2st( 'N', 'N', 'U', 0, -1, a, 1, d, e,
1002  $ c, 1, w, 1, info )
1003  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
1004  infot = 7
1005  CALL zhetrd_hb2st( 'N', 'N', 'U', 0, 1, a, 1, d, e,
1006  $ c, 1, w, 1, info )
1007  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
1008  infot = 11
1009  CALL zhetrd_hb2st( 'N', 'N', 'U', 0, 0, a, 1, d, e,
1010  $ c, 0, w, 1, info )
1011  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
1012  infot = 13
1013  CALL zhetrd_hb2st( 'N', 'N', 'U', 0, 0, a, 1, d, e,
1014  $ c, 1, w, 0, info )
1015  CALL chkxer( 'ZHETRD_HB2ST', infot, nout, lerr, ok )
1016  nt = nt + 9
1017 *
1018 * ZHBEVD
1019 *
1020  srnamt = 'ZHBEVD'
1021  infot = 1
1022  CALL zhbevd( '/', 'U', 0, 0, a, 1, x, z, 1, w, 1, rw, 1, iw, 1,
1023  $ info )
1024  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1025  infot = 2
1026  CALL zhbevd( 'N', '/', 0, 0, a, 1, x, z, 1, w, 1, rw, 1, iw, 1,
1027  $ info )
1028  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1029  infot = 3
1030  CALL zhbevd( 'N', 'U', -1, 0, a, 1, x, z, 1, w, 1, rw, 1, iw,
1031  $ 1, info )
1032  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1033  infot = 4
1034  CALL zhbevd( 'N', 'U', 0, -1, a, 1, x, z, 1, w, 1, rw, 1, iw,
1035  $ 1, info )
1036  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1037  infot = 6
1038  CALL zhbevd( 'N', 'U', 2, 1, a, 1, x, z, 1, w, 2, rw, 2, iw, 1,
1039  $ info )
1040  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1041  infot = 9
1042  CALL zhbevd( 'V', 'U', 2, 1, a, 2, x, z, 1, w, 8, rw, 25, iw,
1043  $ 12, info )
1044  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1045  infot = 11
1046  CALL zhbevd( 'N', 'U', 1, 0, a, 1, x, z, 1, w, 0, rw, 1, iw, 1,
1047  $ info )
1048  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1049  infot = 11
1050  CALL zhbevd( 'N', 'U', 2, 1, a, 2, x, z, 2, w, 1, rw, 2, iw, 1,
1051  $ info )
1052  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1053  infot = 11
1054  CALL zhbevd( 'V', 'U', 2, 1, a, 2, x, z, 2, w, 2, rw, 25, iw,
1055  $ 12, info )
1056  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1057  infot = 13
1058  CALL zhbevd( 'N', 'U', 1, 0, a, 1, x, z, 1, w, 1, rw, 0, iw, 1,
1059  $ info )
1060  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1061  infot = 13
1062  CALL zhbevd( 'N', 'U', 2, 1, a, 2, x, z, 2, w, 2, rw, 1, iw, 1,
1063  $ info )
1064  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1065  infot = 13
1066  CALL zhbevd( 'V', 'U', 2, 1, a, 2, x, z, 2, w, 8, rw, 2, iw,
1067  $ 12, info )
1068  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1069  infot = 15
1070  CALL zhbevd( 'N', 'U', 1, 0, a, 1, x, z, 1, w, 1, rw, 1, iw, 0,
1071  $ info )
1072  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1073  infot = 15
1074  CALL zhbevd( 'N', 'U', 2, 1, a, 2, x, z, 2, w, 2, rw, 2, iw, 0,
1075  $ info )
1076  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1077  infot = 15
1078  CALL zhbevd( 'V', 'U', 2, 1, a, 2, x, z, 2, w, 8, rw, 25, iw,
1079  $ 2, info )
1080  CALL chkxer( 'ZHBEVD', infot, nout, lerr, ok )
1081  nt = nt + 15
1082 *
1083 * ZHBEVD_2STAGE
1084 *
1085  srnamt = 'ZHBEVD_2STAGE'
1086  infot = 1
1087  CALL zhbevd_2stage( '/', 'U', 0, 0, a, 1, x, z, 1,
1088  $ w, 1, rw, 1, iw, 1, info )
1089  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1090  infot = 1
1091  CALL zhbevd_2stage( 'V', 'U', 0, 0, a, 1, x, z, 1,
1092  $ w, 1, rw, 1, iw, 1, info )
1093  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1094  infot = 2
1095  CALL zhbevd_2stage( 'N', '/', 0, 0, a, 1, x, z, 1,
1096  $ w, 1, rw, 1, iw, 1, info )
1097  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1098  infot = 3
1099  CALL zhbevd_2stage( 'N', 'U', -1, 0, a, 1, x, z, 1,
1100  $ w, 1, rw, 1, iw, 1, info )
1101  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1102  infot = 4
1103  CALL zhbevd_2stage( 'N', 'U', 0, -1, a, 1, x, z, 1,
1104  $ w, 1, rw, 1, iw, 1, info )
1105  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1106  infot = 6
1107  CALL zhbevd_2stage( 'N', 'U', 2, 1, a, 1, x, z, 1,
1108  $ w, 2, rw, 2, iw, 1, info )
1109  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1110  infot = 9
1111  CALL zhbevd_2stage( 'N', 'U', 2, 1, a, 2, x, z, 0,
1112  $ w, 8, rw, 25, iw, 12, info )
1113  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1114  infot = 11
1115  CALL zhbevd_2stage( 'N', 'U', 1, 0, a, 1, x, z, 1,
1116  $ w, 0, rw, 1, iw, 1, info )
1117  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1118  infot = 11
1119  CALL zhbevd_2stage( 'N', 'U', 2, 1, a, 2, x, z, 2,
1120  $ w, 1, rw, 2, iw, 1, info )
1121  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1122 * INFOT = 11
1123 * CALL ZHBEVD_2STAGE( 'V', 'U', 2, 1, A, 2, X, Z, 2,
1124 * $ W, 2, RW, 25, IW, 12, INFO )
1125 * CALL CHKXER( 'ZHBEVD_2STAGE', INFOT, NOUT, LERR, OK )
1126  infot = 13
1127  CALL zhbevd_2stage( 'N', 'U', 1, 0, a, 1, x, z, 1,
1128  $ w, 1, rw, 0, iw, 1, info )
1129  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1130  infot = 13
1131  CALL zhbevd_2stage( 'N', 'U', 2, 1, a, 2, x, z, 2,
1132  $ w, 25, rw, 1, iw, 1, info )
1133  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1134 * INFOT = 13
1135 * CALL ZHBEVD_2STAGE( 'V', 'U', 2, 1, A, 2, X, Z, 2,
1136 * $ W, 25, RW, 2, IW, 12, INFO )
1137 * CALL CHKXER( 'ZHBEVD_2STAGE', INFOT, NOUT, LERR, OK )
1138  infot = 15
1139  CALL zhbevd_2stage( 'N', 'U', 1, 0, a, 1, x, z, 1,
1140  $ w, 1, rw, 1, iw, 0, info )
1141  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1142  infot = 15
1143  CALL zhbevd_2stage( 'N', 'U', 2, 1, a, 2, x, z, 2,
1144  $ w, 25, rw, 2, iw, 0, info )
1145  CALL chkxer( 'ZHBEVD_2STAGE', infot, nout, lerr, ok )
1146 * INFOT = 15
1147 * CALL ZHBEVD_2STAGE( 'V', 'U', 2, 1, A, 2, X, Z, 2,
1148 * $ W, 25, RW, 25, IW, 2, INFO )
1149 * CALL CHKXER( 'ZHBEVD_2STAGE', INFOT, NOUT, LERR, OK )
1150  nt = nt + 13
1151 *
1152 * ZHBEV
1153 *
1154  srnamt = 'ZHBEV '
1155  infot = 1
1156  CALL zhbev( '/', 'U', 0, 0, a, 1, x, z, 1, w, rw, info )
1157  CALL chkxer( 'ZHBEV ', infot, nout, lerr, ok )
1158  infot = 2
1159  CALL zhbev( 'N', '/', 0, 0, a, 1, x, z, 1, w, rw, info )
1160  CALL chkxer( 'ZHBEV ', infot, nout, lerr, ok )
1161  infot = 3
1162  CALL zhbev( 'N', 'U', -1, 0, a, 1, x, z, 1, w, rw, info )
1163  CALL chkxer( 'ZHBEV ', infot, nout, lerr, ok )
1164  infot = 4
1165  CALL zhbev( 'N', 'U', 0, -1, a, 1, x, z, 1, w, rw, info )
1166  CALL chkxer( 'ZHBEV ', infot, nout, lerr, ok )
1167  infot = 6
1168  CALL zhbev( 'N', 'U', 2, 1, a, 1, x, z, 1, w, rw, info )
1169  CALL chkxer( 'ZHBEV ', infot, nout, lerr, ok )
1170  infot = 9
1171  CALL zhbev( 'V', 'U', 2, 0, a, 1, x, z, 1, w, rw, info )
1172  CALL chkxer( 'ZHBEV ', infot, nout, lerr, ok )
1173  nt = nt + 6
1174 *
1175 * ZHBEV_2STAGE
1176 *
1177  srnamt = 'ZHBEV_2STAGE '
1178  infot = 1
1179  CALL zhbev_2stage( '/', 'U', 0, 0, a, 1, x,
1180  $ z, 1, w, 0, rw, info )
1181  CALL chkxer( 'ZHBEV_2STAGE ', infot, nout, lerr, ok )
1182  infot = 1
1183  CALL zhbev_2stage( 'V', 'U', 0, 0, a, 1, x,
1184  $ z, 1, w, 0, rw, info )
1185  CALL chkxer( 'ZHBEV_2STAGE ', infot, nout, lerr, ok )
1186  infot = 2
1187  CALL zhbev_2stage( 'N', '/', 0, 0, a, 1, x,
1188  $ z, 1, w, 0, rw, info )
1189  CALL chkxer( 'ZHBEV_2STAGE ', infot, nout, lerr, ok )
1190  infot = 3
1191  CALL zhbev_2stage( 'N', 'U', -1, 0, a, 1, x,
1192  $ z, 1, w, 0, rw, info )
1193  CALL chkxer( 'ZHBEV_2STAGE ', infot, nout, lerr, ok )
1194  infot = 4
1195  CALL zhbev_2stage( 'N', 'U', 0, -1, a, 1, x,
1196  $ z, 1, w, 0, rw, info )
1197  CALL chkxer( 'ZHBEV_2STAGE ', infot, nout, lerr, ok )
1198  infot = 6
1199  CALL zhbev_2stage( 'N', 'U', 2, 1, a, 1, x,
1200  $ z, 1, w, 0, rw, info )
1201  CALL chkxer( 'ZHBEV_2STAGE ', infot, nout, lerr, ok )
1202  infot = 9
1203  CALL zhbev_2stage( 'N', 'U', 2, 0, a, 1, x,
1204  $ z, 0, w, 0, rw, info )
1205  CALL chkxer( 'ZHBEV_2STAGE ', infot, nout, lerr, ok )
1206  infot = 11
1207  CALL zhbev_2stage( 'N', 'U', 2, 0, a, 1, x,
1208  $ z, 1, w, 0, rw, info )
1209  CALL chkxer( 'ZHBEV_2STAGE ', infot, nout, lerr, ok )
1210  nt = nt + 8
1211 *
1212 * ZHBEVX
1213 *
1214  srnamt = 'ZHBEVX'
1215  infot = 1
1216  CALL zhbevx( '/', 'A', 'U', 0, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1217  $ 0, 0.0d0, m, x, z, 1, w, rw, iw, i3, info )
1218  CALL chkxer( 'ZHBEVX', infot, nout, lerr, ok )
1219  infot = 2
1220  CALL zhbevx( 'V', '/', 'U', 0, 0, a, 1, q, 1, 0.0d0, 1.0d0, 1,
1221  $ 0, 0.0d0, m, x, z, 1, w, rw, iw, i3, info )
1222  CALL chkxer( 'ZHBEVX', infot, nout, lerr, ok )
1223  infot = 3
1224  CALL zhbevx( 'V', 'A', '/', 0, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1225  $ 0, 0.0d0, m, x, z, 1, w, rw, iw, i3, info )
1226  infot = 4
1227  CALL zhbevx( 'V', 'A', 'U', -1, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1228  $ 0, 0.0d0, m, x, z, 1, w, rw, iw, i3, info )
1229  CALL chkxer( 'ZHBEVX', infot, nout, lerr, ok )
1230  infot = 5
1231  CALL zhbevx( 'V', 'A', 'U', 0, -1, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1232  $ 0, 0.0d0, m, x, z, 1, w, rw, iw, i3, info )
1233  CALL chkxer( 'ZHBEVX', infot, nout, lerr, ok )
1234  infot = 7
1235  CALL zhbevx( 'V', 'A', 'U', 2, 1, a, 1, q, 2, 0.0d0, 0.0d0, 0,
1236  $ 0, 0.0d0, m, x, z, 2, w, rw, iw, i3, info )
1237  CALL chkxer( 'ZHBEVX', infot, nout, lerr, ok )
1238  infot = 9
1239  CALL zhbevx( 'V', 'A', 'U', 2, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1240  $ 0, 0.0d0, m, x, z, 2, w, rw, iw, i3, info )
1241  CALL chkxer( 'ZHBEVX', infot, nout, lerr, ok )
1242  infot = 11
1243  CALL zhbevx( 'V', 'V', 'U', 1, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1244  $ 0, 0.0d0, m, x, z, 1, w, rw, iw, i3, info )
1245  CALL chkxer( 'ZHBEVX', infot, nout, lerr, ok )
1246  infot = 12
1247  CALL zhbevx( 'V', 'I', 'U', 1, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1248  $ 0, 0.0d0, m, x, z, 1, w, rw, iw, i3, info )
1249  CALL chkxer( 'ZHBEVX', infot, nout, lerr, ok )
1250  infot = 13
1251  CALL zhbevx( 'V', 'I', 'U', 1, 0, a, 1, q, 1, 0.0d0, 0.0d0, 1,
1252  $ 2, 0.0d0, m, x, z, 1, w, rw, iw, i3, info )
1253  CALL chkxer( 'ZHBEVX', infot, nout, lerr, ok )
1254  infot = 18
1255  CALL zhbevx( 'V', 'A', 'U', 2, 0, a, 1, q, 2, 0.0d0, 0.0d0, 0,
1256  $ 0, 0.0d0, m, x, z, 1, w, rw, iw, i3, info )
1257  CALL chkxer( 'ZHBEVX', infot, nout, lerr, ok )
1258  nt = nt + 11
1259 *
1260 * ZHBEVX_2STAGE
1261 *
1262  srnamt = 'ZHBEVX_2STAGE'
1263  infot = 1
1264  CALL zhbevx_2stage( '/', 'A', 'U', 0, 0, a, 1, q, 1,
1265  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
1266  $ m, x, z, 1, w, 0, rw, iw, i3, info )
1267  infot = 1
1268  CALL zhbevx_2stage( 'V', 'A', 'U', 0, 0, a, 1, q, 1,
1269  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
1270  $ m, x, z, 1, w, 0, rw, iw, i3, info )
1271  CALL chkxer( 'ZHBEVX_2STAGE', infot, nout, lerr, ok )
1272  infot = 2
1273  CALL zhbevx_2stage( 'N', '/', 'U', 0, 0, a, 1, q, 1,
1274  $ 0.0d0, 1.0d0, 1, 0, 0.0d0,
1275  $ m, x, z, 1, w, 0, rw, iw, i3, info )
1276  CALL chkxer( 'ZHBEVX_2STAGE', infot, nout, lerr, ok )
1277  infot = 3
1278  CALL zhbevx_2stage( 'N', 'A', '/', 0, 0, a, 1, q, 1,
1279  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
1280  $ m, x, z, 1, w, 0, rw, iw, i3, info )
1281  infot = 4
1282  CALL zhbevx_2stage( 'N', 'A', 'U', -1, 0, a, 1, q, 1,
1283  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
1284  $ m, x, z, 1, w, 0, rw, iw, i3, info )
1285  CALL chkxer( 'ZHBEVX_2STAGE', infot, nout, lerr, ok )
1286  infot = 5
1287  CALL zhbevx_2stage( 'N', 'A', 'U', 0, -1, a, 1, q, 1,
1288  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
1289  $ m, x, z, 1, w, 0, rw, iw, i3, info )
1290  CALL chkxer( 'ZHBEVX_2STAGE', infot, nout, lerr, ok )
1291  infot = 7
1292  CALL zhbevx_2stage( 'N', 'A', 'U', 2, 1, a, 1, q, 2,
1293  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
1294  $ m, x, z, 2, w, 0, rw, iw, i3, info )
1295  CALL chkxer( 'ZHBEVX_2STAGE', infot, nout, lerr, ok )
1296 * INFOT = 9
1297 * CALL ZHBEVX_2STAGE( 'V', 'A', 'U', 2, 0, A, 1, Q, 1,
1298 * $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
1299 * $ M, X, Z, 2, W, 0, RW, IW, I3, INFO )
1300 * CALL CHKXER( 'ZHBEVX_2STAGE', INFOT, NOUT, LERR, OK )
1301  infot = 11
1302  CALL zhbevx_2stage( 'N', 'V', 'U', 1, 0, a, 1, q, 1,
1303  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
1304  $ m, x, z, 1, w, 0, rw, iw, i3, info )
1305  CALL chkxer( 'ZHBEVX_2STAGE', infot, nout, lerr, ok )
1306  infot = 12
1307  CALL zhbevx_2stage( 'N', 'I', 'U', 1, 0, a, 1, q, 1,
1308  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
1309  $ m, x, z, 1, w, 0, rw, iw, i3, info )
1310  CALL chkxer( 'ZHBEVX_2STAGE', infot, nout, lerr, ok )
1311  infot = 13
1312  CALL zhbevx_2stage( 'N', 'I', 'U', 1, 0, a, 1, q, 1,
1313  $ 0.0d0, 0.0d0, 1, 2, 0.0d0,
1314  $ m, x, z, 1, w, 0, rw, iw, i3, info )
1315  CALL chkxer( 'ZHBEVX_2STAGE', infot, nout, lerr, ok )
1316  infot = 18
1317  CALL zhbevx_2stage( 'N', 'A', 'U', 2, 0, a, 1, q, 2,
1318  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
1319  $ m, x, z, 0, w, 0, rw, iw, i3, info )
1320  CALL chkxer( 'ZHBEVX_2STAGE', infot, nout, lerr, ok )
1321  infot = 20
1322  CALL zhbevx_2stage( 'N', 'A', 'U', 2, 0, a, 1, q, 2,
1323  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
1324  $ m, x, z, 1, w, 0, rw, iw, i3, info )
1325  CALL chkxer( 'ZHBEVX_2STAGE', infot, nout, lerr, ok )
1326  nt = nt + 12
1327  END IF
1328 *
1329 * Print a summary line.
1330 *
1331  IF( ok ) THEN
1332  WRITE( nout, fmt = 9999 )path, nt
1333  ELSE
1334  WRITE( nout, fmt = 9998 )path
1335  END IF
1336 *
1337  9999 FORMAT( 1x, a3, ' routines passed the tests of the error exits',
1338  $ ' (', i3, ' tests done)' )
1339  9998 FORMAT( ' *** ', a3, ' routines failed the tests of the error ',
1340  $ 'exits ***' )
1341 *
1342  RETURN
1343 *
1344 * End of ZERRST
1345 *
1346  END
zheev
subroutine zheev(JOBZ, UPLO, N, A, LDA, W, WORK, LWORK, RWORK, INFO)
ZHEEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for HE matrices
Definition: zheev.f:142
zstein
subroutine zstein(N, D, E, M, W, IBLOCK, ISPLIT, Z, LDZ, WORK, IWORK, IFAIL, INFO)
ZSTEIN
Definition: zstein.f:184
zheevr_2stage
subroutine zheevr_2stage(JOBZ, RANGE, UPLO, N, A, LDA, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, ISUPPZ, WORK, LWORK, RWORK, LRWORK, IWORK, LIWORK, INFO)
ZHEEVR_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for HE mat...
Definition: zheevr_2stage.f:408
zhpev
subroutine zhpev(JOBZ, UPLO, N, AP, W, Z, LDZ, WORK, RWORK, INFO)
ZHPEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrices
Definition: zhpev.f:140
zhbevx
subroutine zhbevx(JOBZ, RANGE, UPLO, N, KD, AB, LDAB, Q, LDQ, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, RWORK, IWORK, IFAIL, INFO)
ZHBEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
Definition: zhbevx.f:269
zungtr
subroutine zungtr(UPLO, N, A, LDA, TAU, WORK, LWORK, INFO)
ZUNGTR
Definition: zungtr.f:125
zpteqr
subroutine zpteqr(COMPZ, N, D, E, Z, LDZ, WORK, INFO)
ZPTEQR
Definition: zpteqr.f:147
zhbevd
subroutine zhbevd(JOBZ, UPLO, N, KD, AB, LDAB, W, Z, LDZ, WORK, LWORK, RWORK, LRWORK, IWORK, LIWORK, INFO)
ZHBEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
Definition: zhbevd.f:217
zheevr
subroutine zheevr(JOBZ, RANGE, UPLO, N, A, LDA, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, ISUPPZ, WORK, LWORK, RWORK, LRWORK, IWORK, LIWORK, INFO)
ZHEEVR computes the eigenvalues and, optionally, the left and/or right eigenvectors for HE matrices
Definition: zheevr.f:359
zsteqr
subroutine zsteqr(COMPZ, N, D, E, Z, LDZ, WORK, INFO)
ZSTEQR
Definition: zsteqr.f:134
chkxer
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
Definition: cblat2.f:3199
zhbev_2stage
subroutine zhbev_2stage(JOBZ, UPLO, N, KD, AB, LDAB, W, Z, LDZ, WORK, LWORK, RWORK, INFO)
ZHBEV_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER m...
Definition: zhbev_2stage.f:213
zhetrd_he2hb
subroutine zhetrd_he2hb(UPLO, N, KD, A, LDA, AB, LDAB, TAU, WORK, LWORK, INFO)
ZHETRD_HE2HB
Definition: zhetrd_he2hb.f:245
zerrst
subroutine zerrst(PATH, NUNIT)
ZERRST
Definition: zerrst.f:64
zhetrd_hb2st
subroutine zhetrd_hb2st(STAGE1, VECT, UPLO, N, KD, AB, LDAB, D, E, HOUS, LHOUS, WORK, LWORK, INFO)
ZHETRD_HB2ST reduces a complex Hermitian band matrix A to real symmetric tridiagonal form T
Definition: zhetrd_hb2st.F:232
zhbevx_2stage
subroutine zhbevx_2stage(JOBZ, RANGE, UPLO, N, KD, AB, LDAB, Q, LDQ, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, LWORK, RWORK, IWORK, IFAIL, INFO)
ZHBEVX_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER ...
Definition: zhbevx_2stage.f:329
zupmtr
subroutine zupmtr(SIDE, UPLO, TRANS, M, N, AP, TAU, C, LDC, WORK, INFO)
ZUPMTR
Definition: zupmtr.f:152
zhbtrd
subroutine zhbtrd(VECT, UPLO, N, KD, AB, LDAB, D, E, Q, LDQ, WORK, INFO)
ZHBTRD
Definition: zhbtrd.f:165
zhetrd
subroutine zhetrd(UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO)
ZHETRD
Definition: zhetrd.f:194
zheevd
subroutine zheevd(JOBZ, UPLO, N, A, LDA, W, WORK, LWORK, RWORK, LRWORK, IWORK, LIWORK, INFO)
ZHEEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for HE matrices
Definition: zheevd.f:207
zstedc
subroutine zstedc(COMPZ, N, D, E, Z, LDZ, WORK, LWORK, RWORK, LRWORK, IWORK, LIWORK, INFO)
ZSTEDC
Definition: zstedc.f:214
zheevd_2stage
subroutine zheevd_2stage(JOBZ, UPLO, N, A, LDA, W, WORK, LWORK, RWORK, LRWORK, IWORK, LIWORK, INFO)
ZHEEVD_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for HE mat...
Definition: zheevd_2stage.f:255
zunmtr
subroutine zunmtr(SIDE, UPLO, TRANS, M, N, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
ZUNMTR
Definition: zunmtr.f:173
zheevx
subroutine zheevx(JOBZ, RANGE, UPLO, N, A, LDA, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, LWORK, RWORK, IWORK, IFAIL, INFO)
ZHEEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for HE matrices
Definition: zheevx.f:261
zhptrd
subroutine zhptrd(UPLO, N, AP, D, E, TAU, INFO)
ZHPTRD
Definition: zhptrd.f:153
zheevx_2stage
subroutine zheevx_2stage(JOBZ, RANGE, UPLO, N, A, LDA, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, LWORK, RWORK, IWORK, IFAIL, INFO)
ZHEEVX_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for HE mat...
Definition: zheevx_2stage.f:308
zhpevx
subroutine zhpevx(JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, RWORK, IWORK, IFAIL, INFO)
ZHPEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
Definition: zhpevx.f:242
zhbev
subroutine zhbev(JOBZ, UPLO, N, KD, AB, LDAB, W, Z, LDZ, WORK, RWORK, INFO)
ZHBEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrices
Definition: zhbev.f:154
zupgtr
subroutine zupgtr(UPLO, N, AP, TAU, Q, LDQ, WORK, INFO)
ZUPGTR
Definition: zupgtr.f:116
zhbevd_2stage
subroutine zhbevd_2stage(JOBZ, UPLO, N, KD, AB, LDAB, W, Z, LDZ, WORK, LWORK, RWORK, LRWORK, IWORK, LIWORK, INFO)
ZHBEVD_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER ...
Definition: zhbevd_2stage.f:262
zheev_2stage
subroutine zheev_2stage(JOBZ, UPLO, N, A, LDA, W, WORK, LWORK, RWORK, INFO)
ZHEEV_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for HE matr...
Definition: zheev_2stage.f:191
zhetrd_2stage
subroutine zhetrd_2stage(VECT, UPLO, N, A, LDA, D, E, TAU, HOUS2, LHOUS2, WORK, LWORK, INFO)
ZHETRD_2STAGE
Definition: zhetrd_2stage.f:226
zhpevd
subroutine zhpevd(JOBZ, UPLO, N, AP, W, Z, LDZ, WORK, LWORK, RWORK, LRWORK, IWORK, LIWORK, INFO)
ZHPEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
Definition: zhpevd.f:202