56 SUBROUTINE zerrvx( PATH, NUNIT )
77 INTEGER I, INFO, J, N_ERR_BNDS, NPARAMS
78 DOUBLE PRECISION RCOND, RPVGRW, BERR
82 DOUBLE PRECISION C( nmax ), R( nmax ), R1( nmax ), R2( nmax ),
83 $ rf( nmax ), rw( nmax ), err_bnds_n( nmax, 3 ),
84 $ err_bnds_c( nmax, 3 ), params( 1 )
85 COMPLEX*16 A( nmax, nmax ), AF( nmax, nmax ), B( nmax ),
86 $ w( 2*nmax ), x( nmax )
106 COMMON / infoc / infot, nout, ok, lerr
107 COMMON / srnamc / srnamt
110 INTRINSIC dble, dcmplx
115 WRITE( nout, fmt = * )
122 a( i, j ) = dcmplx( 1.d0 / dble( i+j ),
123 $ -1.d0 / dble( i+j ) )
124 af( i, j ) = dcmplx( 1.d0 / dble( i+j ),
125 $ -1.d0 / dble( i+j ) )
139 IF( lsamen( 2, c2,
'GE' ) )
THEN
145 CALL zgesv( -1, 0, a, 1, ip, b, 1, info )
146 CALL chkxer(
'ZGESV ', infot, nout, lerr, ok )
148 CALL zgesv( 0, -1, a, 1, ip, b, 1, info )
149 CALL chkxer(
'ZGESV ', infot, nout, lerr, ok )
151 CALL zgesv( 2, 1, a, 1, ip, b, 2, info )
152 CALL chkxer(
'ZGESV ', infot, nout, lerr, ok )
154 CALL zgesv( 2, 1, a, 2, ip, b, 1, info )
155 CALL chkxer(
'ZGESV ', infot, nout, lerr, ok )
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 )
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 )
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 )
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 )
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 )
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 )
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 )
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 )
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 )
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 )
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 )
214 CALL zgesvxx(
'/',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b,
215 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
216 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
218 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
220 CALL zgesvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, c, b,
221 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
222 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
224 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
226 CALL zgesvxx(
'N',
'N', -1, 0, a, 1, af, 1, ip, eq, r, c, b,
227 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
228 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
230 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
232 CALL zgesvxx(
'N',
'N', 0, -1, a, 1, af, 1, ip, eq, r, c, b,
233 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
234 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
236 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
238 CALL zgesvxx(
'N',
'N', 2, 1, a, 1, af, 2, ip, eq, r, c, b,
239 $ 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
240 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
242 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
244 CALL zgesvxx(
'N',
'N', 2, 1, a, 2, af, 1, ip, eq, r, c, b,
245 $ 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
246 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
248 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
251 CALL zgesvxx(
'F',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b,
252 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
253 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
255 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
258 CALL zgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b,
259 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
260 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
262 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
265 CALL zgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b,
266 $ 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
267 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
269 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
271 CALL zgesvxx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b,
272 $ 1, x, 2, rcond, rpvgrw, berr, n_err_bnds,
273 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
275 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
277 CALL zgesvxx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b,
278 $ 2, x, 1, rcond, rpvgrw, berr, n_err_bnds,
279 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
281 CALL chkxer(
'ZGESVXX', infot, nout, lerr, ok )
283 ELSE IF( lsamen( 2, c2,
'GB' ) )
THEN
289 CALL zgbsv( -1, 0, 0, 0, a, 1, ip, b, 1, info )
290 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
292 CALL zgbsv( 1, -1, 0, 0, a, 1, ip, b, 1, info )
293 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
295 CALL zgbsv( 1, 0, -1, 0, a, 1, ip, b, 1, info )
296 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
298 CALL zgbsv( 0, 0, 0, -1, a, 1, ip, b, 1, info )
299 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
301 CALL zgbsv( 1, 1, 1, 0, a, 3, ip, b, 1, info )
302 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
304 CALL zgbsv( 2, 0, 0, 0, a, 1, ip, b, 1, info )
305 CALL chkxer(
'ZGBSV ', infot, nout, lerr, ok )
311 CALL zgbsvx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
312 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
313 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
315 CALL zgbsvx(
'N',
'/', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
316 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
317 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
319 CALL zgbsvx(
'N',
'N', -1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
320 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
321 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
323 CALL zgbsvx(
'N',
'N', 1, -1, 0, 0, a, 1, af, 1, ip, eq, r, c,
324 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
325 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
327 CALL zgbsvx(
'N',
'N', 1, 0, -1, 0, a, 1, af, 1, ip, eq, r, c,
328 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
329 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
331 CALL zgbsvx(
'N',
'N', 0, 0, 0, -1, a, 1, af, 1, ip, eq, r, c,
332 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
333 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
335 CALL zgbsvx(
'N',
'N', 1, 1, 1, 0, a, 2, af, 4, ip, eq, r, c,
336 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
337 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
339 CALL zgbsvx(
'N',
'N', 1, 1, 1, 0, a, 3, af, 3, ip, eq, r, c,
340 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
341 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
344 CALL zgbsvx(
'F',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
345 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
346 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
349 CALL zgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
350 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
351 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
354 CALL zgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
355 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
356 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
358 CALL zgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
359 $ b, 1, x, 2, rcond, r1, r2, w, rw, info )
360 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
362 CALL zgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
363 $ b, 2, x, 1, rcond, r1, r2, w, rw, info )
364 CALL chkxer(
'ZGBSVX', infot, nout, lerr, ok )
372 CALL zgbsvxx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
373 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
374 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
376 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
378 CALL zgbsvxx(
'N',
'/', 0, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
379 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
380 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
382 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
384 CALL zgbsvxx(
'N',
'N', -1, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
385 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
386 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
388 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
390 CALL zgbsvxx(
'N',
'N', 2, -1, 1, 0, a, 1, af, 1, ip, eq,
391 $ r, c, b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
392 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
394 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
396 CALL zgbsvxx(
'N',
'N', 2, 1, -1, 0, a, 1, af, 1, ip, eq,
397 $ r, c, b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
398 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
400 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
402 CALL zgbsvxx(
'N',
'N', 0, 1, 1, -1, a, 1, af, 1, ip, eq, r, c,
403 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
404 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
406 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
408 CALL zgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 2, af, 2, ip, eq, r, c,
409 $ b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
410 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
412 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
414 CALL zgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 3, ip, eq, r, c,
415 $ b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
416 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
418 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
421 CALL zgbsvxx(
'F',
'N', 0, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
422 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
423 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
425 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
428 CALL zgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
429 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
430 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
432 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
435 CALL zgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
436 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
437 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
439 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
441 CALL zgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
442 $ b, 1, x, 2, rcond, rpvgrw, berr, n_err_bnds,
443 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
445 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
447 CALL zgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
448 $ b, 2, x, 1, rcond, rpvgrw, berr, n_err_bnds,
449 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
451 CALL chkxer(
'ZGBSVXX', infot, nout, lerr, ok )
453 ELSE IF( lsamen( 2, c2,
'GT' ) )
THEN
459 CALL zgtsv( -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
461 CALL chkxer(
'ZGTSV ', infot, nout, lerr, ok )
463 CALL zgtsv( 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
465 CALL chkxer(
'ZGTSV ', infot, nout, lerr, ok )
467 CALL zgtsv( 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1, info )
468 CALL chkxer(
'ZGTSV ', infot, nout, lerr, ok )
474 CALL zgtsvx(
'/',
'N', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
475 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
476 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
477 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
479 CALL zgtsvx(
'N',
'/', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
480 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
481 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
482 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
484 CALL zgtsvx(
'N',
'N', -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
485 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
486 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
487 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
489 CALL zgtsvx(
'N',
'N', 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
490 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
491 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
492 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
494 CALL zgtsvx(
'N',
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
495 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
496 $ ip, b, 1, x, 2, rcond, r1, r2, w, rw, info )
497 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
499 CALL zgtsvx(
'N',
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
500 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
501 $ ip, b, 2, x, 1, rcond, r1, r2, w, rw, info )
502 CALL chkxer(
'ZGTSVX', infot, nout, lerr, ok )
504 ELSE IF( lsamen( 2, c2,
'HR' ) )
THEN
508 srnamt =
'ZHESV_ROOK'
510 CALL zhesv_rook(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
511 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
513 CALL zhesv_rook(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
514 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
516 CALL zhesv_rook(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
517 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
519 CALL zhesv_rook(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
520 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
522 ELSE IF( lsamen( 2, c2,
'PO' ) )
THEN
528 CALL zposv(
'/', 0, 0, a, 1, b, 1, info )
529 CALL chkxer(
'ZPOSV ', infot, nout, lerr, ok )
531 CALL zposv(
'U', -1, 0, a, 1, b, 1, info )
532 CALL chkxer(
'ZPOSV ', infot, nout, lerr, ok )
534 CALL zposv(
'U', 0, -1, a, 1, b, 1, info )
535 CALL chkxer(
'ZPOSV ', infot, nout, lerr, ok )
537 CALL zposv(
'U', 2, 0, a, 1, b, 2, info )
538 CALL chkxer(
'ZPOSV ', infot, nout, lerr, ok )
540 CALL zposv(
'U', 2, 0, a, 2, b, 1, info )
541 CALL chkxer(
'ZPOSV ', infot, nout, lerr, ok )
547 CALL zposvx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
548 $ rcond, r1, r2, w, rw, info )
549 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
551 CALL zposvx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
552 $ rcond, r1, r2, w, rw, info )
553 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
555 CALL zposvx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
556 $ rcond, r1, r2, w, rw, info )
557 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
559 CALL zposvx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
560 $ rcond, r1, r2, w, rw, info )
561 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
563 CALL zposvx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
564 $ rcond, r1, r2, w, rw, info )
565 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
567 CALL zposvx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
568 $ rcond, r1, r2, w, rw, info )
569 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
572 CALL zposvx(
'F',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
573 $ rcond, r1, r2, w, rw, info )
574 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
577 CALL zposvx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
578 $ rcond, r1, r2, w, rw, info )
579 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
581 CALL zposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 1, x, 2,
582 $ rcond, r1, r2, w, rw, info )
583 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
585 CALL zposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
586 $ rcond, r1, r2, w, rw, info )
587 CALL chkxer(
'ZPOSVX', infot, nout, lerr, ok )
595 CALL zposvxx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
596 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
597 $ err_bnds_c, nparams, params, w, rw, info )
598 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
600 CALL zposvxx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
601 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
602 $ err_bnds_c, nparams, params, w, rw, info )
603 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
605 CALL zposvxx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
606 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
607 $ err_bnds_c, nparams, params, w, rw, info )
608 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
610 CALL zposvxx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
611 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
612 $ err_bnds_c, nparams, params, w, rw, info )
613 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
615 CALL zposvxx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
616 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
617 $ err_bnds_c, nparams, params, w, rw, info )
618 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
620 CALL zposvxx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
621 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
622 $ err_bnds_c, nparams, params, w, rw, info )
623 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
626 CALL zposvxx(
'F',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
627 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
628 $ err_bnds_c, nparams, params, w, rw, info )
629 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
632 CALL zposvxx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
633 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
634 $ err_bnds_c, nparams, params, w, rw, info )
635 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
637 CALL zposvxx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 1, x, 2,
638 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
639 $ err_bnds_c, nparams, params, w, rw, info )
640 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
642 CALL zposvxx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
643 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
644 $ err_bnds_c, nparams, params, w, rw, info )
645 CALL chkxer(
'ZPOSVXX', infot, nout, lerr, ok )
647 ELSE IF( lsamen( 2, c2,
'PP' ) )
THEN
653 CALL zppsv(
'/', 0, 0, a, b, 1, info )
654 CALL chkxer(
'ZPPSV ', infot, nout, lerr, ok )
656 CALL zppsv(
'U', -1, 0, a, b, 1, info )
657 CALL chkxer(
'ZPPSV ', infot, nout, lerr, ok )
659 CALL zppsv(
'U', 0, -1, a, b, 1, info )
660 CALL chkxer(
'ZPPSV ', infot, nout, lerr, ok )
662 CALL zppsv(
'U', 2, 0, a, b, 1, info )
663 CALL chkxer(
'ZPPSV ', infot, nout, lerr, ok )
669 CALL zppsvx(
'/',
'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
670 $ r1, r2, w, rw, info )
671 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
673 CALL zppsvx(
'N',
'/', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
674 $ r1, r2, w, rw, info )
675 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
677 CALL zppsvx(
'N',
'U', -1, 0, a, af, eq, c, b, 1, x, 1, rcond,
678 $ r1, r2, w, rw, info )
679 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
681 CALL zppsvx(
'N',
'U', 0, -1, a, af, eq, c, b, 1, x, 1, rcond,
682 $ r1, r2, w, rw, info )
683 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
686 CALL zppsvx(
'F',
'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
687 $ r1, r2, w, rw, info )
688 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
691 CALL zppsvx(
'F',
'U', 1, 0, a, af, eq, c, b, 1, x, 1, rcond,
692 $ r1, r2, w, rw, info )
693 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
695 CALL zppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 1, x, 2, rcond,
696 $ r1, r2, w, rw, info )
697 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
699 CALL zppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 2, x, 1, rcond,
700 $ r1, r2, w, rw, info )
701 CALL chkxer(
'ZPPSVX', infot, nout, lerr, ok )
703 ELSE IF( lsamen( 2, c2,
'PB' ) )
THEN
709 CALL zpbsv(
'/', 0, 0, 0, a, 1, b, 1, info )
710 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
712 CALL zpbsv(
'U', -1, 0, 0, a, 1, b, 1, info )
713 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
715 CALL zpbsv(
'U', 1, -1, 0, a, 1, b, 1, info )
716 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
718 CALL zpbsv(
'U', 0, 0, -1, a, 1, b, 1, info )
719 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
721 CALL zpbsv(
'U', 1, 1, 0, a, 1, b, 2, info )
722 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
724 CALL zpbsv(
'U', 2, 0, 0, a, 1, b, 1, info )
725 CALL chkxer(
'ZPBSV ', infot, nout, lerr, ok )
731 CALL zpbsvx(
'/',
'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
732 $ rcond, r1, r2, w, rw, info )
733 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
735 CALL zpbsvx(
'N',
'/', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
736 $ rcond, r1, r2, w, rw, info )
737 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
739 CALL zpbsvx(
'N',
'U', -1, 0, 0, a, 1, af, 1, eq, c, b, 1, x,
740 $ 1, rcond, r1, r2, w, rw, info )
741 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
743 CALL zpbsvx(
'N',
'U', 1, -1, 0, a, 1, af, 1, eq, c, b, 1, x,
744 $ 1, rcond, r1, r2, w, rw, info )
745 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
747 CALL zpbsvx(
'N',
'U', 0, 0, -1, a, 1, af, 1, eq, c, b, 1, x,
748 $ 1, rcond, r1, r2, w, rw, info )
749 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
751 CALL zpbsvx(
'N',
'U', 1, 1, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
752 $ rcond, r1, r2, w, rw, info )
753 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
755 CALL zpbsvx(
'N',
'U', 1, 1, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
756 $ rcond, r1, r2, w, rw, info )
757 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
760 CALL zpbsvx(
'F',
'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
761 $ rcond, r1, r2, w, rw, info )
762 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
765 CALL zpbsvx(
'F',
'U', 1, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
766 $ rcond, r1, r2, w, rw, info )
767 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
769 CALL zpbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 2,
770 $ rcond, r1, r2, w, rw, info )
771 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
773 CALL zpbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 2, x, 1,
774 $ rcond, r1, r2, w, rw, info )
775 CALL chkxer(
'ZPBSVX', infot, nout, lerr, ok )
777 ELSE IF( lsamen( 2, c2,
'PT' ) )
THEN
783 CALL zptsv( -1, 0, r, a( 1, 1 ), b, 1, info )
784 CALL chkxer(
'ZPTSV ', infot, nout, lerr, ok )
786 CALL zptsv( 0, -1, r, a( 1, 1 ), b, 1, info )
787 CALL chkxer(
'ZPTSV ', infot, nout, lerr, ok )
789 CALL zptsv( 2, 0, r, a( 1, 1 ), b, 1, info )
790 CALL chkxer(
'ZPTSV ', infot, nout, lerr, ok )
796 CALL zptsvx(
'/', 0, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
797 $ 1, rcond, r1, r2, w, rw, info )
798 CALL chkxer(
'ZPTSVX', infot, nout, lerr, ok )
800 CALL zptsvx(
'N', -1, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
801 $ 1, rcond, r1, r2, w, rw, info )
802 CALL chkxer(
'ZPTSVX', infot, nout, lerr, ok )
804 CALL zptsvx(
'N', 0, -1, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
805 $ 1, rcond, r1, r2, w, rw, info )
806 CALL chkxer(
'ZPTSVX', infot, nout, lerr, ok )
808 CALL zptsvx(
'N', 2, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
809 $ 2, rcond, r1, r2, w, rw, info )
810 CALL chkxer(
'ZPTSVX', infot, nout, lerr, ok )
812 CALL zptsvx(
'N', 2, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 2, x,
813 $ 1, rcond, r1, r2, w, rw, info )
814 CALL chkxer(
'ZPTSVX', infot, nout, lerr, ok )
816 ELSE IF( lsamen( 2, c2,
'HE' ) )
THEN
822 CALL zhesv(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
823 CALL chkxer(
'ZHESV ', infot, nout, lerr, ok )
825 CALL zhesv(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
826 CALL chkxer(
'ZHESV ', infot, nout, lerr, ok )
828 CALL zhesv(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
829 CALL chkxer(
'ZHESV ', infot, nout, lerr, ok )
831 CALL zhesv(
'U', 2, 0, a, 1, ip, b, 2, w, 1, info )
832 CALL chkxer(
'ZHESV ', infot, nout, lerr, ok )
834 CALL zhesv(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
835 CALL chkxer(
'ZHESV ', infot, nout, lerr, ok )
841 CALL zhesvx(
'/',
'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
842 $ rcond, r1, r2, w, 1, rw, info )
843 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
845 CALL zhesvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
846 $ rcond, r1, r2, w, 1, rw, info )
847 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
849 CALL zhesvx(
'N',
'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
850 $ rcond, r1, r2, w, 1, rw, info )
851 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
853 CALL zhesvx(
'N',
'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
854 $ rcond, r1, r2, w, 1, rw, info )
855 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
857 CALL zhesvx(
'N',
'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
858 $ rcond, r1, r2, w, 4, rw, info )
859 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
861 CALL zhesvx(
'N',
'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
862 $ rcond, r1, r2, w, 4, rw, info )
863 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
865 CALL zhesvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
866 $ rcond, r1, r2, w, 4, rw, info )
867 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
869 CALL zhesvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
870 $ rcond, r1, r2, w, 4, rw, info )
871 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
873 CALL zhesvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
874 $ rcond, r1, r2, w, 3, rw, info )
875 CALL chkxer(
'ZHESVX', infot, nout, lerr, ok )
883 CALL zhesvxx(
'/',
'U', 0, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
884 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
885 $ err_bnds_c, nparams, params, w, rw, info )
886 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
888 CALL zhesvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
889 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
890 $ err_bnds_c, nparams, params, w, rw, info )
891 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
893 CALL zhesvxx(
'N',
'U', -1, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
894 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
895 $ err_bnds_c, nparams, params, w, rw, info )
896 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
898 CALL zhesvxx(
'N',
'U', 0, -1, a, 1, af, 1, ip, eq, c, b, 1, x,
899 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
900 $ err_bnds_c, nparams, params, w, rw, info )
901 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
903 CALL zhesvxx(
'N',
'U', 2, 0, a, 1, af, 2, ip, eq, c, b, 2, x,
904 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
905 $ err_bnds_c, nparams, params, w, rw, info )
906 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
908 CALL zhesvxx(
'N',
'U', 2, 0, a, 2, af, 1, ip, eq, c, b, 2, x,
909 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
910 $ err_bnds_c, nparams, params, w, rw, info )
911 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
914 CALL zhesvxx(
'F',
'U', 0, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
915 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
916 $ err_bnds_c, nparams, params, w, rw, info )
917 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
920 CALL zhesvxx(
'F',
'U', 1, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
921 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
922 $ err_bnds_c, nparams, params, w, rw, info )
923 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
925 CALL zhesvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, c, b, 1, x,
926 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
927 $ err_bnds_c, nparams, params, w, rw, info )
928 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
930 CALL zhesvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, c, b, 2, x,
931 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
932 $ err_bnds_c, nparams, params, w, rw, info )
933 CALL chkxer(
'ZHESVXX', infot, nout, lerr, ok )
935 ELSE IF( lsamen( 2, c2,
'HR' ) )
THEN
939 srnamt =
'ZHESV_ROOK'
941 CALL zhesv_rook(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
942 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
944 CALL zhesv_rook(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
945 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
947 CALL zhesv_rook(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
948 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
950 CALL zhesv_rook(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
951 CALL chkxer(
'ZHESV_ROOK', infot, nout, lerr, ok )
953 ELSE IF( lsamen( 2, c2,
'HP' ) )
THEN
959 CALL zhpsv(
'/', 0, 0, a, ip, b, 1, info )
960 CALL chkxer(
'ZHPSV ', infot, nout, lerr, ok )
962 CALL zhpsv(
'U', -1, 0, a, ip, b, 1, info )
963 CALL chkxer(
'ZHPSV ', infot, nout, lerr, ok )
965 CALL zhpsv(
'U', 0, -1, a, ip, b, 1, info )
966 CALL chkxer(
'ZHPSV ', infot, nout, lerr, ok )
968 CALL zhpsv(
'U', 2, 0, a, ip, b, 1, info )
969 CALL chkxer(
'ZHPSV ', infot, nout, lerr, ok )
975 CALL zhpsvx(
'/',
'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
977 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
979 CALL zhpsvx(
'N',
'/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
981 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
983 CALL zhpsvx(
'N',
'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
985 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
987 CALL zhpsvx(
'N',
'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
989 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
991 CALL zhpsvx(
'N',
'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
993 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
995 CALL zhpsvx(
'N',
'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
997 CALL chkxer(
'ZHPSVX', infot, nout, lerr, ok )
999 ELSE IF( lsamen( 2, c2,
'SY' ) )
THEN
1005 CALL zsysv(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
1006 CALL chkxer(
'ZSYSV ', infot, nout, lerr, ok )
1008 CALL zsysv(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
1009 CALL chkxer(
'ZSYSV ', infot, nout, lerr, ok )
1011 CALL zsysv(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
1012 CALL chkxer(
'ZSYSV ', infot, nout, lerr, ok )
1014 CALL zsysv(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
1015 CALL chkxer(
'ZSYSV ', infot, nout, lerr, ok )
1021 CALL zsysvx(
'/',
'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
1022 $ rcond, r1, r2, w, 1, rw, info )
1023 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1025 CALL zsysvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
1026 $ rcond, r1, r2, w, 1, rw, info )
1027 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1029 CALL zsysvx(
'N',
'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
1030 $ rcond, r1, r2, w, 1, rw, info )
1031 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1033 CALL zsysvx(
'N',
'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
1034 $ rcond, r1, r2, w, 1, rw, info )
1035 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1037 CALL zsysvx(
'N',
'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
1038 $ rcond, r1, r2, w, 4, rw, info )
1039 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1041 CALL zsysvx(
'N',
'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
1042 $ rcond, r1, r2, w, 4, rw, info )
1043 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1045 CALL zsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
1046 $ rcond, r1, r2, w, 4, rw, info )
1047 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1049 CALL zsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
1050 $ rcond, r1, r2, w, 4, rw, info )
1051 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1053 CALL zsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
1054 $ rcond, r1, r2, w, 3, rw, info )
1055 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1064 CALL zsysvxx(
'/',
'U', 0, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
1065 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1066 $ err_bnds_c, nparams, params, w, rw, info )
1067 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1069 CALL zsysvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
1070 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1071 $ err_bnds_c, nparams, params, w, rw, info )
1072 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1074 CALL zsysvxx(
'N',
'U', -1, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
1075 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1076 $ err_bnds_c, nparams, params, w, rw, info )
1077 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1080 CALL zsysvxx(
'N',
'U', 0, -1, a, 1, af, 1, ip, eq, r, b, 1, x,
1081 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1082 $ err_bnds_c, nparams, params, w, rw, info )
1083 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1086 CALL zsysvxx(
'N',
'U', 2, 0, a, 1, af, 2, ip, eq, r, b, 2, x,
1087 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1088 $ err_bnds_c, nparams, params, w, rw, info )
1089 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1091 CALL zsysvxx(
'N',
'U', 2, 0, a, 2, af, 1, ip, eq, r, b, 2, x,
1092 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1093 $ err_bnds_c, nparams, params, w, rw, info )
1094 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1097 CALL zsysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 1, x,
1098 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1099 $ err_bnds_c, nparams, params, w, rw, info )
1100 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1102 CALL zsysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 2, x,
1103 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1104 $ err_bnds_c, nparams, params, w, rw, info )
1105 CALL chkxer(
'ZSYSVXX', infot, nout, lerr, ok )
1106 CALL chkxer(
'ZSYSVX', infot, nout, lerr, ok )
1108 ELSE IF( lsamen( 2, c2,
'SR' ) )
THEN
1112 srnamt =
'ZSYSV_ROOK'
1114 CALL zsysv_rook(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
1115 CALL chkxer(
'ZSYSV_ROOK', infot, nout, lerr, ok )
1117 CALL zsysv_rook(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
1118 CALL chkxer(
'ZSYSV_ROOK', infot, nout, lerr, ok )
1120 CALL zsysv_rook(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
1121 CALL chkxer(
'ZSYSV_ROOK', infot, nout, lerr, ok )
1123 CALL zsysv_rook(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
1124 CALL chkxer(
'ZSYSV_ROOK', infot, nout, lerr, ok )
1126 ELSE IF( lsamen( 2, c2,
'SP' ) )
THEN
1132 CALL zspsv(
'/', 0, 0, a, ip, b, 1, info )
1133 CALL chkxer(
'ZSPSV ', infot, nout, lerr, ok )
1135 CALL zspsv(
'U', -1, 0, a, ip, b, 1, info )
1136 CALL chkxer(
'ZSPSV ', infot, nout, lerr, ok )
1138 CALL zspsv(
'U', 0, -1, a, ip, b, 1, info )
1139 CALL chkxer(
'ZSPSV ', infot, nout, lerr, ok )
1141 CALL zspsv(
'U', 2, 0, a, ip, b, 1, info )
1142 CALL chkxer(
'ZSPSV ', infot, nout, lerr, ok )
1148 CALL zspsvx(
'/',
'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1150 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1152 CALL zspsvx(
'N',
'/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1154 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1156 CALL zspsvx(
'N',
'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1158 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1160 CALL zspsvx(
'N',
'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
1162 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1164 CALL zspsvx(
'N',
'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
1166 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1168 CALL zspsvx(
'N',
'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
1170 CALL chkxer(
'ZSPSVX', infot, nout, lerr, ok )
1176 WRITE( nout, fmt = 9999 )path
1178 WRITE( nout, fmt = 9998 )path
1181 9999
FORMAT( 1x, a3,
' drivers passed the tests of the error exits' )
1182 9998
FORMAT(
' *** ', a3,
' drivers failed the tests of the error ',
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 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 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 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 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 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 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 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 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 zppsv(UPLO, N, NRHS, AP, B, LDB, INFO)
ZPPSV computes the solution to system of linear equations A * X = B for OTHER 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 chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
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 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 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 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 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) ...
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 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 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 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 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 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 zgtsv(N, NRHS, DL, D, DU, B, LDB, INFO)
ZGTSV computes the solution to system of linear equations A * X = B for GT 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 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 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 zptsv(N, NRHS, D, E, B, LDB, INFO)
ZPTSV computes the solution to system of linear equations A * X = B for PT 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 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 ...