LAPACK  3.9.0
LAPACK: Linear Algebra PACKage
zerrvx.f
Go to the documentation of this file.
1 *> \brief \b ZERRVX
2 *
3 * =========== DOCUMENTATION ===========
4 *
5 * Online html documentation available at
6 * http://www.netlib.org/lapack/explore-html/
7 *
8 * Definition:
9 * ===========
10 *
11 * SUBROUTINE ZERRVX( PATH, NUNIT )
12 *
13 * .. Scalar Arguments ..
14 * CHARACTER*3 PATH
15 * INTEGER NUNIT
16 * ..
17 *
18 *
19 *> \par Purpose:
20 * =============
21 *>
22 *> \verbatim
23 *>
24 *> ZERRVX tests the error exits for the COMPLEX*16 driver routines
25 *> for solving linear systems of equations.
26 *> \endverbatim
27 *
28 * Arguments:
29 * ==========
30 *
31 *> \param[in] PATH
32 *> \verbatim
33 *> PATH is CHARACTER*3
34 *> The LAPACK path name for the routines to be tested.
35 *> \endverbatim
36 *>
37 *> \param[in] NUNIT
38 *> \verbatim
39 *> NUNIT is INTEGER
40 *> The unit number for output.
41 *> \endverbatim
42 *
43 * Authors:
44 * ========
45 *
46 *> \author Univ. of Tennessee
47 *> \author Univ. of California Berkeley
48 *> \author Univ. of Colorado Denver
49 *> \author NAG Ltd.
50 *
51 *> \date November 2017
52 *
53 *> \ingroup complex16_lin
54 *
55 * =====================================================================
56  SUBROUTINE zerrvx( PATH, NUNIT )
57 *
58 * -- LAPACK test routine (version 3.8.0) --
59 * -- LAPACK is a software package provided by Univ. of Tennessee, --
60 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
61 * November 2017
62 *
63 * .. Scalar Arguments ..
64  CHARACTER*3 PATH
65  INTEGER NUNIT
66 * ..
67 *
68 * =====================================================================
69 *
70 * .. Parameters ..
71  INTEGER NMAX
72  parameter( nmax = 4 )
73 * ..
74 * .. Local Scalars ..
75  CHARACTER EQ
76  CHARACTER*2 C2
77  INTEGER I, INFO, J
78  DOUBLE PRECISION RCOND
79 * ..
80 * .. Local Arrays ..
81  INTEGER IP( NMAX )
82  DOUBLE PRECISION C( NMAX ), R( NMAX ), R1( NMAX ), R2( NMAX ),
83  $ RF( NMAX ), RW( NMAX )
84  COMPLEX*16 A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
85  $ E( NMAX ), W( 2*NMAX ), X( NMAX )
86 * ..
87 * .. External Functions ..
88  LOGICAL LSAMEN
89  EXTERNAL lsamen
90 * ..
91 * .. External Subroutines ..
92  EXTERNAL chkxer, zgbsv, zgbsvx, zgesv, zgesvx, zgtsv,
93  $ zgtsvx, zhesv, zhesv_rk, zhesv_rook, zhesvx,
94  $ zhpsv, zhpsvx, zpbsv, zpbsvx, zposv, zposvx,
95  $ zppsv, zppsvx, zptsv, zptsvx, zspsv, zspsvx,
96  $ zsysv, zsysv_aa, zsysv_rk, zsysv_rook,
97  $ zsysvx, zhesv_aa_2stage
98 * ..
99 * .. Scalars in Common ..
100  LOGICAL LERR, OK
101  CHARACTER*32 SRNAMT
102  INTEGER INFOT, NOUT
103 * ..
104 * .. Common blocks ..
105  COMMON / infoc / infot, nout, ok, lerr
106  COMMON / srnamc / srnamt
107 * ..
108 * .. Intrinsic Functions ..
109  INTRINSIC dble, dcmplx
110 * ..
111 * .. Executable Statements ..
112 *
113  nout = nunit
114  WRITE( nout, fmt = * )
115  c2 = path( 2: 3 )
116 *
117 * Set the variables to innocuous values.
118 *
119  DO 20 j = 1, nmax
120  DO 10 i = 1, nmax
121  a( i, j ) = dcmplx( 1.d0 / dble( i+j ),
122  $ -1.d0 / dble( i+j ) )
123  af( i, j ) = dcmplx( 1.d0 / dble( i+j ),
124  $ -1.d0 / dble( i+j ) )
125  10 CONTINUE
126  b( j ) = 0.d0
127  e( j ) = 0.d0
128  r1( j ) = 0.d0
129  r2( j ) = 0.d0
130  w( j ) = 0.d0
131  x( j ) = 0.d0
132  c( j ) = 0.d0
133  r( j ) = 0.d0
134  ip( j ) = j
135  20 CONTINUE
136  eq = ' '
137  ok = .true.
138 *
139  IF( lsamen( 2, c2, 'GE' ) ) THEN
140 *
141 * ZGESV
142 *
143  srnamt = 'ZGESV '
144  infot = 1
145  CALL zgesv( -1, 0, a, 1, ip, b, 1, info )
146  CALL chkxer( 'ZGESV ', infot, nout, lerr, ok )
147  infot = 2
148  CALL zgesv( 0, -1, a, 1, ip, b, 1, info )
149  CALL chkxer( 'ZGESV ', infot, nout, lerr, ok )
150  infot = 4
151  CALL zgesv( 2, 1, a, 1, ip, b, 2, info )
152  CALL chkxer( 'ZGESV ', infot, nout, lerr, ok )
153  infot = 7
154  CALL zgesv( 2, 1, a, 2, ip, b, 1, info )
155  CALL chkxer( 'ZGESV ', infot, nout, lerr, ok )
156 *
157 * ZGESVX
158 *
159  srnamt = 'ZGESVX'
160  infot = 1
161  CALL zgesvx( '/', 'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
162  $ x, 1, rcond, r1, r2, w, rw, info )
163  CALL chkxer( 'ZGESVX', infot, nout, lerr, ok )
164  infot = 2
165  CALL zgesvx( 'N', '/', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
166  $ x, 1, rcond, r1, r2, w, rw, info )
167  CALL chkxer( 'ZGESVX', infot, nout, lerr, ok )
168  infot = 3
169  CALL zgesvx( 'N', 'N', -1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
170  $ x, 1, rcond, r1, r2, w, rw, info )
171  CALL chkxer( 'ZGESVX', infot, nout, lerr, ok )
172  infot = 4
173  CALL zgesvx( 'N', 'N', 0, -1, a, 1, af, 1, ip, eq, r, c, b, 1,
174  $ x, 1, rcond, r1, r2, w, rw, info )
175  CALL chkxer( 'ZGESVX', infot, nout, lerr, ok )
176  infot = 6
177  CALL zgesvx( 'N', 'N', 2, 1, a, 1, af, 2, ip, eq, r, c, b, 2,
178  $ x, 2, rcond, r1, r2, w, rw, info )
179  CALL chkxer( 'ZGESVX', infot, nout, lerr, ok )
180  infot = 8
181  CALL zgesvx( 'N', 'N', 2, 1, a, 2, af, 1, ip, eq, r, c, b, 2,
182  $ x, 2, rcond, r1, r2, w, rw, info )
183  CALL chkxer( 'ZGESVX', infot, nout, lerr, ok )
184  infot = 10
185  eq = '/'
186  CALL zgesvx( 'F', 'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
187  $ x, 1, rcond, r1, r2, w, rw, info )
188  CALL chkxer( 'ZGESVX', infot, nout, lerr, ok )
189  infot = 11
190  eq = 'R'
191  CALL zgesvx( 'F', 'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
192  $ x, 1, rcond, r1, r2, w, rw, info )
193  CALL chkxer( 'ZGESVX', infot, nout, lerr, ok )
194  infot = 12
195  eq = 'C'
196  CALL zgesvx( 'F', 'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
197  $ x, 1, rcond, r1, r2, w, rw, info )
198  CALL chkxer( 'ZGESVX', infot, nout, lerr, ok )
199  infot = 14
200  CALL zgesvx( 'N', 'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 1,
201  $ x, 2, rcond, r1, r2, w, rw, info )
202  CALL chkxer( 'ZGESVX', infot, nout, lerr, ok )
203  infot = 16
204  CALL zgesvx( 'N', 'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 2,
205  $ x, 1, rcond, r1, r2, w, rw, info )
206  CALL chkxer( 'ZGESVX', infot, nout, lerr, ok )
207 *
208  ELSE IF( lsamen( 2, c2, 'GB' ) ) THEN
209 *
210 * ZGBSV
211 *
212  srnamt = 'ZGBSV '
213  infot = 1
214  CALL zgbsv( -1, 0, 0, 0, a, 1, ip, b, 1, info )
215  CALL chkxer( 'ZGBSV ', infot, nout, lerr, ok )
216  infot = 2
217  CALL zgbsv( 1, -1, 0, 0, a, 1, ip, b, 1, info )
218  CALL chkxer( 'ZGBSV ', infot, nout, lerr, ok )
219  infot = 3
220  CALL zgbsv( 1, 0, -1, 0, a, 1, ip, b, 1, info )
221  CALL chkxer( 'ZGBSV ', infot, nout, lerr, ok )
222  infot = 4
223  CALL zgbsv( 0, 0, 0, -1, a, 1, ip, b, 1, info )
224  CALL chkxer( 'ZGBSV ', infot, nout, lerr, ok )
225  infot = 6
226  CALL zgbsv( 1, 1, 1, 0, a, 3, ip, b, 1, info )
227  CALL chkxer( 'ZGBSV ', infot, nout, lerr, ok )
228  infot = 9
229  CALL zgbsv( 2, 0, 0, 0, a, 1, ip, b, 1, info )
230  CALL chkxer( 'ZGBSV ', infot, nout, lerr, ok )
231 *
232 * ZGBSVX
233 *
234  srnamt = 'ZGBSVX'
235  infot = 1
236  CALL zgbsvx( '/', 'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
237  $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
238  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
239  infot = 2
240  CALL zgbsvx( 'N', '/', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
241  $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
242  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
243  infot = 3
244  CALL zgbsvx( 'N', 'N', -1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
245  $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
246  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
247  infot = 4
248  CALL zgbsvx( 'N', 'N', 1, -1, 0, 0, a, 1, af, 1, ip, eq, r, c,
249  $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
250  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
251  infot = 5
252  CALL zgbsvx( 'N', 'N', 1, 0, -1, 0, a, 1, af, 1, ip, eq, r, c,
253  $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
254  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
255  infot = 6
256  CALL zgbsvx( 'N', 'N', 0, 0, 0, -1, a, 1, af, 1, ip, eq, r, c,
257  $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
258  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
259  infot = 8
260  CALL zgbsvx( 'N', 'N', 1, 1, 1, 0, a, 2, af, 4, ip, eq, r, c,
261  $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
262  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
263  infot = 10
264  CALL zgbsvx( 'N', 'N', 1, 1, 1, 0, a, 3, af, 3, ip, eq, r, c,
265  $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
266  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
267  infot = 12
268  eq = '/'
269  CALL zgbsvx( 'F', 'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
270  $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
271  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
272  infot = 13
273  eq = 'R'
274  CALL zgbsvx( 'F', 'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
275  $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
276  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
277  infot = 14
278  eq = 'C'
279  CALL zgbsvx( 'F', 'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
280  $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
281  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
282  infot = 16
283  CALL zgbsvx( 'N', 'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
284  $ b, 1, x, 2, rcond, r1, r2, w, rw, info )
285  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
286  infot = 18
287  CALL zgbsvx( 'N', 'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
288  $ b, 2, x, 1, rcond, r1, r2, w, rw, info )
289  CALL chkxer( 'ZGBSVX', infot, nout, lerr, ok )
290 *
291  ELSE IF( lsamen( 2, c2, 'GT' ) ) THEN
292 *
293 * ZGTSV
294 *
295  srnamt = 'ZGTSV '
296  infot = 1
297  CALL zgtsv( -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
298  $ info )
299  CALL chkxer( 'ZGTSV ', infot, nout, lerr, ok )
300  infot = 2
301  CALL zgtsv( 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
302  $ info )
303  CALL chkxer( 'ZGTSV ', infot, nout, lerr, ok )
304  infot = 7
305  CALL zgtsv( 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1, info )
306  CALL chkxer( 'ZGTSV ', infot, nout, lerr, ok )
307 *
308 * ZGTSVX
309 *
310  srnamt = 'ZGTSVX'
311  infot = 1
312  CALL zgtsvx( '/', 'N', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
313  $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
314  $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
315  CALL chkxer( 'ZGTSVX', infot, nout, lerr, ok )
316  infot = 2
317  CALL zgtsvx( 'N', '/', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
318  $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
319  $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
320  CALL chkxer( 'ZGTSVX', infot, nout, lerr, ok )
321  infot = 3
322  CALL zgtsvx( 'N', 'N', -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
323  $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
324  $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
325  CALL chkxer( 'ZGTSVX', infot, nout, lerr, ok )
326  infot = 4
327  CALL zgtsvx( 'N', 'N', 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
328  $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
329  $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
330  CALL chkxer( 'ZGTSVX', infot, nout, lerr, ok )
331  infot = 14
332  CALL zgtsvx( 'N', 'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
333  $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
334  $ ip, b, 1, x, 2, rcond, r1, r2, w, rw, info )
335  CALL chkxer( 'ZGTSVX', infot, nout, lerr, ok )
336  infot = 16
337  CALL zgtsvx( 'N', 'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
338  $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
339  $ ip, b, 2, x, 1, rcond, r1, r2, w, rw, info )
340  CALL chkxer( 'ZGTSVX', infot, nout, lerr, ok )
341 *
342  ELSE IF( lsamen( 2, c2, 'PO' ) ) THEN
343 *
344 * ZPOSV
345 *
346  srnamt = 'ZPOSV '
347  infot = 1
348  CALL zposv( '/', 0, 0, a, 1, b, 1, info )
349  CALL chkxer( 'ZPOSV ', infot, nout, lerr, ok )
350  infot = 2
351  CALL zposv( 'U', -1, 0, a, 1, b, 1, info )
352  CALL chkxer( 'ZPOSV ', infot, nout, lerr, ok )
353  infot = 3
354  CALL zposv( 'U', 0, -1, a, 1, b, 1, info )
355  CALL chkxer( 'ZPOSV ', infot, nout, lerr, ok )
356  infot = 5
357  CALL zposv( 'U', 2, 0, a, 1, b, 2, info )
358  CALL chkxer( 'ZPOSV ', infot, nout, lerr, ok )
359  infot = 7
360  CALL zposv( 'U', 2, 0, a, 2, b, 1, info )
361  CALL chkxer( 'ZPOSV ', infot, nout, lerr, ok )
362 *
363 * ZPOSVX
364 *
365  srnamt = 'ZPOSVX'
366  infot = 1
367  CALL zposvx( '/', 'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
368  $ rcond, r1, r2, w, rw, info )
369  CALL chkxer( 'ZPOSVX', infot, nout, lerr, ok )
370  infot = 2
371  CALL zposvx( 'N', '/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
372  $ rcond, r1, r2, w, rw, info )
373  CALL chkxer( 'ZPOSVX', infot, nout, lerr, ok )
374  infot = 3
375  CALL zposvx( 'N', 'U', -1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
376  $ rcond, r1, r2, w, rw, info )
377  CALL chkxer( 'ZPOSVX', infot, nout, lerr, ok )
378  infot = 4
379  CALL zposvx( 'N', 'U', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
380  $ rcond, r1, r2, w, rw, info )
381  CALL chkxer( 'ZPOSVX', infot, nout, lerr, ok )
382  infot = 6
383  CALL zposvx( 'N', 'U', 2, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
384  $ rcond, r1, r2, w, rw, info )
385  CALL chkxer( 'ZPOSVX', infot, nout, lerr, ok )
386  infot = 8
387  CALL zposvx( 'N', 'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
388  $ rcond, r1, r2, w, rw, info )
389  CALL chkxer( 'ZPOSVX', infot, nout, lerr, ok )
390  infot = 9
391  eq = '/'
392  CALL zposvx( 'F', 'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
393  $ rcond, r1, r2, w, rw, info )
394  CALL chkxer( 'ZPOSVX', infot, nout, lerr, ok )
395  infot = 10
396  eq = 'Y'
397  CALL zposvx( 'F', 'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
398  $ rcond, r1, r2, w, rw, info )
399  CALL chkxer( 'ZPOSVX', infot, nout, lerr, ok )
400  infot = 12
401  CALL zposvx( 'N', 'U', 2, 0, a, 2, af, 2, eq, c, b, 1, x, 2,
402  $ rcond, r1, r2, w, rw, info )
403  CALL chkxer( 'ZPOSVX', infot, nout, lerr, ok )
404  infot = 14
405  CALL zposvx( 'N', 'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
406  $ rcond, r1, r2, w, rw, info )
407  CALL chkxer( 'ZPOSVX', infot, nout, lerr, ok )
408 *
409  ELSE IF( lsamen( 2, c2, 'PP' ) ) THEN
410 *
411 * ZPPSV
412 *
413  srnamt = 'ZPPSV '
414  infot = 1
415  CALL zppsv( '/', 0, 0, a, b, 1, info )
416  CALL chkxer( 'ZPPSV ', infot, nout, lerr, ok )
417  infot = 2
418  CALL zppsv( 'U', -1, 0, a, b, 1, info )
419  CALL chkxer( 'ZPPSV ', infot, nout, lerr, ok )
420  infot = 3
421  CALL zppsv( 'U', 0, -1, a, b, 1, info )
422  CALL chkxer( 'ZPPSV ', infot, nout, lerr, ok )
423  infot = 6
424  CALL zppsv( 'U', 2, 0, a, b, 1, info )
425  CALL chkxer( 'ZPPSV ', infot, nout, lerr, ok )
426 *
427 * ZPPSVX
428 *
429  srnamt = 'ZPPSVX'
430  infot = 1
431  CALL zppsvx( '/', 'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
432  $ r1, r2, w, rw, info )
433  CALL chkxer( 'ZPPSVX', infot, nout, lerr, ok )
434  infot = 2
435  CALL zppsvx( 'N', '/', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
436  $ r1, r2, w, rw, info )
437  CALL chkxer( 'ZPPSVX', infot, nout, lerr, ok )
438  infot = 3
439  CALL zppsvx( 'N', 'U', -1, 0, a, af, eq, c, b, 1, x, 1, rcond,
440  $ r1, r2, w, rw, info )
441  CALL chkxer( 'ZPPSVX', infot, nout, lerr, ok )
442  infot = 4
443  CALL zppsvx( 'N', 'U', 0, -1, a, af, eq, c, b, 1, x, 1, rcond,
444  $ r1, r2, w, rw, info )
445  CALL chkxer( 'ZPPSVX', infot, nout, lerr, ok )
446  infot = 7
447  eq = '/'
448  CALL zppsvx( 'F', 'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
449  $ r1, r2, w, rw, info )
450  CALL chkxer( 'ZPPSVX', infot, nout, lerr, ok )
451  infot = 8
452  eq = 'Y'
453  CALL zppsvx( 'F', 'U', 1, 0, a, af, eq, c, b, 1, x, 1, rcond,
454  $ r1, r2, w, rw, info )
455  CALL chkxer( 'ZPPSVX', infot, nout, lerr, ok )
456  infot = 10
457  CALL zppsvx( 'N', 'U', 2, 0, a, af, eq, c, b, 1, x, 2, rcond,
458  $ r1, r2, w, rw, info )
459  CALL chkxer( 'ZPPSVX', infot, nout, lerr, ok )
460  infot = 12
461  CALL zppsvx( 'N', 'U', 2, 0, a, af, eq, c, b, 2, x, 1, rcond,
462  $ r1, r2, w, rw, info )
463  CALL chkxer( 'ZPPSVX', infot, nout, lerr, ok )
464 *
465  ELSE IF( lsamen( 2, c2, 'PB' ) ) THEN
466 *
467 * ZPBSV
468 *
469  srnamt = 'ZPBSV '
470  infot = 1
471  CALL zpbsv( '/', 0, 0, 0, a, 1, b, 1, info )
472  CALL chkxer( 'ZPBSV ', infot, nout, lerr, ok )
473  infot = 2
474  CALL zpbsv( 'U', -1, 0, 0, a, 1, b, 1, info )
475  CALL chkxer( 'ZPBSV ', infot, nout, lerr, ok )
476  infot = 3
477  CALL zpbsv( 'U', 1, -1, 0, a, 1, b, 1, info )
478  CALL chkxer( 'ZPBSV ', infot, nout, lerr, ok )
479  infot = 4
480  CALL zpbsv( 'U', 0, 0, -1, a, 1, b, 1, info )
481  CALL chkxer( 'ZPBSV ', infot, nout, lerr, ok )
482  infot = 6
483  CALL zpbsv( 'U', 1, 1, 0, a, 1, b, 2, info )
484  CALL chkxer( 'ZPBSV ', infot, nout, lerr, ok )
485  infot = 8
486  CALL zpbsv( 'U', 2, 0, 0, a, 1, b, 1, info )
487  CALL chkxer( 'ZPBSV ', infot, nout, lerr, ok )
488 *
489 * ZPBSVX
490 *
491  srnamt = 'ZPBSVX'
492  infot = 1
493  CALL zpbsvx( '/', 'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
494  $ rcond, r1, r2, w, rw, info )
495  CALL chkxer( 'ZPBSVX', infot, nout, lerr, ok )
496  infot = 2
497  CALL zpbsvx( 'N', '/', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
498  $ rcond, r1, r2, w, rw, info )
499  CALL chkxer( 'ZPBSVX', infot, nout, lerr, ok )
500  infot = 3
501  CALL zpbsvx( 'N', 'U', -1, 0, 0, a, 1, af, 1, eq, c, b, 1, x,
502  $ 1, rcond, r1, r2, w, rw, info )
503  CALL chkxer( 'ZPBSVX', infot, nout, lerr, ok )
504  infot = 4
505  CALL zpbsvx( 'N', 'U', 1, -1, 0, a, 1, af, 1, eq, c, b, 1, x,
506  $ 1, rcond, r1, r2, w, rw, info )
507  CALL chkxer( 'ZPBSVX', infot, nout, lerr, ok )
508  infot = 5
509  CALL zpbsvx( 'N', 'U', 0, 0, -1, a, 1, af, 1, eq, c, b, 1, x,
510  $ 1, rcond, r1, r2, w, rw, info )
511  CALL chkxer( 'ZPBSVX', infot, nout, lerr, ok )
512  infot = 7
513  CALL zpbsvx( 'N', 'U', 1, 1, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
514  $ rcond, r1, r2, w, rw, info )
515  CALL chkxer( 'ZPBSVX', infot, nout, lerr, ok )
516  infot = 9
517  CALL zpbsvx( 'N', 'U', 1, 1, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
518  $ rcond, r1, r2, w, rw, info )
519  CALL chkxer( 'ZPBSVX', infot, nout, lerr, ok )
520  infot = 10
521  eq = '/'
522  CALL zpbsvx( 'F', 'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
523  $ rcond, r1, r2, w, rw, info )
524  CALL chkxer( 'ZPBSVX', infot, nout, lerr, ok )
525  infot = 11
526  eq = 'Y'
527  CALL zpbsvx( 'F', 'U', 1, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
528  $ rcond, r1, r2, w, rw, info )
529  CALL chkxer( 'ZPBSVX', infot, nout, lerr, ok )
530  infot = 13
531  CALL zpbsvx( 'N', 'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 2,
532  $ rcond, r1, r2, w, rw, info )
533  CALL chkxer( 'ZPBSVX', infot, nout, lerr, ok )
534  infot = 15
535  CALL zpbsvx( 'N', 'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 2, x, 1,
536  $ rcond, r1, r2, w, rw, info )
537  CALL chkxer( 'ZPBSVX', infot, nout, lerr, ok )
538 *
539  ELSE IF( lsamen( 2, c2, 'PT' ) ) THEN
540 *
541 * ZPTSV
542 *
543  srnamt = 'ZPTSV '
544  infot = 1
545  CALL zptsv( -1, 0, r, a( 1, 1 ), b, 1, info )
546  CALL chkxer( 'ZPTSV ', infot, nout, lerr, ok )
547  infot = 2
548  CALL zptsv( 0, -1, r, a( 1, 1 ), b, 1, info )
549  CALL chkxer( 'ZPTSV ', infot, nout, lerr, ok )
550  infot = 6
551  CALL zptsv( 2, 0, r, a( 1, 1 ), b, 1, info )
552  CALL chkxer( 'ZPTSV ', infot, nout, lerr, ok )
553 *
554 * ZPTSVX
555 *
556  srnamt = 'ZPTSVX'
557  infot = 1
558  CALL zptsvx( '/', 0, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
559  $ 1, rcond, r1, r2, w, rw, info )
560  CALL chkxer( 'ZPTSVX', infot, nout, lerr, ok )
561  infot = 2
562  CALL zptsvx( 'N', -1, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
563  $ 1, rcond, r1, r2, w, rw, info )
564  CALL chkxer( 'ZPTSVX', infot, nout, lerr, ok )
565  infot = 3
566  CALL zptsvx( 'N', 0, -1, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
567  $ 1, rcond, r1, r2, w, rw, info )
568  CALL chkxer( 'ZPTSVX', infot, nout, lerr, ok )
569  infot = 9
570  CALL zptsvx( 'N', 2, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
571  $ 2, rcond, r1, r2, w, rw, info )
572  CALL chkxer( 'ZPTSVX', infot, nout, lerr, ok )
573  infot = 11
574  CALL zptsvx( 'N', 2, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 2, x,
575  $ 1, rcond, r1, r2, w, rw, info )
576  CALL chkxer( 'ZPTSVX', infot, nout, lerr, ok )
577 *
578  ELSE IF( lsamen( 2, c2, 'HE' ) ) THEN
579 *
580 * ZHESV
581 *
582  srnamt = 'ZHESV '
583  infot = 1
584  CALL zhesv( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
585  CALL chkxer( 'ZHESV ', infot, nout, lerr, ok )
586  infot = 2
587  CALL zhesv( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
588  CALL chkxer( 'ZHESV ', infot, nout, lerr, ok )
589  infot = 3
590  CALL zhesv( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
591  CALL chkxer( 'ZHESV ', infot, nout, lerr, ok )
592  infot = 5
593  CALL zhesv( 'U', 2, 0, a, 1, ip, b, 2, w, 1, info )
594  CALL chkxer( 'ZHESV ', infot, nout, lerr, ok )
595  infot = 8
596  CALL zhesv( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
597  CALL chkxer( 'ZHESV ', infot, nout, lerr, ok )
598  infot = 10
599  CALL zhesv( 'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
600  CALL chkxer( 'ZHESV ', infot, nout, lerr, ok )
601  infot = 10
602  CALL zhesv( 'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
603  CALL chkxer( 'ZHESV ', infot, nout, lerr, ok )
604 *
605 * ZHESVX
606 *
607  srnamt = 'ZHESVX'
608  infot = 1
609  CALL zhesvx( '/', 'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
610  $ rcond, r1, r2, w, 1, rw, info )
611  CALL chkxer( 'ZHESVX', infot, nout, lerr, ok )
612  infot = 2
613  CALL zhesvx( 'N', '/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
614  $ rcond, r1, r2, w, 1, rw, info )
615  CALL chkxer( 'ZHESVX', infot, nout, lerr, ok )
616  infot = 3
617  CALL zhesvx( 'N', 'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
618  $ rcond, r1, r2, w, 1, rw, info )
619  CALL chkxer( 'ZHESVX', infot, nout, lerr, ok )
620  infot = 4
621  CALL zhesvx( 'N', 'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
622  $ rcond, r1, r2, w, 1, rw, info )
623  CALL chkxer( 'ZHESVX', infot, nout, lerr, ok )
624  infot = 6
625  CALL zhesvx( 'N', 'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
626  $ rcond, r1, r2, w, 4, rw, info )
627  CALL chkxer( 'ZHESVX', infot, nout, lerr, ok )
628  infot = 8
629  CALL zhesvx( 'N', 'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
630  $ rcond, r1, r2, w, 4, rw, info )
631  CALL chkxer( 'ZHESVX', infot, nout, lerr, ok )
632  infot = 11
633  CALL zhesvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
634  $ rcond, r1, r2, w, 4, rw, info )
635  CALL chkxer( 'ZHESVX', infot, nout, lerr, ok )
636  infot = 13
637  CALL zhesvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
638  $ rcond, r1, r2, w, 4, rw, info )
639  CALL chkxer( 'ZHESVX', infot, nout, lerr, ok )
640  infot = 18
641  CALL zhesvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
642  $ rcond, r1, r2, w, 3, rw, info )
643  CALL chkxer( 'ZHESVX', infot, nout, lerr, ok )
644 *
645  ELSE IF( lsamen( 2, c2, 'HR' ) ) THEN
646 *
647 * ZHESV_ROOK
648 *
649  srnamt = 'ZHESV_ROOK'
650  infot = 1
651  CALL zhesv_rook( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
652  CALL chkxer( 'ZHESV_ROOK', infot, nout, lerr, ok )
653  infot = 2
654  CALL zhesv_rook( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
655  CALL chkxer( 'ZHESV_ROOK', infot, nout, lerr, ok )
656  infot = 3
657  CALL zhesv_rook( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
658  CALL chkxer( 'ZHESV_ROOK', infot, nout, lerr, ok )
659  infot = 8
660  CALL zhesv_rook( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
661  CALL chkxer( 'ZHESV_ROOK', infot, nout, lerr, ok )
662  infot = 10
663  CALL zhesv_rook( 'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
664  CALL chkxer( 'ZHESV_ROOK', infot, nout, lerr, ok )
665  infot = 10
666  CALL zhesv_rook( 'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
667  CALL chkxer( 'ZHESV_ROOK', infot, nout, lerr, ok )
668 *
669  ELSE IF( lsamen( 2, c2, 'HK' ) ) THEN
670 *
671 * ZSYSV_RK
672 *
673 * Test error exits of the driver that uses factorization
674 * of a Hermitian indefinite matrix with rook
675 * (bounded Bunch-Kaufman) pivoting with the new storage
676 * format for factors L ( or U) and D.
677 *
678 * L (or U) is stored in A, diagonal of D is stored on the
679 * diagonal of A, subdiagonal of D is stored in a separate array E.
680 *
681  srnamt = 'ZHESV_RK'
682  infot = 1
683  CALL zhesv_rk( '/', 0, 0, a, 1, e, ip, b, 1, w, 1, info )
684  CALL chkxer( 'ZHESV_RK', infot, nout, lerr, ok )
685  infot = 2
686  CALL zhesv_rk( 'U', -1, 0, a, 1, e, ip, b, 1, w, 1, info )
687  CALL chkxer( 'ZHESV_RK', infot, nout, lerr, ok )
688  infot = 3
689  CALL zhesv_rk( 'U', 0, -1, a, 1, e, ip, b, 1, w, 1, info )
690  CALL chkxer( 'ZHESV_RK', infot, nout, lerr, ok )
691  infot = 5
692  CALL zhesv_rk( 'U', 2, 0, a, 1, e, ip, b, 2, w, 1, info )
693  CALL chkxer( 'ZHESV_RK', infot, nout, lerr, ok )
694  infot = 9
695  CALL zhesv_rk( 'U', 2, 0, a, 2, e, ip, b, 1, w, 1, info )
696  CALL chkxer( 'ZHESV_RK', infot, nout, lerr, ok )
697  infot = 11
698  CALL zhesv_rk( 'U', 0, 0, a, 1, e, ip, b, 1, w, 0, info )
699  CALL chkxer( 'ZHESV_RK', infot, nout, lerr, ok )
700  infot = 11
701  CALL zhesv_rk( 'U', 0, 0, a, 1, e, ip, b, 1, w, -2, info )
702  CALL chkxer( 'ZHESV_RK', infot, nout, lerr, ok )
703 *
704  ELSE IF( lsamen( 2, c2, 'HA' ) ) THEN
705 *
706 * ZHESV_AA
707 *
708  srnamt = 'ZHESV_AA'
709  infot = 1
710  CALL zhesv_aa( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
711  CALL chkxer( 'ZHESV_AA', infot, nout, lerr, ok )
712  infot = 2
713  CALL zhesv_aa( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
714  CALL chkxer( 'ZHESV_AA', infot, nout, lerr, ok )
715  infot = 3
716  CALL zhesv_aa( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
717  CALL chkxer( 'ZHESV_AA', infot, nout, lerr, ok )
718  infot = 8
719  CALL zhesv_aa( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
720  CALL chkxer( 'ZHESV_AA', infot, nout, lerr, ok )
721 *
722  ELSE IF( lsamen( 2, c2, 'H2' ) ) THEN
723 *
724 * ZHESV_AASEN_2STAGE
725 *
726  srnamt = 'ZHESV_AA_2STAGE'
727  infot = 1
728  CALL zhesv_aa_2stage( '/', 0, 0, a, 1, a, 1, ip, ip, b, 1,
729  $ w, 1, info )
730  CALL chkxer( 'ZHESV_AA_2STAGE', infot, nout, lerr, ok )
731  infot = 2
732  CALL zhesv_aa_2stage( 'U', -1, 0, a, 1, a, 1, ip, ip, b, 1,
733  $ w, 1, info )
734  CALL chkxer( 'ZHESV_AA_2STAGE', infot, nout, lerr, ok )
735  infot = 3
736  CALL zhesv_aa_2stage( 'U', 0, -1, a, 1, a, 1, ip, ip, b, 1,
737  $ w, 1, info )
738  CALL chkxer( 'ZHESV_AA_2STAGE', infot, nout, lerr, ok )
739  infot = 5
740  CALL zhesv_aa_2stage( 'U', 2, 1, a, 1, a, 1, ip, ip, b, 1,
741  $ w, 1, info )
742  CALL chkxer( 'ZHESV_AA_2STAGE', infot, nout, lerr, ok )
743  infot = 11
744  CALL zhesv_aa_2stage( 'U', 2, 1, a, 2, a, 8, ip, ip, b, 1,
745  $ w, 1, info )
746  CALL chkxer( 'ZHESV_AA_2STAGE', infot, nout, lerr, ok )
747  infot = 7
748  CALL zhesv_aa_2stage( 'U', 2, 1, a, 2, a, 1, ip, ip, b, 2,
749  $ w, 1, info )
750  CALL chkxer( 'ZHESV_AA_2STAGE', infot, nout, lerr, ok )
751 *
752  ELSE IF( lsamen( 2, c2, 'S2' ) ) THEN
753 *
754 * ZSYSV_AASEN_2STAGE
755 *
756  srnamt = 'ZSYSV_AA_2STAGE'
757  infot = 1
758  CALL zsysv_aa_2stage( '/', 0, 0, a, 1, a, 1, ip, ip, b, 1,
759  $ w, 1, info )
760  CALL chkxer( 'ZSYSV_AA_2STAGE', infot, nout, lerr, ok )
761  infot = 2
762  CALL zsysv_aa_2stage( 'U', -1, 0, a, 1, a, 1, ip, ip, b, 1,
763  $ w, 1, info )
764  CALL chkxer( 'ZSYSV_AA_2STAGE', infot, nout, lerr, ok )
765  infot = 3
766  CALL zsysv_aa_2stage( 'U', 0, -1, a, 1, a, 1, ip, ip, b, 1,
767  $ w, 1, info )
768  CALL chkxer( 'ZSYSV_AA_2STAGE', infot, nout, lerr, ok )
769  infot = 5
770  CALL zsysv_aa_2stage( 'U', 2, 1, a, 1, a, 1, ip, ip, b, 1,
771  $ w, 1, info )
772  CALL chkxer( 'ZSYSV_AA_2STAGE', infot, nout, lerr, ok )
773  infot = 11
774  CALL zsysv_aa_2stage( 'U', 2, 1, a, 2, a, 8, ip, ip, b, 1,
775  $ w, 1, info )
776  CALL chkxer( 'ZSYSV_AA_2STAGE', infot, nout, lerr, ok )
777  infot = 7
778  CALL zsysv_aa_2stage( 'U', 2, 1, a, 2, a, 1, ip, ip, b, 2,
779  $ w, 1, info )
780  CALL chkxer( 'ZSYSV_AA_2STAGE', infot, nout, lerr, ok )
781 **
782  ELSE IF( lsamen( 2, c2, 'HP' ) ) THEN
783 *
784 * ZHPSV
785 *
786  srnamt = 'ZHPSV '
787  infot = 1
788  CALL zhpsv( '/', 0, 0, a, ip, b, 1, info )
789  CALL chkxer( 'ZHPSV ', infot, nout, lerr, ok )
790  infot = 2
791  CALL zhpsv( 'U', -1, 0, a, ip, b, 1, info )
792  CALL chkxer( 'ZHPSV ', infot, nout, lerr, ok )
793  infot = 3
794  CALL zhpsv( 'U', 0, -1, a, ip, b, 1, info )
795  CALL chkxer( 'ZHPSV ', infot, nout, lerr, ok )
796  infot = 7
797  CALL zhpsv( 'U', 2, 0, a, ip, b, 1, info )
798  CALL chkxer( 'ZHPSV ', infot, nout, lerr, ok )
799 *
800 * ZHPSVX
801 *
802  srnamt = 'ZHPSVX'
803  infot = 1
804  CALL zhpsvx( '/', 'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
805  $ r2, w, rw, info )
806  CALL chkxer( 'ZHPSVX', infot, nout, lerr, ok )
807  infot = 2
808  CALL zhpsvx( 'N', '/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
809  $ r2, w, rw, info )
810  CALL chkxer( 'ZHPSVX', infot, nout, lerr, ok )
811  infot = 3
812  CALL zhpsvx( 'N', 'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
813  $ r2, w, rw, info )
814  CALL chkxer( 'ZHPSVX', infot, nout, lerr, ok )
815  infot = 4
816  CALL zhpsvx( 'N', 'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
817  $ r2, w, rw, info )
818  CALL chkxer( 'ZHPSVX', infot, nout, lerr, ok )
819  infot = 9
820  CALL zhpsvx( 'N', 'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
821  $ r2, w, rw, info )
822  CALL chkxer( 'ZHPSVX', infot, nout, lerr, ok )
823  infot = 11
824  CALL zhpsvx( 'N', 'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
825  $ r2, w, rw, info )
826  CALL chkxer( 'ZHPSVX', infot, nout, lerr, ok )
827 *
828  ELSE IF( lsamen( 2, c2, 'SY' ) ) THEN
829 *
830 * ZSYSV
831 *
832  srnamt = 'ZSYSV '
833  infot = 1
834  CALL zsysv( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
835  CALL chkxer( 'ZSYSV ', infot, nout, lerr, ok )
836  infot = 2
837  CALL zsysv( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
838  CALL chkxer( 'ZSYSV ', infot, nout, lerr, ok )
839  infot = 3
840  CALL zsysv( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
841  CALL chkxer( 'ZSYSV ', infot, nout, lerr, ok )
842  infot = 8
843  CALL zsysv( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
844  CALL chkxer( 'ZSYSV ', infot, nout, lerr, ok )
845  infot = 10
846  CALL zsysv( 'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
847  CALL chkxer( 'ZSYSV ', infot, nout, lerr, ok )
848  infot = 10
849  CALL zsysv( 'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
850  CALL chkxer( 'ZSYSV ', infot, nout, lerr, ok )
851 *
852 * ZSYSVX
853 *
854  srnamt = 'ZSYSVX'
855  infot = 1
856  CALL zsysvx( '/', 'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
857  $ rcond, r1, r2, w, 1, rw, info )
858  CALL chkxer( 'ZSYSVX', infot, nout, lerr, ok )
859  infot = 2
860  CALL zsysvx( 'N', '/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
861  $ rcond, r1, r2, w, 1, rw, info )
862  CALL chkxer( 'ZSYSVX', infot, nout, lerr, ok )
863  infot = 3
864  CALL zsysvx( 'N', 'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
865  $ rcond, r1, r2, w, 1, rw, info )
866  CALL chkxer( 'ZSYSVX', infot, nout, lerr, ok )
867  infot = 4
868  CALL zsysvx( 'N', 'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
869  $ rcond, r1, r2, w, 1, rw, info )
870  CALL chkxer( 'ZSYSVX', infot, nout, lerr, ok )
871  infot = 6
872  CALL zsysvx( 'N', 'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
873  $ rcond, r1, r2, w, 4, rw, info )
874  CALL chkxer( 'ZSYSVX', infot, nout, lerr, ok )
875  infot = 8
876  CALL zsysvx( 'N', 'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
877  $ rcond, r1, r2, w, 4, rw, info )
878  CALL chkxer( 'ZSYSVX', infot, nout, lerr, ok )
879  infot = 11
880  CALL zsysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
881  $ rcond, r1, r2, w, 4, rw, info )
882  CALL chkxer( 'ZSYSVX', infot, nout, lerr, ok )
883  infot = 13
884  CALL zsysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
885  $ rcond, r1, r2, w, 4, rw, info )
886  CALL chkxer( 'ZSYSVX', infot, nout, lerr, ok )
887  infot = 18
888  CALL zsysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
889  $ rcond, r1, r2, w, 3, rw, info )
890  CALL chkxer( 'ZSYSVX', infot, nout, lerr, ok )
891 *
892  ELSE IF( lsamen( 2, c2, 'SR' ) ) THEN
893 *
894 * ZSYSV_ROOK
895 *
896  srnamt = 'ZSYSV_ROOK'
897  infot = 1
898  CALL zsysv_rook( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
899  CALL chkxer( 'ZSYSV_ROOK', infot, nout, lerr, ok )
900  infot = 2
901  CALL zsysv_rook( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
902  CALL chkxer( 'ZSYSV_ROOK', infot, nout, lerr, ok )
903  infot = 3
904  CALL zsysv_rook( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
905  CALL chkxer( 'ZSYSV_ROOK', infot, nout, lerr, ok )
906  infot = 8
907  CALL zsysv_rook( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
908  CALL chkxer( 'ZSYSV_ROOK', infot, nout, lerr, ok )
909  infot = 10
910  CALL zsysv_rook( 'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
911  CALL chkxer( 'ZSYSV_ROOK', infot, nout, lerr, ok )
912  infot = 10
913  CALL zsysv_rook( 'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
914 *
915  ELSE IF( lsamen( 2, c2, 'SK' ) ) THEN
916 *
917 * ZSYSV_RK
918 *
919 * Test error exits of the driver that uses factorization
920 * of a symmetric indefinite matrix with rook
921 * (bounded Bunch-Kaufman) pivoting with the new storage
922 * format for factors L ( or U) and D.
923 *
924 * L (or U) is stored in A, diagonal of D is stored on the
925 * diagonal of A, subdiagonal of D is stored in a separate array E.
926 *
927  srnamt = 'ZSYSV_RK'
928  infot = 1
929  CALL zsysv_rk( '/', 0, 0, a, 1, e, ip, b, 1, w, 1, info )
930  CALL chkxer( 'ZSYSV_RK', infot, nout, lerr, ok )
931  infot = 2
932  CALL zsysv_rk( 'U', -1, 0, a, 1, e, ip, b, 1, w, 1, info )
933  CALL chkxer( 'ZSYSV_RK', infot, nout, lerr, ok )
934  infot = 3
935  CALL zsysv_rk( 'U', 0, -1, a, 1, e, ip, b, 1, w, 1, info )
936  CALL chkxer( 'ZSYSV_RK', infot, nout, lerr, ok )
937  infot = 5
938  CALL zsysv_rk( 'U', 2, 0, a, 1, e, ip, b, 2, w, 1, info )
939  CALL chkxer( 'ZSYSV_RK', infot, nout, lerr, ok )
940  infot = 9
941  CALL zsysv_rk( 'U', 2, 0, a, 2, e, ip, b, 1, w, 1, info )
942  CALL chkxer( 'ZSYSV_RK', infot, nout, lerr, ok )
943  infot = 11
944  CALL zsysv_rk( 'U', 0, 0, a, 1, e, ip, b, 1, w, 0, info )
945  CALL chkxer( 'ZSYSV_RK', infot, nout, lerr, ok )
946  infot = 11
947  CALL zsysv_rk( 'U', 0, 0, a, 1, e, ip, b, 1, w, -2, info )
948  CALL chkxer( 'ZSYSV_RK', infot, nout, lerr, ok )
949 *
950  ELSE IF( lsamen( 2, c2, 'SP' ) ) THEN
951 *
952 * ZSPSV
953 *
954  srnamt = 'ZSPSV '
955  infot = 1
956  CALL zspsv( '/', 0, 0, a, ip, b, 1, info )
957  CALL chkxer( 'ZSPSV ', infot, nout, lerr, ok )
958  infot = 2
959  CALL zspsv( 'U', -1, 0, a, ip, b, 1, info )
960  CALL chkxer( 'ZSPSV ', infot, nout, lerr, ok )
961  infot = 3
962  CALL zspsv( 'U', 0, -1, a, ip, b, 1, info )
963  CALL chkxer( 'ZSPSV ', infot, nout, lerr, ok )
964  infot = 7
965  CALL zspsv( 'U', 2, 0, a, ip, b, 1, info )
966  CALL chkxer( 'ZSPSV ', infot, nout, lerr, ok )
967 *
968 * ZSPSVX
969 *
970  srnamt = 'ZSPSVX'
971  infot = 1
972  CALL zspsvx( '/', 'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
973  $ r2, w, rw, info )
974  CALL chkxer( 'ZSPSVX', infot, nout, lerr, ok )
975  infot = 2
976  CALL zspsvx( 'N', '/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
977  $ r2, w, rw, info )
978  CALL chkxer( 'ZSPSVX', infot, nout, lerr, ok )
979  infot = 3
980  CALL zspsvx( 'N', 'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
981  $ r2, w, rw, info )
982  CALL chkxer( 'ZSPSVX', infot, nout, lerr, ok )
983  infot = 4
984  CALL zspsvx( 'N', 'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
985  $ r2, w, rw, info )
986  CALL chkxer( 'ZSPSVX', infot, nout, lerr, ok )
987  infot = 9
988  CALL zspsvx( 'N', 'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
989  $ r2, w, rw, info )
990  CALL chkxer( 'ZSPSVX', infot, nout, lerr, ok )
991  infot = 11
992  CALL zspsvx( 'N', 'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
993  $ r2, w, rw, info )
994  CALL chkxer( 'ZSPSVX', infot, nout, lerr, ok )
995  END IF
996 *
997 * Print a summary line.
998 *
999  IF( ok ) THEN
1000  WRITE( nout, fmt = 9999 )path
1001  ELSE
1002  WRITE( nout, fmt = 9998 )path
1003  END IF
1004 *
1005  9999 FORMAT( 1x, a3, ' drivers passed the tests of the error exits' )
1006  9998 FORMAT( ' *** ', a3, ' drivers failed the tests of the error ',
1007  $ 'exits ***' )
1008 *
1009  RETURN
1010 *
1011 * End of ZERRVX
1012 *
1013  END
zgesv
subroutine zgesv(N, NRHS, A, LDA, IPIV, B, LDB, INFO)
ZGESV computes the solution to system of linear equations A * X = B for GE matrices (simple driver)
Definition: zgesv.f:124
zposvx
subroutine zposvx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, EQUED, S, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
ZPOSVX computes the solution to system of linear equations A * X = B for PO matrices
Definition: zposvx.f:308
zerrvx
subroutine zerrvx(PATH, NUNIT)
ZERRVX
Definition: zerrvx.f:57
zhpsv
subroutine zhpsv(UPLO, N, NRHS, AP, IPIV, B, LDB, INFO)
ZHPSV computes the solution to system of linear equations A * X = B for OTHER matrices
Definition: zhpsv.f:164
zhesv_aa_2stage
subroutine zhesv_aa_2stage(UPLO, N, NRHS, A, LDA, TB, LTB, IPIV, IPIV2, B, LDB, WORK, LWORK, INFO)
ZHESV_AA_2STAGE computes the solution to system of linear equations A * X = B for HE matrices
Definition: zhesv_aa_2stage.f:189
zptsvx
subroutine zptsvx(FACT, N, NRHS, D, E, DF, EF, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
ZPTSVX computes the solution to system of linear equations A * X = B for PT matrices
Definition: zptsvx.f:236
zhesv_aa
subroutine zhesv_aa(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
ZHESV_AA computes the solution to system of linear equations A * X = B for HE matrices
Definition: zhesv_aa.f:164
zgbsv
subroutine zgbsv(N, KL, KU, NRHS, AB, LDAB, IPIV, B, LDB, INFO)
ZGBSV computes the solution to system of linear equations A * X = B for GB matrices (simple driver)
Definition: zgbsv.f:164
zhesv
subroutine zhesv(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
ZHESV computes the solution to system of linear equations A * X = B for HE matrices
Definition: zhesv.f:173
zhesv_rk
subroutine zhesv_rk(UPLO, N, NRHS, A, LDA, E, IPIV, B, LDB, WORK, LWORK, INFO)
ZHESV_RK computes the solution to system of linear equations A * X = B for SY matrices
Definition: zhesv_rk.f:230
zsysv_aa
subroutine zsysv_aa(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
ZSYSV_AA computes the solution to system of linear equations A * X = B for SY matrices
Definition: zsysv_aa.f:164
zsysv
subroutine zsysv(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
ZSYSV computes the solution to system of linear equations A * X = B for SY matrices
Definition: zsysv.f:173
zgbsvx
subroutine zgbsvx(FACT, TRANS, N, KL, KU, NRHS, AB, LDAB, AFB, LDAFB, IPIV, EQUED, R, C, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
ZGBSVX computes the solution to system of linear equations A * X = B for GB matrices
Definition: zgbsvx.f:372
zgesvx
subroutine zgesvx(FACT, TRANS, N, NRHS, A, LDA, AF, LDAF, IPIV, EQUED, R, C, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
ZGESVX computes the solution to system of linear equations A * X = B for GE matrices
Definition: zgesvx.f:352
zhpsvx
subroutine zhpsvx(FACT, UPLO, N, NRHS, AP, AFP, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
ZHPSVX computes the solution to system of linear equations A * X = B for OTHER matrices
Definition: zhpsvx.f:279
zgtsv
subroutine zgtsv(N, NRHS, DL, D, DU, B, LDB, INFO)
ZGTSV computes the solution to system of linear equations A * X = B for GT matrices
Definition: zgtsv.f:126
chkxer
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
Definition: cblat2.f:3199
zhesvx
subroutine zhesvx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, LWORK, RWORK, INFO)
ZHESVX computes the solution to system of linear equations A * X = B for HE matrices
Definition: zhesvx.f:287
zgtsvx
subroutine zgtsvx(FACT, TRANS, N, NRHS, DL, D, DU, DLF, DF, DUF, DU2, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
ZGTSVX computes the solution to system of linear equations A * X = B for GT matrices
Definition: zgtsvx.f:296
zhesv_rook
subroutine zhesv_rook(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
ZHESV_ROOK computes the solution to a system of linear equations A * X = B for HE matrices using the ...
Definition: zhesv_rook.f:207
zspsv
subroutine zspsv(UPLO, N, NRHS, AP, IPIV, B, LDB, INFO)
ZSPSV computes the solution to system of linear equations A * X = B for OTHER matrices
Definition: zspsv.f:164
zsysvx
subroutine zsysvx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, LWORK, RWORK, INFO)
ZSYSVX computes the solution to system of linear equations A * X = B for SY matrices
Definition: zsysvx.f:287
zpbsv
subroutine zpbsv(UPLO, N, KD, NRHS, AB, LDAB, B, LDB, INFO)
ZPBSV computes the solution to system of linear equations A * X = B for OTHER matrices
Definition: zpbsv.f:166
zppsvx
subroutine zppsvx(FACT, UPLO, N, NRHS, AP, AFP, EQUED, S, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
ZPPSVX computes the solution to system of linear equations A * X = B for OTHER matrices
Definition: zppsvx.f:313
zppsv
subroutine zppsv(UPLO, N, NRHS, AP, B, LDB, INFO)
ZPPSV computes the solution to system of linear equations A * X = B for OTHER matrices
Definition: zppsv.f:146
zsysv_aa_2stage
subroutine zsysv_aa_2stage(UPLO, N, NRHS, A, LDA, TB, LTB, IPIV, IPIV2, B, LDB, WORK, LWORK, INFO)
ZSYSV_AA_2STAGE computes the solution to system of linear equations A * X = B for SY matrices
Definition: zsysv_aa_2stage.f:188
zsysv_rook
subroutine zsysv_rook(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
ZSYSV_ROOK computes the solution to system of linear equations A * X = B for SY matrices
Definition: zsysv_rook.f:206
zpbsvx
subroutine zpbsvx(FACT, UPLO, N, KD, NRHS, AB, LDAB, AFB, LDAFB, EQUED, S, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
ZPBSVX computes the solution to system of linear equations A * X = B for OTHER matrices
Definition: zpbsvx.f:344
zptsv
subroutine zptsv(N, NRHS, D, E, B, LDB, INFO)
ZPTSV computes the solution to system of linear equations A * X = B for PT matrices
Definition: zptsv.f:117
zposv
subroutine zposv(UPLO, N, NRHS, A, LDA, B, LDB, INFO)
ZPOSV computes the solution to system of linear equations A * X = B for PO matrices
Definition: zposv.f:132
zsysv_rk
subroutine zsysv_rk(UPLO, N, NRHS, A, LDA, E, IPIV, B, LDB, WORK, LWORK, INFO)
ZSYSV_RK computes the solution to system of linear equations A * X = B for SY matrices
Definition: zsysv_rk.f:230
zspsvx
subroutine zspsvx(FACT, UPLO, N, NRHS, AP, AFP, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
ZSPSVX computes the solution to system of linear equations A * X = B for OTHER matrices
Definition: zspsvx.f:279