56 SUBROUTINE cerrvx( PATH, NUNIT )
77 INTEGER I, INFO, J, N_ERR_BNDS, NPARAMS
78 REAL RCOND, RPVGRW, BERR
82 REAL 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 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 cmplx, real
115 WRITE( nout, fmt = * )
122 a( i, j ) = cmplx( 1. /
REAL( I+J ), -1. /
REAL( I+J ) )
123 af( i, j ) = cmplx( 1. /
REAL( I+J ), -1. /
REAL( I+J ) )
137 IF( lsamen( 2, c2,
'GE' ) )
THEN
143 CALL cgesv( -1, 0, a, 1, ip, b, 1, info )
144 CALL chkxer(
'CGESV ', infot, nout, lerr, ok )
146 CALL cgesv( 0, -1, a, 1, ip, b, 1, info )
147 CALL chkxer(
'CGESV ', infot, nout, lerr, ok )
149 CALL cgesv( 2, 1, a, 1, ip, b, 2, info )
150 CALL chkxer(
'CGESV ', infot, nout, lerr, ok )
152 CALL cgesv( 2, 1, a, 2, ip, b, 1, info )
153 CALL chkxer(
'CGESV ', infot, nout, lerr, ok )
159 CALL cgesvx(
'/',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
160 $ x, 1, rcond, r1, r2, w, rw, info )
161 CALL chkxer(
'CGESVX', infot, nout, lerr, ok )
163 CALL cgesvx(
'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(
'CGESVX', infot, nout, lerr, ok )
167 CALL cgesvx(
'N',
'N', -1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
168 $ x, 1, rcond, r1, r2, w, rw, info )
169 CALL chkxer(
'CGESVX', infot, nout, lerr, ok )
171 CALL cgesvx(
'N',
'N', 0, -1, a, 1, af, 1, ip, eq, r, c, b, 1,
172 $ x, 1, rcond, r1, r2, w, rw, info )
173 CALL chkxer(
'CGESVX', infot, nout, lerr, ok )
175 CALL cgesvx(
'N',
'N', 2, 1, a, 1, af, 2, ip, eq, r, c, b, 2,
176 $ x, 2, rcond, r1, r2, w, rw, info )
177 CALL chkxer(
'CGESVX', infot, nout, lerr, ok )
179 CALL cgesvx(
'N',
'N', 2, 1, a, 2, af, 1, ip, eq, r, c, b, 2,
180 $ x, 2, rcond, r1, r2, w, rw, info )
181 CALL chkxer(
'CGESVX', infot, nout, lerr, ok )
184 CALL cgesvx(
'F',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
185 $ x, 1, rcond, r1, r2, w, rw, info )
186 CALL chkxer(
'CGESVX', infot, nout, lerr, ok )
189 CALL cgesvx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
190 $ x, 1, rcond, r1, r2, w, rw, info )
191 CALL chkxer(
'CGESVX', infot, nout, lerr, ok )
194 CALL cgesvx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
195 $ x, 1, rcond, r1, r2, w, rw, info )
196 CALL chkxer(
'CGESVX', infot, nout, lerr, ok )
198 CALL cgesvx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 1,
199 $ x, 2, rcond, r1, r2, w, rw, info )
200 CALL chkxer(
'CGESVX', infot, nout, lerr, ok )
202 CALL cgesvx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 2,
203 $ x, 1, rcond, r1, r2, w, rw, info )
204 CALL chkxer(
'CGESVX', infot, nout, lerr, ok )
212 CALL cgesvxx(
'/',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
213 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
214 $ err_bnds_c, nparams, params, w, rw, info )
215 CALL chkxer(
'CGESVXX', infot, nout, lerr, ok )
217 CALL cgesvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
218 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
219 $ err_bnds_c, nparams, params, w, rw, info )
220 CALL chkxer(
'CGESVXX', infot, nout, lerr, ok )
222 CALL cgesvxx(
'N',
'N', -1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
223 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
224 $ err_bnds_c, nparams, params, w, rw, info )
225 CALL chkxer(
'CGESVXX', infot, nout, lerr, ok )
227 CALL cgesvxx(
'N',
'N', 0, -1, a, 1, af, 1, ip, eq, r, c, b, 1,
228 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
229 $ err_bnds_c, nparams, params, w, rw, info )
230 CALL chkxer(
'CGESVXX', infot, nout, lerr, ok )
232 CALL cgesvxx(
'N',
'N', 2, 1, a, 1, af, 2, ip, eq, r, c, b, 2,
233 $ x, 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
234 $ err_bnds_c, nparams, params, w, rw, info )
235 CALL chkxer(
'CGESVXX', infot, nout, lerr, ok )
237 CALL cgesvxx(
'N',
'N', 2, 1, a, 2, af, 1, ip, eq, r, c, b, 2,
238 $ x, 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
239 $ err_bnds_c, nparams, params, w, rw, info )
240 CALL chkxer(
'CGESVXX', infot, nout, lerr, ok )
243 CALL cgesvxx(
'F',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
244 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
245 $ err_bnds_c, nparams, params, w, rw, info )
246 CALL chkxer(
'CGESVXX', infot, nout, lerr, ok )
249 CALL cgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
250 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
251 $ err_bnds_c, nparams, params, w, rw, info )
252 CALL chkxer(
'CGESVXX', infot, nout, lerr, ok )
255 CALL cgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
256 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
257 $ err_bnds_c, nparams, params, w, rw, info )
258 CALL chkxer(
'CGESVXX', infot, nout, lerr, ok )
260 CALL cgesvxx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 1,
261 $ x, 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
262 $ err_bnds_c, nparams, params, w, rw, info )
263 CALL chkxer(
'CGESVXX', infot, nout, lerr, ok )
265 CALL cgesvxx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 2,
266 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
267 $ err_bnds_c, nparams, params, w, rw, info )
268 CALL chkxer(
'CGESVXX', infot, nout, lerr, ok )
270 ELSE IF( lsamen( 2, c2,
'GB' ) )
THEN
276 CALL cgbsv( -1, 0, 0, 0, a, 1, ip, b, 1, info )
277 CALL chkxer(
'CGBSV ', infot, nout, lerr, ok )
279 CALL cgbsv( 1, -1, 0, 0, a, 1, ip, b, 1, info )
280 CALL chkxer(
'CGBSV ', infot, nout, lerr, ok )
282 CALL cgbsv( 1, 0, -1, 0, a, 1, ip, b, 1, info )
283 CALL chkxer(
'CGBSV ', infot, nout, lerr, ok )
285 CALL cgbsv( 0, 0, 0, -1, a, 1, ip, b, 1, info )
286 CALL chkxer(
'CGBSV ', infot, nout, lerr, ok )
288 CALL cgbsv( 1, 1, 1, 0, a, 3, ip, b, 1, info )
289 CALL chkxer(
'CGBSV ', infot, nout, lerr, ok )
291 CALL cgbsv( 2, 0, 0, 0, a, 1, ip, b, 1, info )
292 CALL chkxer(
'CGBSV ', infot, nout, lerr, ok )
298 CALL cgbsvx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
299 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
300 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
302 CALL cgbsvx(
'N',
'/', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
303 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
304 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
306 CALL cgbsvx(
'N',
'N', -1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
307 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
308 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
310 CALL cgbsvx(
'N',
'N', 1, -1, 0, 0, a, 1, af, 1, ip, eq, r, c,
311 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
312 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
314 CALL cgbsvx(
'N',
'N', 1, 0, -1, 0, a, 1, af, 1, ip, eq, r, c,
315 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
316 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
318 CALL cgbsvx(
'N',
'N', 0, 0, 0, -1, a, 1, af, 1, ip, eq, r, c,
319 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
320 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
322 CALL cgbsvx(
'N',
'N', 1, 1, 1, 0, a, 2, af, 4, ip, eq, r, c,
323 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
324 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
326 CALL cgbsvx(
'N',
'N', 1, 1, 1, 0, a, 3, af, 3, ip, eq, r, c,
327 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
328 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
331 CALL cgbsvx(
'F',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
332 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
333 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
336 CALL cgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
337 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
338 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
341 CALL cgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
342 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
343 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
345 CALL cgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
346 $ b, 1, x, 2, rcond, r1, r2, w, rw, info )
347 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
349 CALL cgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
350 $ b, 2, x, 1, rcond, r1, r2, w, rw, info )
351 CALL chkxer(
'CGBSVX', infot, nout, lerr, ok )
359 CALL cgbsvxx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
360 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
361 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
363 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
365 CALL cgbsvxx(
'N',
'/', 0, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
366 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
367 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
369 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
371 CALL cgbsvxx(
'N',
'N', -1, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
372 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
373 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
375 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
377 CALL cgbsvxx(
'N',
'N', 2, -1, 1, 0, a, 1, af, 1, ip, eq,
378 $ r, c, b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
379 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
381 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
383 CALL cgbsvxx(
'N',
'N', 2, 1, -1, 0, a, 1, af, 1, ip, eq,
384 $ r, c, b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
385 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
387 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
389 CALL cgbsvxx(
'N',
'N', 0, 1, 1, -1, a, 1, af, 1, ip, eq, r, c,
390 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
391 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
393 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
395 CALL cgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 2, af, 2, ip, eq, r, c,
396 $ b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
397 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
399 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
401 CALL cgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 3, ip, eq, r, c,
402 $ b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
403 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
405 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
408 CALL cgbsvxx(
'F',
'N', 0, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
409 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
410 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
412 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
415 CALL cgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
416 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
417 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
419 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
422 CALL cgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
423 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
424 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
426 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
428 CALL cgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
429 $ b, 1, x, 2, rcond, rpvgrw, berr, n_err_bnds,
430 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
432 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
434 CALL cgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
435 $ b, 2, x, 1, rcond, rpvgrw, berr, n_err_bnds,
436 $ err_bnds_n, err_bnds_c, nparams, params, w, rw,
438 CALL chkxer(
'CGBSVXX', infot, nout, lerr, ok )
440 ELSE IF( lsamen( 2, c2,
'GT' ) )
THEN
446 CALL cgtsv( -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
448 CALL chkxer(
'CGTSV ', infot, nout, lerr, ok )
450 CALL cgtsv( 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
452 CALL chkxer(
'CGTSV ', infot, nout, lerr, ok )
454 CALL cgtsv( 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1, info )
455 CALL chkxer(
'CGTSV ', infot, nout, lerr, ok )
461 CALL cgtsvx(
'/',
'N', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
462 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
463 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
464 CALL chkxer(
'CGTSVX', infot, nout, lerr, ok )
466 CALL cgtsvx(
'N',
'/', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
467 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
468 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
469 CALL chkxer(
'CGTSVX', infot, nout, lerr, ok )
471 CALL cgtsvx(
'N',
'N', -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
472 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
473 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
474 CALL chkxer(
'CGTSVX', infot, nout, lerr, ok )
476 CALL cgtsvx(
'N',
'N', 0, -1, 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(
'CGTSVX', infot, nout, lerr, ok )
481 CALL cgtsvx(
'N',
'N', 2, 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, 2, rcond, r1, r2, w, rw, info )
484 CALL chkxer(
'CGTSVX', infot, nout, lerr, ok )
486 CALL cgtsvx(
'N',
'N', 2, 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, 2, x, 1, rcond, r1, r2, w, rw, info )
489 CALL chkxer(
'CGTSVX', infot, nout, lerr, ok )
491 ELSE IF( lsamen( 2, c2,
'PO' ) )
THEN
497 CALL cposv(
'/', 0, 0, a, 1, b, 1, info )
498 CALL chkxer(
'CPOSV ', infot, nout, lerr, ok )
500 CALL cposv(
'U', -1, 0, a, 1, b, 1, info )
501 CALL chkxer(
'CPOSV ', infot, nout, lerr, ok )
503 CALL cposv(
'U', 0, -1, a, 1, b, 1, info )
504 CALL chkxer(
'CPOSV ', infot, nout, lerr, ok )
506 CALL cposv(
'U', 2, 0, a, 1, b, 2, info )
507 CALL chkxer(
'CPOSV ', infot, nout, lerr, ok )
509 CALL cposv(
'U', 2, 0, a, 2, b, 1, info )
510 CALL chkxer(
'CPOSV ', infot, nout, lerr, ok )
516 CALL cposvx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
517 $ rcond, r1, r2, w, rw, info )
518 CALL chkxer(
'CPOSVX', infot, nout, lerr, ok )
520 CALL cposvx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
521 $ rcond, r1, r2, w, rw, info )
522 CALL chkxer(
'CPOSVX', infot, nout, lerr, ok )
524 CALL cposvx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
525 $ rcond, r1, r2, w, rw, info )
526 CALL chkxer(
'CPOSVX', infot, nout, lerr, ok )
528 CALL cposvx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
529 $ rcond, r1, r2, w, rw, info )
530 CALL chkxer(
'CPOSVX', infot, nout, lerr, ok )
532 CALL cposvx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
533 $ rcond, r1, r2, w, rw, info )
534 CALL chkxer(
'CPOSVX', infot, nout, lerr, ok )
536 CALL cposvx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
537 $ rcond, r1, r2, w, rw, info )
538 CALL chkxer(
'CPOSVX', infot, nout, lerr, ok )
541 CALL cposvx(
'F',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
542 $ rcond, r1, r2, w, rw, info )
543 CALL chkxer(
'CPOSVX', infot, nout, lerr, ok )
546 CALL cposvx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
547 $ rcond, r1, r2, w, rw, info )
548 CALL chkxer(
'CPOSVX', infot, nout, lerr, ok )
550 CALL cposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 1, x, 2,
551 $ rcond, r1, r2, w, rw, info )
552 CALL chkxer(
'CPOSVX', infot, nout, lerr, ok )
554 CALL cposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
555 $ rcond, r1, r2, w, rw, info )
556 CALL chkxer(
'CPOSVX', infot, nout, lerr, ok )
564 CALL cposvxx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
565 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
566 $ err_bnds_c, nparams, params, w, rw, info )
567 CALL chkxer(
'CPOSVXX', infot, nout, lerr, ok )
569 CALL cposvxx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
570 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
571 $ err_bnds_c, nparams, params, w, rw, info )
572 CALL chkxer(
'CPOSVXX', infot, nout, lerr, ok )
574 CALL cposvxx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
575 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
576 $ err_bnds_c, nparams, params, w, rw, info )
577 CALL chkxer(
'CPOSVXX', infot, nout, lerr, ok )
579 CALL cposvxx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
580 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
581 $ err_bnds_c, nparams, params, w, rw, info )
582 CALL chkxer(
'CPOSVXX', infot, nout, lerr, ok )
584 CALL cposvxx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
585 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
586 $ err_bnds_c, nparams, params, w, rw, info )
587 CALL chkxer(
'CPOSVXX', infot, nout, lerr, ok )
589 CALL cposvxx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
590 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
591 $ err_bnds_c, nparams, params, w, rw, info )
592 CALL chkxer(
'CPOSVXX', infot, nout, lerr, ok )
595 CALL cposvxx(
'F',
'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(
'CPOSVXX', infot, nout, lerr, ok )
601 CALL cposvxx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
602 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
603 $ err_bnds_c, nparams, params, w, rw, info )
604 CALL chkxer(
'CPOSVXX', infot, nout, lerr, ok )
606 CALL cposvxx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 1, x, 2,
607 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
608 $ err_bnds_c, nparams, params, w, rw, info )
609 CALL chkxer(
'CPOSVXX', infot, nout, lerr, ok )
611 CALL cposvxx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
612 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
613 $ err_bnds_c, nparams, params, w, rw, info )
614 CALL chkxer(
'CPOSVXX', infot, nout, lerr, ok )
616 ELSE IF( lsamen( 2, c2,
'PP' ) )
THEN
622 CALL cppsv(
'/', 0, 0, a, b, 1, info )
623 CALL chkxer(
'CPPSV ', infot, nout, lerr, ok )
625 CALL cppsv(
'U', -1, 0, a, b, 1, info )
626 CALL chkxer(
'CPPSV ', infot, nout, lerr, ok )
628 CALL cppsv(
'U', 0, -1, a, b, 1, info )
629 CALL chkxer(
'CPPSV ', infot, nout, lerr, ok )
631 CALL cppsv(
'U', 2, 0, a, b, 1, info )
632 CALL chkxer(
'CPPSV ', infot, nout, lerr, ok )
638 CALL cppsvx(
'/',
'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
639 $ r1, r2, w, rw, info )
640 CALL chkxer(
'CPPSVX', infot, nout, lerr, ok )
642 CALL cppsvx(
'N',
'/', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
643 $ r1, r2, w, rw, info )
644 CALL chkxer(
'CPPSVX', infot, nout, lerr, ok )
646 CALL cppsvx(
'N',
'U', -1, 0, a, af, eq, c, b, 1, x, 1, rcond,
647 $ r1, r2, w, rw, info )
648 CALL chkxer(
'CPPSVX', infot, nout, lerr, ok )
650 CALL cppsvx(
'N',
'U', 0, -1, a, af, eq, c, b, 1, x, 1, rcond,
651 $ r1, r2, w, rw, info )
652 CALL chkxer(
'CPPSVX', infot, nout, lerr, ok )
655 CALL cppsvx(
'F',
'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
656 $ r1, r2, w, rw, info )
657 CALL chkxer(
'CPPSVX', infot, nout, lerr, ok )
660 CALL cppsvx(
'F',
'U', 1, 0, a, af, eq, c, b, 1, x, 1, rcond,
661 $ r1, r2, w, rw, info )
662 CALL chkxer(
'CPPSVX', infot, nout, lerr, ok )
664 CALL cppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 1, x, 2, rcond,
665 $ r1, r2, w, rw, info )
666 CALL chkxer(
'CPPSVX', infot, nout, lerr, ok )
668 CALL cppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 2, x, 1, rcond,
669 $ r1, r2, w, rw, info )
670 CALL chkxer(
'CPPSVX', infot, nout, lerr, ok )
672 ELSE IF( lsamen( 2, c2,
'PB' ) )
THEN
678 CALL cpbsv(
'/', 0, 0, 0, a, 1, b, 1, info )
679 CALL chkxer(
'CPBSV ', infot, nout, lerr, ok )
681 CALL cpbsv(
'U', -1, 0, 0, a, 1, b, 1, info )
682 CALL chkxer(
'CPBSV ', infot, nout, lerr, ok )
684 CALL cpbsv(
'U', 1, -1, 0, a, 1, b, 1, info )
685 CALL chkxer(
'CPBSV ', infot, nout, lerr, ok )
687 CALL cpbsv(
'U', 0, 0, -1, a, 1, b, 1, info )
688 CALL chkxer(
'CPBSV ', infot, nout, lerr, ok )
690 CALL cpbsv(
'U', 1, 1, 0, a, 1, b, 2, info )
691 CALL chkxer(
'CPBSV ', infot, nout, lerr, ok )
693 CALL cpbsv(
'U', 2, 0, 0, a, 1, b, 1, info )
694 CALL chkxer(
'CPBSV ', infot, nout, lerr, ok )
700 CALL cpbsvx(
'/',
'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
701 $ rcond, r1, r2, w, rw, info )
702 CALL chkxer(
'CPBSVX', infot, nout, lerr, ok )
704 CALL cpbsvx(
'N',
'/', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
705 $ rcond, r1, r2, w, rw, info )
706 CALL chkxer(
'CPBSVX', infot, nout, lerr, ok )
708 CALL cpbsvx(
'N',
'U', -1, 0, 0, a, 1, af, 1, eq, c, b, 1, x,
709 $ 1, rcond, r1, r2, w, rw, info )
710 CALL chkxer(
'CPBSVX', infot, nout, lerr, ok )
712 CALL cpbsvx(
'N',
'U', 1, -1, 0, a, 1, af, 1, eq, c, b, 1, x,
713 $ 1, rcond, r1, r2, w, rw, info )
714 CALL chkxer(
'CPBSVX', infot, nout, lerr, ok )
716 CALL cpbsvx(
'N',
'U', 0, 0, -1, a, 1, af, 1, eq, c, b, 1, x,
717 $ 1, rcond, r1, r2, w, rw, info )
718 CALL chkxer(
'CPBSVX', infot, nout, lerr, ok )
720 CALL cpbsvx(
'N',
'U', 1, 1, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
721 $ rcond, r1, r2, w, rw, info )
722 CALL chkxer(
'CPBSVX', infot, nout, lerr, ok )
724 CALL cpbsvx(
'N',
'U', 1, 1, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
725 $ rcond, r1, r2, w, rw, info )
726 CALL chkxer(
'CPBSVX', infot, nout, lerr, ok )
729 CALL cpbsvx(
'F',
'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
730 $ rcond, r1, r2, w, rw, info )
731 CALL chkxer(
'CPBSVX', infot, nout, lerr, ok )
734 CALL cpbsvx(
'F',
'U', 1, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
735 $ rcond, r1, r2, w, rw, info )
736 CALL chkxer(
'CPBSVX', infot, nout, lerr, ok )
738 CALL cpbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 2,
739 $ rcond, r1, r2, w, rw, info )
740 CALL chkxer(
'CPBSVX', infot, nout, lerr, ok )
742 CALL cpbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 2, x, 1,
743 $ rcond, r1, r2, w, rw, info )
744 CALL chkxer(
'CPBSVX', infot, nout, lerr, ok )
746 ELSE IF( lsamen( 2, c2,
'PT' ) )
THEN
752 CALL cptsv( -1, 0, r, a( 1, 1 ), b, 1, info )
753 CALL chkxer(
'CPTSV ', infot, nout, lerr, ok )
755 CALL cptsv( 0, -1, r, a( 1, 1 ), b, 1, info )
756 CALL chkxer(
'CPTSV ', infot, nout, lerr, ok )
758 CALL cptsv( 2, 0, r, a( 1, 1 ), b, 1, info )
759 CALL chkxer(
'CPTSV ', infot, nout, lerr, ok )
765 CALL cptsvx(
'/', 0, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
766 $ 1, rcond, r1, r2, w, rw, info )
767 CALL chkxer(
'CPTSVX', infot, nout, lerr, ok )
769 CALL cptsvx(
'N', -1, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
770 $ 1, rcond, r1, r2, w, rw, info )
771 CALL chkxer(
'CPTSVX', infot, nout, lerr, ok )
773 CALL cptsvx(
'N', 0, -1, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
774 $ 1, rcond, r1, r2, w, rw, info )
775 CALL chkxer(
'CPTSVX', infot, nout, lerr, ok )
777 CALL cptsvx(
'N', 2, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
778 $ 2, rcond, r1, r2, w, rw, info )
779 CALL chkxer(
'CPTSVX', infot, nout, lerr, ok )
781 CALL cptsvx(
'N', 2, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 2, x,
782 $ 1, rcond, r1, r2, w, rw, info )
783 CALL chkxer(
'CPTSVX', infot, nout, lerr, ok )
785 ELSE IF( lsamen( 2, c2,
'HE' ) )
THEN
791 CALL chesv(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
792 CALL chkxer(
'CHESV ', infot, nout, lerr, ok )
794 CALL chesv(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
795 CALL chkxer(
'CHESV ', infot, nout, lerr, ok )
797 CALL chesv(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
798 CALL chkxer(
'CHESV ', infot, nout, lerr, ok )
800 CALL chesv(
'U', 2, 0, a, 1, ip, b, 2, w, 1, info )
801 CALL chkxer(
'CHESV ', infot, nout, lerr, ok )
803 CALL chesv(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
804 CALL chkxer(
'CHESV ', infot, nout, lerr, ok )
810 CALL chesvx(
'/',
'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
811 $ rcond, r1, r2, w, 1, rw, info )
812 CALL chkxer(
'CHESVX', infot, nout, lerr, ok )
814 CALL chesvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
815 $ rcond, r1, r2, w, 1, rw, info )
816 CALL chkxer(
'CHESVX', infot, nout, lerr, ok )
818 CALL chesvx(
'N',
'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
819 $ rcond, r1, r2, w, 1, rw, info )
820 CALL chkxer(
'CHESVX', infot, nout, lerr, ok )
822 CALL chesvx(
'N',
'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
823 $ rcond, r1, r2, w, 1, rw, info )
824 CALL chkxer(
'CHESVX', infot, nout, lerr, ok )
826 CALL chesvx(
'N',
'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
827 $ rcond, r1, r2, w, 4, rw, info )
828 CALL chkxer(
'CHESVX', infot, nout, lerr, ok )
830 CALL chesvx(
'N',
'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
831 $ rcond, r1, r2, w, 4, rw, info )
832 CALL chkxer(
'CHESVX', infot, nout, lerr, ok )
834 CALL chesvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
835 $ rcond, r1, r2, w, 4, rw, info )
836 CALL chkxer(
'CHESVX', infot, nout, lerr, ok )
838 CALL chesvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
839 $ rcond, r1, r2, w, 4, rw, info )
840 CALL chkxer(
'CHESVX', infot, nout, lerr, ok )
842 CALL chesvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
843 $ rcond, r1, r2, w, 3, rw, info )
844 CALL chkxer(
'CHESVX', infot, nout, lerr, ok )
852 CALL chesvxx(
'/',
'U', 0, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
853 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
854 $ err_bnds_c, nparams, params, w, rw, info )
855 CALL chkxer(
'CHESVXX', infot, nout, lerr, ok )
857 CALL chesvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
858 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
859 $ err_bnds_c, nparams, params, w, rw, info )
860 CALL chkxer(
'CHESVXX', infot, nout, lerr, ok )
862 CALL chesvxx(
'N',
'U', -1, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
863 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
864 $ err_bnds_c, nparams, params, w, rw, info )
865 CALL chkxer(
'CHESVXX', infot, nout, lerr, ok )
867 CALL chesvxx(
'N',
'U', 0, -1, a, 1, af, 1, ip, eq, c, b, 1, x,
868 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
869 $ err_bnds_c, nparams, params, w, rw, info )
870 CALL chkxer(
'CHESVXX', infot, nout, lerr, ok )
872 CALL chesvxx(
'N',
'U', 2, 0, a, 1, af, 2, ip, eq, c, b, 2, x,
873 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
874 $ err_bnds_c, nparams, params, w, rw, info )
875 CALL chkxer(
'CHESVXX', infot, nout, lerr, ok )
877 CALL chesvxx(
'N',
'U', 2, 0, a, 2, af, 1, ip, eq, c, b, 2, x,
878 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
879 $ err_bnds_c, nparams, params, w, rw, info )
880 CALL chkxer(
'CHESVXX', infot, nout, lerr, ok )
883 CALL chesvxx(
'F',
'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(
'CHESVXX', infot, nout, lerr, ok )
889 CALL chesvxx(
'F',
'U', 1, 0, a, 1, af, 1, ip, eq, c, b, 1, x,
890 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
891 $ err_bnds_c, nparams, params, w, rw, info )
892 CALL chkxer(
'CHESVXX', infot, nout, lerr, ok )
894 CALL chesvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, c, b, 1, x,
895 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
896 $ err_bnds_c, nparams, params, w, rw, info )
897 CALL chkxer(
'CHESVXX', infot, nout, lerr, ok )
899 CALL chesvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, c, b, 2, x,
900 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
901 $ err_bnds_c, nparams, params, w, rw, info )
902 CALL chkxer(
'CHESVXX', infot, nout, lerr, ok )
904 ELSE IF( lsamen( 2, c2,
'HR' ) )
THEN
908 srnamt =
'CHESV_ROOK'
910 CALL chesv_rook(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
911 CALL chkxer(
'CHESV_ROOK', infot, nout, lerr, ok )
913 CALL chesv_rook(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
914 CALL chkxer(
'CHESV_ROOK', infot, nout, lerr, ok )
916 CALL chesv_rook(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
917 CALL chkxer(
'CHESV_ROOK', infot, nout, lerr, ok )
919 CALL chesv_rook(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
920 CALL chkxer(
'CHESV_ROOK', infot, nout, lerr, ok )
922 ELSE IF( lsamen( 2, c2,
'HP' ) )
THEN
928 CALL chpsv(
'/', 0, 0, a, ip, b, 1, info )
929 CALL chkxer(
'CHPSV ', infot, nout, lerr, ok )
931 CALL chpsv(
'U', -1, 0, a, ip, b, 1, info )
932 CALL chkxer(
'CHPSV ', infot, nout, lerr, ok )
934 CALL chpsv(
'U', 0, -1, a, ip, b, 1, info )
935 CALL chkxer(
'CHPSV ', infot, nout, lerr, ok )
937 CALL chpsv(
'U', 2, 0, a, ip, b, 1, info )
938 CALL chkxer(
'CHPSV ', infot, nout, lerr, ok )
944 CALL chpsvx(
'/',
'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
946 CALL chkxer(
'CHPSVX', infot, nout, lerr, ok )
948 CALL chpsvx(
'N',
'/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
950 CALL chkxer(
'CHPSVX', infot, nout, lerr, ok )
952 CALL chpsvx(
'N',
'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
954 CALL chkxer(
'CHPSVX', infot, nout, lerr, ok )
956 CALL chpsvx(
'N',
'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
958 CALL chkxer(
'CHPSVX', infot, nout, lerr, ok )
960 CALL chpsvx(
'N',
'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
962 CALL chkxer(
'CHPSVX', infot, nout, lerr, ok )
964 CALL chpsvx(
'N',
'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
966 CALL chkxer(
'CHPSVX', infot, nout, lerr, ok )
968 ELSE IF( lsamen( 2, c2,
'SY' ) )
THEN
974 CALL csysv(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
975 CALL chkxer(
'CSYSV ', infot, nout, lerr, ok )
977 CALL csysv(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
978 CALL chkxer(
'CSYSV ', infot, nout, lerr, ok )
980 CALL csysv(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
981 CALL chkxer(
'CSYSV ', infot, nout, lerr, ok )
983 CALL csysv(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
984 CALL chkxer(
'CSYSV ', infot, nout, lerr, ok )
990 CALL csysvx(
'/',
'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
991 $ rcond, r1, r2, w, 1, rw, info )
992 CALL chkxer(
'CSYSVX', infot, nout, lerr, ok )
994 CALL csysvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
995 $ rcond, r1, r2, w, 1, rw, info )
996 CALL chkxer(
'CSYSVX', infot, nout, lerr, ok )
998 CALL csysvx(
'N',
'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
999 $ rcond, r1, r2, w, 1, rw, info )
1000 CALL chkxer(
'CSYSVX', infot, nout, lerr, ok )
1002 CALL csysvx(
'N',
'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
1003 $ rcond, r1, r2, w, 1, rw, info )
1004 CALL chkxer(
'CSYSVX', infot, nout, lerr, ok )
1006 CALL csysvx(
'N',
'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
1007 $ rcond, r1, r2, w, 4, rw, info )
1008 CALL chkxer(
'CSYSVX', infot, nout, lerr, ok )
1010 CALL csysvx(
'N',
'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
1011 $ rcond, r1, r2, w, 4, rw, info )
1012 CALL chkxer(
'CSYSVX', infot, nout, lerr, ok )
1014 CALL csysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
1015 $ rcond, r1, r2, w, 4, rw, info )
1016 CALL chkxer(
'CSYSVX', infot, nout, lerr, ok )
1018 CALL csysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
1019 $ rcond, r1, r2, w, 4, rw, info )
1020 CALL chkxer(
'CSYSVX', infot, nout, lerr, ok )
1022 CALL csysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
1023 $ rcond, r1, r2, w, 3, rw, info )
1024 CALL chkxer(
'CSYSVX', infot, nout, lerr, ok )
1033 CALL csysvxx(
'/',
'U', 0, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
1034 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1035 $ err_bnds_c, nparams, params, w, rw, info )
1036 CALL chkxer(
'CSYSVXX', infot, nout, lerr, ok )
1038 CALL csysvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
1039 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1040 $ err_bnds_c, nparams, params, w, rw, info )
1041 CALL chkxer(
'CSYSVXX', infot, nout, lerr, ok )
1043 CALL csysvxx(
'N',
'U', -1, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
1044 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1045 $ err_bnds_c, nparams, params, w, rw, info )
1046 CALL chkxer(
'CSYSVXX', infot, nout, lerr, ok )
1049 CALL csysvxx(
'N',
'U', 0, -1, a, 1, af, 1, ip, eq, r, b, 1, x,
1050 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1051 $ err_bnds_c, nparams, params, w, rw, info )
1052 CALL chkxer(
'CSYSVXX', infot, nout, lerr, ok )
1055 CALL csysvxx(
'N',
'U', 2, 0, a, 1, af, 2, ip, eq, r, b, 2, x,
1056 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1057 $ err_bnds_c, nparams, params, w, rw, info )
1058 CALL chkxer(
'CSYSVXX', infot, nout, lerr, ok )
1060 CALL csysvxx(
'N',
'U', 2, 0, a, 2, af, 1, ip, eq, r, b, 2, x,
1061 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1062 $ err_bnds_c, nparams, params, w, rw, info )
1063 CALL chkxer(
'CSYSVXX', infot, nout, lerr, ok )
1066 CALL csysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 1, x,
1067 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1068 $ err_bnds_c, nparams, params, w, rw, info )
1069 CALL chkxer(
'CSYSVXX', infot, nout, lerr, ok )
1071 CALL csysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 2, x,
1072 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
1073 $ err_bnds_c, nparams, params, w, rw, info )
1074 CALL chkxer(
'CSYSVXX', infot, nout, lerr, ok )
1076 ELSE IF( lsamen( 2, c2,
'SR' ) )
THEN
1080 srnamt =
'CSYSV_ROOK'
1082 CALL csysv_rook(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
1083 CALL chkxer(
'CSYSV_ROOK', infot, nout, lerr, ok )
1085 CALL csysv_rook(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
1086 CALL chkxer(
'CSYSV_ROOK', infot, nout, lerr, ok )
1088 CALL csysv_rook(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
1089 CALL chkxer(
'CSYSV_ROOK', infot, nout, lerr, ok )
1091 CALL csysv_rook(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
1092 CALL chkxer(
'CSYSV_ROOK', infot, nout, lerr, ok )
1094 ELSE IF( lsamen( 2, c2,
'SP' ) )
THEN
1100 CALL cspsv(
'/', 0, 0, a, ip, b, 1, info )
1101 CALL chkxer(
'CSPSV ', infot, nout, lerr, ok )
1103 CALL cspsv(
'U', -1, 0, a, ip, b, 1, info )
1104 CALL chkxer(
'CSPSV ', infot, nout, lerr, ok )
1106 CALL cspsv(
'U', 0, -1, a, ip, b, 1, info )
1107 CALL chkxer(
'CSPSV ', infot, nout, lerr, ok )
1109 CALL cspsv(
'U', 2, 0, a, ip, b, 1, info )
1110 CALL chkxer(
'CSPSV ', infot, nout, lerr, ok )
1116 CALL cspsvx(
'/',
'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1118 CALL chkxer(
'CSPSVX', infot, nout, lerr, ok )
1120 CALL cspsvx(
'N',
'/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1122 CALL chkxer(
'CSPSVX', infot, nout, lerr, ok )
1124 CALL cspsvx(
'N',
'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1126 CALL chkxer(
'CSPSVX', infot, nout, lerr, ok )
1128 CALL cspsvx(
'N',
'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
1130 CALL chkxer(
'CSPSVX', infot, nout, lerr, ok )
1132 CALL cspsvx(
'N',
'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
1134 CALL chkxer(
'CSPSVX', infot, nout, lerr, ok )
1136 CALL cspsvx(
'N',
'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
1138 CALL chkxer(
'CSPSVX', infot, nout, lerr, ok )
1144 WRITE( nout, fmt = 9999 )path
1146 WRITE( nout, fmt = 9998 )path
1149 9999
FORMAT( 1x, a3,
' drivers passed the tests of the error exits' )
1150 9998
FORMAT(
' *** ', a3,
' drivers failed the tests of the error ',
subroutine chpsv(UPLO, N, NRHS, AP, IPIV, B, LDB, INFO)
CHPSV computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine cposvx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, EQUED, S, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
CPOSVX computes the solution to system of linear equations A * X = B for PO matrices ...
subroutine cgtsvx(FACT, TRANS, N, NRHS, DL, D, DU, DLF, DF, DUF, DU2, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
CGTSVX computes the solution to system of linear equations A * X = B for GT matrices ...
subroutine cptsv(N, NRHS, D, E, B, LDB, INFO)
CPTSV computes the solution to system of linear equations A * X = B for PT matrices ...
subroutine cposvxx(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)
CPOSVXX computes the solution to system of linear equations A * X = B for PO matrices ...
subroutine cspsvx(FACT, UPLO, N, NRHS, AP, AFP, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
CSPSVX computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine csysv(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
CSYSV computes the solution to system of linear equations A * X = B for SY matrices ...
subroutine chesvxx(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)
CHESVXX computes the solution to system of linear equations A * X = B for HE matrices ...
subroutine cgbsvx(FACT, TRANS, N, KL, KU, NRHS, AB, LDAB, AFB, LDAFB, IPIV, EQUED, R, C, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
CGBSVX computes the solution to system of linear equations A * X = B for GB matrices ...
subroutine chesv_rook(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
CHESV_ROOK computes the solution to a system of linear equations A * X = B for HE matrices using the ...
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine cspsv(UPLO, N, NRHS, AP, IPIV, B, LDB, INFO)
CSPSV computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine cgesvxx(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)
CGESVXX computes the solution to system of linear equations A * X = B for GE matrices ...
subroutine csysvxx(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)
CSYSVXX computes the solution to system of linear equations A * X = B for SY matrices ...
subroutine chesv(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
CHESV computes the solution to system of linear equations A * X = B for HE matrices ...
subroutine cgbsvxx(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)
CGBSVXX computes the solution to system of linear equations A * X = B for GB matrices ...
subroutine cgesvx(FACT, TRANS, N, NRHS, A, LDA, AF, LDAF, IPIV, EQUED, R, C, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
CGESVX computes the solution to system of linear equations A * X = B for GE matrices ...
subroutine csysv_rook(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
CSYSV_ROOK computes the solution to system of linear equations A * X = B for SY matrices ...
subroutine cgbsv(N, KL, KU, NRHS, AB, LDAB, IPIV, B, LDB, INFO)
CGBSV computes the solution to system of linear equations A * X = B for GB matrices (simple driver) ...
subroutine cpbsvx(FACT, UPLO, N, KD, NRHS, AB, LDAB, AFB, LDAFB, EQUED, S, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
CPBSVX computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine cppsv(UPLO, N, NRHS, AP, B, LDB, INFO)
CPPSV computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine cgtsv(N, NRHS, DL, D, DU, B, LDB, INFO)
CGTSV computes the solution to system of linear equations A * X = B for GT matrices ...
subroutine csysvx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, LWORK, RWORK, INFO)
CSYSVX computes the solution to system of linear equations A * X = B for SY matrices ...
subroutine cgesv(N, NRHS, A, LDA, IPIV, B, LDB, INFO)
CGESV computes the solution to system of linear equations A * X = B for GE matrices (simple driver) ...
subroutine chpsvx(FACT, UPLO, N, NRHS, AP, AFP, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
CHPSVX computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine cptsvx(FACT, N, NRHS, D, E, DF, EF, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
CPTSVX computes the solution to system of linear equations A * X = B for PT matrices ...
subroutine cerrvx(PATH, NUNIT)
CERRVX
subroutine cposv(UPLO, N, NRHS, A, LDA, B, LDB, INFO)
CPOSV computes the solution to system of linear equations A * X = B for PO matrices ...
subroutine chesvx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, LWORK, RWORK, INFO)
CHESVX computes the solution to system of linear equations A * X = B for HE matrices ...
subroutine cpbsv(UPLO, N, KD, NRHS, AB, LDAB, B, LDB, INFO)
CPBSV computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine cppsvx(FACT, UPLO, N, NRHS, AP, AFP, EQUED, S, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, RWORK, INFO)
CPPSVX computes the solution to system of linear equations A * X = B for OTHER matrices ...