76 parameter( nmax = 3, liw = 12*nmax, lw = 20*nmax )
80 INTEGER I, INFO, J, M, N, NSPLIT, NT
83 INTEGER I1( NMAX ), I2( NMAX ), I3( NMAX ), IW( LIW )
84 REAL A( NMAX, NMAX ), C( NMAX, NMAX ), D( NMAX ),
85 $ E( NMAX ), Q( NMAX, NMAX ), R( NMAX ),
86 $ TAU( NMAX ), W( LW ), X( NMAX ),
110 COMMON / infoc / infot, nout, ok, lerr
111 COMMON / srnamc / srnamt
119 WRITE( nout, fmt = * )
126 a( i, j ) = 1. / real( i+j )
141 IF( lsamen( 2, c2,
'ST' ) )
THEN
147 CALL ssytrd(
'/', 0, a, 1, d, e, tau, w, 1, info )
148 CALL chkxer(
'SSYTRD', infot, nout, lerr, ok )
150 CALL ssytrd(
'U', -1, a, 1, d, e, tau, w, 1, info )
151 CALL chkxer(
'SSYTRD', infot, nout, lerr, ok )
153 CALL ssytrd(
'U', 2, a, 1, d, e, tau, w, 1, info )
154 CALL chkxer(
'SSYTRD', infot, nout, lerr, ok )
156 CALL ssytrd(
'U', 0, a, 1, d, e, tau, w, 0, info )
157 CALL chkxer(
'SSYTRD', infot, nout, lerr, ok )
164 CALL ssytd2(
'/', 0, a, 1, d, e, tau, info )
165 CALL chkxer(
'SSYTD2', infot, nout, lerr, ok )
167 CALL ssytd2(
'U', -1, a, 1, d, e, tau, info )
168 CALL chkxer(
'SSYTD2', infot, nout, lerr, ok )
170 CALL ssytd2(
'U', 2, a, 1, d, e, tau, info )
171 CALL chkxer(
'SSYTD2', infot, nout, lerr, ok )
176 srnamt =
'SSYTRD_2STAGE'
180 CALL chkxer(
'SSYTRD_2STAGE', infot, nout, lerr, ok )
184 CALL chkxer(
'SSYTRD_2STAGE', infot, nout, lerr, ok )
188 CALL chkxer(
'SSYTRD_2STAGE', infot, nout, lerr, ok )
192 CALL chkxer(
'SSYTRD_2STAGE', infot, nout, lerr, ok )
196 CALL chkxer(
'SSYTRD_2STAGE', infot, nout, lerr, ok )
200 CALL chkxer(
'SSYTRD_2STAGE', infot, nout, lerr, ok )
204 CALL chkxer(
'SSYTRD_2STAGE', infot, nout, lerr, ok )
209 srnamt =
'SSYTRD_SY2SB'
211 CALL ssytrd_sy2sb(
'/', 0, 0, a, 1, c, 1, tau, w, 1, info )
212 CALL chkxer(
'SSYTRD_SY2SB', infot, nout, lerr, ok )
214 CALL ssytrd_sy2sb(
'U', -1, 0, a, 1, c, 1, tau, w, 1, info )
215 CALL chkxer(
'SSYTRD_SY2SB', infot, nout, lerr, ok )
217 CALL ssytrd_sy2sb(
'U', 0, -1, a, 1, c, 1, tau, w, 1, info )
218 CALL chkxer(
'SSYTRD_SY2SB', infot, nout, lerr, ok )
220 CALL ssytrd_sy2sb(
'U', 2, 0, a, 1, c, 1, tau, w, 1, info )
221 CALL chkxer(
'SSYTRD_SY2SB', infot, nout, lerr, ok )
223 CALL ssytrd_sy2sb(
'U', 0, 2, a, 1, c, 1, tau, w, 1, info )
224 CALL chkxer(
'SSYTRD_SY2SB', infot, nout, lerr, ok )
226 CALL ssytrd_sy2sb(
'U', 0, 0, a, 1, c, 1, tau, w, 0, info )
227 CALL chkxer(
'SSYTRD_SY2SB', infot, nout, lerr, ok )
232 srnamt =
'SSYTRD_SB2ST'
236 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
240 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
244 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
248 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
252 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
256 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
260 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
264 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
268 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
275 CALL sorgtr(
'/', 0, a, 1, tau, w, 1, info )
276 CALL chkxer(
'SORGTR', infot, nout, lerr, ok )
278 CALL sorgtr(
'U', -1, a, 1, tau, w, 1, info )
279 CALL chkxer(
'SORGTR', infot, nout, lerr, ok )
281 CALL sorgtr(
'U', 2, a, 1, tau, w, 1, info )
282 CALL chkxer(
'SORGTR', infot, nout, lerr, ok )
284 CALL sorgtr(
'U', 3, a, 3, tau, w, 1, info )
285 CALL chkxer(
'SORGTR', infot, nout, lerr, ok )
292 CALL sormtr(
'/',
'U',
'N', 0, 0, a, 1, tau, c, 1, w, 1, info )
293 CALL chkxer(
'SORMTR', infot, nout, lerr, ok )
295 CALL sormtr(
'L',
'/',
'N', 0, 0, a, 1, tau, c, 1, w, 1, info )
296 CALL chkxer(
'SORMTR', infot, nout, lerr, ok )
298 CALL sormtr(
'L',
'U',
'/', 0, 0, a, 1, tau, c, 1, w, 1, info )
299 CALL chkxer(
'SORMTR', infot, nout, lerr, ok )
301 CALL sormtr(
'L',
'U',
'N', -1, 0, a, 1, tau, c, 1, w, 1,
303 CALL chkxer(
'SORMTR', infot, nout, lerr, ok )
305 CALL sormtr(
'L',
'U',
'N', 0, -1, a, 1, tau, c, 1, w, 1,
307 CALL chkxer(
'SORMTR', infot, nout, lerr, ok )
309 CALL sormtr(
'L',
'U',
'N', 2, 0, a, 1, tau, c, 2, w, 1, info )
310 CALL chkxer(
'SORMTR', infot, nout, lerr, ok )
312 CALL sormtr(
'R',
'U',
'N', 0, 2, a, 1, tau, c, 1, w, 1, info )
313 CALL chkxer(
'SORMTR', infot, nout, lerr, ok )
315 CALL sormtr(
'L',
'U',
'N', 2, 0, a, 2, tau, c, 1, w, 1, info )
316 CALL chkxer(
'SORMTR', infot, nout, lerr, ok )
318 CALL sormtr(
'L',
'U',
'N', 0, 2, a, 1, tau, c, 1, w, 1, info )
319 CALL chkxer(
'SORMTR', infot, nout, lerr, ok )
321 CALL sormtr(
'R',
'U',
'N', 2, 0, a, 1, tau, c, 2, w, 1, info )
322 CALL chkxer(
'SORMTR', infot, nout, lerr, ok )
329 CALL ssptrd(
'/', 0, a, d, e, tau, info )
330 CALL chkxer(
'SSPTRD', infot, nout, lerr, ok )
332 CALL ssptrd(
'U', -1, a, d, e, tau, info )
333 CALL chkxer(
'SSPTRD', infot, nout, lerr, ok )
340 CALL sopgtr(
'/', 0, a, tau, z, 1, w, info )
341 CALL chkxer(
'SOPGTR', infot, nout, lerr, ok )
343 CALL sopgtr(
'U', -1, a, tau, z, 1, w, info )
344 CALL chkxer(
'SOPGTR', infot, nout, lerr, ok )
346 CALL sopgtr(
'U', 2, a, tau, z, 1, w, info )
347 CALL chkxer(
'SOPGTR', infot, nout, lerr, ok )
354 CALL sopmtr(
'/',
'U',
'N', 0, 0, a, tau, c, 1, w, info )
355 CALL chkxer(
'SOPMTR', infot, nout, lerr, ok )
357 CALL sopmtr(
'L',
'/',
'N', 0, 0, a, tau, c, 1, w, info )
358 CALL chkxer(
'SOPMTR', infot, nout, lerr, ok )
360 CALL sopmtr(
'L',
'U',
'/', 0, 0, a, tau, c, 1, w, info )
361 CALL chkxer(
'SOPMTR', infot, nout, lerr, ok )
363 CALL sopmtr(
'L',
'U',
'N', -1, 0, a, tau, c, 1, w, info )
364 CALL chkxer(
'SOPMTR', infot, nout, lerr, ok )
366 CALL sopmtr(
'L',
'U',
'N', 0, -1, a, tau, c, 1, w, info )
367 CALL chkxer(
'SOPMTR', infot, nout, lerr, ok )
369 CALL sopmtr(
'L',
'U',
'N', 2, 0, a, tau, c, 1, w, info )
370 CALL chkxer(
'SOPMTR', infot, nout, lerr, ok )
377 CALL spteqr(
'/', 0, d, e, z, 1, w, info )
378 CALL chkxer(
'SPTEQR', infot, nout, lerr, ok )
380 CALL spteqr(
'N', -1, d, e, z, 1, w, info )
381 CALL chkxer(
'SPTEQR', infot, nout, lerr, ok )
383 CALL spteqr(
'V', 2, d, e, z, 1, w, info )
384 CALL chkxer(
'SPTEQR', infot, nout, lerr, ok )
391 CALL sstebz(
'/',
'E', 0, 0.0, 1.0, 1, 0, 0.0, d, e, m, nsplit,
392 $ x, i1, i2, w, iw, info )
393 CALL chkxer(
'SSTEBZ', infot, nout, lerr, ok )
395 CALL sstebz(
'A',
'/', 0, 0.0, 0.0, 0, 0, 0.0, d, e, m, nsplit,
396 $ x, i1, i2, w, iw, info )
397 CALL chkxer(
'SSTEBZ', infot, nout, lerr, ok )
399 CALL sstebz(
'A',
'E', -1, 0.0, 0.0, 0, 0, 0.0, d, e, m,
400 $ nsplit, x, i1, i2, w, iw, info )
401 CALL chkxer(
'SSTEBZ', infot, nout, lerr, ok )
403 CALL sstebz(
'V',
'E', 0, 0.0, 0.0, 0, 0, 0.0, d, e, m, nsplit,
404 $ x, i1, i2, w, iw, info )
405 CALL chkxer(
'SSTEBZ', infot, nout, lerr, ok )
407 CALL sstebz(
'I',
'E', 0, 0.0, 0.0, 0, 0, 0.0, d, e, m, nsplit,
408 $ x, i1, i2, w, iw, info )
409 CALL chkxer(
'SSTEBZ', infot, nout, lerr, ok )
411 CALL sstebz(
'I',
'E', 1, 0.0, 0.0, 2, 1, 0.0, d, e, m, nsplit,
412 $ x, i1, i2, w, iw, info )
413 CALL chkxer(
'SSTEBZ', infot, nout, lerr, ok )
415 CALL sstebz(
'I',
'E', 1, 0.0, 0.0, 1, 0, 0.0, d, e, m, nsplit,
416 $ x, i1, i2, w, iw, info )
417 CALL chkxer(
'SSTEBZ', infot, nout, lerr, ok )
419 CALL sstebz(
'I',
'E', 1, 0.0, 0.0, 1, 2, 0.0, d, e, m, nsplit,
420 $ x, i1, i2, w, iw, info )
421 CALL chkxer(
'SSTEBZ', infot, nout, lerr, ok )
428 CALL sstein( -1, d, e, 0, x, i1, i2, z, 1, w, iw, i3, info )
429 CALL chkxer(
'SSTEIN', infot, nout, lerr, ok )
431 CALL sstein( 0, d, e, -1, x, i1, i2, z, 1, w, iw, i3, info )
432 CALL chkxer(
'SSTEIN', infot, nout, lerr, ok )
434 CALL sstein( 0, d, e, 1, x, i1, i2, z, 1, w, iw, i3, info )
435 CALL chkxer(
'SSTEIN', infot, nout, lerr, ok )
437 CALL sstein( 2, d, e, 0, x, i1, i2, z, 1, w, iw, i3, info )
438 CALL chkxer(
'SSTEIN', infot, nout, lerr, ok )
445 CALL ssteqr(
'/', 0, d, e, z, 1, w, info )
446 CALL chkxer(
'SSTEQR', infot, nout, lerr, ok )
448 CALL ssteqr(
'N', -1, d, e, z, 1, w, info )
449 CALL chkxer(
'SSTEQR', infot, nout, lerr, ok )
451 CALL ssteqr(
'V', 2, d, e, z, 1, w, info )
452 CALL chkxer(
'SSTEQR', infot, nout, lerr, ok )
459 CALL ssterf( -1, d, e, info )
460 CALL chkxer(
'SSTERF', infot, nout, lerr, ok )
467 CALL sstedc(
'/', 0, d, e, z, 1, w, 1, iw, 1, info )
468 CALL chkxer(
'SSTEDC', infot, nout, lerr, ok )
470 CALL sstedc(
'N', -1, d, e, z, 1, w, 1, iw, 1, info )
471 CALL chkxer(
'SSTEDC', infot, nout, lerr, ok )
473 CALL sstedc(
'V', 2, d, e, z, 1, w, 23, iw, 28, info )
474 CALL chkxer(
'SSTEDC', infot, nout, lerr, ok )
476 CALL sstedc(
'N', 1, d, e, z, 1, w, 0, iw, 1, info )
477 CALL chkxer(
'SSTEDC', infot, nout, lerr, ok )
479 CALL sstedc(
'I', 2, d, e, z, 2, w, 0, iw, 12, info )
480 CALL chkxer(
'SSTEDC', infot, nout, lerr, ok )
482 CALL sstedc(
'V', 2, d, e, z, 2, w, 0, iw, 28, info )
483 CALL chkxer(
'SSTEDC', infot, nout, lerr, ok )
485 CALL sstedc(
'N', 1, d, e, z, 1, w, 1, iw, 0, info )
486 CALL chkxer(
'SSTEDC', infot, nout, lerr, ok )
488 CALL sstedc(
'I', 2, d, e, z, 2, w, 19, iw, 0, info )
489 CALL chkxer(
'SSTEDC', infot, nout, lerr, ok )
491 CALL sstedc(
'V', 2, d, e, z, 2, w, 23, iw, 0, info )
492 CALL chkxer(
'SSTEDC', infot, nout, lerr, ok )
499 CALL sstevd(
'/', 0, d, e, z, 1, w, 1, iw, 1, info )
500 CALL chkxer(
'SSTEVD', infot, nout, lerr, ok )
502 CALL sstevd(
'N', -1, d, e, z, 1, w, 1, iw, 1, info )
503 CALL chkxer(
'SSTEVD', infot, nout, lerr, ok )
505 CALL sstevd(
'V', 2, d, e, z, 1, w, 19, iw, 12, info )
506 CALL chkxer(
'SSTEVD', infot, nout, lerr, ok )
508 CALL sstevd(
'N', 1, d, e, z, 1, w, 0, iw, 1, info )
509 CALL chkxer(
'SSTEVD', infot, nout, lerr, ok )
511 CALL sstevd(
'V', 2, d, e, z, 2, w, 12, iw, 12, info )
512 CALL chkxer(
'SSTEVD', infot, nout, lerr, ok )
514 CALL sstevd(
'N', 0, d, e, z, 1, w, 1, iw, 0, info )
515 CALL chkxer(
'SSTEVD', infot, nout, lerr, ok )
517 CALL sstevd(
'V', 2, d, e, z, 2, w, 19, iw, 11, info )
518 CALL chkxer(
'SSTEVD', infot, nout, lerr, ok )
525 CALL sstev(
'/', 0, d, e, z, 1, w, info )
526 CALL chkxer(
'SSTEV ', infot, nout, lerr, ok )
528 CALL sstev(
'N', -1, d, e, z, 1, w, info )
529 CALL chkxer(
'SSTEV ', infot, nout, lerr, ok )
531 CALL sstev(
'V', 2, d, e, z, 1, w, info )
532 CALL chkxer(
'SSTEV ', infot, nout, lerr, ok )
539 CALL sstevx(
'/',
'A', 0, d, e, 0.0, 0.0, 0, 0, 0.0, m, x, z,
540 $ 1, w, iw, i3, info )
541 CALL chkxer(
'SSTEVX', infot, nout, lerr, ok )
543 CALL sstevx(
'N',
'/', 0, d, e, 0.0, 1.0, 1, 0, 0.0, m, x, z,
544 $ 1, w, iw, i3, info )
545 CALL chkxer(
'SSTEVX', infot, nout, lerr, ok )
547 CALL sstevx(
'N',
'A', -1, d, e, 0.0, 0.0, 0, 0, 0.0, m, x, z,
548 $ 1, w, iw, i3, info )
549 CALL chkxer(
'SSTEVX', infot, nout, lerr, ok )
551 CALL sstevx(
'N',
'V', 1, d, e, 0.0, 0.0, 0, 0, 0.0, m, x, z,
552 $ 1, w, iw, i3, info )
553 CALL chkxer(
'SSTEVX', infot, nout, lerr, ok )
555 CALL sstevx(
'N',
'I', 1, d, e, 0.0, 0.0, 0, 0, 0.0, m, x, z,
556 $ 1, w, iw, i3, info )
557 CALL chkxer(
'SSTEVX', infot, nout, lerr, ok )
559 CALL sstevx(
'N',
'I', 1, d, e, 0.0, 0.0, 2, 1, 0.0, m, x, z,
560 $ 1, w, iw, i3, info )
561 CALL chkxer(
'SSTEVX', infot, nout, lerr, ok )
563 CALL sstevx(
'N',
'I', 2, d, e, 0.0, 0.0, 2, 1, 0.0, m, x, z,
564 $ 1, w, iw, i3, info )
565 CALL chkxer(
'SSTEVX', infot, nout, lerr, ok )
567 CALL sstevx(
'N',
'I', 1, d, e, 0.0, 0.0, 1, 2, 0.0, m, x, z,
568 $ 1, w, iw, i3, info )
569 CALL chkxer(
'SSTEVX', infot, nout, lerr, ok )
571 CALL sstevx(
'V',
'A', 2, d, e, 0.0, 0.0, 0, 0, 0.0, m, x, z,
572 $ 1, w, iw, i3, info )
573 CALL chkxer(
'SSTEVX', infot, nout, lerr, ok )
581 CALL sstevr(
'/',
'A', 0, d, e, 0.0, 0.0, 1, 1, 0.0, m, r, z,
582 $ 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
583 CALL chkxer(
'SSTEVR', infot, nout, lerr, ok )
585 CALL sstevr(
'V',
'/', 0, d, e, 0.0, 0.0, 1, 1, 0.0, m, r, z,
586 $ 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
587 CALL chkxer(
'SSTEVR', infot, nout, lerr, ok )
589 CALL sstevr(
'V',
'A', -1, d, e, 0.0, 0.0, 1, 1, 0.0, m, r, z,
590 $ 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
591 CALL chkxer(
'SSTEVR', infot, nout, lerr, ok )
593 CALL sstevr(
'V',
'V', 1, d, e, 0.0, 0.0, 1, 1, 0.0, m, r, z,
594 $ 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
595 CALL chkxer(
'SSTEVR', infot, nout, lerr, ok )
597 CALL sstevr(
'V',
'I', 1, d, e, 0.0, 0.0, 0, 1, 0.0, m, w, z,
598 $ 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
599 CALL chkxer(
'SSTEVR', infot, nout, lerr, ok )
602 CALL sstevr(
'V',
'I', 2, d, e, 0.0, 0.0, 2, 1, 0.0, m, w, z,
603 $ 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
604 CALL chkxer(
'SSTEVR', infot, nout, lerr, ok )
607 CALL sstevr(
'V',
'I', 1, d, e, 0.0, 0.0, 1, 1, 0.0, m, w, z,
608 $ 0, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
609 CALL chkxer(
'SSTEVR', infot, nout, lerr, ok )
611 CALL sstevr(
'V',
'I', 1, d, e, 0.0, 0.0, 1, 1, 0.0, m, w, z,
612 $ 1, iw, x, 20*n-1, iw( 2*n+1 ), 10*n, info )
613 CALL chkxer(
'SSTEVR', infot, nout, lerr, ok )
615 CALL sstevr(
'V',
'I', 1, d, e, 0.0, 0.0, 1, 1, 0.0, m, w, z,
616 $ 1, iw, x, 20*n, iw( 2*n+1 ), 10*n-1, info )
617 CALL chkxer(
'SSTEVR', infot, nout, lerr, ok )
624 CALL ssyevd(
'/',
'U', 0, a, 1, x, w, 1, iw, 1, info )
625 CALL chkxer(
'SSYEVD', infot, nout, lerr, ok )
627 CALL ssyevd(
'N',
'/', 0, a, 1, x, w, 1, iw, 1, info )
628 CALL chkxer(
'SSYEVD', infot, nout, lerr, ok )
630 CALL ssyevd(
'N',
'U', -1, a, 1, x, w, 1, iw, 1, info )
631 CALL chkxer(
'SSYEVD', infot, nout, lerr, ok )
633 CALL ssyevd(
'N',
'U', 2, a, 1, x, w, 3, iw, 1, info )
634 CALL chkxer(
'SSYEVD', infot, nout, lerr, ok )
636 CALL ssyevd(
'N',
'U', 1, a, 1, x, w, 0, iw, 1, info )
637 CALL chkxer(
'SSYEVD', infot, nout, lerr, ok )
639 CALL ssyevd(
'N',
'U', 2, a, 2, x, w, 4, iw, 1, info )
640 CALL chkxer(
'SSYEVD', infot, nout, lerr, ok )
642 CALL ssyevd(
'V',
'U', 2, a, 2, x, w, 20, iw, 12, info )
643 CALL chkxer(
'SSYEVD', infot, nout, lerr, ok )
645 CALL ssyevd(
'N',
'U', 1, a, 1, x, w, 1, iw, 0, info )
646 CALL chkxer(
'SSYEVD', infot, nout, lerr, ok )
648 CALL ssyevd(
'N',
'U', 2, a, 2, x, w, 5, iw, 0, info )
649 CALL chkxer(
'SSYEVD', infot, nout, lerr, ok )
651 CALL ssyevd(
'V',
'U', 2, a, 2, x, w, 27, iw, 11, info )
652 CALL chkxer(
'SSYEVD', infot, nout, lerr, ok )
657 srnamt =
'SSYEVD_2STAGE'
659 CALL ssyevd_2stage(
'/',
'U', 0, a, 1, x, w, 1, iw, 1, info )
660 CALL chkxer(
'SSYEVD_2STAGE', infot, nout, lerr, ok )
662 CALL ssyevd_2stage(
'V',
'U', 0, a, 1, x, w, 1, iw, 1, info )
663 CALL chkxer(
'SSYEVD_2STAGE', infot, nout, lerr, ok )
665 CALL ssyevd_2stage(
'N',
'/', 0, a, 1, x, w, 1, iw, 1, info )
666 CALL chkxer(
'SSYEVD_2STAGE', infot, nout, lerr, ok )
668 CALL ssyevd_2stage(
'N',
'U', -1, a, 1, x, w, 1, iw, 1, info )
669 CALL chkxer(
'SSYEVD_2STAGE', infot, nout, lerr, ok )
671 CALL ssyevd_2stage(
'N',
'U', 2, a, 1, x, w, 3, iw, 1, info )
672 CALL chkxer(
'SSYEVD_2STAGE', infot, nout, lerr, ok )
674 CALL ssyevd_2stage(
'N',
'U', 1, a, 1, x, w, 0, iw, 1, info )
675 CALL chkxer(
'SSYEVD_2STAGE', infot, nout, lerr, ok )
677 CALL ssyevd_2stage(
'N',
'U', 2, a, 2, x, w, 4, iw, 1, info )
678 CALL chkxer(
'SSYEVD_2STAGE', infot, nout, lerr, ok )
683 CALL ssyevd_2stage(
'N',
'U', 1, a, 1, x, w, 1, iw, 0, info )
684 CALL chkxer(
'SSYEVD_2STAGE', infot, nout, lerr, ok )
686 CALL ssyevd_2stage(
'N',
'U', 2, a, 2, x, w, 25, iw, 0, info )
687 CALL chkxer(
'SSYEVD_2STAGE', infot, nout, lerr, ok )
698 CALL ssyevr(
'/',
'A',
'U', 0, a, 1, 0.0, 0.0, 1, 1, 0.0, m, r,
699 $ z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
700 CALL chkxer(
'SSYEVR', infot, nout, lerr, ok )
702 CALL ssyevr(
'V',
'/',
'U', 0, a, 1, 0.0, 0.0, 1, 1, 0.0, m, r,
703 $ z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
704 CALL chkxer(
'SSYEVR', infot, nout, lerr, ok )
706 CALL ssyevr(
'V',
'A',
'/', -1, a, 1, 0.0, 0.0, 1, 1, 0.0, m,
707 $ r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
708 CALL chkxer(
'SSYEVR', infot, nout, lerr, ok )
710 CALL ssyevr(
'V',
'A',
'U', -1, a, 1, 0.0, 0.0, 1, 1, 0.0, m,
711 $ r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
712 CALL chkxer(
'SSYEVR', infot, nout, lerr, ok )
714 CALL ssyevr(
'V',
'A',
'U', 2, a, 1, 0.0, 0.0, 1, 1, 0.0, m, r,
715 $ z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
716 CALL chkxer(
'SSYEVR', infot, nout, lerr, ok )
718 CALL ssyevr(
'V',
'V',
'U', 1, a, 1, 0.0e0, 0.0e0, 1, 1, 0.0,
719 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
720 CALL chkxer(
'SSYEVR', infot, nout, lerr, ok )
722 CALL ssyevr(
'V',
'I',
'U', 1, a, 1, 0.0e0, 0.0e0, 0, 1, 0.0,
723 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
724 CALL chkxer(
'SSYEVR', infot, nout, lerr, ok )
727 CALL ssyevr(
'V',
'I',
'U', 2, a, 2, 0.0e0, 0.0e0, 2, 1, 0.0,
728 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
729 CALL chkxer(
'SSYEVR', infot, nout, lerr, ok )
731 CALL ssyevr(
'V',
'I',
'U', 1, a, 1, 0.0e0, 0.0e0, 1, 1, 0.0,
732 $ m, r, z, 0, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
733 CALL chkxer(
'SSYEVR', infot, nout, lerr, ok )
735 CALL ssyevr(
'V',
'I',
'U', 1, a, 1, 0.0e0, 0.0e0, 1, 1, 0.0,
736 $ m, r, z, 1, iw, q, 26*n-1, iw( 2*n+1 ), 10*n,
738 CALL chkxer(
'SSYEVR', infot, nout, lerr, ok )
740 CALL ssyevr(
'V',
'I',
'U', 1, a, 1, 0.0e0, 0.0e0, 1, 1, 0.0,
741 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n-1,
743 CALL chkxer(
'SSYEVR', infot, nout, lerr, ok )
748 srnamt =
'SSYEVR_2STAGE'
752 $ 0.0e0, 0.0e0, 1, 1, 0.0e0,
753 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
754 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
757 $ 0.0e0, 0.0e0, 1, 1, 0.0e0,
758 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
759 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
762 $ 0.0e0, 0.0e0, 1, 1, 0.0e0,
763 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
764 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
767 $ 0.0e0, 0.0e0, 1, 1, 0.0e0,
768 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
769 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
772 $ 0.0e0, 0.0e0, 1, 1, 0.0e0,
773 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
774 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
777 $ 0.0e0, 0.0e0, 1, 1, 0.0e0,
778 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
779 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
782 $ 0.0e0, 0.0e0, 1, 1, 0.0e0,
783 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
784 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
787 $ 0.0e0, 0.0e0, 0, 1, 0.0e0,
788 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
789 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
792 $ 0.0e0, 0.0e0, 2, 1, 0.0e0,
793 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
794 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
797 $ 0.0e0, 0.0e0, 1, 1, 0.0e0,
798 $ m, r, z, 0, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
799 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
802 $ 0.0e0, 0.0e0, 1, 1, 0.0e0,
803 $ m, r, z, 1, iw, q, 0, iw( 2*n+1 ), 10*n,
805 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
808 $ 0.0e0, 0.0e0, 1, 1, 0.0e0,
809 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 0,
811 CALL chkxer(
'SSYEVR_2STAGE', infot, nout, lerr, ok )
818 CALL ssyev(
'/',
'U', 0, a, 1, x, w, 1, info )
819 CALL chkxer(
'SSYEV ', infot, nout, lerr, ok )
821 CALL ssyev(
'N',
'/', 0, a, 1, x, w, 1, info )
822 CALL chkxer(
'SSYEV ', infot, nout, lerr, ok )
824 CALL ssyev(
'N',
'U', -1, a, 1, x, w, 1, info )
825 CALL chkxer(
'SSYEV ', infot, nout, lerr, ok )
827 CALL ssyev(
'N',
'U', 2, a, 1, x, w, 3, info )
828 CALL chkxer(
'SSYEV ', infot, nout, lerr, ok )
830 CALL ssyev(
'N',
'U', 1, a, 1, x, w, 1, info )
831 CALL chkxer(
'SSYEV ', infot, nout, lerr, ok )
836 srnamt =
'SSYEV_2STAGE '
839 CALL chkxer(
'SSYEV_2STAGE ', infot, nout, lerr, ok )
842 CALL chkxer(
'SSYEV_2STAGE ', infot, nout, lerr, ok )
845 CALL chkxer(
'SSYEV_2STAGE ', infot, nout, lerr, ok )
848 CALL chkxer(
'SSYEV_2STAGE ', infot, nout, lerr, ok )
851 CALL chkxer(
'SSYEV_2STAGE ', infot, nout, lerr, ok )
854 CALL chkxer(
'SSYEV_2STAGE ', infot, nout, lerr, ok )
861 CALL ssyevx(
'/',
'A',
'U', 0, a, 1, 0.0, 0.0, 0, 0, 0.0, m, x,
862 $ z, 1, w, 1, iw, i3, info )
863 CALL chkxer(
'SSYEVX', infot, nout, lerr, ok )
865 CALL ssyevx(
'N',
'/',
'U', 0, a, 1, 0.0, 1.0, 1, 0, 0.0, m, x,
866 $ z, 1, w, 1, iw, i3, info )
867 CALL chkxer(
'SSYEVX', infot, nout, lerr, ok )
869 CALL ssyevx(
'N',
'A',
'/', 0, a, 1, 0.0, 0.0, 0, 0, 0.0, m, x,
870 $ z, 1, w, 1, iw, i3, info )
872 CALL ssyevx(
'N',
'A',
'U', -1, a, 1, 0.0, 0.0, 0, 0, 0.0, m,
873 $ x, z, 1, w, 1, iw, i3, info )
874 CALL chkxer(
'SSYEVX', infot, nout, lerr, ok )
876 CALL ssyevx(
'N',
'A',
'U', 2, a, 1, 0.0, 0.0, 0, 0, 0.0, m, x,
877 $ z, 1, w, 16, iw, i3, info )
878 CALL chkxer(
'SSYEVX', infot, nout, lerr, ok )
880 CALL ssyevx(
'N',
'V',
'U', 1, a, 1, 0.0, 0.0, 0, 0, 0.0, m, x,
881 $ z, 1, w, 8, iw, i3, info )
882 CALL chkxer(
'SSYEVX', infot, nout, lerr, ok )
884 CALL ssyevx(
'N',
'I',
'U', 1, a, 1, 0.0, 0.0, 0, 0, 0.0, m, x,
885 $ z, 1, w, 8, iw, i3, info )
886 CALL chkxer(
'SSYEVX', infot, nout, lerr, ok )
888 CALL ssyevx(
'N',
'I',
'U', 1, a, 1, 0.0, 0.0, 2, 1, 0.0, m, x,
889 $ z, 1, w, 8, iw, i3, info )
890 CALL chkxer(
'SSYEVX', infot, nout, lerr, ok )
892 CALL ssyevx(
'N',
'I',
'U', 2, a, 2, 0.0, 0.0, 2, 1, 0.0, m, x,
893 $ z, 1, w, 16, iw, i3, info )
894 CALL chkxer(
'SSYEVX', infot, nout, lerr, ok )
896 CALL ssyevx(
'N',
'I',
'U', 1, a, 1, 0.0, 0.0, 1, 2, 0.0, m, x,
897 $ z, 1, w, 8, iw, i3, info )
898 CALL chkxer(
'SSYEVX', infot, nout, lerr, ok )
900 CALL ssyevx(
'V',
'A',
'U', 2, a, 2, 0.0, 0.0, 0, 0, 0.0, m, x,
901 $ z, 1, w, 16, iw, i3, info )
902 CALL chkxer(
'SSYEVX', infot, nout, lerr, ok )
904 CALL ssyevx(
'V',
'A',
'U', 1, a, 1, 0.0, 0.0, 0, 0, 0.0, m, x,
905 $ z, 1, w, 0, iw, i3, info )
906 CALL chkxer(
'SSYEVX', infot, nout, lerr, ok )
911 srnamt =
'SSYEVX_2STAGE'
914 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
915 $ m, x, z, 1, w, 1, iw, i3, info )
916 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
919 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
920 $ m, x, z, 1, w, 1, iw, i3, info )
921 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
924 $ 0.0e0, 1.0e0, 1, 0, 0.0e0,
925 $ m, x, z, 1, w, 1, iw, i3, info )
926 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
929 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
930 $ m, x, z, 1, w, 1, iw, i3, info )
933 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
934 $ m, x, z, 1, w, 1, iw, i3, info )
935 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
938 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
939 $ m, x, z, 1, w, 16, iw, i3, info )
940 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
943 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
944 $ m, x, z, 1, w, 8, iw, i3, info )
945 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
948 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
949 $ m, x, z, 1, w, 8, iw, i3, info )
950 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
953 $ 0.0e0, 0.0e0, 2, 1, 0.0e0,
954 $ m, x, z, 1, w, 8, iw, i3, info )
955 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
958 $ 0.0e0, 0.0e0, 2, 1, 0.0e0,
959 $ m, x, z, 1, w, 16, iw, i3, info )
960 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
963 $ 0.0e0, 0.0e0, 1, 2, 0.0e0,
964 $ m, x, z, 1, w, 8, iw, i3, info )
965 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
968 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
969 $ m, x, z, 0, w, 16, iw, i3, info )
970 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
973 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
974 $ m, x, z, 1, w, 0, iw, i3, info )
975 CALL chkxer(
'SSYEVX_2STAGE', infot, nout, lerr, ok )
982 CALL sspevd(
'/',
'U', 0, a, x, z, 1, w, 1, iw, 1, info )
983 CALL chkxer(
'SSPEVD', infot, nout, lerr, ok )
985 CALL sspevd(
'N',
'/', 0, a, x, z, 1, w, 1, iw, 1, info )
986 CALL chkxer(
'SSPEVD', infot, nout, lerr, ok )
988 CALL sspevd(
'N',
'U', -1, a, x, z, 1, w, 1, iw, 1, info )
989 CALL chkxer(
'SSPEVD', infot, nout, lerr, ok )
991 CALL sspevd(
'V',
'U', 2, a, x, z, 1, w, 23, iw, 12, info )
992 CALL chkxer(
'SSPEVD', infot, nout, lerr, ok )
994 CALL sspevd(
'N',
'U', 1, a, x, z, 1, w, 0, iw, 1, info )
995 CALL chkxer(
'SSPEVD', infot, nout, lerr, ok )
997 CALL sspevd(
'N',
'U', 2, a, x, z, 1, w, 3, iw, 1, info )
998 CALL chkxer(
'SSPEVD', infot, nout, lerr, ok )
1000 CALL sspevd(
'V',
'U', 2, a, x, z, 2, w, 16, iw, 12, info )
1001 CALL chkxer(
'SSPEVD', infot, nout, lerr, ok )
1003 CALL sspevd(
'N',
'U', 1, a, x, z, 1, w, 1, iw, 0, info )
1004 CALL chkxer(
'SSPEVD', infot, nout, lerr, ok )
1006 CALL sspevd(
'N',
'U', 2, a, x, z, 1, w, 4, iw, 0, info )
1007 CALL chkxer(
'SSPEVD', infot, nout, lerr, ok )
1009 CALL sspevd(
'V',
'U', 2, a, x, z, 2, w, 23, iw, 11, info )
1010 CALL chkxer(
'SSPEVD', infot, nout, lerr, ok )
1017 CALL sspev(
'/',
'U', 0, a, w, z, 1, x, info )
1018 CALL chkxer(
'SSPEV ', infot, nout, lerr, ok )
1020 CALL sspev(
'N',
'/', 0, a, w, z, 1, x, info )
1021 CALL chkxer(
'SSPEV ', infot, nout, lerr, ok )
1023 CALL sspev(
'N',
'U', -1, a, w, z, 1, x, info )
1024 CALL chkxer(
'SSPEV ', infot, nout, lerr, ok )
1026 CALL sspev(
'V',
'U', 2, a, w, z, 1, x, info )
1027 CALL chkxer(
'SSPEV ', infot, nout, lerr, ok )
1034 CALL sspevx(
'/',
'A',
'U', 0, a, 0.0, 0.0, 0, 0, 0.0, m, x, z,
1035 $ 1, w, iw, i3, info )
1036 CALL chkxer(
'SSPEVX', infot, nout, lerr, ok )
1038 CALL sspevx(
'N',
'/',
'U', 0, a, 0.0, 0.0, 0, 0, 0.0, m, x, z,
1039 $ 1, w, iw, i3, info )
1040 CALL chkxer(
'SSPEVX', infot, nout, lerr, ok )
1042 CALL sspevx(
'N',
'A',
'/', 0, a, 0.0, 0.0, 0, 0, 0.0, m, x, z,
1043 $ 1, w, iw, i3, info )
1045 CALL sspevx(
'N',
'A',
'U', -1, a, 0.0, 0.0, 0, 0, 0.0, m, x,
1046 $ z, 1, w, iw, i3, info )
1047 CALL chkxer(
'SSPEVX', infot, nout, lerr, ok )
1049 CALL sspevx(
'N',
'V',
'U', 1, a, 0.0, 0.0, 0, 0, 0.0, m, x, z,
1050 $ 1, w, iw, i3, info )
1051 CALL chkxer(
'SSPEVX', infot, nout, lerr, ok )
1053 CALL sspevx(
'N',
'I',
'U', 1, a, 0.0, 0.0, 0, 0, 0.0, m, x, z,
1054 $ 1, w, iw, i3, info )
1055 CALL chkxer(
'SSPEVX', infot, nout, lerr, ok )
1057 CALL sspevx(
'N',
'I',
'U', 1, a, 0.0, 0.0, 2, 1, 0.0, m, x, z,
1058 $ 1, w, iw, i3, info )
1059 CALL chkxer(
'SSPEVX', infot, nout, lerr, ok )
1061 CALL sspevx(
'N',
'I',
'U', 2, a, 0.0, 0.0, 2, 1, 0.0, m, x, z,
1062 $ 1, w, iw, i3, info )
1063 CALL chkxer(
'SSPEVX', infot, nout, lerr, ok )
1065 CALL sspevx(
'N',
'I',
'U', 1, a, 0.0, 0.0, 1, 2, 0.0, m, x, z,
1066 $ 1, w, iw, i3, info )
1067 CALL chkxer(
'SSPEVX', infot, nout, lerr, ok )
1069 CALL sspevx(
'V',
'A',
'U', 2, a, 0.0, 0.0, 0, 0, 0.0, m, x, z,
1070 $ 1, w, iw, i3, info )
1071 CALL chkxer(
'SSPEVX', infot, nout, lerr, ok )
1076 ELSE IF( lsamen( 2, c2,
'SB' ) )
THEN
1082 CALL ssbtrd(
'/',
'U', 0, 0, a, 1, d, e, z, 1, w, info )
1083 CALL chkxer(
'SSBTRD', infot, nout, lerr, ok )
1085 CALL ssbtrd(
'N',
'/', 0, 0, a, 1, d, e, z, 1, w, info )
1086 CALL chkxer(
'SSBTRD', infot, nout, lerr, ok )
1088 CALL ssbtrd(
'N',
'U', -1, 0, a, 1, d, e, z, 1, w, info )
1089 CALL chkxer(
'SSBTRD', infot, nout, lerr, ok )
1091 CALL ssbtrd(
'N',
'U', 0, -1, a, 1, d, e, z, 1, w, info )
1092 CALL chkxer(
'SSBTRD', infot, nout, lerr, ok )
1094 CALL ssbtrd(
'N',
'U', 1, 1, a, 1, d, e, z, 1, w, info )
1095 CALL chkxer(
'SSBTRD', infot, nout, lerr, ok )
1097 CALL ssbtrd(
'V',
'U', 2, 0, a, 1, d, e, z, 1, w, info )
1098 CALL chkxer(
'SSBTRD', infot, nout, lerr, ok )
1103 srnamt =
'SSYTRD_SB2ST'
1106 $ c, 1, w, 1, info )
1107 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
1110 $ c, 1, w, 1, info )
1111 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
1114 $ c, 1, w, 1, info )
1115 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
1118 $ c, 1, w, 1, info )
1119 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
1122 $ c, 1, w, 1, info )
1123 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
1126 $ c, 1, w, 1, info )
1127 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
1130 $ c, 1, w, 1, info )
1131 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
1134 $ c, 0, w, 1, info )
1135 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
1138 $ c, 1, w, 0, info )
1139 CALL chkxer(
'SSYTRD_SB2ST', infot, nout, lerr, ok )
1146 CALL ssbevd(
'/',
'U', 0, 0, a, 1, x, z, 1, w, 1, iw, 1, info )
1147 CALL chkxer(
'SSBEVD', infot, nout, lerr, ok )
1149 CALL ssbevd(
'N',
'/', 0, 0, a, 1, x, z, 1, w, 1, iw, 1, info )
1150 CALL chkxer(
'SSBEVD', infot, nout, lerr, ok )
1152 CALL ssbevd(
'N',
'U', -1, 0, a, 1, x, z, 1, w, 1, iw, 1,
1154 CALL chkxer(
'SSBEVD', infot, nout, lerr, ok )
1156 CALL ssbevd(
'N',
'U', 0, -1, a, 1, x, z, 1, w, 1, iw, 1,
1158 CALL chkxer(
'SSBEVD', infot, nout, lerr, ok )
1160 CALL ssbevd(
'N',
'U', 2, 1, a, 1, x, z, 1, w, 4, iw, 1, info )
1161 CALL chkxer(
'SSBEVD', infot, nout, lerr, ok )
1163 CALL ssbevd(
'V',
'U', 2, 1, a, 2, x, z, 1, w, 25, iw, 12,
1165 CALL chkxer(
'SSBEVD', infot, nout, lerr, ok )
1167 CALL ssbevd(
'N',
'U', 1, 0, a, 1, x, z, 1, w, 0, iw, 1, info )
1168 CALL chkxer(
'SSBEVD', infot, nout, lerr, ok )
1170 CALL ssbevd(
'N',
'U', 2, 0, a, 1, x, z, 1, w, 3, iw, 1, info )
1171 CALL chkxer(
'SSBEVD', infot, nout, lerr, ok )
1173 CALL ssbevd(
'V',
'U', 2, 0, a, 1, x, z, 2, w, 18, iw, 12,
1175 CALL chkxer(
'SSBEVD', infot, nout, lerr, ok )
1177 CALL ssbevd(
'N',
'U', 1, 0, a, 1, x, z, 1, w, 1, iw, 0, info )
1178 CALL chkxer(
'SSBEVD', infot, nout, lerr, ok )
1180 CALL ssbevd(
'V',
'U', 2, 0, a, 1, x, z, 2, w, 25, iw, 11,
1182 CALL chkxer(
'SSBEVD', infot, nout, lerr, ok )
1187 srnamt =
'SSBEVD_2STAGE'
1189 CALL ssbevd_2stage(
'/',
'U', 0, 0, a, 1, x, z, 1, w,
1191 CALL chkxer(
'SSBEVD_2STAGE', infot, nout, lerr, ok )
1193 CALL ssbevd_2stage(
'V',
'U', 0, 0, a, 1, x, z, 1, w,
1195 CALL chkxer(
'SSBEVD_2STAGE', infot, nout, lerr, ok )
1197 CALL ssbevd_2stage(
'N',
'/', 0, 0, a, 1, x, z, 1, w,
1199 CALL chkxer(
'SSBEVD_2STAGE', infot, nout, lerr, ok )
1201 CALL ssbevd_2stage(
'N',
'U', -1, 0, a, 1, x, z, 1, w,
1203 CALL chkxer(
'SSBEVD_2STAGE', infot, nout, lerr, ok )
1205 CALL ssbevd_2stage(
'N',
'U', 0, -1, a, 1, x, z, 1, w,
1207 CALL chkxer(
'SSBEVD_2STAGE', infot, nout, lerr, ok )
1209 CALL ssbevd_2stage(
'N',
'U', 2, 1, a, 1, x, z, 1, w,
1211 CALL chkxer(
'SSBEVD_2STAGE', infot, nout, lerr, ok )
1217 CALL ssbevd_2stage(
'N',
'U', 1, 0, a, 1, x, z, 1, w,
1219 CALL chkxer(
'SSBEVD_2STAGE', infot, nout, lerr, ok )
1221 CALL ssbevd_2stage(
'N',
'U', 2, 0, a, 1, x, z, 1, w,
1223 CALL chkxer(
'SSBEVD_2STAGE', infot, nout, lerr, ok )
1229 CALL ssbevd_2stage(
'N',
'U', 1, 0, a, 1, x, z, 1, w,
1231 CALL chkxer(
'SSBEVD_2STAGE', infot, nout, lerr, ok )
1243 CALL ssbev(
'/',
'U', 0, 0, a, 1, x, z, 1, w, info )
1244 CALL chkxer(
'SSBEV ', infot, nout, lerr, ok )
1246 CALL ssbev(
'N',
'/', 0, 0, a, 1, x, z, 1, w, info )
1247 CALL chkxer(
'SSBEV ', infot, nout, lerr, ok )
1249 CALL ssbev(
'N',
'U', -1, 0, a, 1, x, z, 1, w, info )
1250 CALL chkxer(
'SSBEV ', infot, nout, lerr, ok )
1252 CALL ssbev(
'N',
'U', 0, -1, a, 1, x, z, 1, w, info )
1253 CALL chkxer(
'SSBEV ', infot, nout, lerr, ok )
1255 CALL ssbev(
'N',
'U', 2, 1, a, 1, x, z, 1, w, info )
1256 CALL chkxer(
'SSBEV ', infot, nout, lerr, ok )
1258 CALL ssbev(
'V',
'U', 2, 0, a, 1, x, z, 1, w, info )
1259 CALL chkxer(
'SSBEV ', infot, nout, lerr, ok )
1264 srnamt =
'SSBEV_2STAGE '
1266 CALL ssbev_2stage(
'/',
'U', 0, 0, a, 1, x, z, 1, w, 0, info )
1267 CALL chkxer(
'SSBEV_2STAGE ', infot, nout, lerr, ok )
1269 CALL ssbev_2stage(
'V',
'U', 0, 0, a, 1, x, z, 1, w, 0, info )
1270 CALL chkxer(
'SSBEV_2STAGE ', infot, nout, lerr, ok )
1272 CALL ssbev_2stage(
'N',
'/', 0, 0, a, 1, x, z, 1, w, 0, info )
1273 CALL chkxer(
'SSBEV_2STAGE ', infot, nout, lerr, ok )
1275 CALL ssbev_2stage(
'N',
'U', -1, 0, a, 1, x, z, 1, w, 0, info )
1276 CALL chkxer(
'SSBEV_2STAGE ', infot, nout, lerr, ok )
1278 CALL ssbev_2stage(
'N',
'U', 0, -1, a, 1, x, z, 1, w, 0, info )
1279 CALL chkxer(
'SSBEV_2STAGE ', infot, nout, lerr, ok )
1281 CALL ssbev_2stage(
'N',
'U', 2, 1, a, 1, x, z, 1, w, 0, info )
1282 CALL chkxer(
'SSBEV_2STAGE ', infot, nout, lerr, ok )
1284 CALL ssbev_2stage(
'N',
'U', 2, 0, a, 1, x, z, 0, w, 0, info )
1285 CALL chkxer(
'SSBEV_2STAGE ', infot, nout, lerr, ok )
1287 CALL ssbev_2stage(
'N',
'U', 0, 0, a, 1, x, z, 1, w, 0, info )
1288 CALL chkxer(
'SSBEV_2STAGE ', infot, nout, lerr, ok )
1295 CALL ssbevx(
'/',
'A',
'U', 0, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1296 $ 0.0, m, x, z, 1, w, iw, i3, info )
1297 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1299 CALL ssbevx(
'N',
'/',
'U', 0, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1300 $ 0.0, m, x, z, 1, w, iw, i3, info )
1301 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1303 CALL ssbevx(
'N',
'A',
'/', 0, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1304 $ 0.0, m, x, z, 1, w, iw, i3, info )
1305 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1307 CALL ssbevx(
'N',
'A',
'U', -1, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1308 $ 0.0, m, x, z, 1, w, iw, i3, info )
1309 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1311 CALL ssbevx(
'N',
'A',
'U', 0, -1, a, 1, q, 1, 0.0, 0.0, 0, 0,
1312 $ 0.0, m, x, z, 1, w, iw, i3, info )
1313 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1315 CALL ssbevx(
'N',
'A',
'U', 2, 1, a, 1, q, 1, 0.0, 0.0, 0, 0,
1316 $ 0.0, m, x, z, 1, w, iw, i3, info )
1317 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1319 CALL ssbevx(
'V',
'A',
'U', 2, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1320 $ 0.0, m, x, z, 2, w, iw, i3, info )
1321 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1323 CALL ssbevx(
'N',
'V',
'U', 1, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1324 $ 0.0, m, x, z, 1, w, iw, i3, info )
1325 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1327 CALL ssbevx(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1328 $ 0.0, m, x, z, 1, w, iw, i3, info )
1329 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1331 CALL ssbevx(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0, 0.0, 2, 1,
1332 $ 0.0, m, x, z, 1, w, iw, i3, info )
1333 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1335 CALL ssbevx(
'N',
'I',
'U', 2, 0, a, 1, q, 1, 0.0, 0.0, 2, 1,
1336 $ 0.0, m, x, z, 1, w, iw, i3, info )
1337 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1339 CALL ssbevx(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0, 0.0, 1, 2,
1340 $ 0.0, m, x, z, 1, w, iw, i3, info )
1341 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1343 CALL ssbevx(
'V',
'A',
'U', 2, 0, a, 1, q, 2, 0.0, 0.0, 0, 0,
1344 $ 0.0, m, x, z, 1, w, iw, i3, info )
1345 CALL chkxer(
'SSBEVX', infot, nout, lerr, ok )
1350 srnamt =
'SSBEVX_2STAGE'
1352 CALL ssbevx_2stage(
'/',
'A',
'U', 0, 0, a, 1, q, 1, 0.0e0,
1353 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1354 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1356 CALL ssbevx_2stage(
'V',
'A',
'U', 0, 0, a, 1, q, 1, 0.0e0,
1357 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1358 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1360 CALL ssbevx_2stage(
'N',
'/',
'U', 0, 0, a, 1, q, 1, 0.0e0,
1361 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1362 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1364 CALL ssbevx_2stage(
'N',
'A',
'/', 0, 0, a, 1, q, 1, 0.0e0,
1365 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1366 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1368 CALL ssbevx_2stage(
'N',
'A',
'U', -1, 0, a, 1, q, 1, 0.0e0,
1369 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1370 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1372 CALL ssbevx_2stage(
'N',
'A',
'U', 0, -1, a, 1, q, 1, 0.0e0,
1373 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1374 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1376 CALL ssbevx_2stage(
'N',
'A',
'U', 2, 1, a, 1, q, 1, 0.0e0,
1377 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1378 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1384 CALL ssbevx_2stage(
'N',
'V',
'U', 1, 0, a, 1, q, 1, 0.0e0,
1385 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1386 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1388 CALL ssbevx_2stage(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0e0,
1389 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1390 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1392 CALL ssbevx_2stage(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0e0,
1393 $ 0.0e0, 2, 1, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1394 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1396 CALL ssbevx_2stage(
'N',
'I',
'U', 2, 0, a, 1, q, 1, 0.0e0,
1397 $ 0.0e0, 2, 1, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1398 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1400 CALL ssbevx_2stage(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0e0,
1401 $ 0.0e0, 1, 2, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1402 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1408 CALL ssbevx_2stage(
'N',
'A',
'U', 0, 0, a, 1, q, 1, 0.0e0,
1409 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1410 CALL chkxer(
'SSBEVX_2STAGE', infot, nout, lerr, ok )
1418 WRITE( nout, fmt = 9999 )path, nt
1420 WRITE( nout, fmt = 9998 )path
1423 9999
FORMAT( 1x, a3,
' routines passed the tests of the error exits',
1424 $
' (', i3,
' tests done)' )
1425 9998
FORMAT(
' *** ', a3,
' routines failed the tests of the error ',
subroutine chkxer(srnamt, infot, nout, lerr, ok)
subroutine ssbev_2stage(jobz, uplo, n, kd, ab, ldab, w, z, ldz, work, lwork, info)
SSBEV_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER m...
subroutine ssbev(jobz, uplo, n, kd, ab, ldab, w, z, ldz, work, info)
SSBEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrices
subroutine ssbevd_2stage(jobz, uplo, n, kd, ab, ldab, w, z, ldz, work, lwork, iwork, liwork, info)
SSBEVD_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER ...
subroutine ssbevd(jobz, uplo, n, kd, ab, ldab, w, z, ldz, work, lwork, iwork, liwork, info)
SSBEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine ssbevx_2stage(jobz, range, uplo, n, kd, ab, ldab, q, ldq, vl, vu, il, iu, abstol, m, w, z, ldz, work, lwork, iwork, ifail, info)
SSBEVX_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER ...
subroutine ssbevx(jobz, range, uplo, n, kd, ab, ldab, q, ldq, vl, vu, il, iu, abstol, m, w, z, ldz, work, iwork, ifail, info)
SSBEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine ssbtrd(vect, uplo, n, kd, ab, ldab, d, e, q, ldq, work, info)
SSBTRD
subroutine ssyev_2stage(jobz, uplo, n, a, lda, w, work, lwork, info)
SSYEV_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matr...
subroutine ssyev(jobz, uplo, n, a, lda, w, work, lwork, info)
SSYEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
subroutine ssyevd_2stage(jobz, uplo, n, a, lda, w, work, lwork, iwork, liwork, info)
SSYEVD_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY mat...
subroutine ssyevd(jobz, uplo, n, a, lda, w, work, lwork, iwork, liwork, info)
SSYEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
subroutine ssyevr_2stage(jobz, range, uplo, n, a, lda, vl, vu, il, iu, abstol, m, w, z, ldz, isuppz, work, lwork, iwork, liwork, info)
SSYEVR_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY mat...
subroutine ssyevr(jobz, range, uplo, n, a, lda, vl, vu, il, iu, abstol, m, w, z, ldz, isuppz, work, lwork, iwork, liwork, info)
SSYEVR computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
subroutine ssyevx_2stage(jobz, range, uplo, n, a, lda, vl, vu, il, iu, abstol, m, w, z, ldz, work, lwork, iwork, ifail, info)
SSYEVX_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY mat...
subroutine ssyevx(jobz, range, uplo, n, a, lda, vl, vu, il, iu, abstol, m, w, z, ldz, work, lwork, iwork, ifail, info)
SSYEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
subroutine ssytd2(uplo, n, a, lda, d, e, tau, info)
SSYTD2 reduces a symmetric matrix to real symmetric tridiagonal form by an orthogonal similarity tran...
subroutine ssytrd_2stage(vect, uplo, n, a, lda, d, e, tau, hous2, lhous2, work, lwork, info)
SSYTRD_2STAGE
subroutine ssytrd_sb2st(stage1, vect, uplo, n, kd, ab, ldab, d, e, hous, lhous, work, lwork, info)
SSYTRD_SB2ST reduces a real symmetric band matrix A to real symmetric tridiagonal form T
subroutine ssytrd_sy2sb(uplo, n, kd, a, lda, ab, ldab, tau, work, lwork, info)
SSYTRD_SY2SB
subroutine ssytrd(uplo, n, a, lda, d, e, tau, work, lwork, info)
SSYTRD
subroutine sspev(jobz, uplo, n, ap, w, z, ldz, work, info)
SSPEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrices
subroutine sspevd(jobz, uplo, n, ap, w, z, ldz, work, lwork, iwork, liwork, info)
SSPEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine sspevx(jobz, range, uplo, n, ap, vl, vu, il, iu, abstol, m, w, z, ldz, work, iwork, ifail, info)
SSPEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine ssptrd(uplo, n, ap, d, e, tau, info)
SSPTRD
subroutine spteqr(compz, n, d, e, z, ldz, work, info)
SPTEQR
subroutine sstebz(range, order, n, vl, vu, il, iu, abstol, d, e, m, nsplit, w, iblock, isplit, work, iwork, info)
SSTEBZ
subroutine sstedc(compz, n, d, e, z, ldz, work, lwork, iwork, liwork, info)
SSTEDC
subroutine sstein(n, d, e, m, w, iblock, isplit, z, ldz, work, iwork, ifail, info)
SSTEIN
subroutine ssteqr(compz, n, d, e, z, ldz, work, info)
SSTEQR
subroutine ssterf(n, d, e, info)
SSTERF
subroutine sstev(jobz, n, d, e, z, ldz, work, info)
SSTEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrices
subroutine sstevd(jobz, n, d, e, z, ldz, work, lwork, iwork, liwork, info)
SSTEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine sstevr(jobz, range, n, d, e, vl, vu, il, iu, abstol, m, w, z, ldz, isuppz, work, lwork, iwork, liwork, info)
SSTEVR computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine sstevx(jobz, range, n, d, e, vl, vu, il, iu, abstol, m, w, z, ldz, work, iwork, ifail, info)
SSTEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine sorgtr(uplo, n, a, lda, tau, work, lwork, info)
SORGTR
subroutine sormtr(side, uplo, trans, m, n, a, lda, tau, c, ldc, work, lwork, info)
SORMTR
subroutine sopgtr(uplo, n, ap, tau, q, ldq, work, info)
SOPGTR
subroutine sopmtr(side, uplo, trans, m, n, ap, tau, c, ldc, work, info)
SOPMTR
subroutine serrst(path, nunit)
SERRST