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