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