56 SUBROUTINE zerrvx( PATH, NUNIT )
72 parameter ( nmax = 4 )
74 parameter ( one = 1.0d+0 )
79 INTEGER i, info, j, n_err_bnds, nparams
80 DOUBLE PRECISION rcond, rpvgrw, berr
84 DOUBLE PRECISION c( nmax ), r( nmax ), r1( nmax ), r2( nmax ),
85 $ rf( nmax ), rw( nmax ), err_bnds_n( nmax, 3 ),
86 $ err_bnds_c( nmax, 3 ), params( 1 )
87 COMPLEX*16 a( nmax, nmax ), af( nmax, nmax ), b( nmax ),
88 $ w( 2*nmax ), x( nmax )
108 COMMON / infoc / infot, nout, ok, lerr
109 COMMON / srnamc / srnamt
112 INTRINSIC dble, dcmplx
117 WRITE( nout, fmt = * )
124 a( i, j ) = dcmplx( 1.d0 / dble( i+j ),
125 $ -1.d0 / dble( i+j ) )
126 af( i, j ) = dcmplx( 1.d0 / dble( i+j ),
127 $ -1.d0 / dble( i+j ) )
141 IF(
lsamen( 2, c2,
'GE' ) )
THEN
147 CALL zgesv( -1, 0, a, 1, ip, b, 1, info )
148 CALL chkxer(
'ZGESV ', infot, nout, lerr, ok )
150 CALL zgesv( 0, -1, a, 1, ip, b, 1, info )
151 CALL chkxer(
'ZGESV ', infot, nout, lerr, ok )
153 CALL zgesv( 2, 1, a, 1, ip, b, 2, info )
154 CALL chkxer(
'ZGESV ', infot, nout, lerr, ok )
156 CALL zgesv( 2, 1, a, 2, ip, b, 1, info )
157 CALL chkxer(
'ZGESV ', infot, nout, lerr, ok )
163 CALL zgesvx(
'/',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
164 $ x, 1, rcond, r1, r2, w, rw, info )
165 CALL chkxer(
'ZGESVX', infot, nout, lerr, ok )
167 CALL zgesvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
168 $ x, 1, rcond, r1, r2, w, rw, info )
169 CALL chkxer(
'ZGESVX', infot, nout, lerr, ok )
171 CALL zgesvx(
'N',
'N', -1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
172 $ x, 1, rcond, r1, r2, w, rw, info )
173 CALL chkxer(
'ZGESVX', infot, nout, lerr, ok )
175 CALL zgesvx(
'N',
'N', 0, -1, a, 1, af, 1, ip, eq, r, c, b, 1,
176 $ x, 1, rcond, r1, r2, w, rw, info )
177 CALL chkxer(
'ZGESVX', infot, nout, lerr, ok )
179 CALL zgesvx(
'N',
'N', 2, 1, a, 1, af, 2, ip, eq, r, c, b, 2,
180 $ x, 2, rcond, r1, r2, w, rw, info )
181 CALL chkxer(
'ZGESVX', infot, nout, lerr, ok )
183 CALL zgesvx(
'N',
'N', 2, 1, a, 2, af, 1, ip, eq, r, c, b, 2,
184 $ x, 2, rcond, r1, r2, w, rw, info )
185 CALL chkxer(
'ZGESVX', infot, nout, lerr, ok )
188 CALL zgesvx(
'F',
'N', 0, 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 )
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 )
198 CALL zgesvx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
199 $ x, 1, rcond, r1, r2, w, rw, info )
200 CALL chkxer(
'ZGESVX', infot, nout, lerr, ok )
202 CALL zgesvx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 1,
203 $ x, 2, rcond, r1, r2, w, rw, info )
204 CALL chkxer(
'ZGESVX', infot, nout, lerr, ok )
206 CALL zgesvx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 2,
207 $ x, 1, rcond, r1, r2, w, rw, info )
208 CALL chkxer(
'ZGESVX', infot, nout, lerr, ok )
216 CALL zgesvxx(
'/',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b,
217 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
218 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
220 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
222 CALL zgesvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, c, b,
223 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
224 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
226 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
228 CALL zgesvxx(
'N',
'N', -1, 0, a, 1, af, 1, ip, eq, r, c, b,
229 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
230 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
232 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
234 CALL zgesvxx(
'N',
'N', 0, -1, a, 1, af, 1, ip, eq, r, c, b,
235 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
236 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
238 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
240 CALL zgesvxx(
'N',
'N', 2, 1, a, 1, af, 2, ip, eq, r, c, b,
241 $ 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
242 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
244 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
246 CALL zgesvxx(
'N',
'N', 2, 1, a, 2, af, 1, ip, eq, r, c, b,
247 $ 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
248 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
250 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
253 CALL zgesvxx(
'F',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b,
254 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
255 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
257 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
260 CALL zgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b,
261 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
262 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
264 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
267 CALL zgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b,
268 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
269 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
271 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
273 CALL zgesvxx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b,
274 $ 1, x, 2, rcond, rpvgrw, berr, n_err_bnds,
275 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
277 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
279 CALL zgesvxx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b,
280 $ 2, x, 1, rcond, rpvgrw, berr, n_err_bnds,
281 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
283 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
285 ELSE IF(
lsamen( 2, c2,
'GB' ) )
THEN
291 CALL zgbsv( -1, 0, 0, 0, a, 1, ip, b, 1, info )
292 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
294 CALL zgbsv( 1, -1, 0, 0, a, 1, ip, b, 1, info )
295 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
297 CALL zgbsv( 1, 0, -1, 0, a, 1, ip, b, 1, info )
298 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
300 CALL zgbsv( 0, 0, 0, -1, a, 1, ip, b, 1, info )
301 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
303 CALL zgbsv( 1, 1, 1, 0, a, 3, ip, b, 1, info )
304 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
306 CALL zgbsv( 2, 0, 0, 0, a, 1, ip, b, 1, info )
307 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
313 CALL zgbsvx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
314 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
315 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
317 CALL zgbsvx(
'N',
'/', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
318 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
319 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
321 CALL zgbsvx(
'N',
'N', -1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
322 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
323 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
325 CALL zgbsvx(
'N',
'N', 1, -1, 0, 0, a, 1, af, 1, ip, eq, r, c,
326 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
327 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
329 CALL zgbsvx(
'N',
'N', 1, 0, -1, 0, a, 1, af, 1, ip, eq, r, c,
330 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
331 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
333 CALL zgbsvx(
'N',
'N', 0, 0, 0, -1, a, 1, af, 1, ip, eq, r, c,
334 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
335 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
337 CALL zgbsvx(
'N',
'N', 1, 1, 1, 0, a, 2, af, 4, ip, eq, r, c,
338 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
339 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
341 CALL zgbsvx(
'N',
'N', 1, 1, 1, 0, a, 3, af, 3, ip, eq, r, c,
342 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
343 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
346 CALL zgbsvx(
'F',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
347 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
348 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
351 CALL zgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
352 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
353 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
356 CALL zgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
357 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
358 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
360 CALL zgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
361 $ b, 1, x, 2, rcond, r1, r2, w, rw, info )
362 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
364 CALL zgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
365 $ b, 2, x, 1, rcond, r1, r2, w, rw, info )
366 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
374 CALL zgbsvxx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
375 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
376 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
378 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
380 CALL zgbsvxx(
'N',
'/', 0, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
381 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
382 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
384 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
386 CALL zgbsvxx(
'N',
'N', -1, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
387 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
388 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
390 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
392 CALL zgbsvxx(
'N',
'N', 2, -1, 1, 0, a, 1, af, 1, ip, eq,
393 $ r, c, b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
394 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
396 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
398 CALL zgbsvxx(
'N',
'N', 2, 1, -1, 0, a, 1, af, 1, ip, eq,
399 $ r, c, b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
400 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
402 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
404 CALL zgbsvxx(
'N',
'N', 0, 1, 1, -1, a, 1, af, 1, ip, eq, r, c,
405 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
406 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
408 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
410 CALL zgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 2, af, 2, ip, eq, r, c,
411 $ b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
412 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
414 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
416 CALL zgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 3, ip, eq, r, c,
417 $ b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
418 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
420 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
423 CALL zgbsvxx(
'F',
'N', 0, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
424 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
425 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
427 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
430 CALL zgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
431 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
432 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
434 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
437 CALL zgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
438 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
439 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
441 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
443 CALL zgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
444 $ b, 1, x, 2, rcond, rpvgrw, berr, n_err_bnds,
445 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
447 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
449 CALL zgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
450 $ b, 2, x, 1, rcond, rpvgrw, berr, n_err_bnds,
451 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
453 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
455 ELSE IF(
lsamen( 2, c2,
'GT' ) )
THEN
461 CALL zgtsv( -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
463 CALL chkxer(
'ZGTSV ', infot, nout, lerr, ok )
465 CALL zgtsv( 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
467 CALL chkxer(
'ZGTSV ', infot, nout, lerr, ok )
469 CALL zgtsv( 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1, info )
470 CALL chkxer(
'ZGTSV ', infot, nout, lerr, ok )
476 CALL zgtsvx(
'/',
'N', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
477 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
478 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
479 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
481 CALL zgtsvx(
'N',
'/', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
482 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
483 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
484 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
486 CALL zgtsvx(
'N',
'N', -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
487 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
488 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
489 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
491 CALL zgtsvx(
'N',
'N', 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
492 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
493 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
494 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
496 CALL zgtsvx(
'N',
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
497 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
498 $ ip, b, 1, x, 2, rcond, r1, r2, w, rw, info )
499 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
501 CALL zgtsvx(
'N',
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
502 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
503 $ ip, b, 2, x, 1, rcond, r1, r2, w, rw, info )
504 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
506 ELSE IF(
lsamen( 2, c2,
'HR' ) )
THEN
510 srnamt =
'ZHESV_ROOK'
512 CALL zhesv_rook(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
513 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
515 CALL zhesv_rook(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
516 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
518 CALL zhesv_rook(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
519 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
521 CALL zhesv_rook(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
522 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
524 ELSE IF(
lsamen( 2, c2,
'PO' ) )
THEN
530 CALL zposv(
'/', 0, 0, a, 1, b, 1, info )
531 CALL chkxer(
'ZPOSV ', infot, nout, lerr, ok )
533 CALL zposv(
'U', -1, 0, a, 1, b, 1, info )
534 CALL chkxer(
'ZPOSV ', infot, nout, lerr, ok )
536 CALL zposv(
'U', 0, -1, a, 1, b, 1, info )
537 CALL chkxer(
'ZPOSV ', infot, nout, lerr, ok )
539 CALL zposv(
'U', 2, 0, a, 1, b, 2, info )
540 CALL chkxer(
'ZPOSV ', infot, nout, lerr, ok )
542 CALL zposv(
'U', 2, 0, a, 2, b, 1, info )
543 CALL chkxer(
'ZPOSV ', infot, nout, lerr, ok )
549 CALL zposvx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
550 $ rcond, r1, r2, w, rw, info )
551 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
553 CALL zposvx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
554 $ rcond, r1, r2, w, rw, info )
555 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
557 CALL zposvx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
558 $ rcond, r1, r2, w, rw, info )
559 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
561 CALL zposvx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
562 $ rcond, r1, r2, w, rw, info )
563 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
565 CALL zposvx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
566 $ rcond, r1, r2, w, rw, info )
567 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
569 CALL zposvx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
570 $ rcond, r1, r2, w, rw, info )
571 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
574 CALL zposvx(
'F',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
575 $ rcond, r1, r2, w, rw, info )
576 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
579 CALL zposvx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
580 $ rcond, r1, r2, w, rw, info )
581 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
583 CALL zposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 1, x, 2,
584 $ rcond, r1, r2, w, rw, info )
585 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
587 CALL zposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
588 $ rcond, r1, r2, w, rw, info )
589 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
597 CALL zposvxx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
598 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
599 $ err_bnds_c, nparams, params, w, rw, info )
600 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
602 CALL zposvxx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
603 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
604 $ err_bnds_c, nparams, params, w, rw, info )
605 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
607 CALL zposvxx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
608 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
609 $ err_bnds_c, nparams, params, w, rw, info )
610 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
612 CALL zposvxx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
613 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
614 $ err_bnds_c, nparams, params, w, rw, info )
615 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
617 CALL zposvxx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
618 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
619 $ err_bnds_c, nparams, params, w, rw, info )
620 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
622 CALL zposvxx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
623 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
624 $ err_bnds_c, nparams, params, w, rw, info )
625 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
628 CALL zposvxx(
'F',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
629 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
630 $ err_bnds_c, nparams, params, w, rw, info )
631 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
634 CALL zposvxx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
635 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
636 $ err_bnds_c, nparams, params, w, rw, info )
637 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
639 CALL zposvxx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 1, x, 2,
640 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
641 $ err_bnds_c, nparams, params, w, rw, info )
642 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
644 CALL zposvxx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
645 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
646 $ err_bnds_c, nparams, params, w, rw, info )
647 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
649 ELSE IF(
lsamen( 2, c2,
'PP' ) )
THEN
655 CALL zppsv(
'/', 0, 0, a, b, 1, info )
656 CALL chkxer(
'ZPPSV ', infot, nout, lerr, ok )
658 CALL zppsv(
'U', -1, 0, a, b, 1, info )
659 CALL chkxer(
'ZPPSV ', infot, nout, lerr, ok )
661 CALL zppsv(
'U', 0, -1, a, b, 1, info )
662 CALL chkxer(
'ZPPSV ', infot, nout, lerr, ok )
664 CALL zppsv(
'U', 2, 0, a, b, 1, info )
665 CALL chkxer(
'ZPPSV ', infot, nout, lerr, ok )
671 CALL zppsvx(
'/',
'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
672 $ r1, r2, w, rw, info )
673 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
675 CALL zppsvx(
'N',
'/', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
676 $ r1, r2, w, rw, info )
677 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
679 CALL zppsvx(
'N',
'U', -1, 0, a, af, eq, c, b, 1, x, 1, rcond,
680 $ r1, r2, w, rw, info )
681 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
683 CALL zppsvx(
'N',
'U', 0, -1, a, af, eq, c, b, 1, x, 1, rcond,
684 $ r1, r2, w, rw, info )
685 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
688 CALL zppsvx(
'F',
'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
689 $ r1, r2, w, rw, info )
690 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
693 CALL zppsvx(
'F',
'U', 1, 0, a, af, eq, c, b, 1, x, 1, rcond,
694 $ r1, r2, w, rw, info )
695 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
697 CALL zppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 1, x, 2, rcond,
698 $ r1, r2, w, rw, info )
699 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
701 CALL zppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 2, x, 1, rcond,
702 $ r1, r2, w, rw, info )
703 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
705 ELSE IF(
lsamen( 2, c2,
'PB' ) )
THEN
711 CALL zpbsv(
'/', 0, 0, 0, a, 1, b, 1, info )
712 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
714 CALL zpbsv(
'U', -1, 0, 0, a, 1, b, 1, info )
715 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
717 CALL zpbsv(
'U', 1, -1, 0, a, 1, b, 1, info )
718 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
720 CALL zpbsv(
'U', 0, 0, -1, a, 1, b, 1, info )
721 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
723 CALL zpbsv(
'U', 1, 1, 0, a, 1, b, 2, info )
724 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
726 CALL zpbsv(
'U', 2, 0, 0, a, 1, b, 1, info )
727 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
733 CALL zpbsvx(
'/',
'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
734 $ rcond, r1, r2, w, rw, info )
735 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
737 CALL zpbsvx(
'N',
'/', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
738 $ rcond, r1, r2, w, rw, info )
739 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
741 CALL zpbsvx(
'N',
'U', -1, 0, 0, a, 1, af, 1, eq, c, b, 1, x,
742 $ 1, rcond, r1, r2, w, rw, info )
743 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
745 CALL zpbsvx(
'N',
'U', 1, -1, 0, a, 1, af, 1, eq, c, b, 1, x,
746 $ 1, rcond, r1, r2, w, rw, info )
747 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
749 CALL zpbsvx(
'N',
'U', 0, 0, -1, a, 1, af, 1, eq, c, b, 1, x,
750 $ 1, rcond, r1, r2, w, rw, info )
751 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
753 CALL zpbsvx(
'N',
'U', 1, 1, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
754 $ rcond, r1, r2, w, rw, info )
755 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
757 CALL zpbsvx(
'N',
'U', 1, 1, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
758 $ rcond, r1, r2, w, rw, info )
759 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
762 CALL zpbsvx(
'F',
'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
763 $ rcond, r1, r2, w, rw, info )
764 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
767 CALL zpbsvx(
'F',
'U', 1, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
768 $ rcond, r1, r2, w, rw, info )
769 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
771 CALL zpbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 2,
772 $ rcond, r1, r2, w, rw, info )
773 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
775 CALL zpbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 2, x, 1,
776 $ rcond, r1, r2, w, rw, info )
777 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
779 ELSE IF(
lsamen( 2, c2,
'PT' ) )
THEN
785 CALL zptsv( -1, 0, r, a( 1, 1 ), b, 1, info )
786 CALL chkxer(
'ZPTSV ', infot, nout, lerr, ok )
788 CALL zptsv( 0, -1, r, a( 1, 1 ), b, 1, info )
789 CALL chkxer(
'ZPTSV ', infot, nout, lerr, ok )
791 CALL zptsv( 2, 0, r, a( 1, 1 ), b, 1, info )
792 CALL chkxer(
'ZPTSV ', infot, nout, lerr, ok )
798 CALL zptsvx(
'/', 0, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
799 $ 1, rcond, r1, r2, w, rw, info )
800 CALL chkxer(
'ZPTSVX', infot, nout, lerr, ok )
802 CALL zptsvx(
'N', -1, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
803 $ 1, rcond, r1, r2, w, rw, info )
804 CALL chkxer(
'ZPTSVX', infot, nout, lerr, ok )
806 CALL zptsvx(
'N', 0, -1, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
807 $ 1, rcond, r1, r2, w, rw, info )
808 CALL chkxer(
'ZPTSVX', infot, nout, lerr, ok )
810 CALL zptsvx(
'N', 2, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
811 $ 2, rcond, r1, r2, w, rw, info )
812 CALL chkxer(
'ZPTSVX', infot, nout, lerr, ok )
814 CALL zptsvx(
'N', 2, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 2, x,
815 $ 1, rcond, r1, r2, w, rw, info )
816 CALL chkxer(
'ZPTSVX', infot, nout, lerr, ok )
818 ELSE IF(
lsamen( 2, c2,
'HE' ) )
THEN
824 CALL zhesv(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
825 CALL chkxer(
'ZHESV ', infot, nout, lerr, ok )
827 CALL zhesv(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
828 CALL chkxer(
'ZHESV ', infot, nout, lerr, ok )
830 CALL zhesv(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
831 CALL chkxer(
'ZHESV ', infot, nout, lerr, ok )
833 CALL zhesv(
'U', 2, 0, a, 1, ip, b, 2, w, 1, info )
834 CALL chkxer(
'ZHESV ', infot, nout, lerr, ok )
836 CALL zhesv(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
837 CALL chkxer(
'ZHESV ', infot, nout, lerr, ok )
843 CALL zhesvx(
'/',
'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
844 $ rcond, r1, r2, w, 1, rw, info )
845 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
847 CALL zhesvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
848 $ rcond, r1, r2, w, 1, rw, info )
849 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
851 CALL zhesvx(
'N',
'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
852 $ rcond, r1, r2, w, 1, rw, info )
853 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
855 CALL zhesvx(
'N',
'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
856 $ rcond, r1, r2, w, 1, rw, info )
857 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
859 CALL zhesvx(
'N',
'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
860 $ rcond, r1, r2, w, 4, rw, info )
861 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
863 CALL zhesvx(
'N',
'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
864 $ rcond, r1, r2, w, 4, rw, info )
865 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
867 CALL zhesvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
868 $ rcond, r1, r2, w, 4, rw, info )
869 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
871 CALL zhesvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
872 $ rcond, r1, r2, w, 4, rw, info )
873 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
875 CALL zhesvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
876 $ rcond, r1, r2, w, 3, rw, info )
877 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
885 CALL zhesvxx(
'/',
'U', 0, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
886 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
887 $ err_bnds_c, nparams, params, w, rw, info )
888 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
890 CALL zhesvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
891 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
892 $ err_bnds_c, nparams, params, w, rw, info )
893 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
895 CALL zhesvxx(
'N',
'U', -1, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
896 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
897 $ err_bnds_c, nparams, params, w, rw, info )
898 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
900 CALL zhesvxx(
'N',
'U', 0, -1, a, 1, af, 1, ip, eq, c, b, 1, x,
901 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
902 $ err_bnds_c, nparams, params, w, rw, info )
903 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
905 CALL zhesvxx(
'N',
'U', 2, 0, a, 1, af, 2, ip, eq, c, b, 2, x,
906 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
907 $ err_bnds_c, nparams, params, w, rw, info )
908 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
910 CALL zhesvxx(
'N',
'U', 2, 0, a, 2, af, 1, ip, eq, c, b, 2, x,
911 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
912 $ err_bnds_c, nparams, params, w, rw, info )
913 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
916 CALL zhesvxx(
'F',
'U', 0, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
917 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
918 $ err_bnds_c, nparams, params, w, rw, info )
919 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
922 CALL zhesvxx(
'F',
'U', 1, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
923 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
924 $ err_bnds_c, nparams, params, w, rw, info )
925 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
927 CALL zhesvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, c, b, 1, x,
928 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
929 $ err_bnds_c, nparams, params, w, rw, info )
930 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
932 CALL zhesvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, c, b, 2, x,
933 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
934 $ err_bnds_c, nparams, params, w, rw, info )
935 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
937 ELSE IF(
lsamen( 2, c2,
'HR' ) )
THEN
941 srnamt =
'ZHESV_ROOK'
943 CALL zhesv_rook(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
944 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
946 CALL zhesv_rook(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
947 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
949 CALL zhesv_rook(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
950 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
952 CALL zhesv_rook(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
953 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
955 ELSE IF(
lsamen( 2, c2,
'HP' ) )
THEN
961 CALL zhpsv(
'/', 0, 0, a, ip, b, 1, info )
962 CALL chkxer(
'ZHPSV ', infot, nout, lerr, ok )
964 CALL zhpsv(
'U', -1, 0, a, ip, b, 1, info )
965 CALL chkxer(
'ZHPSV ', infot, nout, lerr, ok )
967 CALL zhpsv(
'U', 0, -1, a, ip, b, 1, info )
968 CALL chkxer(
'ZHPSV ', infot, nout, lerr, ok )
970 CALL zhpsv(
'U', 2, 0, a, ip, b, 1, info )
971 CALL chkxer(
'ZHPSV ', infot, nout, lerr, ok )
977 CALL zhpsvx(
'/',
'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
979 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
981 CALL zhpsvx(
'N',
'/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
983 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
985 CALL zhpsvx(
'N',
'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
987 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
989 CALL zhpsvx(
'N',
'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
991 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
993 CALL zhpsvx(
'N',
'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
995 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
997 CALL zhpsvx(
'N',
'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
999 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
1001 ELSE IF(
lsamen( 2, c2,
'SY' ) )
THEN
1007 CALL zsysv(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
1008 CALL chkxer(
'ZSYSV ', infot, nout, lerr, ok )
1010 CALL zsysv(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
1011 CALL chkxer(
'ZSYSV ', infot, nout, lerr, ok )
1013 CALL zsysv(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
1014 CALL chkxer(
'ZSYSV ', infot, nout, lerr, ok )
1016 CALL zsysv(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
1017 CALL chkxer(
'ZSYSV ', infot, nout, lerr, ok )
1023 CALL zsysvx(
'/',
'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
1024 $ rcond, r1, r2, w, 1, rw, info )
1025 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1027 CALL zsysvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
1028 $ rcond, r1, r2, w, 1, rw, info )
1029 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1031 CALL zsysvx(
'N',
'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
1032 $ rcond, r1, r2, w, 1, rw, info )
1033 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1035 CALL zsysvx(
'N',
'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
1036 $ rcond, r1, r2, w, 1, rw, info )
1037 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1039 CALL zsysvx(
'N',
'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
1040 $ rcond, r1, r2, w, 4, rw, info )
1041 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1043 CALL zsysvx(
'N',
'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
1044 $ rcond, r1, r2, w, 4, rw, info )
1045 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1047 CALL zsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
1048 $ rcond, r1, r2, w, 4, rw, info )
1049 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1051 CALL zsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
1052 $ rcond, r1, r2, w, 4, rw, info )
1053 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1055 CALL zsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
1056 $ rcond, r1, r2, w, 3, rw, info )
1057 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1066 CALL zsysvxx(
'/',
'U', 0, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
1067 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1068 $ err_bnds_c, nparams, params, w, rw, info )
1069 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1071 CALL zsysvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
1072 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1073 $ err_bnds_c, nparams, params, w, rw, info )
1074 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1076 CALL zsysvxx(
'N',
'U', -1, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
1077 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1078 $ err_bnds_c, nparams, params, w, rw, info )
1079 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1082 CALL zsysvxx(
'N',
'U', 0, -1, a, 1, af, 1, ip, eq, r, b, 1, x,
1083 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1084 $ err_bnds_c, nparams, params, w, rw, info )
1085 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1088 CALL zsysvxx(
'N',
'U', 2, 0, a, 1, af, 2, ip, eq, r, b, 2, x,
1089 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1090 $ err_bnds_c, nparams, params, w, rw, info )
1091 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1093 CALL zsysvxx(
'N',
'U', 2, 0, a, 2, af, 1, ip, eq, r, b, 2, x,
1094 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1095 $ err_bnds_c, nparams, params, w, rw, info )
1096 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1098 CALL zsysvxx(
'F',
'U', 2, 0, a, 2, af, 2, ip,
'A', r, b, 2, x,
1099 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1100 $ err_bnds_c, nparams, params, w, rw, info )
1101 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1104 CALL zsysvxx(
'F',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 2, x,
1105 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1106 $ err_bnds_c, nparams, params, w, rw, info )
1107 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1111 CALL zsysvxx(
'F',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 2, x,
1112 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1113 $ err_bnds_c, nparams, params, w, rw, info )
1114 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1117 CALL zsysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 1, x,
1118 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1119 $ err_bnds_c, nparams, params, w, rw, info )
1120 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1122 CALL zsysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 2, x,
1123 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1124 $ err_bnds_c, nparams, params, w, rw, info )
1125 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1127 ELSE IF(
lsamen( 2, c2,
'SR' ) )
THEN
1131 srnamt =
'ZSYSV_ROOK'
1133 CALL zsysv_rook(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
1134 CALL chkxer(
'ZSYSV_ROOK', infot, nout, lerr, ok )
1136 CALL zsysv_rook(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
1137 CALL chkxer(
'ZSYSV_ROOK', infot, nout, lerr, ok )
1139 CALL zsysv_rook(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
1140 CALL chkxer(
'ZSYSV_ROOK', infot, nout, lerr, ok )
1142 CALL zsysv_rook(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
1143 CALL chkxer(
'ZSYSV_ROOK', infot, nout, lerr, ok )
1145 ELSE IF(
lsamen( 2, c2,
'SP' ) )
THEN
1151 CALL zspsv(
'/', 0, 0, a, ip, b, 1, info )
1152 CALL chkxer(
'ZSPSV ', infot, nout, lerr, ok )
1154 CALL zspsv(
'U', -1, 0, a, ip, b, 1, info )
1155 CALL chkxer(
'ZSPSV ', infot, nout, lerr, ok )
1157 CALL zspsv(
'U', 0, -1, a, ip, b, 1, info )
1158 CALL chkxer(
'ZSPSV ', infot, nout, lerr, ok )
1160 CALL zspsv(
'U', 2, 0, a, ip, b, 1, info )
1161 CALL chkxer(
'ZSPSV ', infot, nout, lerr, ok )
1167 CALL zspsvx(
'/',
'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1169 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1171 CALL zspsvx(
'N',
'/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1173 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1175 CALL zspsvx(
'N',
'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1177 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1179 CALL zspsvx(
'N',
'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
1181 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1183 CALL zspsvx(
'N',
'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
1185 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1187 CALL zspsvx(
'N',
'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
1189 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1195 WRITE( nout, fmt = 9999 )path
1197 WRITE( nout, fmt = 9998 )path
1200 9999
FORMAT( 1x, a3,
' drivers passed the tests of the error exits' )
1201 9998
FORMAT(
' *** ', a3,
' drivers failed the tests of the error ',
subroutine zgesvxx(FACT, TRANS, N, NRHS, A, LDA, AF, LDAF, IPIV, EQUED, R, C, B, LDB, X, LDX, RCOND, RPVGRW, BERR, N_ERR_BNDS, ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, RWORK, INFO)
ZGESVXX computes the solution to system of linear equations A * X = B for GE matrices ...
subroutine zhesvxx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, IPIV, EQUED, S, B, LDB, X, LDX, RCOND, RPVGRW, BERR, N_ERR_BNDS, ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, RWORK, INFO)
ZHESVXX computes the solution to system of linear equations A * X = B for HE matrices ...
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) ...
logical function lsamen(N, CA, CB)
LSAMEN
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 ...
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 ...
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 ...
subroutine zgbsvxx(FACT, TRANS, N, KL, KU, NRHS, AB, LDAB, AFB, LDAFB, IPIV, EQUED, R, C, B, LDB, X, LDX, RCOND, RPVGRW, BERR, N_ERR_BNDS, ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, RWORK, INFO)
ZGBSVXX computes the solution to system of linear equations A * X = B for GB matrices ...
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 ...
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) ...
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 ...
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 ...
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 ...
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine zsysvxx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, IPIV, EQUED, S, B, LDB, X, LDX, RCOND, RPVGRW, BERR, N_ERR_BNDS, ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, RWORK, INFO)
ZSYSVXX computes the solution to system of linear equations A * X = B for SY matrices ...
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 ...
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 ...
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 ...
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 ...
subroutine zppsv(UPLO, N, NRHS, AP, B, LDB, INFO)
ZPPSV computes the solution to system of linear equations A * X = B for OTHER matrices ...
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 ...
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 ...
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 ...
subroutine zerrvx(PATH, NUNIT)
ZERRVX
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 ...
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 ...
subroutine zposvxx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, EQUED, S, B, LDB, X, LDX, RCOND, RPVGRW, BERR, N_ERR_BNDS, ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, RWORK, INFO)
ZPOSVXX computes the solution to system of linear equations A * X = B for PO matrices ...
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 ...
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 ...
subroutine zptsv(N, NRHS, D, E, B, LDB, INFO)
ZPTSV computes the solution to system of linear equations A * X = B for PT matrices ...
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 ...
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 ...