81 INTEGER IP( NMAX ), IW( NMAX )
82 REAL A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
83 $ C( NMAX ), E( NMAX ), R( NMAX ), R1( NMAX ),
84 $ R2( NMAX ), W( 2*NMAX ), X( NMAX )
103 COMMON / infoc / infot, nout, ok, lerr
104 COMMON / srnamc / srnamt
112 WRITE( nout, fmt = * )
119 a( i, j ) = 1. / real( i+j )
120 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 )
204 ELSE IF(
lsamen( 2, c2,
'GB' ) )
THEN
210 CALL sgbsv( -1, 0, 0, 0, a, 1, ip, b, 1, info )
211 CALL chkxer(
'SGBSV ', infot, nout, lerr, ok )
213 CALL sgbsv( 1, -1, 0, 0, a, 1, ip, b, 1, info )
214 CALL chkxer(
'SGBSV ', infot, nout, lerr, ok )
216 CALL sgbsv( 1, 0, -1, 0, a, 1, ip, b, 1, info )
217 CALL chkxer(
'SGBSV ', infot, nout, lerr, ok )
219 CALL sgbsv( 0, 0, 0, -1, a, 1, ip, b, 1, info )
220 CALL chkxer(
'SGBSV ', infot, nout, lerr, ok )
222 CALL sgbsv( 1, 1, 1, 0, a, 3, ip, b, 1, info )
223 CALL chkxer(
'SGBSV ', infot, nout, lerr, ok )
225 CALL sgbsv( 2, 0, 0, 0, a, 1, ip, b, 1, info )
226 CALL chkxer(
'SGBSV ', infot, nout, lerr, ok )
232 CALL sgbsvx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
233 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
234 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
236 CALL sgbsvx(
'N',
'/', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
237 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
238 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
240 CALL sgbsvx(
'N',
'N', -1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
241 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
242 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
244 CALL sgbsvx(
'N',
'N', 1, -1, 0, 0, a, 1, af, 1, ip, eq, r, c,
245 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
246 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
248 CALL sgbsvx(
'N',
'N', 1, 0, -1, 0, a, 1, af, 1, ip, eq, r, c,
249 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
250 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
252 CALL sgbsvx(
'N',
'N', 0, 0, 0, -1, a, 1, af, 1, ip, eq, r, c,
253 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
254 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
256 CALL sgbsvx(
'N',
'N', 1, 1, 1, 0, a, 2, af, 4, ip, eq, r, c,
257 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
258 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
260 CALL sgbsvx(
'N',
'N', 1, 1, 1, 0, a, 3, af, 3, ip, eq, r, c,
261 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
262 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
265 CALL sgbsvx(
'F',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
266 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
267 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
270 CALL sgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
271 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
272 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
275 CALL sgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
276 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
277 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
279 CALL sgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
280 $ b, 1, x, 2, rcond, r1, r2, w, iw, info )
281 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
283 CALL sgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
284 $ b, 2, x, 1, rcond, r1, r2, w, iw, info )
285 CALL chkxer(
'SGBSVX', infot, nout, lerr, ok )
287 ELSE IF(
lsamen( 2, c2,
'GT' ) )
THEN
293 CALL sgtsv( -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
295 CALL chkxer(
'SGTSV ', infot, nout, lerr, ok )
297 CALL sgtsv( 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
299 CALL chkxer(
'SGTSV ', infot, nout, lerr, ok )
301 CALL sgtsv( 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1, info )
302 CALL chkxer(
'SGTSV ', infot, nout, lerr, ok )
308 CALL sgtsvx(
'/',
'N', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
309 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
310 $ ip, b, 1, x, 1, rcond, r1, r2, w, iw, info )
311 CALL chkxer(
'SGTSVX', infot, nout, lerr, ok )
313 CALL sgtsvx(
'N',
'/', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
314 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
315 $ ip, b, 1, x, 1, rcond, r1, r2, w, iw, info )
316 CALL chkxer(
'SGTSVX', infot, nout, lerr, ok )
318 CALL sgtsvx(
'N',
'N', -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
319 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
320 $ ip, b, 1, x, 1, rcond, r1, r2, w, iw, info )
321 CALL chkxer(
'SGTSVX', infot, nout, lerr, ok )
323 CALL sgtsvx(
'N',
'N', 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
324 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
325 $ ip, b, 1, x, 1, rcond, r1, r2, w, iw, info )
326 CALL chkxer(
'SGTSVX', infot, nout, lerr, ok )
328 CALL sgtsvx(
'N',
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
329 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
330 $ ip, b, 1, x, 2, rcond, r1, r2, w, iw, info )
331 CALL chkxer(
'SGTSVX', infot, nout, lerr, ok )
333 CALL sgtsvx(
'N',
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
334 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
335 $ ip, b, 2, x, 1, rcond, r1, r2, w, iw, info )
336 CALL chkxer(
'SGTSVX', infot, nout, lerr, ok )
338 ELSE IF(
lsamen( 2, c2,
'PO' ) )
THEN
344 CALL sposv(
'/', 0, 0, a, 1, b, 1, info )
345 CALL chkxer(
'SPOSV ', infot, nout, lerr, ok )
347 CALL sposv(
'U', -1, 0, a, 1, b, 1, info )
348 CALL chkxer(
'SPOSV ', infot, nout, lerr, ok )
350 CALL sposv(
'U', 0, -1, a, 1, b, 1, info )
351 CALL chkxer(
'SPOSV ', infot, nout, lerr, ok )
353 CALL sposv(
'U', 2, 0, a, 1, b, 2, info )
354 CALL chkxer(
'SPOSV ', infot, nout, lerr, ok )
356 CALL sposv(
'U', 2, 0, a, 2, b, 1, info )
357 CALL chkxer(
'SPOSV ', infot, nout, lerr, ok )
363 CALL sposvx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
364 $ rcond, r1, r2, w, iw, info )
365 CALL chkxer(
'SPOSVX', infot, nout, lerr, ok )
367 CALL sposvx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
368 $ rcond, r1, r2, w, iw, info )
369 CALL chkxer(
'SPOSVX', infot, nout, lerr, ok )
371 CALL sposvx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
372 $ rcond, r1, r2, w, iw, info )
373 CALL chkxer(
'SPOSVX', infot, nout, lerr, ok )
375 CALL sposvx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
376 $ rcond, r1, r2, w, iw, info )
377 CALL chkxer(
'SPOSVX', infot, nout, lerr, ok )
379 CALL sposvx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
380 $ rcond, r1, r2, w, iw, info )
381 CALL chkxer(
'SPOSVX', infot, nout, lerr, ok )
383 CALL sposvx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
384 $ rcond, r1, r2, w, iw, info )
385 CALL chkxer(
'SPOSVX', infot, nout, lerr, ok )
388 CALL sposvx(
'F',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
389 $ rcond, r1, r2, w, iw, info )
390 CALL chkxer(
'SPOSVX', infot, nout, lerr, ok )
393 CALL sposvx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
394 $ rcond, r1, r2, w, iw, info )
395 CALL chkxer(
'SPOSVX', infot, nout, lerr, ok )
397 CALL sposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 1, x, 2,
398 $ rcond, r1, r2, w, iw, info )
399 CALL chkxer(
'SPOSVX', infot, nout, lerr, ok )
401 CALL sposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
402 $ rcond, r1, r2, w, iw, info )
403 CALL chkxer(
'SPOSVX', infot, nout, lerr, ok )
405 ELSE IF(
lsamen( 2, c2,
'PP' ) )
THEN
411 CALL sppsv(
'/', 0, 0, a, b, 1, info )
412 CALL chkxer(
'SPPSV ', infot, nout, lerr, ok )
414 CALL sppsv(
'U', -1, 0, a, b, 1, info )
415 CALL chkxer(
'SPPSV ', infot, nout, lerr, ok )
417 CALL sppsv(
'U', 0, -1, a, b, 1, info )
418 CALL chkxer(
'SPPSV ', infot, nout, lerr, ok )
420 CALL sppsv(
'U', 2, 0, a, b, 1, info )
421 CALL chkxer(
'SPPSV ', infot, nout, lerr, ok )
427 CALL sppsvx(
'/',
'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
428 $ r1, r2, w, iw, info )
429 CALL chkxer(
'SPPSVX', infot, nout, lerr, ok )
431 CALL sppsvx(
'N',
'/', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
432 $ r1, r2, w, iw, info )
433 CALL chkxer(
'SPPSVX', infot, nout, lerr, ok )
435 CALL sppsvx(
'N',
'U', -1, 0, a, af, eq, c, b, 1, x, 1, rcond,
436 $ r1, r2, w, iw, info )
437 CALL chkxer(
'SPPSVX', infot, nout, lerr, ok )
439 CALL sppsvx(
'N',
'U', 0, -1, a, af, eq, c, b, 1, x, 1, rcond,
440 $ r1, r2, w, iw, info )
441 CALL chkxer(
'SPPSVX', infot, nout, lerr, ok )
444 CALL sppsvx(
'F',
'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
445 $ r1, r2, w, iw, info )
446 CALL chkxer(
'SPPSVX', infot, nout, lerr, ok )
449 CALL sppsvx(
'F',
'U', 1, 0, a, af, eq, c, b, 1, x, 1, rcond,
450 $ r1, r2, w, iw, info )
451 CALL chkxer(
'SPPSVX', infot, nout, lerr, ok )
453 CALL sppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 1, x, 2, rcond,
454 $ r1, r2, w, iw, info )
455 CALL chkxer(
'SPPSVX', infot, nout, lerr, ok )
457 CALL sppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 2, x, 1, rcond,
458 $ r1, r2, w, iw, info )
459 CALL chkxer(
'SPPSVX', infot, nout, lerr, ok )
461 ELSE IF(
lsamen( 2, c2,
'PB' ) )
THEN
467 CALL spbsv(
'/', 0, 0, 0, a, 1, b, 1, info )
468 CALL chkxer(
'SPBSV ', infot, nout, lerr, ok )
470 CALL spbsv(
'U', -1, 0, 0, a, 1, b, 1, info )
471 CALL chkxer(
'SPBSV ', infot, nout, lerr, ok )
473 CALL spbsv(
'U', 1, -1, 0, a, 1, b, 1, info )
474 CALL chkxer(
'SPBSV ', infot, nout, lerr, ok )
476 CALL spbsv(
'U', 0, 0, -1, a, 1, b, 1, info )
477 CALL chkxer(
'SPBSV ', infot, nout, lerr, ok )
479 CALL spbsv(
'U', 1, 1, 0, a, 1, b, 2, info )
480 CALL chkxer(
'SPBSV ', infot, nout, lerr, ok )
482 CALL spbsv(
'U', 2, 0, 0, a, 1, b, 1, info )
483 CALL chkxer(
'SPBSV ', infot, nout, lerr, ok )
489 CALL spbsvx(
'/',
'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
490 $ rcond, r1, r2, w, iw, info )
491 CALL chkxer(
'SPBSVX', infot, nout, lerr, ok )
493 CALL spbsvx(
'N',
'/', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
494 $ rcond, r1, r2, w, iw, info )
495 CALL chkxer(
'SPBSVX', infot, nout, lerr, ok )
497 CALL spbsvx(
'N',
'U', -1, 0, 0, a, 1, af, 1, eq, c, b, 1, x,
498 $ 1, rcond, r1, r2, w, iw, info )
499 CALL chkxer(
'SPBSVX', infot, nout, lerr, ok )
501 CALL spbsvx(
'N',
'U', 1, -1, 0, a, 1, af, 1, eq, c, b, 1, x,
502 $ 1, rcond, r1, r2, w, iw, info )
503 CALL chkxer(
'SPBSVX', infot, nout, lerr, ok )
505 CALL spbsvx(
'N',
'U', 0, 0, -1, a, 1, af, 1, eq, c, b, 1, x,
506 $ 1, rcond, r1, r2, w, iw, info )
507 CALL chkxer(
'SPBSVX', infot, nout, lerr, ok )
509 CALL spbsvx(
'N',
'U', 1, 1, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
510 $ rcond, r1, r2, w, iw, info )
511 CALL chkxer(
'SPBSVX', infot, nout, lerr, ok )
513 CALL spbsvx(
'N',
'U', 1, 1, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
514 $ rcond, r1, r2, w, iw, info )
515 CALL chkxer(
'SPBSVX', infot, nout, lerr, ok )
518 CALL spbsvx(
'F',
'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
519 $ rcond, r1, r2, w, iw, info )
520 CALL chkxer(
'SPBSVX', infot, nout, lerr, ok )
523 CALL spbsvx(
'F',
'U', 1, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
524 $ rcond, r1, r2, w, iw, info )
525 CALL chkxer(
'SPBSVX', infot, nout, lerr, ok )
527 CALL spbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 2,
528 $ rcond, r1, r2, w, iw, info )
529 CALL chkxer(
'SPBSVX', infot, nout, lerr, ok )
531 CALL spbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 2, x, 1,
532 $ rcond, r1, r2, w, iw, info )
533 CALL chkxer(
'SPBSVX', infot, nout, lerr, ok )
535 ELSE IF(
lsamen( 2, c2,
'PT' ) )
THEN
541 CALL sptsv( -1, 0, a( 1, 1 ), a( 1, 2 ), b, 1, info )
542 CALL chkxer(
'SPTSV ', infot, nout, lerr, ok )
544 CALL sptsv( 0, -1, a( 1, 1 ), a( 1, 2 ), b, 1, info )
545 CALL chkxer(
'SPTSV ', infot, nout, lerr, ok )
547 CALL sptsv( 2, 0, a( 1, 1 ), a( 1, 2 ), b, 1, info )
548 CALL chkxer(
'SPTSV ', infot, nout, lerr, ok )
554 CALL sptsvx(
'/', 0, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
555 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
556 CALL chkxer(
'SPTSVX', infot, nout, lerr, ok )
558 CALL sptsvx(
'N', -1, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
559 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
560 CALL chkxer(
'SPTSVX', infot, nout, lerr, ok )
562 CALL sptsvx(
'N', 0, -1, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
563 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
564 CALL chkxer(
'SPTSVX', infot, nout, lerr, ok )
566 CALL sptsvx(
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
567 $ af( 1, 2 ), b, 1, x, 2, rcond, r1, r2, w, info )
568 CALL chkxer(
'SPTSVX', infot, nout, lerr, ok )
570 CALL sptsvx(
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
571 $ af( 1, 2 ), b, 2, x, 1, rcond, r1, r2, w, info )
572 CALL chkxer(
'SPTSVX', infot, nout, lerr, ok )
574 ELSE IF(
lsamen( 2, c2,
'SY' ) )
THEN
580 CALL ssysv(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
581 CALL chkxer(
'SSYSV ', infot, nout, lerr, ok )
583 CALL ssysv(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
584 CALL chkxer(
'SSYSV ', infot, nout, lerr, ok )
586 CALL ssysv(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
587 CALL chkxer(
'SSYSV ', infot, nout, lerr, ok )
589 CALL ssysv(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
590 CALL chkxer(
'SSYSV ', infot, nout, lerr, ok )
592 CALL ssysv(
'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
593 CALL chkxer(
'SSYSV ', infot, nout, lerr, ok )
595 CALL ssysv(
'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
596 CALL chkxer(
'SSYSV ', infot, nout, lerr, ok )
602 CALL ssysvx(
'/',
'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
603 $ rcond, r1, r2, w, 1, iw, info )
604 CALL chkxer(
'SSYSVX', infot, nout, lerr, ok )
606 CALL ssysvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
607 $ rcond, r1, r2, w, 1, iw, info )
608 CALL chkxer(
'SSYSVX', infot, nout, lerr, ok )
610 CALL ssysvx(
'N',
'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
611 $ rcond, r1, r2, w, 1, iw, info )
612 CALL chkxer(
'SSYSVX', infot, nout, lerr, ok )
614 CALL ssysvx(
'N',
'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
615 $ rcond, r1, r2, w, 1, iw, info )
616 CALL chkxer(
'SSYSVX', infot, nout, lerr, ok )
618 CALL ssysvx(
'N',
'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
619 $ rcond, r1, r2, w, 4, iw, info )
620 CALL chkxer(
'SSYSVX', infot, nout, lerr, ok )
622 CALL ssysvx(
'N',
'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
623 $ rcond, r1, r2, w, 4, iw, info )
624 CALL chkxer(
'SSYSVX', infot, nout, lerr, ok )
626 CALL ssysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
627 $ rcond, r1, r2, w, 4, iw, info )
628 CALL chkxer(
'SSYSVX', infot, nout, lerr, ok )
630 CALL ssysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
631 $ rcond, r1, r2, w, 4, iw, info )
632 CALL chkxer(
'SSYSVX', infot, nout, lerr, ok )
634 CALL ssysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
635 $ rcond, r1, r2, w, 3, iw, info )
636 CALL chkxer(
'SSYSVX', infot, nout, lerr, ok )
639 ELSE IF(
lsamen( 2, c2,
'SR' ) )
THEN
643 srnamt =
'SSYSV_ROOK'
645 CALL ssysv_rook(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
646 CALL chkxer(
'SSYSV_ROOK', infot, nout, lerr, ok )
648 CALL ssysv_rook(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
649 CALL chkxer(
'SSYSV_ROOK', infot, nout, lerr, ok )
651 CALL ssysv_rook(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
652 CALL chkxer(
'SSYSV_ROOK', infot, nout, lerr, ok )
654 CALL ssysv_rook(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
655 CALL chkxer(
'SSYSV_ROOK', infot, nout, lerr, ok )
657 CALL ssysv_rook(
'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
658 CALL chkxer(
'SSYSV_ROOK', infot, nout, lerr, ok )
660 CALL ssysv_rook(
'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
661 CALL chkxer(
'SSYSV_ROOK', infot, nout, lerr, ok )
663 ELSE IF(
lsamen( 2, c2,
'SK' ) )
THEN
677 CALL ssysv_rk(
'/', 0, 0, a, 1, e, ip, b, 1, w, 1, info )
678 CALL chkxer(
'SSYSV_RK', infot, nout, lerr, ok )
680 CALL ssysv_rk(
'U', -1, 0, a, 1, e, ip, b, 1, w, 1, info )
681 CALL chkxer(
'SSYSV_RK', infot, nout, lerr, ok )
683 CALL ssysv_rk(
'U', 0, -1, a, 1, e, ip, b, 1, w, 1, info )
684 CALL chkxer(
'SSYSV_RK', infot, nout, lerr, ok )
686 CALL ssysv_rk(
'U', 2, 0, a, 1, e, ip, b, 2, w, 1, info )
687 CALL chkxer(
'SSYSV_RK', infot, nout, lerr, ok )
689 CALL ssysv_rk(
'U', 2, 0, a, 2, e, ip, b, 1, w, 1, info )
690 CALL chkxer(
'SSYSV_RK', infot, nout, lerr, ok )
692 CALL ssysv_rk(
'U', 0, 0, a, 1, e, ip, b, 1, w, 0, info )
693 CALL chkxer(
'SSYSV_RK', infot, nout, lerr, ok )
695 CALL ssysv_rk(
'U', 0, 0, a, 1, e, ip, b, 1, w, -2, info )
696 CALL chkxer(
'SSYSV_RK', infot, nout, lerr, ok )
698 ELSE IF(
lsamen( 2, c2,
'SA' ) )
THEN
704 CALL ssysv_aa(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
705 CALL chkxer(
'SSYSV_AA', infot, nout, lerr, ok )
707 CALL ssysv_aa(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
708 CALL chkxer(
'SSYSV_AA', infot, nout, lerr, ok )
710 CALL ssysv_aa(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
711 CALL chkxer(
'SSYSV_AA', infot, nout, lerr, ok )
713 CALL ssysv_aa(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
714 CALL chkxer(
'SSYSV_AA', infot, nout, lerr, ok )
716 ELSE IF(
lsamen( 2, c2,
'S2' ) )
THEN
720 srnamt =
'SSYSV_AA_2STAGE'
722 CALL ssysv_aa_2stage(
'/', 0, 0, a, 1, a, 1, ip, ip, b, 1,
724 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
726 CALL ssysv_aa_2stage(
'U', -1, 0, a, 1, a, 1, ip, ip, b, 1,
728 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
730 CALL ssysv_aa_2stage(
'U', 0, -1, a, 1, a, 1, ip, ip, b, 1,
732 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
734 CALL ssysv_aa_2stage(
'U', 2, 1, a, 1, a, 1, ip, ip, b, 1,
736 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
738 CALL ssysv_aa_2stage(
'U', 2, 1, a, 2, a, 8, ip, ip, b, 1,
740 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
742 CALL ssysv_aa_2stage(
'U', 2, 1, a, 2, a, 1, ip, ip, b, 2,
744 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
746 ELSE IF(
lsamen( 2, c2,
'SP' ) )
THEN
752 CALL sspsv(
'/', 0, 0, a, ip, b, 1, info )
753 CALL chkxer(
'SSPSV ', infot, nout, lerr, ok )
755 CALL sspsv(
'U', -1, 0, a, ip, b, 1, info )
756 CALL chkxer(
'SSPSV ', infot, nout, lerr, ok )
758 CALL sspsv(
'U', 0, -1, a, ip, b, 1, info )
759 CALL chkxer(
'SSPSV ', infot, nout, lerr, ok )
761 CALL sspsv(
'U', 2, 0, a, ip, b, 1, info )
762 CALL chkxer(
'SSPSV ', infot, nout, lerr, ok )
768 CALL sspsvx(
'/',
'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
770 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
772 CALL sspsvx(
'N',
'/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
774 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
776 CALL sspsvx(
'N',
'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
778 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
780 CALL sspsvx(
'N',
'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
782 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
784 CALL sspsvx(
'N',
'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
786 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
788 CALL sspsvx(
'N',
'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
790 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
796 WRITE( nout, fmt = 9999 )path
798 WRITE( nout, fmt = 9998 )path
801 9999
FORMAT( 1x, a3,
' drivers passed the tests of the error exits' )
802 9998
FORMAT(
' *** ', a3,
' drivers failed the tests of the error ',