77 INTEGER i, info,
j, n_err_bnds, nparams
78 REAL rcond, rpvgrw, berr
81 INTEGER ip( nmax ), iw( nmax )
82 REAL a( nmax, nmax ), af( nmax, nmax ),
b( nmax ),
83 $ c( nmax ), r( nmax ), r1( nmax ), r2( nmax ),
84 $ w( 2*nmax ), x( nmax ), err_bnds_n( nmax, 3 ),
85 $ err_bnds_c( nmax, 3 ), params( 1 )
104 COMMON / infoc / infot, nout, ok, lerr
105 COMMON / srnamc / srnamt
113 WRITE( nout, fmt = * )
120 a( i,
j ) = 1. /
REAL( i+
j )
121 af( i,
j ) = 1. /
REAL( i+
j )
135 IF(
lsamen( 2, c2,
'GE' ) )
THEN
141 CALL
sgesv( -1, 0, a, 1, ip,
b, 1, info )
142 CALL
chkxer(
'SGESV ', infot, nout, lerr, ok )
144 CALL
sgesv( 0, -1, a, 1, ip,
b, 1, info )
145 CALL
chkxer(
'SGESV ', infot, nout, lerr, ok )
147 CALL
sgesv( 2, 1, a, 1, ip,
b, 2, info )
148 CALL
chkxer(
'SGESV ', infot, nout, lerr, ok )
150 CALL
sgesv( 2, 1, a, 2, ip,
b, 1, info )
151 CALL
chkxer(
'SGESV ', infot, nout, lerr, ok )
157 CALL
sgesvx(
'/',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
158 $ x, 1, rcond, r1, r2, w, iw, info )
159 CALL
chkxer(
'SGESVX', infot, nout, lerr, ok )
161 CALL
sgesvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
162 $ x, 1, rcond, r1, r2, w, iw, info )
163 CALL
chkxer(
'SGESVX', infot, nout, lerr, ok )
165 CALL
sgesvx(
'N',
'N', -1, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
166 $ x, 1, rcond, r1, r2, w, iw, info )
167 CALL
chkxer(
'SGESVX', infot, nout, lerr, ok )
169 CALL
sgesvx(
'N',
'N', 0, -1, a, 1, af, 1, ip, eq, r, c,
b, 1,
170 $ x, 1, rcond, r1, r2, w, iw, info )
171 CALL
chkxer(
'SGESVX', infot, nout, lerr, ok )
173 CALL
sgesvx(
'N',
'N', 2, 1, a, 1, af, 2, ip, eq, r, c,
b, 2,
174 $ x, 2, rcond, r1, r2, w, iw, info )
175 CALL
chkxer(
'SGESVX', infot, nout, lerr, ok )
177 CALL
sgesvx(
'N',
'N', 2, 1, a, 2, af, 1, ip, eq, r, c,
b, 2,
178 $ x, 2, rcond, r1, r2, w, iw, info )
179 CALL
chkxer(
'SGESVX', infot, nout, lerr, ok )
182 CALL
sgesvx(
'F',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
183 $ x, 1, rcond, r1, r2, w, iw, info )
184 CALL
chkxer(
'SGESVX', infot, nout, lerr, ok )
187 CALL
sgesvx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
188 $ x, 1, rcond, r1, r2, w, iw, info )
189 CALL
chkxer(
'SGESVX', infot, nout, lerr, ok )
192 CALL
sgesvx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
193 $ x, 1, rcond, r1, r2, w, iw, info )
194 CALL
chkxer(
'SGESVX', infot, nout, lerr, ok )
196 CALL
sgesvx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c,
b, 1,
197 $ x, 2, rcond, r1, r2, w, iw, info )
198 CALL
chkxer(
'SGESVX', infot, nout, lerr, ok )
200 CALL
sgesvx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c,
b, 2,
201 $ x, 1, rcond, r1, r2, w, iw, info )
202 CALL
chkxer(
'SGESVX', infot, nout, lerr, ok )
210 CALL
sgesvxx(
'/',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
211 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
212 $ err_bnds_c, nparams, params, w, iw, info )
213 CALL
chkxer(
'SGESVXX', infot, nout, lerr, ok )
215 CALL
sgesvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
216 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
217 $ err_bnds_c, nparams, params, w, iw, info )
218 CALL
chkxer(
'SGESVXX', infot, nout, lerr, ok )
220 CALL
sgesvxx(
'N',
'N', -1, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
221 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
222 $ err_bnds_c, nparams, params, w, iw, info )
223 CALL
chkxer(
'SGESVXX', infot, nout, lerr, ok )
225 CALL
sgesvxx(
'N',
'N', 0, -1, a, 1, af, 1, ip, eq, r, c,
b, 1,
226 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
227 $ err_bnds_c, nparams, params, w, iw, info )
228 CALL
chkxer(
'SGESVXX', infot, nout, lerr, ok )
230 CALL
sgesvxx(
'N',
'N', 2, 1, a, 1, af, 2, ip, eq, r, c,
b, 2,
231 $ x, 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
232 $ err_bnds_c, nparams, params, w, iw, info )
233 CALL
chkxer(
'SGESVXX', infot, nout, lerr, ok )
235 CALL
sgesvxx(
'N',
'N', 2, 1, a, 2, af, 1, ip, eq, r, c,
b, 2,
236 $ x, 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
237 $ err_bnds_c, nparams, params, w, iw, info )
238 CALL
chkxer(
'SGESVXX', infot, nout, lerr, ok )
241 CALL
sgesvxx(
'F',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
242 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
243 $ err_bnds_c, nparams, params, w, iw, info )
244 CALL
chkxer(
'SGESVXX', infot, nout, lerr, ok )
247 CALL
sgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
248 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
249 $ err_bnds_c, nparams, params, w, iw, info )
250 CALL
chkxer(
'SGESVXX', infot, nout, lerr, ok )
253 CALL
sgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c,
b, 1,
254 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
255 $ err_bnds_c, nparams, params, w, iw, info )
256 CALL
chkxer(
'SGESVXX', infot, nout, lerr, ok )
258 CALL
sgesvxx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c,
b, 1,
259 $ x, 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
260 $ err_bnds_c, nparams, params, w, iw, info )
261 CALL
chkxer(
'SGESVXX', infot, nout, lerr, ok )
263 CALL
sgesvxx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c,
b, 2,
264 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
265 $ err_bnds_c, nparams, params, w, iw, info )
266 CALL
chkxer(
'SGESVXX', infot, nout, lerr, ok )
268 ELSE IF(
lsamen( 2, c2,
'GB' ) )
THEN
274 CALL
sgbsv( -1, 0, 0, 0, a, 1, ip,
b, 1, info )
275 CALL
chkxer(
'SGBSV ', infot, nout, lerr, ok )
277 CALL
sgbsv( 1, -1, 0, 0, a, 1, ip,
b, 1, info )
278 CALL
chkxer(
'SGBSV ', infot, nout, lerr, ok )
280 CALL
sgbsv( 1, 0, -1, 0, a, 1, ip,
b, 1, info )
281 CALL
chkxer(
'SGBSV ', infot, nout, lerr, ok )
283 CALL
sgbsv( 0, 0, 0, -1, a, 1, ip,
b, 1, info )
284 CALL
chkxer(
'SGBSV ', infot, nout, lerr, ok )
286 CALL
sgbsv( 1, 1, 1, 0, a, 3, ip,
b, 1, info )
287 CALL
chkxer(
'SGBSV ', infot, nout, lerr, ok )
289 CALL
sgbsv( 2, 0, 0, 0, a, 1, ip,
b, 1, info )
290 CALL
chkxer(
'SGBSV ', infot, nout, lerr, ok )
296 CALL
sgbsvx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
297 $
b, 1, x, 1, rcond, r1, r2, w, iw, info )
298 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
300 CALL
sgbsvx(
'N',
'/', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
301 $
b, 1, x, 1, rcond, r1, r2, w, iw, info )
302 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
304 CALL
sgbsvx(
'N',
'N', -1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
305 $
b, 1, x, 1, rcond, r1, r2, w, iw, info )
306 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
308 CALL
sgbsvx(
'N',
'N', 1, -1, 0, 0, a, 1, af, 1, ip, eq, r, c,
309 $
b, 1, x, 1, rcond, r1, r2, w, iw, info )
310 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
312 CALL
sgbsvx(
'N',
'N', 1, 0, -1, 0, a, 1, af, 1, ip, eq, r, c,
313 $
b, 1, x, 1, rcond, r1, r2, w, iw, info )
314 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
316 CALL
sgbsvx(
'N',
'N', 0, 0, 0, -1, a, 1, af, 1, ip, eq, r, c,
317 $
b, 1, x, 1, rcond, r1, r2, w, iw, info )
318 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
320 CALL
sgbsvx(
'N',
'N', 1, 1, 1, 0, a, 2, af, 4, ip, eq, r, c,
321 $
b, 1, x, 1, rcond, r1, r2, w, iw, info )
322 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
324 CALL
sgbsvx(
'N',
'N', 1, 1, 1, 0, a, 3, af, 3, ip, eq, r, c,
325 $
b, 1, x, 1, rcond, r1, r2, w, iw, info )
326 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
329 CALL
sgbsvx(
'F',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
330 $
b, 1, x, 1, rcond, r1, r2, w, iw, info )
331 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
334 CALL
sgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
335 $
b, 1, x, 1, rcond, r1, r2, w, iw, info )
336 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
339 CALL
sgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
340 $
b, 1, x, 1, rcond, r1, r2, w, iw, info )
341 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
343 CALL
sgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
344 $
b, 1, x, 2, rcond, r1, r2, w, iw, info )
345 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
347 CALL
sgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
348 $
b, 2, x, 1, rcond, r1, r2, w, iw, info )
349 CALL
chkxer(
'SGBSVX', infot, nout, lerr, ok )
357 CALL
sgbsvxx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
358 $
b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
359 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
361 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
363 CALL
sgbsvxx(
'N',
'/', 0, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
364 $
b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
365 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
367 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
369 CALL
sgbsvxx(
'N',
'N', -1, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
370 $
b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
371 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
373 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
375 CALL
sgbsvxx(
'N',
'N', 2, -1, 1, 0, a, 1, af, 1, ip, eq,
376 $ r, c,
b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
377 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
379 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
381 CALL
sgbsvxx(
'N',
'N', 2, 1, -1, 0, a, 1, af, 1, ip, eq,
382 $ r, c,
b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
383 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
385 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
387 CALL
sgbsvxx(
'N',
'N', 0, 1, 1, -1, a, 1, af, 1, ip, eq, r, c,
388 $
b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
389 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
391 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
393 CALL
sgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 2, af, 2, ip, eq, r, c,
394 $
b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
395 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
397 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
399 CALL
sgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 3, ip, eq, r, c,
400 $
b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
401 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
403 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
406 CALL
sgbsvxx(
'F',
'N', 0, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
407 $
b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
408 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
410 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
413 CALL
sgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
414 $
b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
415 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
417 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
420 CALL
sgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
421 $
b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
422 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
424 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
426 CALL
sgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
427 $
b, 1, x, 2, rcond, rpvgrw, berr, n_err_bnds,
428 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
430 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
432 CALL
sgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
433 $
b, 2, x, 1, rcond, rpvgrw, berr, n_err_bnds,
434 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
436 CALL
chkxer(
'SGBSVXX', infot, nout, lerr, ok )
438 ELSE IF(
lsamen( 2, c2,
'GT' ) )
THEN
444 CALL
sgtsv( -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
b, 1,
446 CALL
chkxer(
'SGTSV ', infot, nout, lerr, ok )
448 CALL
sgtsv( 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
b, 1,
450 CALL
chkxer(
'SGTSV ', infot, nout, lerr, ok )
452 CALL
sgtsv( 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
b, 1, info )
453 CALL
chkxer(
'SGTSV ', infot, nout, lerr, ok )
459 CALL
sgtsvx(
'/',
'N', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
460 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
461 $ ip,
b, 1, x, 1, rcond, r1, r2, w, iw, info )
462 CALL
chkxer(
'SGTSVX', infot, nout, lerr, ok )
464 CALL
sgtsvx(
'N',
'/', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
465 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
466 $ ip,
b, 1, x, 1, rcond, r1, r2, w, iw, info )
467 CALL
chkxer(
'SGTSVX', infot, nout, lerr, ok )
469 CALL
sgtsvx(
'N',
'N', -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
470 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
471 $ ip,
b, 1, x, 1, rcond, r1, r2, w, iw, info )
472 CALL
chkxer(
'SGTSVX', infot, nout, lerr, ok )
474 CALL
sgtsvx(
'N',
'N', 0, -1, 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, iw, info )
477 CALL
chkxer(
'SGTSVX', infot, nout, lerr, ok )
479 CALL
sgtsvx(
'N',
'N', 2, 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, 2, rcond, r1, r2, w, iw, info )
482 CALL
chkxer(
'SGTSVX', infot, nout, lerr, ok )
484 CALL
sgtsvx(
'N',
'N', 2, 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, 2, x, 1, rcond, r1, r2, w, iw, info )
487 CALL
chkxer(
'SGTSVX', infot, nout, lerr, ok )
489 ELSE IF(
lsamen( 2, c2,
'PO' ) )
THEN
495 CALL
sposv(
'/', 0, 0, a, 1,
b, 1, info )
496 CALL
chkxer(
'SPOSV ', infot, nout, lerr, ok )
498 CALL
sposv(
'U', -1, 0, a, 1,
b, 1, info )
499 CALL
chkxer(
'SPOSV ', infot, nout, lerr, ok )
501 CALL
sposv(
'U', 0, -1, a, 1,
b, 1, info )
502 CALL
chkxer(
'SPOSV ', infot, nout, lerr, ok )
504 CALL
sposv(
'U', 2, 0, a, 1,
b, 2, info )
505 CALL
chkxer(
'SPOSV ', infot, nout, lerr, ok )
507 CALL
sposv(
'U', 2, 0, a, 2,
b, 1, info )
508 CALL
chkxer(
'SPOSV ', infot, nout, lerr, ok )
514 CALL
sposvx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
515 $ rcond, r1, r2, w, iw, info )
516 CALL
chkxer(
'SPOSVX', infot, nout, lerr, ok )
518 CALL
sposvx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
519 $ rcond, r1, r2, w, iw, info )
520 CALL
chkxer(
'SPOSVX', infot, nout, lerr, ok )
522 CALL
sposvx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
523 $ rcond, r1, r2, w, iw, info )
524 CALL
chkxer(
'SPOSVX', infot, nout, lerr, ok )
526 CALL
sposvx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c,
b, 1, x, 1,
527 $ rcond, r1, r2, w, iw, info )
528 CALL
chkxer(
'SPOSVX', infot, nout, lerr, ok )
530 CALL
sposvx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c,
b, 2, x, 2,
531 $ rcond, r1, r2, w, iw, info )
532 CALL
chkxer(
'SPOSVX', infot, nout, lerr, ok )
534 CALL
sposvx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c,
b, 2, x, 2,
535 $ rcond, r1, r2, w, iw, info )
536 CALL
chkxer(
'SPOSVX', infot, nout, lerr, ok )
539 CALL
sposvx(
'F',
'U', 0, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
540 $ rcond, r1, r2, w, iw, info )
541 CALL
chkxer(
'SPOSVX', infot, nout, lerr, ok )
544 CALL
sposvx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
545 $ rcond, r1, r2, w, iw, info )
546 CALL
chkxer(
'SPOSVX', infot, nout, lerr, ok )
548 CALL
sposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c,
b, 1, x, 2,
549 $ rcond, r1, r2, w, iw, info )
550 CALL
chkxer(
'SPOSVX', infot, nout, lerr, ok )
552 CALL
sposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c,
b, 2, x, 1,
553 $ rcond, r1, r2, w, iw, info )
554 CALL
chkxer(
'SPOSVX', infot, nout, lerr, ok )
562 CALL
sposvxx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
563 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
564 $ err_bnds_c, nparams, params, w, iw, info )
565 CALL
chkxer(
'SPOSVXX', infot, nout, lerr, ok )
567 CALL
sposvxx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
568 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
569 $ err_bnds_c, nparams, params, w, iw, info )
570 CALL
chkxer(
'SPOSVXX', infot, nout, lerr, ok )
572 CALL
sposvxx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
573 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
574 $ err_bnds_c, nparams, params, w, iw, info )
575 CALL
chkxer(
'SPOSVXX', infot, nout, lerr, ok )
577 CALL
sposvxx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c,
b, 1, x, 1,
578 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
579 $ err_bnds_c, nparams, params, w, iw, info )
580 CALL
chkxer(
'SPOSVXX', infot, nout, lerr, ok )
582 CALL
sposvxx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c,
b, 2, x, 2,
583 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
584 $ err_bnds_c, nparams, params, w, iw, info )
585 CALL
chkxer(
'SPOSVXX', infot, nout, lerr, ok )
587 CALL
sposvxx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c,
b, 2, x, 2,
588 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
589 $ err_bnds_c, nparams, params, w, iw, info )
590 CALL
chkxer(
'SPOSVXX', infot, nout, lerr, ok )
593 CALL
sposvxx(
'F',
'U', 0, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
594 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
595 $ err_bnds_c, nparams, params, w, iw, info )
596 CALL
chkxer(
'SPOSVXX', infot, nout, lerr, ok )
599 CALL
sposvxx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
600 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
601 $ err_bnds_c, nparams, params, w, iw, info )
602 CALL
chkxer(
'SPOSVXX', infot, nout, lerr, ok )
604 CALL
sposvxx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c,
b, 1, x, 2,
605 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
606 $ err_bnds_c, nparams, params, w, iw, info )
607 CALL
chkxer(
'SPOSVXX', infot, nout, lerr, ok )
609 CALL
sposvxx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c,
b, 2, x, 1,
610 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
611 $ err_bnds_c, nparams, params, w, iw, info )
612 CALL
chkxer(
'SPOSVXX', infot, nout, lerr, ok )
614 ELSE IF(
lsamen( 2, c2,
'PP' ) )
THEN
620 CALL
sppsv(
'/', 0, 0, a,
b, 1, info )
621 CALL
chkxer(
'SPPSV ', infot, nout, lerr, ok )
623 CALL
sppsv(
'U', -1, 0, a,
b, 1, info )
624 CALL
chkxer(
'SPPSV ', infot, nout, lerr, ok )
626 CALL
sppsv(
'U', 0, -1, a,
b, 1, info )
627 CALL
chkxer(
'SPPSV ', infot, nout, lerr, ok )
629 CALL
sppsv(
'U', 2, 0, a,
b, 1, info )
630 CALL
chkxer(
'SPPSV ', infot, nout, lerr, ok )
636 CALL
sppsvx(
'/',
'U', 0, 0, a, af, eq, c,
b, 1, x, 1, rcond,
637 $ r1, r2, w, iw, info )
638 CALL
chkxer(
'SPPSVX', infot, nout, lerr, ok )
640 CALL
sppsvx(
'N',
'/', 0, 0, a, af, eq, c,
b, 1, x, 1, rcond,
641 $ r1, r2, w, iw, info )
642 CALL
chkxer(
'SPPSVX', infot, nout, lerr, ok )
644 CALL
sppsvx(
'N',
'U', -1, 0, a, af, eq, c,
b, 1, x, 1, rcond,
645 $ r1, r2, w, iw, info )
646 CALL
chkxer(
'SPPSVX', infot, nout, lerr, ok )
648 CALL
sppsvx(
'N',
'U', 0, -1, a, af, eq, c,
b, 1, x, 1, rcond,
649 $ r1, r2, w, iw, info )
650 CALL
chkxer(
'SPPSVX', infot, nout, lerr, ok )
653 CALL
sppsvx(
'F',
'U', 0, 0, a, af, eq, c,
b, 1, x, 1, rcond,
654 $ r1, r2, w, iw, info )
655 CALL
chkxer(
'SPPSVX', infot, nout, lerr, ok )
658 CALL
sppsvx(
'F',
'U', 1, 0, a, af, eq, c,
b, 1, x, 1, rcond,
659 $ r1, r2, w, iw, info )
660 CALL
chkxer(
'SPPSVX', infot, nout, lerr, ok )
662 CALL
sppsvx(
'N',
'U', 2, 0, a, af, eq, c,
b, 1, x, 2, rcond,
663 $ r1, r2, w, iw, info )
664 CALL
chkxer(
'SPPSVX', infot, nout, lerr, ok )
666 CALL
sppsvx(
'N',
'U', 2, 0, a, af, eq, c,
b, 2, x, 1, rcond,
667 $ r1, r2, w, iw, info )
668 CALL
chkxer(
'SPPSVX', infot, nout, lerr, ok )
670 ELSE IF(
lsamen( 2, c2,
'PB' ) )
THEN
676 CALL
spbsv(
'/', 0, 0, 0, a, 1,
b, 1, info )
677 CALL
chkxer(
'SPBSV ', infot, nout, lerr, ok )
679 CALL
spbsv(
'U', -1, 0, 0, a, 1,
b, 1, info )
680 CALL
chkxer(
'SPBSV ', infot, nout, lerr, ok )
682 CALL
spbsv(
'U', 1, -1, 0, a, 1,
b, 1, info )
683 CALL
chkxer(
'SPBSV ', infot, nout, lerr, ok )
685 CALL
spbsv(
'U', 0, 0, -1, a, 1,
b, 1, info )
686 CALL
chkxer(
'SPBSV ', infot, nout, lerr, ok )
688 CALL
spbsv(
'U', 1, 1, 0, a, 1,
b, 2, info )
689 CALL
chkxer(
'SPBSV ', infot, nout, lerr, ok )
691 CALL
spbsv(
'U', 2, 0, 0, a, 1,
b, 1, info )
692 CALL
chkxer(
'SPBSV ', infot, nout, lerr, ok )
698 CALL
spbsvx(
'/',
'U', 0, 0, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
699 $ rcond, r1, r2, w, iw, info )
700 CALL
chkxer(
'SPBSVX', infot, nout, lerr, ok )
702 CALL
spbsvx(
'N',
'/', 0, 0, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
703 $ rcond, r1, r2, w, iw, info )
704 CALL
chkxer(
'SPBSVX', infot, nout, lerr, ok )
706 CALL
spbsvx(
'N',
'U', -1, 0, 0, a, 1, af, 1, eq, c,
b, 1, x,
707 $ 1, rcond, r1, r2, w, iw, info )
708 CALL
chkxer(
'SPBSVX', infot, nout, lerr, ok )
710 CALL
spbsvx(
'N',
'U', 1, -1, 0, a, 1, af, 1, eq, c,
b, 1, x,
711 $ 1, rcond, r1, r2, w, iw, info )
712 CALL
chkxer(
'SPBSVX', infot, nout, lerr, ok )
714 CALL
spbsvx(
'N',
'U', 0, 0, -1, a, 1, af, 1, eq, c,
b, 1, x,
715 $ 1, rcond, r1, r2, w, iw, info )
716 CALL
chkxer(
'SPBSVX', infot, nout, lerr, ok )
718 CALL
spbsvx(
'N',
'U', 1, 1, 0, a, 1, af, 2, eq, c,
b, 2, x, 2,
719 $ rcond, r1, r2, w, iw, info )
720 CALL
chkxer(
'SPBSVX', infot, nout, lerr, ok )
722 CALL
spbsvx(
'N',
'U', 1, 1, 0, a, 2, af, 1, eq, c,
b, 2, x, 2,
723 $ rcond, r1, r2, w, iw, info )
724 CALL
chkxer(
'SPBSVX', infot, nout, lerr, ok )
727 CALL
spbsvx(
'F',
'U', 0, 0, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
728 $ rcond, r1, r2, w, iw, info )
729 CALL
chkxer(
'SPBSVX', infot, nout, lerr, ok )
732 CALL
spbsvx(
'F',
'U', 1, 0, 0, a, 1, af, 1, eq, c,
b, 1, x, 1,
733 $ rcond, r1, r2, w, iw, info )
734 CALL
chkxer(
'SPBSVX', infot, nout, lerr, ok )
736 CALL
spbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c,
b, 1, x, 2,
737 $ rcond, r1, r2, w, iw, info )
738 CALL
chkxer(
'SPBSVX', infot, nout, lerr, ok )
740 CALL
spbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c,
b, 2, x, 1,
741 $ rcond, r1, r2, w, iw, info )
742 CALL
chkxer(
'SPBSVX', infot, nout, lerr, ok )
744 ELSE IF(
lsamen( 2, c2,
'PT' ) )
THEN
750 CALL
sptsv( -1, 0, a( 1, 1 ), a( 1, 2 ),
b, 1, info )
751 CALL
chkxer(
'SPTSV ', infot, nout, lerr, ok )
753 CALL
sptsv( 0, -1, a( 1, 1 ), a( 1, 2 ),
b, 1, info )
754 CALL
chkxer(
'SPTSV ', infot, nout, lerr, ok )
756 CALL
sptsv( 2, 0, a( 1, 1 ), a( 1, 2 ),
b, 1, info )
757 CALL
chkxer(
'SPTSV ', infot, nout, lerr, ok )
763 CALL
sptsvx(
'/', 0, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
764 $ af( 1, 2 ),
b, 1, x, 1, rcond, r1, r2, w, info )
765 CALL
chkxer(
'SPTSVX', infot, nout, lerr, ok )
767 CALL
sptsvx(
'N', -1, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
768 $ af( 1, 2 ),
b, 1, x, 1, rcond, r1, r2, w, info )
769 CALL
chkxer(
'SPTSVX', infot, nout, lerr, ok )
771 CALL
sptsvx(
'N', 0, -1, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
772 $ af( 1, 2 ),
b, 1, x, 1, rcond, r1, r2, w, info )
773 CALL
chkxer(
'SPTSVX', infot, nout, lerr, ok )
775 CALL
sptsvx(
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
776 $ af( 1, 2 ),
b, 1, x, 2, rcond, r1, r2, w, info )
777 CALL
chkxer(
'SPTSVX', infot, nout, lerr, ok )
779 CALL
sptsvx(
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
780 $ af( 1, 2 ),
b, 2, x, 1, rcond, r1, r2, w, info )
781 CALL
chkxer(
'SPTSVX', infot, nout, lerr, ok )
783 ELSE IF(
lsamen( 2, c2,
'SY' ) )
THEN
789 CALL
ssysv(
'/', 0, 0, a, 1, ip,
b, 1, w, 1, info )
790 CALL
chkxer(
'SSYSV ', infot, nout, lerr, ok )
792 CALL
ssysv(
'U', -1, 0, a, 1, ip,
b, 1, w, 1, info )
793 CALL
chkxer(
'SSYSV ', infot, nout, lerr, ok )
795 CALL
ssysv(
'U', 0, -1, a, 1, ip,
b, 1, w, 1, info )
796 CALL
chkxer(
'SSYSV ', infot, nout, lerr, ok )
798 CALL
ssysv(
'U', 2, 0, a, 2, ip,
b, 1, w, 1, info )
799 CALL
chkxer(
'SSYSV ', infot, nout, lerr, ok )
805 CALL
ssysvx(
'/',
'U', 0, 0, a, 1, af, 1, ip,
b, 1, x, 1,
806 $ rcond, r1, r2, w, 1, iw, info )
807 CALL
chkxer(
'SSYSVX', infot, nout, lerr, ok )
809 CALL
ssysvx(
'N',
'/', 0, 0, a, 1, af, 1, ip,
b, 1, x, 1,
810 $ rcond, r1, r2, w, 1, iw, info )
811 CALL
chkxer(
'SSYSVX', infot, nout, lerr, ok )
813 CALL
ssysvx(
'N',
'U', -1, 0, a, 1, af, 1, ip,
b, 1, x, 1,
814 $ rcond, r1, r2, w, 1, iw, info )
815 CALL
chkxer(
'SSYSVX', infot, nout, lerr, ok )
817 CALL
ssysvx(
'N',
'U', 0, -1, a, 1, af, 1, ip,
b, 1, x, 1,
818 $ rcond, r1, r2, w, 1, iw, info )
819 CALL
chkxer(
'SSYSVX', infot, nout, lerr, ok )
821 CALL
ssysvx(
'N',
'U', 2, 0, a, 1, af, 2, ip,
b, 2, x, 2,
822 $ rcond, r1, r2, w, 4, iw, info )
823 CALL
chkxer(
'SSYSVX', infot, nout, lerr, ok )
825 CALL
ssysvx(
'N',
'U', 2, 0, a, 2, af, 1, ip,
b, 2, x, 2,
826 $ rcond, r1, r2, w, 4, iw, info )
827 CALL
chkxer(
'SSYSVX', infot, nout, lerr, ok )
829 CALL
ssysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip,
b, 1, x, 2,
830 $ rcond, r1, r2, w, 4, iw, info )
831 CALL
chkxer(
'SSYSVX', infot, nout, lerr, ok )
833 CALL
ssysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip,
b, 2, x, 1,
834 $ rcond, r1, r2, w, 4, iw, info )
835 CALL
chkxer(
'SSYSVX', infot, nout, lerr, ok )
837 CALL
ssysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip,
b, 2, x, 2,
838 $ rcond, r1, r2, w, 3, iw, info )
839 CALL
chkxer(
'SSYSVX', infot, nout, lerr, ok )
848 CALL
ssysvxx(
'/',
'U', 0, 0, a, 1, af, 1, ip, eq, r,
b, 1, x,
849 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
850 $ err_bnds_c, nparams, params, w, iw, info )
851 CALL
chkxer(
'SSYSVXX', infot, nout, lerr, ok )
853 CALL
ssysvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r,
b, 1, x,
854 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
855 $ err_bnds_c, nparams, params, w, iw, info )
856 CALL
chkxer(
'SSYSVXX', infot, nout, lerr, ok )
858 CALL
ssysvxx(
'N',
'U', -1, 0, a, 1, af, 1, ip, eq, r,
b, 1, x,
859 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
860 $ err_bnds_c, nparams, params, w, iw, info )
861 CALL
chkxer(
'SSYSVXX', infot, nout, lerr, ok )
864 CALL
ssysvxx(
'N',
'U', 0, -1, a, 1, af, 1, ip, eq, r,
b, 1, x,
865 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
866 $ err_bnds_c, nparams, params, w, iw, info )
867 CALL
chkxer(
'SSYSVXX', infot, nout, lerr, ok )
870 CALL
ssysvxx(
'N',
'U', 2, 0, a, 1, af, 2, ip, eq, r,
b, 2, x,
871 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
872 $ err_bnds_c, nparams, params, w, iw, info )
873 CALL
chkxer(
'SSYSVXX', infot, nout, lerr, ok )
875 CALL
ssysvxx(
'N',
'U', 2, 0, a, 2, af, 1, ip, eq, r,
b, 2, x,
876 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
877 $ err_bnds_c, nparams, params, w, iw, info )
878 CALL
chkxer(
'SSYSVXX', infot, nout, lerr, ok )
881 CALL
ssysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r,
b, 1, x,
882 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
883 $ err_bnds_c, nparams, params, w, iw, info )
884 CALL
chkxer(
'SSYSVXX', infot, nout, lerr, ok )
886 CALL
ssysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r,
b, 2, x,
887 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
888 $ err_bnds_c, nparams, params, w, iw, info )
889 CALL
chkxer(
'SSYSVXX', infot, nout, lerr, ok )
893 srnamt =
'SSYSV_ROOK'
895 CALL
ssysv_rook(
'/', 0, 0, a, 1, ip,
b, 1, w, 1, info )
896 CALL
chkxer(
'SSYSV_ROOK', infot, nout, lerr, ok )
898 CALL
ssysv_rook(
'U', -1, 0, a, 1, ip,
b, 1, w, 1, info )
899 CALL
chkxer(
'SSYSV_ROOK', infot, nout, lerr, ok )
901 CALL
ssysv_rook(
'U', 0, -1, a, 1, ip,
b, 1, w, 1, info )
902 CALL
chkxer(
'SSYSV_ROOK', infot, nout, lerr, ok )
904 CALL
ssysv_rook(
'U', 2, 0, a, 2, ip,
b, 1, w, 1, info )
905 CALL
chkxer(
'SSYSV_ROOK', infot, nout, lerr, ok )
907 ELSE IF(
lsamen( 2, c2,
'SP' ) )
THEN
913 CALL
sspsv(
'/', 0, 0, a, ip,
b, 1, info )
914 CALL
chkxer(
'SSPSV ', infot, nout, lerr, ok )
916 CALL
sspsv(
'U', -1, 0, a, ip,
b, 1, info )
917 CALL
chkxer(
'SSPSV ', infot, nout, lerr, ok )
919 CALL
sspsv(
'U', 0, -1, a, ip,
b, 1, info )
920 CALL
chkxer(
'SSPSV ', infot, nout, lerr, ok )
922 CALL
sspsv(
'U', 2, 0, a, ip,
b, 1, info )
923 CALL
chkxer(
'SSPSV ', infot, nout, lerr, ok )
929 CALL
sspsvx(
'/',
'U', 0, 0, a, af, ip,
b, 1, x, 1, rcond, r1,
931 CALL
chkxer(
'SSPSVX', infot, nout, lerr, ok )
933 CALL
sspsvx(
'N',
'/', 0, 0, a, af, ip,
b, 1, x, 1, rcond, r1,
935 CALL
chkxer(
'SSPSVX', infot, nout, lerr, ok )
937 CALL
sspsvx(
'N',
'U', -1, 0, a, af, ip,
b, 1, x, 1, rcond, r1,
939 CALL
chkxer(
'SSPSVX', infot, nout, lerr, ok )
941 CALL
sspsvx(
'N',
'U', 0, -1, a, af, ip,
b, 1, x, 1, rcond, r1,
943 CALL
chkxer(
'SSPSVX', infot, nout, lerr, ok )
945 CALL
sspsvx(
'N',
'U', 2, 0, a, af, ip,
b, 1, x, 2, rcond, r1,
947 CALL
chkxer(
'SSPSVX', infot, nout, lerr, ok )
949 CALL
sspsvx(
'N',
'U', 2, 0, a, af, ip,
b, 2, x, 1, rcond, r1,
951 CALL
chkxer(
'SSPSVX', infot, nout, lerr, ok )
957 WRITE( nout, fmt = 9999 )path
959 WRITE( nout, fmt = 9998 )path
962 9999
FORMAT( 1x, a3,
' drivers passed the tests of the error exits' )
963 9998
FORMAT(
' *** ', a3,
' drivers failed the tests of the error ',