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( 1
x, 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)
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
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 ...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
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 ...
LOGICAL function lsamen(N, CA, CB)
LSAMEN
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 ...