78 INTEGER IP( NMAX ), IW( NMAX )
79 REAL A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
80 $ C( NMAX ), E( NMAX ), R( NMAX ), R1( NMAX ),
81 $ R2( NMAX ), W( 2*NMAX ), X( NMAX )
100 COMMON / infoc / infot, nout, ok, lerr
101 COMMON / srnamc / srnamt
109 WRITE( nout, fmt = * )
116 a( i, j ) = 1. / real( i+j )
117 af( i, j ) = 1. / real( i+j )
132 IF( lsamen( 2, c2,
'GE' ) )
THEN
138 CALL sgesv( -1, 0, a, 1, ip, b, 1, info )
139 CALL chkxer(
'SGESV ', infot, nout, lerr, ok )
141 CALL sgesv( 0, -1, a, 1, ip, b, 1, info )
144 CALL SGESV( 2, 1, A, 1, IP, B, 2, INFO )
145 CALL CHKXER( 'sgesv ', INFOT, NOUT, LERR, OK )
147 CALL SGESV( 2, 1, A, 2, IP, B, 1, INFO )
148 CALL CHKXER( 'sgesv ', INFOT, NOUT, LERR, OK )
154 CALL SGESVX( '/
', 'n
', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
155 $ X, 1, RCOND, R1, R2, W, IW, INFO )
156 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
158 CALL SGESVX( 'n
', '/
', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
159 $ X, 1, RCOND, R1, R2, W, IW, INFO )
160 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
162 CALL SGESVX( 'n
', 'n
', -1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
163 $ X, 1, RCOND, R1, R2, W, IW, INFO )
164 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
166 CALL SGESVX( 'n
', 'n
', 0, -1, A, 1, AF, 1, IP, EQ, R, C, B, 1,
167 $ X, 1, RCOND, R1, R2, W, IW, INFO )
168 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
170 CALL SGESVX( 'n
', 'n
', 2, 1, A, 1, AF, 2, IP, EQ, R, C, B, 2,
171 $ X, 2, RCOND, R1, R2, W, IW, INFO )
172 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
174 CALL SGESVX( 'n
', 'n
', 2, 1, A, 2, AF, 1, IP, EQ, R, C, B, 2,
175 $ X, 2, RCOND, R1, R2, W, IW, INFO )
176 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
179 CALL SGESVX( 'f
', 'n
', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
180 $ X, 1, RCOND, R1, R2, W, IW, INFO )
181 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
184 CALL SGESVX( 'f
', 'n
', 1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
185 $ X, 1, RCOND, R1, R2, W, IW, INFO )
186 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
189 CALL SGESVX( 'f
', 'n
', 1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
190 $ X, 1, RCOND, R1, R2, W, IW, INFO )
191 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
193 CALL SGESVX( 'n
', 'n
', 2, 1, A, 2, AF, 2, IP, EQ, R, C, B, 1,
194 $ X, 2, RCOND, R1, R2, W, IW, INFO )
195 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
197 CALL SGESVX( 'n
', 'n
', 2, 1, A, 2, AF, 2, IP, EQ, R, C, B, 2,
198 $ X, 1, RCOND, R1, R2, W, IW, INFO )
199 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
201 ELSE IF( LSAMEN( 2, C2, 'gb
' ) ) THEN
207 CALL SGBSV( -1, 0, 0, 0, A, 1, IP, B, 1, INFO )
208 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
210 CALL SGBSV( 1, -1, 0, 0, A, 1, IP, B, 1, INFO )
211 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
213 CALL SGBSV( 1, 0, -1, 0, A, 1, IP, B, 1, INFO )
214 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
216 CALL SGBSV( 0, 0, 0, -1, A, 1, IP, B, 1, INFO )
217 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
219 CALL SGBSV( 1, 1, 1, 0, A, 3, IP, B, 1, INFO )
220 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
222 CALL SGBSV( 2, 0, 0, 0, A, 1, IP, B, 1, INFO )
223 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
229 CALL SGBSVX( '/
', 'n
', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
230 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
231 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
233 CALL SGBSVX( 'n
', '/
', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
234 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
235 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
237 CALL SGBSVX( 'n
', 'n
', -1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
238 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
239 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
241 CALL SGBSVX( 'n
', 'n
', 1, -1, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
242 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
243 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
245 CALL SGBSVX( 'n
', 'n
', 1, 0, -1, 0, A, 1, AF, 1, IP, EQ, R, C,
246 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
247 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
249 CALL SGBSVX( 'n
', 'n
', 0, 0, 0, -1, A, 1, AF, 1, IP, EQ, R, C,
250 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
251 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
253 CALL SGBSVX( 'n
', 'n
', 1, 1, 1, 0, A, 2, AF, 4, IP, EQ, R, C,
254 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
255 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
257 CALL SGBSVX( 'n
', 'n
', 1, 1, 1, 0, A, 3, AF, 3, IP, EQ, R, C,
258 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
259 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
262 CALL SGBSVX( 'f
', 'n
', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
263 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
264 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
267 CALL SGBSVX( 'f
', 'n
', 1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
268 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
269 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
272 CALL SGBSVX( 'f
', 'n
', 1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
273 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
274 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
276 CALL SGBSVX( 'n
', 'n
', 2, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
277 $ B, 1, X, 2, RCOND, R1, R2, W, IW, INFO )
278 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
280 CALL SGBSVX( 'n
', 'n
', 2, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
281 $ B, 2, X, 1, RCOND, R1, R2, W, IW, INFO )
282 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
284 ELSE IF( LSAMEN( 2, C2, 'gt
' ) ) THEN
290 CALL SGTSV( -1, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
292 CALL CHKXER( 'sgtsv ', INFOT, NOUT, LERR, OK )
294 CALL SGTSV( 0, -1, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
296 CALL CHKXER( 'sgtsv ', INFOT, NOUT, LERR, OK )
298 CALL SGTSV( 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1, INFO )
299 CALL CHKXER( 'sgtsv ', INFOT, NOUT, LERR, OK )
305 CALL SGTSVX( '/
', 'n
', 0, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
306 $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
307 $ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
308 CALL CHKXER( 'sgtsvx', INFOT, NOUT, LERR, OK )
310 CALL SGTSVX( 'n
', '/
', 0, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
311 $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
312 $ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
313 CALL CHKXER( 'sgtsvx', INFOT, NOUT, LERR, OK )
315 CALL SGTSVX( 'n
', 'n
', -1, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
316 $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
317 $ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
318 CALL CHKXER( 'sgtsvx', INFOT, NOUT, LERR, OK )
320 CALL SGTSVX( 'n
', 'n
', 0, -1, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
321 $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
322 $ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
323 CALL CHKXER( 'sgtsvx', INFOT, NOUT, LERR, OK )
325 CALL SGTSVX( 'n
', 'n
', 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
326 $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
327 $ IP, B, 1, X, 2, RCOND, R1, R2, W, IW, INFO )
328 CALL CHKXER( 'sgtsvx', INFOT, NOUT, LERR, OK )
330 CALL SGTSVX( 'n
', 'n
', 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
331 $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
332 $ IP, B, 2, X, 1, RCOND, R1, R2, W, IW, INFO )
333 CALL CHKXER( 'sgtsvx', INFOT, NOUT, LERR, OK )
335 ELSE IF( LSAMEN( 2, C2, 'po
' ) ) THEN
341 CALL SPOSV( '/
', 0, 0, A, 1, B, 1, INFO )
342 CALL CHKXER( 'sposv ', INFOT, NOUT, LERR, OK )
344 CALL SPOSV( 'u
', -1, 0, A, 1, B, 1, INFO )
345 CALL CHKXER( 'sposv ', INFOT, NOUT, LERR, OK )
347 CALL SPOSV( 'u
', 0, -1, A, 1, B, 1, INFO )
348 CALL CHKXER( 'sposv ', INFOT, NOUT, LERR, OK )
350 CALL SPOSV( 'u
', 2, 0, A, 1, B, 2, INFO )
351 CALL CHKXER( 'sposv ', INFOT, NOUT, LERR, OK )
353 CALL SPOSV( 'u
', 2, 0, A, 2, B, 1, INFO )
354 CALL CHKXER( 'sposv ', INFOT, NOUT, LERR, OK )
360 CALL SPOSVX( '/
', 'u
', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
361 $ RCOND, R1, R2, W, IW, INFO )
362 CALL CHKXER( 'sposvx', infot, nout, lerr, ok )
364 CALL sposvx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
365 $ rcond, r1, r2, w, iw, info )
368 CALL SPOSVX( 'n
', 'u
', -1, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
369 $ RCOND, R1, R2, W, IW, INFO )
370 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
372 CALL SPOSVX( 'n
', 'u
', 0, -1, A, 1, AF, 1, EQ, C, B, 1, X, 1,
373 $ RCOND, R1, R2, W, IW, INFO )
374 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
376 CALL SPOSVX( 'n
', 'u
', 2, 0, A, 1, AF, 2, EQ, C, B, 2, X, 2,
377 $ RCOND, R1, R2, W, IW, INFO )
378 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
380 CALL SPOSVX( 'n
', 'u
', 2, 0, A, 2, AF, 1, EQ, C, B, 2, X, 2,
381 $ RCOND, R1, R2, W, IW, INFO )
382 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
385 CALL SPOSVX( 'f
', 'u
', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
386 $ RCOND, R1, R2, W, IW, INFO )
387 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
390 CALL SPOSVX( 'f
', 'u
', 1, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
391 $ RCOND, R1, R2, W, IW, INFO )
392 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
394 CALL SPOSVX( 'n
', 'u
', 2, 0, A, 2, AF, 2, EQ, C, B, 1, X, 2,
395 $ RCOND, R1, R2, W, IW, INFO )
396 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
398 CALL SPOSVX( 'n
', 'u
', 2, 0, A, 2, AF, 2, EQ, C, B, 2, X, 1,
399 $ RCOND, R1, R2, W, IW, INFO )
400 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
402 ELSE IF( LSAMEN( 2, C2, 'pp
' ) ) THEN
408 CALL SPPSV( '/
', 0, 0, A, B, 1, INFO )
409 CALL CHKXER( 'sppsv ', INFOT, NOUT, LERR, OK )
411 CALL SPPSV( 'u
', -1, 0, A, B, 1, INFO )
412 CALL CHKXER( 'sppsv ', INFOT, NOUT, LERR, OK )
414 CALL SPPSV( 'u
', 0, -1, A, B, 1, INFO )
415 CALL CHKXER( 'sppsv ', INFOT, NOUT, LERR, OK )
417 CALL SPPSV( 'u
', 2, 0, A, B, 1, INFO )
418 CALL CHKXER( 'sppsv ', INFOT, NOUT, LERR, OK )
424 CALL SPPSVX( '/
', 'u
', 0, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
425 $ R1, R2, W, IW, INFO )
426 CALL CHKXER( 'sppsvx', INFOT, NOUT, LERR, OK )
428 CALL SPPSVX( 'n
', '/
', 0, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
429 $ R1, R2, W, IW, INFO )
430 CALL CHKXER( 'sppsvx', INFOT, NOUT, LERR, OK )
432 CALL SPPSVX( 'n
', 'u
', -1, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
433 $ R1, R2, W, IW, INFO )
434 CALL CHKXER( 'sppsvx', INFOT, NOUT, LERR, OK )
436 CALL SPPSVX( 'n
', 'u
', 0, -1, A, AF, EQ, C, B, 1, X, 1, RCOND,
437 $ R1, R2, W, IW, INFO )
438 CALL CHKXER( 'sppsvx', INFOT, NOUT, LERR, OK )
441 CALL SPPSVX( 'f
', 'u
', 0, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
442 $ R1, R2, W, IW, INFO )
443 CALL CHKXER( 'sppsvx', INFOT, NOUT, LERR, OK )
446 CALL SPPSVX( 'f
', 'u
', 1, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
447 $ R1, R2, W, IW, INFO )
448 CALL CHKXER( 'sppsvx', INFOT, NOUT, LERR, OK )
450 CALL SPPSVX( 'n
', 'u
', 2, 0, A, AF, EQ, C, B, 1, X, 2, RCOND,
451 $ R1, R2, W, IW, INFO )
452 CALL CHKXER( 'sppsvx', INFOT, NOUT, LERR, OK )
454 CALL SPPSVX( 'n
', 'u
', 2, 0, A, AF, EQ, C, B, 2, X, 1, RCOND,
455 $ R1, R2, W, IW, INFO )
456 CALL CHKXER( 'sppsvx', INFOT, NOUT, LERR, OK )
458 ELSE IF( LSAMEN( 2, C2, 'pb
' ) ) THEN
464 CALL SPBSV( '/
', 0, 0, 0, A, 1, B, 1, INFO )
465 CALL CHKXER( 'spbsv ', INFOT, NOUT, LERR, OK )
467 CALL SPBSV( 'u
', -1, 0, 0, A, 1, B, 1, INFO )
468 CALL CHKXER( 'spbsv ', INFOT, NOUT, LERR, OK )
470 CALL SPBSV( 'u
', 1, -1, 0, A, 1, B, 1, INFO )
471 CALL CHKXER( 'spbsv ', INFOT, NOUT, LERR, OK )
473 CALL SPBSV( 'u
', 0, 0, -1, A, 1, B, 1, INFO )
474 CALL CHKXER( 'spbsv ', INFOT, NOUT, LERR, OK )
476 CALL SPBSV( 'u
', 1, 1, 0, A, 1, B, 2, INFO )
477 CALL CHKXER( 'spbsv ', INFOT, NOUT, LERR, OK )
479 CALL SPBSV( 'u
', 2, 0, 0, A, 1, B, 1, INFO )
480 CALL CHKXER( 'spbsv ', INFOT, NOUT, LERR, OK )
486 CALL SPBSVX( '/
', 'u
', 0, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
487 $ RCOND, R1, R2, W, IW, INFO )
488 CALL CHKXER( 'spbsvx', INFOT, NOUT, LERR, OK )
490 CALL SPBSVX( 'n
', '/
', 0, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
491 $ RCOND, R1, R2, W, IW, INFO )
492 CALL CHKXER( 'spbsvx', INFOT, NOUT, LERR, OK )
494 CALL SPBSVX( 'n
', 'u
', -1, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X,
495 $ 1, RCOND, R1, R2, W, IW, INFO )
496 CALL CHKXER( 'spbsvx', INFOT, NOUT, LERR, OK )
498 CALL SPBSVX( 'n
', 'u
', 1, -1, 0, A, 1, AF, 1, EQ, C, B, 1, X,
499 $ 1, RCOND, R1, R2, W, IW, INFO )
500 CALL CHKXER( 'spbsvx', INFOT, NOUT, LERR, OK )
502 CALL SPBSVX( 'n
', 'u
', 0, 0, -1, A, 1, AF, 1, EQ, C, B, 1, X,
503 $ 1, RCOND, R1, R2, W, IW, INFO )
504 CALL CHKXER( 'spbsvx', INFOT, NOUT, LERR, OK )
506 CALL SPBSVX( 'n
', 'u
', 1, 1, 0, A, 1, AF, 2, EQ, C, B, 2, X, 2,
507 $ RCOND, R1, R2, W, IW, INFO )
508 CALL CHKXER( 'spbsvx', INFOT, NOUT, LERR, OK )
510 CALL SPBSVX( 'n
', 'u
', 1, 1, 0, A, 2, AF, 1, EQ, C, B, 2, X, 2,
511 $ RCOND, R1, R2, W, IW, INFO )
512 CALL CHKXER( 'spbsvx', INFOT, NOUT, LERR, OK )
515 CALL SPBSVX( 'f
', 'u
', 0, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
516 $ RCOND, R1, R2, W, IW, INFO )
517 CALL CHKXER( 'spbsvx', INFOT, NOUT, LERR, OK )
520 CALL SPBSVX( 'f
', 'u
', 1, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
521 $ RCOND, R1, R2, W, IW, INFO )
522 CALL CHKXER( 'spbsvx', INFOT, NOUT, LERR, OK )
524 CALL SPBSVX( 'n
', 'u
', 2, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 2,
525 $ RCOND, R1, R2, W, IW, INFO )
526 CALL CHKXER( 'spbsvx', INFOT, NOUT, LERR, OK )
528 CALL SPBSVX( 'n
', 'u
', 2, 0, 0, A, 1, AF, 1, EQ, C, B, 2, X, 1,
529 $ RCOND, R1, R2, W, IW, INFO )
530 CALL CHKXER( 'spbsvx', INFOT, NOUT, LERR, OK )
532 ELSE IF( LSAMEN( 2, C2, 'pt
' ) ) THEN
538 CALL SPTSV( -1, 0, A( 1, 1 ), A( 1, 2 ), B, 1, INFO )
539 CALL CHKXER( 'sptsv ', INFOT, NOUT, LERR, OK )
541 CALL SPTSV( 0, -1, A( 1, 1 ), A( 1, 2 ), B, 1, INFO )
542 CALL CHKXER( 'sptsv ', INFOT, NOUT, LERR, OK )
544 CALL SPTSV( 2, 0, A( 1, 1 ), A( 1, 2 ), B, 1, INFO )
545 CALL CHKXER( 'sptsv ', INFOT, NOUT, LERR, OK )
551 CALL SPTSVX( '/
', 0, 0, A( 1, 1 ), A( 1, 2 ), AF( 1, 1 ),
552 $ AF( 1, 2 ), B, 1, X, 1, RCOND, R1, R2, W, INFO )
553 CALL CHKXER( 'sptsvx', INFOT, NOUT, LERR, OK )
555 CALL SPTSVX( 'n
', -1, 0, A( 1, 1 ), A( 1, 2 ), AF( 1, 1 ),
556 $ AF( 1, 2 ), B, 1, X, 1, RCOND, R1, R2, W, INFO )
557 CALL CHKXER( 'sptsvx', INFOT, NOUT, LERR, OK )
559 CALL SPTSVX( 'n
', 0, -1, A( 1, 1 ), A( 1, 2 ), AF( 1, 1 ),
560 $ AF( 1, 2 ), B, 1, X, 1, RCOND, R1, R2, W, INFO )
561 CALL CHKXER( 'sptsvx', INFOT, NOUT, LERR, OK )
563 CALL SPTSVX( 'n
', 2, 0, A( 1, 1 ), A( 1, 2 ), AF( 1, 1 ),
564 $ AF( 1, 2 ), B, 1, X, 2, RCOND, R1, R2, W, INFO )
565 CALL CHKXER( 'sptsvx', INFOT, NOUT, LERR, OK )
567 CALL SPTSVX( 'n
', 2, 0, A( 1, 1 ), A( 1, 2 ), AF( 1, 1 ),
568 $ AF( 1, 2 ), B, 2, X, 1, RCOND, R1, R2, W, INFO )
569 CALL CHKXER( 'sptsvx', INFOT, NOUT, LERR, OK )
571 ELSE IF( LSAMEN( 2, C2, 'sy
' ) ) THEN
577 CALL SSYSV( '/
', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
578 CALL CHKXER( 'ssysv ', INFOT, NOUT, LERR, OK )
580 CALL SSYSV( 'u
', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
581 CALL CHKXER( 'ssysv ', INFOT, NOUT, LERR, OK )
583 CALL SSYSV( 'u
', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
584 CALL CHKXER( 'ssysv ', INFOT, NOUT, LERR, OK )
586 CALL SSYSV( 'u
', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
587 CALL CHKXER( 'ssysv ', INFOT, NOUT, LERR, OK )
589 CALL SSYSV( 'u
', 0, 0, A, 1, IP, B, 1, W, 0, INFO )
590 CALL CHKXER( 'ssysv ', INFOT, NOUT, LERR, OK )
592 CALL SSYSV( 'u
', 0, 0, A, 1, IP, B, 1, W, -2, INFO )
593 CALL CHKXER( 'ssysv ', INFOT, NOUT, LERR, OK )
599 CALL SSYSVX( '/
', 'u
', 0, 0, A, 1, AF, 1, IP, B, 1, X, 1,
600 $ RCOND, R1, R2, W, 1, IW, INFO )
601 CALL CHKXER( 'ssysvx', INFOT, NOUT, LERR, OK )
603 CALL SSYSVX( 'n
', '/
', 0, 0, A, 1, AF, 1, IP, B, 1, X, 1,
604 $ RCOND, R1, R2, W, 1, IW, INFO )
605 CALL CHKXER( 'ssysvx', INFOT, NOUT, LERR, OK )
607 CALL SSYSVX( 'n
', 'u
', -1, 0, A, 1, AF, 1, IP, B, 1, X, 1,
608 $ RCOND, R1, R2, W, 1, IW, INFO )
609 CALL CHKXER( 'ssysvx', INFOT, NOUT, LERR, OK )
611 CALL SSYSVX( 'n
', 'u
', 0, -1, A, 1, AF, 1, IP, B, 1, X, 1,
612 $ RCOND, R1, R2, W, 1, IW, INFO )
613 CALL CHKXER( 'ssysvx', INFOT, NOUT, LERR, OK )
615 CALL SSYSVX( 'n
', 'u
', 2, 0, A, 1, AF, 2, IP, B, 2, X, 2,
616 $ RCOND, R1, R2, W, 4, IW, INFO )
617 CALL CHKXER( 'ssysvx', INFOT, NOUT, LERR, OK )
619 CALL SSYSVX( 'n
', 'u
', 2, 0, A, 2, AF, 1, IP, B, 2, X, 2,
620 $ RCOND, R1, R2, W, 4, IW, INFO )
621 CALL CHKXER( 'ssysvx', INFOT, NOUT, LERR, OK )
623 CALL SSYSVX( 'n
', 'u
', 2, 0, A, 2, AF, 2, IP, B, 1, X, 2,
624 $ RCOND, R1, R2, W, 4, IW, INFO )
625 CALL CHKXER( 'ssysvx', INFOT, NOUT, LERR, OK )
627 CALL SSYSVX( 'n
', 'u
', 2, 0, A, 2, AF, 2, IP, B, 2, X, 1,
628 $ RCOND, R1, R2, W, 4, IW, INFO )
629 CALL CHKXER( 'ssysvx', INFOT, NOUT, LERR, OK )
631 CALL SSYSVX( 'n
', 'u
', 2, 0, A, 2, AF, 2, IP, B, 2, X, 2,
632 $ RCOND, R1, R2, W, 3, IW, INFO )
633 CALL CHKXER( 'ssysvx', INFOT, NOUT, LERR, OK )
636 ELSE IF( LSAMEN( 2, C2, 'sr
' ) ) THEN
642 CALL SSYSV_ROOK( '/
', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
643 CALL CHKXER( 'ssysv_rook', INFOT, NOUT, LERR, OK )
645 CALL SSYSV_ROOK( 'u
', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
646 CALL CHKXER( 'ssysv_rook', INFOT, NOUT, LERR, OK )
648 CALL SSYSV_ROOK( 'u
', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
649 CALL CHKXER( 'ssysv_rook', INFOT, NOUT, LERR, OK )
651 CALL SSYSV_ROOK( 'u
', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
652 CALL CHKXER( 'ssysv_rook', INFOT, NOUT, LERR, OK )
654 CALL SSYSV_ROOK( 'u
', 0, 0, A, 1, IP, B, 1, W, 0, INFO )
655 CALL CHKXER( 'ssysv_rook', INFOT, NOUT, LERR, OK )
657 CALL SSYSV_ROOK( 'u
', 0, 0, A, 1, IP, B, 1, W, -2, INFO )
658 CALL CHKXER( 'ssysv_rook', INFOT, NOUT, LERR, OK )
660 ELSE IF( LSAMEN( 2, C2, 'sk
' ) ) THEN
674 CALL SSYSV_RK( '/
', 0, 0, A, 1, E, IP, B, 1, W, 1, INFO )
675 CALL CHKXER( 'ssysv_rk', INFOT, NOUT, LERR, OK )
677 CALL SSYSV_RK( 'u
', -1, 0, A, 1, E, IP, B, 1, W, 1, INFO )
678 CALL CHKXER( 'ssysv_rk', INFOT, NOUT, LERR, OK )
680 CALL SSYSV_RK( 'u
', 0, -1, A, 1, E, IP, B, 1, W, 1, INFO )
681 CALL CHKXER( 'ssysv_rk', INFOT, NOUT, LERR, OK )
683 CALL SSYSV_RK( 'u
', 2, 0, A, 1, E, IP, B, 2, W, 1, INFO )
684 CALL CHKXER( 'ssysv_rk', INFOT, NOUT, LERR, OK )
686 CALL SSYSV_RK( 'u
', 2, 0, A, 2, E, IP, B, 1, W, 1, INFO )
687 CALL CHKXER( 'ssysv_rk', INFOT, NOUT, LERR, OK )
689 CALL SSYSV_RK( 'u', 0, 0, a, 1, e, ip, b, 1, w, 0, info )
690 CALL chkxer(
'SSYSV_RK', infot, nout, lerr, ok )
692 CALL ssysv_rk(
'U', 0, 0, a, 1, e, ip, b, 1, w, -2, info )
693 CALL chkxer(
'SSYSV_RK', infot, nout, lerr, ok )
695 ELSE IF( lsamen( 2, c2,
'SA' ) )
THEN
701 CALL ssysv_aa(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
702 CALL chkxer(
'SSYSV_AA', infot, nout, lerr, ok )
704 CALL ssysv_aa(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
705 CALL chkxer(
'SSYSV_AA', infot, nout, lerr, ok )
707 CALL ssysv_aa(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
708 CALL chkxer(
'SSYSV_AA', infot, nout, lerr, ok )
710 CALL ssysv_aa(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
711 CALL chkxer(
'SSYSV_AA', infot, nout, lerr, ok )
713 ELSE IF( lsamen( 2, c2,
'S2' ) )
THEN
717 srnamt =
'SSYSV_AA_2STAGE'
719 CALL ssysv_aa_2stage(
'/', 0, 0, a, 1, a, 1, ip, ip, b, 1,
721 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
723 CALL ssysv_aa_2stage(
'U', -1, 0, a, 1, a, 1, ip, ip, b, 1,
725 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
727 CALL ssysv_aa_2stage(
'U', 0, -1, a, 1, a, 1, ip, ip, b, 1,
729 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
731 CALL ssysv_aa_2stage(
'U', 2, 1, a, 1, a, 1, ip, ip, b, 1,
733 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
735 CALL ssysv_aa_2stage(
'U', 2, 1, a, 2, a, 8, ip, ip, b, 1,
737 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
739 CALL ssysv_aa_2stage(
'U', 2, 1, a, 2, a, 1, ip, ip, b, 2,
741 CALL chkxer(
'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
743 ELSE IF( lsamen( 2, c2,
'SP' ) )
THEN
749 CALL sspsv(
'/', 0, 0, a, ip, b, 1, info )
750 CALL chkxer(
'SSPSV ', infot, nout, lerr, ok )
752 CALL sspsv(
'U', -1, 0, a, ip, b, 1, info )
753 CALL chkxer(
'SSPSV ', infot, nout, lerr, ok )
755 CALL sspsv(
'U', 0, -1, a, ip, b, 1, info )
756 CALL chkxer(
'SSPSV ', infot, nout, lerr, ok )
758 CALL sspsv(
'U', 2, 0, a, ip, b, 1, info )
759 CALL chkxer(
'SSPSV ', infot, nout, lerr, ok )
765 CALL sspsvx(
'/',
'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
767 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
769 CALL sspsvx(
'N',
'/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
771 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
773 CALL sspsvx(
'N',
'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
775 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
777 CALL sspsvx(
'N',
'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
779 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
781 CALL sspsvx(
'N',
'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
783 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
785 CALL sspsvx(
'N',
'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
787 CALL chkxer(
'SSPSVX', infot, nout, lerr, ok )
793 WRITE( nout, fmt = 9999 )path
795 WRITE( nout, fmt = 9998 )path
798 9999
FORMAT( 1x, a3,
' drivers passed the tests of the error exits' )
799 9998
FORMAT(
' *** ', a3,
' drivers failed the tests of the error ',