79 INTEGER I, INFO, , M, N, NT
82 INTEGER I1( NMAX ), I2( NMAX ), I3( NMAX ), IW( LIW
83 REAL D( NMAX ), E( NMAX ), R( LW ), RW( LW ),
85 COMPLEX A( NMAX, NMAX ), C( NMAX, NMAX ),
86 $ Q( NMAX, NMAX ), TAU( NMAX ), W( LW ),
109 COMMON / infoc / infot, nout, ok, lerr
110 COMMON / srnamc / srnamt
118 WRITE( nout, fmt = * )
125 a( i, j ) = 1. / real( i+j )
140 IF( lsamen( 2, c2,
'ST' ) )
THEN
146 CALL chetrd( '/
', 0, A, 1, D, E, TAU, W, 1, INFO )
147 CALL CHKXER( 'chetrd', INFOT, NOUT, LERR, OK )
149 CALL CHETRD( 'u
', -1, A, 1, D, E, TAU, W, 1, INFO )
150 CALL CHKXER( 'chetrd', INFOT, NOUT, LERR, OK )
152 CALL CHETRD( 'u
', 2, A, 1, D, E, TAU, W, 1, INFO )
153 CALL CHKXER( 'chetrd', INFOT, NOUT, LERR, OK )
155 CALL CHETRD( 'u
', 0, A, 1, D, E, TAU, W, 0, INFO )
156 CALL CHKXER( 'chetrd', INFOT, NOUT, LERR, OK )
163 CALL CHETRD_2STAGE( '/
', 'u
', 0, A, 1, D, E, TAU,
167 CALL CHETRD_2STAGE( 'h
', 'u
', 0, A, 1, D, E, TAU,
171 CALL CHETRD_2STAGE( 'n
', '/
', 0, A, 1, D, E, TAU,
175 CALL CHETRD_2STAGE( 'n
', 'u
', -1, A, 1, D, E, TAU,
179 CALL CHETRD_2STAGE( 'n
', 'u
', 2, A, 1, D, E, TAU,
183 CALL CHETRD_2STAGE( 'n
', 'u
', 0, A, 1, D, E, TAU,
187 CALL CHETRD_2STAGE( 'n
', 'u
', 0, A, 1, D, E, TAU,
196 CALL CHETRD_HE2HB( '/
', 0, 0, A, 1, C, 1, TAU, W, 1, INFO )
199 CALL CHETRD_HE2HB( 'u
', -1, 0, A, 1, C, 1, TAU, W, 1, INFO )
202 CALL CHETRD_HE2HB( 'u
', 0, -1, A, 1, C, 1, TAU, W, 1, INFO )
205 CALL CHETRD_HE2HB( 'u
', 2, 0, A, 1, C, 1, TAU, W, 1, INFO )
208 CALL CHETRD_HE2HB( 'u
', 0, 2, A, 1, C, 1, TAU, W, 1, INFO )
211 CALL CHETRD_HE2HB( 'u
', 0, 0, A, 1, C, 1, TAU, W, 0, INFO )
221 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
225 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
229 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
233 CALL chkxer(
'CHETRD_HB2ST', infot, nout,
237 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
241 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
245 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
249 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
253 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
260 CALL cungtr(
'/', 0, a, 1, tau, w, 1, info )
261 CALL chkxer(
'CUNGTR', infot, nout, lerr, ok )
263 CALL cungtr(
'U', -1, a, 1, tau, w, 1, info )
264 CALL chkxer(
'CUNGTR', infot, nout, lerr, ok
266 CALL cungtr(
'U', 2, a, 1, tau, w, 1, info )
267 CALL chkxer(
'CUNGTR', infot, nout, lerr, ok )
269 CALL cungtr(
'U', 3, a, 3, tau, w, 1, info )
270 CALL chkxer(
'CUNGTR', infot, nout, lerr, ok )
277 CALL cunmtr(
'/',
'U',
'N', 0, 0, a, 1, tau, c, 1, w, 1, info )
278 CALL chkxer(
'CUNMTR', infot, nout, lerr, ok )
280 CALL cunmtr(
'L',
'/',
'N', 0, 0, a, 1, tau, c, 1, w, 1, info )
281 CALL chkxer(
'CUNMTR', infot, nout, lerr, ok )
283 CALL cunmtr(
'L',
'U',
'/', 0, 0, a, 1, tau, c, 1, w, 1, info )
284 CALL chkxer(
'CUNMTR', infot, nout, lerr, ok )
286 CALL cunmtr(
'L',
'U',
'N', -1, 0, a, 1, tau, c, 1, w, 1,
288 CALL chkxer(
'CUNMTR', infot, nout, lerr, ok )
290 CALL cunmtr(
'L',
'U',
'N', 0, -1, a, 1, tau, c, 1, w, 1,
292 CALL chkxer(
'CUNMTR', infot, nout, lerr, ok )
294 CALL cunmtr(
'L',
'U',
'N', 2, 0, a, 1, tau, c, 2, w, 1, info )
295 CALL chkxer(
'CUNMTR', infot, nout, lerr, ok )
297 CALL cunmtr(
'R',
'U',
'N', 0, 2, a, 1, tau, c, 1, w, 1, info )
298 CALL chkxer(
'CUNMTR', infot, nout, lerr, ok )
300 CALL cunmtr(
'L',
'U',
'N', 2, 0, a, 2, tau, c, 1, w, 1, info )
301 CALL chkxer(
'CUNMTR', infot, nout, lerr, ok )
303 CALL cunmtr(
'L',
'U',
'N', 0, 2, a, 1, tau, c, 1, w, 1, info )
304 CALL chkxer(
'CUNMTR', infot, nout, lerr, ok )
306 CALL cunmtr(
'R',
'U',
'N', 2, 0, a, 1, tau, c, 2, w, 1, info )
307 CALL chkxer(
'CUNMTR', infot, nout, lerr, ok )
314 CALL chptrd(
'/', 0, a, d, e, tau, info )
315 CALL chkxer(
'CHPTRD', infot, nout, lerr, ok )
317 CALL chptrd(
'U', -1, a, d, e, tau, info )
318 CALL chkxer(
'CHPTRD', infot, nout, lerr, ok )
325 CALL cupgtr(
'/', 0, a, tau, z, 1, w, info )
326 CALL chkxer(
'CUPGTR', infot, nout, lerr, ok
328 CALL cupgtr(
'U', -1, a, tau, z, 1, w, info )
329 CALL chkxer(
'CUPGTR', infot, nout, lerr, ok )
331 CALL cupgtr(
'U', 2, a, tau, z, 1, w, info )
332 CALL chkxer(
'CUPGTR', infot, nout, lerr, ok )
339 CALL cupmtr(
'/',
'U',
'N', 0, 0, a, tau, c, 1, w, info )
340 CALL chkxer(
'CUPMTR', infot, nout, lerr, ok )
342 CALL cupmtr(
'L',
'/',
'N', 0, 0, a, tau, c, 1, w, info )
343 CALL chkxer(
'CUPMTR', infot, nout, lerr, ok )
345 CALL cupmtr(
'L',
'U',
'/', 0, 0, a, tau, c, 1, w, info )
346 CALL chkxer(
'CUPMTR', infot, nout, lerr, ok )
348 CALL cupmtr(
'L',
'U',
'N', -1, 0, a, tau, c, 1, w, info )
349 CALL chkxer(
'CUPMTR', infot, nout, lerr, ok )
351 CALL cupmtr(
'L',
'U',
'N', 0, -1, a, tau, c, 1, w, info )
354 CALL CUPMTR( 'l
', 'u
', 'n', 2, 0, a, tau, c, 1, w, info )
355 CALL chkxer(
'CUPMTR', infot, nout, lerr, ok )
362 CALL cpteqr(
'/', 0, d, e, z, 1, rw, info )
363 CALL chkxer(
'CPTEQR', infot, nout, lerr, ok )
365 CALL cpteqr(
'N', -1, d, e, z, 1, rw, info )
366 CALL chkxer(
'CPTEQR', infot, nout, lerr, ok )
368 CALL cpteqr(
'V', 2, d, e, z, 1, rw, info )
369 CALL chkxer(
'CPTEQR', infot, nout, lerr, ok )
376 CALL cstein( -1, d, e, 0, x, i1, i2, z, 1, rw, iw, i3, info )
379 CALL cstein( 0, d, e, -1, x, i1, i2, z, 1, rw, iw, i3, info )
380 CALL chkxer(
'CSTEIN', infot, nout, lerr, ok )
382 CALL cstein( 0, d, e, 1, x, i1, i2, z, 1, rw, iw, i3, info )
383 CALL chkxer(
'CSTEIN', infot, nout, lerr, ok )
385 CALL cstein( 2, d, e, 0, x, i1, i2, z, 1, rw, iw, i3, info )
386 CALL chkxer(
'CSTEIN', infot, nout, lerr, ok )
393 CALL csteqr(
'/', 0, d, e, z, 1, rw, info )
394 CALL chkxer(
'CSTEQR', infot, nout, lerr, ok )
396 CALL csteqr(
'N', -1, d, e, z, 1, rw, info )
397 CALL chkxer(
'CSTEQR', infot, nout, lerr, ok )
399 CALL csteqr(
'V', 2, d, e, z, 1, rw, info )
400 CALL chkxer(
'CSTEQR', infot, nout, lerr, ok )
407 CALL cstedc(
'/', 0, d, e, z, 1, w, 1, rw, 1, iw, 1, info )
408 CALL chkxer(
'CSTEDC', infot, nout, lerr, ok )
410 CALL cstedc(
'N', -1, d, e, z, 1, w, 1, rw, 1, iw, 1, info )
411 CALL chkxer(
'CSTEDC', infot, nout, lerr, ok )
413 CALL cstedc(
'V', 2, d, e, z, 1, w, 4, rw, 23, iw, 28, info )
414 CALL chkxer(
'CSTEDC', infot, nout, lerr, ok )
416 CALL cstedc(
'N', 2, d, e, z, 1, w, 0, rw, 1, iw, 1, info )
417 CALL chkxer(
'CSTEDC', infot, nout, lerr, ok )
419 CALL cstedc(
'V', 2, d, e, z, 2, w, 0, rw, 23, iw, 28, info )
420 CALL chkxer(
'CSTEDC', infot, nout, lerr, ok )
422 CALL cstedc(
'N', 2, d, e, z, 1, w, 1, rw, 0, iw, 1, info )
423 CALL chkxer(
'CSTEDC', infot, nout, lerr, ok )
425 CALL cstedc(
'I', 2, d, e, z, 2, w, 1, rw, 1, iw, 12, info )
426 CALL chkxer(
'CSTEDC', infot, nout, lerr, ok )
428 CALL cstedc(
'V', 2, d, e, z, 2, w, 4, rw, 1, iw, 28, info )
429 CALL chkxer(
'CSTEDC', infot, nout, lerr, ok )
431 CALL cstedc(
'N', 2, d, e, z, 1, w, 1, rw, 1, iw, 0, info )
432 CALL chkxer(
'CSTEDC', infot, nout, lerr, ok )
434 CALL cstedc(
'I', 2, d, e, z, 2, w, 1, rw, 23, iw, 0, info )
437 CALL CSTEDC( 'v
', 2, D, E, Z, 2, W, 4, RW, 23, IW, 0, INFO )
438 CALL CHKXER( 'cstedc', INFOT, NOUT, LERR, OK )
445 CALL CHEEVD( '/
', 'u', 0, a, 1, x, w, 1, rw, 1, iw, 1, info )
446 CALL chkxer(
'CHEEVD', infot, nout, lerr, ok )
448 CALL cheevd(
'N',
'/', 0, a, 1, x, w, 1, rw, 1, iw, 1, info )
449 CALL chkxer(
'CHEEVD', infot, nout, lerr, ok )
451 CALL cheevd(
'N',
'U', -1, a, 1, x, w, 1, rw, 1, iw, 1, info )
452 CALL chkxer(
'CHEEVD', infot, nout, lerr, ok )
454 CALL cheevd(
'N',
'U', 2, a, 1, x, w, 3, rw, 2, iw, 1, info )
455 CALL chkxer(
'CHEEVD', infot, nout, lerr, ok )
457 CALL cheevd(
'N',
'U', 1, a, 1, x, w, 0, rw, 1, iw, 1, info )
458 CALL chkxer(
'CHEEVD', infot, nout, lerr, ok )
460 CALL cheevd(
'N',
'U', 2, a, 2, x, w, 2, rw, 2, iw, 1, info )
461 CALL chkxer(
'CHEEVD', infot, nout, lerr, ok )
463 CALL cheevd(
'V',
'U', 2, a, 2, x, w, 3, rw, 25, iw, 12, info )
464 CALL chkxer(
'CHEEVD', infot, nout, lerr, ok )
466 CALL cheevd( 'n
', 'u
', 1, A, 1, X, W, 1, RW, 0, IW, 1, INFO )
467 CALL CHKXER( 'cheevd', INFOT, NOUT, LERR, OK )
469 CALL CHEEVD( 'n
', 'u
', 2, A, 2, X, W, 3, RW, 1, IW, 1, INFO )
470 CALL CHKXER( 'cheevd', INFOT, NOUT, LERR, OK )
472 CALL CHEEVD( 'v
', 'u
', 2, A, 2, X, W, 8, RW, 18, IW, 12, INFO )
473 CALL CHKXER( 'cheevd', INFOT, NOUT, LERR, OK )
475 CALL CHEEVD( 'n
', 'u
', 1, A, 1, X, W, 1, RW, 1, IW, 0, INFO )
476 CALL CHKXER( 'cheevd', INFOT, NOUT, LERR, OK )
478 CALL CHEEVD( 'v
', 'u
', 2, A, 2, X, W, 8, RW, 25, IW, 11, INFO )
479 CALL CHKXER( 'cheevd', INFOT, NOUT, LERR, OK )
486 CALL CHEEVD_2STAGE( '/
', 'u
', 0, A, 1, X, W, 1,
487 $ RW, 1, IW, 1, INFO )
490 CALL CHEEVD_2STAGE( 'v
', 'u', 0, a, 1, x, w, 1,
491 $ rw, 1, iw, 1, info )
492 CALL chkxer(
'CHEEVD_2STAGE', infot, nout, lerr, ok )
495 $ rw, 1, iw, 1, info )
496 CALL chkxer(
'CHEEVD_2STAGE', infot, nout, lerr, ok )
499 $ rw, 1, iw, 1, info )
500 CALL chkxer(
'CHEEVD_2STAGE', infot, nout, lerr, ok )
503 $ rw, 2, iw, 1, info )
504 CALL chkxer(
'CHEEVD_2STAGE', infot, nout
508 CALL chkxer(
'CHEEVD_2STAGE', infot, nout, lerr, ok )
511 $ rw, 2, iw, 1, info )
512 CALL chkxer(
'CHEEVD_2STAGE', infot, nout, lerr, ok )
519 $ rw, 0, iw, 1, info )
520 CALL chkxer(
'CHEEVD_2STAGE', infot, nout, lerr, ok )
523 $ rw, 1, iw, 1, info )
524 CALL chkxer(
'CHEEVD_2STAGE', infot, nout, lerr, ok )
531 $ rw, 1, iw, 0, info )
532 CALL chkxer(
'CHEEVD_2STAGE', infot, nout, lerr, ok )
543 CALL cheev(
'/',
'U', 0, a, 1, x, w, 1, rw, info )
544 CALL chkxer(
'CHEEV ', infot, nout, lerr, ok )
546 CALL cheev(
'N',
'/', 0, a, 1, x, w, 1, rw, info )
547 CALL chkxer(
'CHEEV ', infot, nout, lerr, ok )
549 CALL cheev(
'N',
'U', -1, a, 1, x, w, 1, rw, info )
550 CALL chkxer(
'CHEEV ', infot, nout, lerr, ok )
552 CALL cheev(
'N',
'U', 2, a, 1, x, w, 3, rw, info )
553 CALL chkxer(
'CHEEV ', infot, nout, lerr, ok )
555 CALL cheev(
'N',
'U', 2, a, 2, x, w, 2, rw, info )
556 CALL chkxer(
'CHEEV ', infot, nout, lerr, ok )
561 srnamt =
'CHEEV_2STAGE '
563 CALL cheev_2stage(
'/',
'U', 0, a, 1, x, w, 1, rw, info )
564 CALL chkxer(
'CHEEV_2STAGE ', infot, nout, lerr, ok )
566 CALL cheev_2stage(
'V',
'U', 0, a, 1, x, w, 1, rw, info )
567 CALL chkxer(
'CHEEV_2STAGE ', infot, nout, lerr, ok )
569 CALL cheev_2stage(
'N',
'/', 0, a, 1, x, w, 1, rw, info )
570 CALL chkxer(
'CHEEV_2STAGE ', infot, nout, lerr, ok )
572 CALL cheev_2stage(
'N',
'U', -1, a, 1, x, w, 1, rw, info )
573 CALL chkxer(
'CHEEV_2STAGE ', infot, nout, lerr, ok )
575 CALL cheev_2stage(
'N',
'U', 2, a, 1, x, w, 3, rw, info )
576 CALL chkxer(
'CHEEV_2STAGE ', infot, nout, lerr, ok )
578 CALL cheev_2stage(
'N',
'U', 2, a, 2, x, w, 2, rw, info )
579 CALL chkxer(
'CHEEV_2STAGE ', infot, nout, lerr, ok )
586 CALL cheevx(
'/',
'A',
'U', 0, a, 1, 0.0, 0.0, 0, 0, 0.0, m, x,
587 $ z, 1, w, 1, rw, iw, i3, info )
588 CALL chkxer(
'CHEEVX', infot, nout, lerr, ok )
590 CALL cheevx(
'V',
'/',
'U', 0, a, 1, 0.0, 1.0, 1, 0, 0.0, m, x,
591 $ z, 1, w, 1, rw, iw, i3, info )
592 CALL chkxer(
'CHEEVX', infot, nout, lerr, ok )
594 CALL cheevx(
'V''A''/', 0, a, 1, 0.0, 0.0, 0, 0, 0.0, m, x,
595 $ z, 1, w, 1, rw, iw, i3, info )
597 CALL cheevx(
'V',
'A',
'U', -1, a, 1, 0.0, 0.0, 0, 0, 0.0, m,
598 $ x, z, 1, w, 1, rw, iw, i3, info )
599 CALL chkxer(
'CHEEVX', infot, nout, lerr, ok )
601 CALL cheevx(
'V',
'A',
'U', 2, a, 1, 0.0, 0.0, 0, 0, 0.0, m, x,
602 $ z, 2, w, 3, rw, iw, i3, info )
605 CALL CHEEVX( 'v
', 'v
', 'u
', 1, A, 1, 0.0, 0.0, 0, 0, 0.0, M, X,
606 $ Z, 1, W, 1, RW, IW, I3, INFO )
607 CALL CHKXER( 'cheevx', INFOT, NOUT, LERR, OK )
609 CALL CHEEVX( 'v
', 'i
', 'u
', 1, A, 1, 0.0, 0.0, 0, 0, 0.0, M, X,
610 $ Z, 1, W, 1, RW, IW, I3, INFO )
611 CALL CHKXER( 'cheevx', INFOT, NOUT, LERR, OK )
613 CALL CHEEVX( 'v
', 'i
', 'u
', 2, A, 2, 0.0, 0.0, 2, 1, 0.0, M, X,
614 $ Z, 2, W, 3, RW, IW, I3, INFO )
615 CALL CHKXER( 'cheevx', INFOT, NOUT, LERR, OK )
617 CALL CHEEVX( 'v
', 'a
', 'u
', 2, A, 2, 0.0, 0.0, 0, 0, 0.0, M, X,
618 $ Z, 1, W, 3, RW, IW, I3, INFO )
619 CALL CHKXER( 'cheevx', INFOT, NOUT, LERR, OK )
621 CALL CHEEVX( 'v
', 'a
', 'u
', 2, A, 2, 0.0, 0.0, 0, 0, 0.0, M, X,
622 $ Z, 2, W, 2, RW, IW, I1, INFO )
623 CALL CHKXER( 'cheevx', INFOT, NOUT, LERR, OK )
630 CALL CHEEVX_2STAGE( '/
', 'a
', 'u
', 0, A, 1,
631 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
632 $ M, X, Z, 1, W, 1, RW, IW, I3, INFO )
635 CALL CHEEVX_2STAGE( 'v
', 'a
', 'u
', 0, A, 1,
636 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
637 $ M, X, Z, 1, W, 1, RW, IW, I3, INFO )
640 CALL CHEEVX_2STAGE( 'n
', '/
', 'u
', 0, A, 1,
641 $ 0.0D0, 1.0D0, 1, 0, 0.0D0,
642 $ M, X, Z, 1, W, 1, RW, IW, I3, INFO )
645 CALL CHEEVX_2STAGE( 'n
', 'a
', '/
', 0, A, 1,
646 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
647 $ M, X, Z, 1, W, 1, RW, IW, I3, INFO )
649 CALL CHEEVX_2STAGE( 'n
', 'a
', 'u
', -1, A, 1,
650 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
651 $ M, X, Z, 1, W, 1, RW, IW, I3, INFO )
654 CALL CHEEVX_2STAGE( 'n
', 'a
', 'u
', 2, A, 1,
655 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
656 $ M, X, Z, 2, W, 3, RW, IW, I3, INFO )
659 CALL CHEEVX_2STAGE( 'n
', 'v
', 'u
', 1, A, 1,
660 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
661 $ M, X, Z, 1, W, 1, RW, IW, I3, INFO )
664 CALL CHEEVX_2STAGE( 'n
', 'i
', 'u
', 1, A, 1,
665 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
666 $ M, X, Z, 1, W, 1, RW, IW, I3, INFO )
669 CALL CHEEVX_2STAGE( 'n
', 'i
', 'u
', 2, A, 2,
670 $ 0.0D0, 0.0D0, 2, 1, 0.0D0,
671 $ M, X, Z, 2, W, 3, RW, IW, I3, INFO )
674 CALL CHEEVX_2STAGE( 'n',
'A',
'U', 2, a, 2,
675 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
676 $ m, x, z, 0, w, 3, rw, iw, i3, info )
677 CALL chkxer(
'CHEEVX_2STAGE', infot, nout, lerr, ok )
680 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
681 $ m, x, z, 2, w, 0, rw, iw, i1, info )
682 CALL chkxer(
'CHEEVX_2STAGE', infot, nout, lerr, ok )
690 CALL CHEEVR( '/
', 'a
', 'u
', 0, A, 1, 0.0, 0.0, 1, 1, 0.0, M, R,
691 $ Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), 10*N,
693 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
695 CALL CHEEVR( 'v
', '/
', 'u
', 0, A, 1, 0.0, 0.0, 1, 1, 0.0, M, R,
696 $ Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), 10*N,
698 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
700 CALL CHEEVR( 'v
', 'a
', '/
', -1, A, 1, 0.0, 0.0, 1, 1, 0.0, M,
701 $ R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), 10*N,
703 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
705 CALL CHEEVR( 'v
', 'a
', 'u
', -1, A, 1, 0.0, 0.0, 1, 1, 0.0, M,
706 $ R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), 10*N,
708 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
710 CALL CHEEVR( 'v
', 'a
', 'u
', 2, A, 1, 0.0, 0.0, 1, 1, 0.0, M, R,
711 $ Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), 10*N,
713 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
715 CALL CHEEVR( 'v
', 'v
', 'u
', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
716 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
718 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
720 CALL CHEEVR( 'v
', 'i
', 'u
', 1, A, 1, 0.0E0, 0.0E0, 0, 1, 0.0,
721 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
723 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
726 CALL CHEEVR( 'v
', 'i
', 'u
', 2, A, 2, 0.0E0, 0.0E0, 2, 1, 0.0,
727 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
729 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
731 CALL CHEEVR( 'v
', 'i
', 'u
', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
732 $ M, R, Z, 0, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
734 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
736 CALL CHEEVR( 'v
', 'i
', 'u
', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
737 $ M, R, Z, 1, IW, Q, 2*N-1, RW, 24*N, IW( 2*N+1 ),
739 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
741 CALL CHEEVR( 'v
', 'i
', 'u
', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
742 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N-1, IW( 2*N-1 ),
744 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
746 CALL CHEEVR( 'v
', 'i
', 'u
', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
747 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW, 10*N-1,
749 CALL CHKXER( 'cheevr', INFOT, NOUT, LERR, OK )
757 CALL CHEEVR_2STAGE( '/
', 'a
', 'u
', 0, A, 1,
758 $ 0.0D0, 0.0D0, 1, 1, 0.0D0,
759 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
763 CALL CHEEVR_2STAGE( 'v
', 'a
', 'u
', 0, A, 1,
764 $ 0.0D0, 0.0D0, 1, 1, 0.0D0,
765 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
769 CALL CHEEVR_2STAGE( 'n
', '/
', 'u
', 0, A, 1,
770 $ 0.0D0, 0.0D0, 1, 1, 0.0D0,
771 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
775 CALL CHEEVR_2STAGE( 'n
', 'a
', '/
', -1, A, 1,
776 $ 0.0D0, 0.0D0, 1, 1, 0.0D0,
777 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N,
778 $ IW( 2*N+1 ), 10*N, INFO )
781 CALL CHEEVR_2STAGE( 'n
', 'a
', 'u
', -1, A, 1,
782 $ 0.0D0, 0.0D0, 1, 1, 0.0D0,
783 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N,
784 $ IW( 2*N+1 ), 10*N, INFO )
787 CALL CHEEVR_2STAGE( 'n
', 'a
', 'u
', 2, A, 1,
788 $ 0.0D0, 0.0D0, 1, 1, 0.0D0,
789 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
793 CALL CHEEVR_2STAGE( 'n
', 'v
', 'u
', 1, A, 1,
794 $ 0.0D0, 0.0D0, 1, 1, 0.0D0,
795 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
799 CALL CHEEVR_2STAGE( 'n
', 'i
', 'u
', 1, A, 1,
800 $ 0.0D0, 0.0D0, 0, 1, 0.0D0,
801 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
805 CALL CHEEVR_2STAGE( 'n
', 'i
', 'u
', 2, A, 2,
806 $ 0.0D0, 0.0D0, 2, 1, 0.0D0,
807 $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
811 CALL CHEEVR_2STAGE( 'n
', 'i
', 'u
', 1, A, 1,
812 $ 0.0D0, 0.0D0, 1, 1, 0.0D0,
813 $ M, R, Z, 0, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ),
817 CALL CHEEVR_2STAGE( 'n
', 'i
', 'u
', 1, A, 1,
818 $ 0.0D0, 0.0D0, 1, 1, 0.0D0,
819 $ M, R, Z, 1, IW, Q, 2*N-1, RW, 24*N, IW( 2*N+1 ),
823 CALL CHEEVR_2STAGE( 'n
', 'i
', 'u
', 1, A, 1,
824 $ 0.0D0, 0.0D0, 1, 1, 0.0D0,
825 $ M, R, Z, 1, IW, Q, 26*N, RW, 24*N-1, IW( 2*N-1 ),
829 CALL CHEEVR_2STAGE( 'n
', 'i
', 'u
', 1, A, 1,
830 $ 0.0D0, 0.0D0, 1, 1, 0.0D0,
831 $ M, R, Z, 1, IW, Q, 26*N, RW, 24*N, IW, 10*N-1,
840 CALL chpevd(
'/',
'U', 0, a, x, z, 1, w, 1, rw, 1, iw, 1,
842 CALL chkxer(
'CHPEVD', infot, nout, lerr, ok )
844 CALL chpevd(
'N',
'/', 0, a, x, z, 1, w, 1, rw, 1, iw, 1,
846 CALL chkxer(
'CHPEVD', infot, nout, lerr, ok )
848 CALL chpevd(
'N',
'U', -1, a, x, z, 1, w, 1, rw, 1, iw, 1,
850 CALL chkxer(
'CHPEVD', infot, nout, lerr, ok )
852 CALL chpevd(
'V',
'U', 2, a, x, z, 1, w, 4, rw, 25, iw, 12,
854 CALL chkxer(
'CHPEVD', infot, nout, lerr, ok )
856 CALL chpevd( 'n
', 'u
', 1, A, X, Z, 1, W, 0, RW, 1, IW, 1,
858 CALL CHKXER( 'chpevd', INFOT, NOUT, LERR, OK )
860 CALL CHPEVD( 'n
', 'u
', 2, A, X, Z, 2, W, 1, RW, 2, IW, 1,
862 CALL CHKXER( 'chpevd', INFOT, NOUT, LERR, OK )
864 CALL CHPEVD( 'v',
'U', 2, a, x, z, 2, w, 2, rw, 25, iw, 12,
868 CALL CHPEVD( 'n
', 'u
', 1, A, X, Z, 1, W, 1, RW, 0, IW, 1,
870 CALL CHKXER( 'chpevd', INFOT, NOUT, LERR, OK )
872 CALL CHPEVD( 'n
', 'u
', 2, A, X, Z, 2, W, 2, RW, 1, IW, 1,
874 CALL CHKXER( 'chpevd', INFOT, NOUT, LERR, OK )
876 CALL CHPEVD( 'v
', 'u
', 2, A, X, Z, 2, W, 4, RW, 18, IW, 12,
878 CALL CHKXER( 'chpevd', INFOT, NOUT, LERR, OK )
880 CALL CHPEVD( 'n
', 'u
', 1, A, X, Z, 1, W, 1, RW, 1, IW, 0,
882 CALL CHKXER( 'chpevd', INFOT, NOUT, LERR, OK )
884 CALL CHPEVD( 'n
', 'u
', 2, A, X, Z, 2, W, 2, RW, 2, IW, 0,
886 CALL CHKXER( 'chpevd', INFOT, NOUT, LERR, OK )
888 CALL CHPEVD( 'v
', 'u
', 2, A, X, Z, 2, W, 4, RW, 25, IW, 2,
890 CALL CHKXER( 'chpevd', INFOT, NOUT, LERR, OK )
897 CALL CHPEV( '/
', 'u
', 0, A, X, Z, 1, W, RW, INFO )
898 CALL CHKXER( 'chpev ', INFOT, NOUT, LERR, OK )
900 CALL CHPEV( 'n
', '/
', 0, A, X, Z, 1, W, RW, INFO )
901 CALL CHKXER( 'chpev ', INFOT, NOUT, LERR, OK )
903 CALL CHPEV( 'n
', 'u
', -1, A, X, Z, 1, W, RW, INFO )
904 CALL CHKXER( 'chpev ', INFOT, NOUT, LERR, OK )
906 CALL CHPEV( 'v
', 'u
', 2, A, X, Z, 1, W, RW, INFO )
907 CALL CHKXER( 'chpev ', INFOT, NOUT, LERR, OK )
914 CALL CHPEVX( '/
', 'a
', 'u
', 0, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
915 $ 1, W, RW, IW, I3, INFO )
916 CALL CHKXER( 'chpevx', INFOT, NOUT, LERR, OK )
918 CALL CHPEVX( 'v
', '/
', 'u
', 0, A, 0.0, 1.0, 1, 0, 0.0, M, X, Z,
919 $ 1, W, RW, IW, I3, INFO )
920 CALL CHKXER( 'chpevx', INFOT, NOUT, LERR, OK )
922 CALL CHPEVX( 'v
', 'a',
'/', 0, a, 0.0, 0.0, 0, 0, 0.0, m, x, z,
923 $ 1, w, rw, iw, i3, info )
924 CALL chkxer(
'CHPEVX', infot, nout, lerr, ok )
926 CALL chpevx(
'V',
'A',
'U', -1, a, 0.0, 0.0, 0, 0, 0.0, m, x,
927 $ z, 1, w, rw, iw, i3, info )
928 CALL chkxer(
'CHPEVX', infot, nout, lerr, ok )
930 CALL chpevx(
'V',
'V',
'U', 1, a, 0.0, 0.0, 0, 0, 0.0, m, x, z,
931 $ 1, w, rw, iw, i3, info )
932 CALL chkxer(
'CHPEVX', infot, nout, lerr, ok )
934 CALL chpevx(
'V',
'I',
'U', 1, a, 0.0, 0.0, 0, 0, 0.0, m, x, z,
935 $ 1, w, rw, iw, i3, info )
936 CALL chkxer(
'CHPEVX', infot, nout, lerr, ok )
938 CALL chpevx(
'V',
'I',
'U', 2, a, 0.0, 0.0, 2, 1, 0.0, m, x, z,
939 $ 2, w, rw, iw, i3, info )
940 CALL chkxer(
'CHPEVX', infot, nout, lerr, ok )
942 CALL chpevx(
'V',
'A',
'U', 2, a, 0.0, 0.0, 0, 0, 0.0, m, x, z,
943 $ 1, w, rw, iw, i3, info )
944 CALL chkxer(
'CHPEVX', infot, nout, lerr, ok )
949 ELSE IF( lsamen( 2, c2,
'HB' ) )
THEN
955 CALL chbtrd(
'/',
'U', 0, 0, a, 1, d, e, z, 1, w, info )
956 CALL chkxer(
'CHBTRD', infot, nout, lerr, ok )
958 CALL chbtrd(
'N',
'/', 0, 0, a, 1, d, e, z, 1, w, info )
959 CALL chkxer(
'CHBTRD', infot, nout, lerr, ok )
961 CALL chbtrd(
'N',
'U', -1, 0, a, 1, d, e, z, 1, w, info )
962 CALL chkxer(
'CHBTRD', infot, nout, lerr, ok )
964 CALL chbtrd(
'N',
'U', 0, -1, a, 1, d, e, z, 1, w, info )
965 CALL chkxer(
'CHBTRD', infot, nout, lerr, ok )
967 CALL chbtrd(
'N',
'U', 1, 1, a, 1, d, e, z, 1, w, info )
968 CALL chkxer(
'CHBTRD', infot, nout, lerr, ok )
970 CALL chbtrd(
'V',
'U', 2, 0, a, 1, d, e, z, 1, w, info )
971 CALL chkxer(
'CHBTRD', infot, nout, lerr, ok )
976 srnamt =
'CHETRD_HB2ST'
982 CALL CHETRD_HB2ST( 'n
', '/
', 'u
', 0, 0, A, 1, D, E,
986 CALL CHETRD_HB2ST( 'n
', 'h
', 'u
', 0, 0, A, 1, D, E,
990 CALL CHETRD_HB2ST( 'n
', 'n
', '/
', 0, 0, A, 1, D, E,
994 CALL CHETRD_HB2ST( 'n
', 'n
', 'u
', -1, 0, A, 1, D, E,
998 CALL CHETRD_HB2ST( 'n',
'N',
'U', 0, -1, a, 1, d, e,
1002 CALL CHETRD_HB2ST( 'n
', 'n
', 'u
', 0, 1, A, 1, D, E,
1003 $ C, 1, W, 1, INFO )
1006 CALL CHETRD_HB2ST( 'n
', 'n
', 'u
', 0, 0, A, 1, D, E,
1007 $ C, 0, W, 1, INFO )
1010 CALL CHETRD_HB2ST( 'n
', 'n
', 'u
', 0, 0, A, 1, D, E,
1011 $ C, 1, W, 0, INFO )
1019 CALL CHBEVD( '/
', 'u
', 0, 0, A, 1, X, Z, 1, W, 1, RW, 1, IW, 1,
1021 CALL CHKXER( 'chbevd', INFOT, NOUT, LERR, OK )
1023 CALL CHBEVD( 'n
', '/
', 0, 0, A, 1, X, Z, 1, W, 1, RW, 1, IW, 1,
1025 CALL CHKXER( 'chbevd', INFOT, NOUT, LERR, OK )
1027 CALL CHBEVD( 'n
', 'u
', -1, 0, A, 1, X, Z, 1, W, 1, RW, 1, IW,
1029 CALL CHKXER( 'chbevd', INFOT, NOUT, LERR, OK )
1031 CALL CHBEVD( 'n
', 'u
', 0, -1, A, 1, X, Z, 1, W, 1, RW, 1, IW,
1033 CALL CHKXER( 'chbevd', INFOT, NOUT, LERR, OK )
1035 CALL CHBEVD( 'n
', 'u
', 2, 1, A, 1, X, Z, 1, W, 2, RW, 2, IW, 1,
1037 CALL CHKXER( 'chbevd', INFOT, NOUT, LERR, OK )
1039 CALL CHBEVD( 'v
', 'u
', 2, 1, A, 2, X, Z, 1, W, 8, RW, 25, IW,
1041 CALL CHKXER( 'chbevd', INFOT, NOUT, LERR, OK )
1043 CALL CHBEVD( 'n
', 'u
', 1, 0, A, 1, X, Z, 1, W, 0, RW, 1, IW, 1,
1045 CALL CHKXER( 'chbevd', INFOT, NOUT, LERR, OK )
1047 CALL CHBEVD( 'n
', 'u
', 2, 1, A, 2, X, Z, 2, W, 1, RW, 2, IW, 1,
1049 CALL CHKXER( 'chbevd', INFOT, NOUT, LERR, OK )
1051 CALL CHBEVD( 'v
', 'u
', 2, 1, A, 2, X, Z, 2, W, 2, RW, 25, IW,
1053 CALL CHKXER( 'chbevd', INFOT, NOUT, LERR, OK )
1055 CALL CHBEVD( 'n
', 'u
', 1, 0, A, 1, X, Z, 1, W, 1, RW, 0, IW, 1,
1057 CALL CHKXER( 'chbevd', infot, nout, lerr, ok )
1059 CALL chbevd(
'N',
'U', 2, 1, a, 2, x, z, 2, w, 2, rw, 1, iw, 1,
1061 CALL chkxer(
'CHBEVD', infot, nout, lerr, ok )
1063 CALL chbevd(
'V',
'U', 2, 1, a, 2, x, z, 2, w, 8, rw, 2, iw,
1065 CALL chkxer(
'CHBEVD', infot, nout, lerr, ok )
1067 CALL chbevd(
'N',
'U', 1, 0, a, 1, x, z, 1, w, 1, rw, 1, iw, 0,
1069 CALL chkxer(
'CHBEVD', infot, nout, lerr, ok )
1071 CALL chbevd(
'N', 'u
', 2, 1, A, 2, X, Z, 2, W, 2, RW, 2, IW, 0,
1073 CALL CHKXER( 'chbevd', INFOT, NOUT, LERR, OK )
1075 CALL CHBEVD( 'v
', 'u
', 2, 1, A, 2, X, Z, 2, W, 8, RW, 25, IW,
1077 CALL CHKXER( 'chbevd', INFOT, NOUT, LERR, OK )
1084 CALL CHBEVD_2STAGE( '/
', 'u
', 0, 0, A, 1, X, Z, 1,
1085 $ W, 1, RW, 1, IW, 1, INFO )
1088 CALL CHBEVD_2STAGE( 'v
', 'u
', 0, 0, A, 1, X, Z, 1,
1089 $ W, 1, RW, 1, IW, 1, INFO )
1092 CALL CHBEVD_2STAGE( 'n
', '/
', 0, 0, A, 1, X, Z, 1,
1093 $ W, 1, RW, 1, IW, 1, INFO )
1096 CALL CHBEVD_2STAGE( 'n
', 'u
', -1, 0, A, 1, X, Z, 1,
1097 $ W, 1, RW, 1, IW, 1, INFO )
1100 CALL CHBEVD_2STAGE( 'n
', 'u
', 0, -1, A, 1, X, Z, 1,
1101 $ W, 1, RW, 1, IW, 1, INFO )
1105 $ W, 2, RW, 2, IW, 1, INFO )
1108 CALL CHBEVD_2STAGE( 'n
', 'u
', 2, 1, A, 2, X, Z, 0,
1109 $ W, 8, RW, 25, IW, 12, INFO )
1112 CALL CHBEVD_2STAGE( 'n
', 'u
', 1, 0, A, 1, X, Z, 1,
1113 $ W, 0, RW, 1, IW, 1, INFO )
1116 CALL CHBEVD_2STAGE( 'n
', 'u
', 2, 1, A, 2, X, Z, 2,
1117 $ W, 1, RW, 2, IW, 1, INFO )
1124 CALL CHBEVD_2STAGE( 'n
', 'u
', 1, 0, A, 1, X, Z, 1,
1125 $ W, 1, RW, 0, IW, 1, INFO )
1128 CALL CHBEVD_2STAGE( 'n
', 'u
', 2, 1, A, 2, X, Z, 2,
1129 $ W, 25, RW, 1, IW, 1, INFO )
1136 CALL CHBEVD_2STAGE( 'n
', 'u
', 1, 0, A, 1, X, Z, 1,
1137 $ W, 1, RW, 1, IW, 0, INFO )
1140 CALL CHBEVD_2STAGE( 'n
', 'u
', 2, 1, A, 2, X, Z, 2,
1141 $ W, 25, RW, 2, IW, 0, INFO )
1153 CALL CHBEV( '/
', 'u
', 0, 0, A, 1, X, Z, 1, W, RW, INFO )
1154 CALL CHKXER( 'chbev ', INFOT, NOUT, LERR, OK )
1156 CALL CHBEV( 'n
', '/
', 0, 0, A, 1, X, Z, 1, W, RW, INFO )
1157 CALL CHKXER( 'chbev ', INFOT, NOUT, LERR, OK )
1159 CALL CHBEV( 'n
', 'u
', -1, 0, A, 1, X, Z, 1, W, RW, INFO )
1160 CALL CHKXER( 'chbev ', INFOT, NOUT, LERR, OK )
1162 CALL CHBEV( 'n
', 'u
', 0, -1, A, 1, X, Z, 1, W, RW, INFO )
1163 CALL CHKXER( 'chbev ', INFOT, NOUT, LERR, OK )
1165 CALL CHBEV( 'n
', 'u
', 2, 1, A, 1, X, Z, 1, W, RW, INFO )
1166 CALL CHKXER( 'chbev ', INFOT, NOUT, LERR, OK )
1168 CALL CHBEV( 'v
', 'u
', 2, 0, A, 1, X, Z, 1, W, RW, INFO )
1169 CALL CHKXER( 'chbev ', INFOT, NOUT, LERR, OK )
1176 CALL CHBEV_2STAGE( '/
', 'u
', 0, 0, A, 1, X,
1177 $ Z, 1, W, 0, RW, INFO )
1180 CALL CHBEV_2STAGE( 'v
', 'u
', 0, 0, A, 1, X,
1181 $ Z, 1, W, 0, RW, INFO )
1184 CALL CHBEV_2STAGE( 'n
', '/
', 0, 0, A, 1, X,
1185 $ Z, 1, W, 0, RW, INFO )
1188 CALL CHBEV_2STAGE( 'n
', 'u
', -1, 0, A, 1, X,
1189 $ Z, 1, W, 0, RW, INFO )
1192 CALL CHBEV_2STAGE( 'n
', 'u
', 0, -1, A, 1, X,
1193 $ Z, 1, W, 0, RW, INFO )
1196 CALL CHBEV_2STAGE( 'n
', 'u
', 2, 1, A, 1, X,
1197 $ Z, 1, W, 0, RW, INFO )
1200 CALL CHBEV_2STAGE( 'n
', 'u
', 2, 0, A, 1, X,
1201 $ Z, 0, W, 0, RW, INFO )
1204 CALL CHBEV_2STAGE( 'n
', 'u
', 2, 0, A, 1, X,
1205 $ Z, 1, W, 0, RW, INFO )
1213 CALL CHBEVX( '/
', 'a
', 'u
', 0, 0, A, 1, Q, 1, 0.0, 0.0, 0, 0,
1214 $ 0.0, M, X, Z, 1, W, RW, IW, I3, INFO )
1215 CALL CHKXER( 'chbevx', infot, nout, lerr, ok )
1217 CALL chbevx(
'V',
'/',
'U', 0, 0, a, 1, q, 1, 0.0, 1.0, 1, 0,
1218 $ 0.0, m, x, z, 1, w, rw, iw, i3, info )
1219 CALL chkxer(
'CHBEVX', infot, nout, lerr, ok )
1221 CALL chbevx(
'V',
'A',
'/', 0, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1222 $ 0.0, m, x, z, 1, w, rw, iw, i3, info )
1224 CALL chbevx(
'V',
'A',
'U', -1, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1225 $ 0.0, m, x, z, 1, w, rw, iw, i3, info )
1226 CALL chkxer(
'CHBEVX', infot, nout, lerr, ok )
1228 CALL chbevx(
'V',
'A',
'U', 0, -1, a, 1, q, 1, 0.0, 0.0, 0, 0,
1229 $ 0.0, m, x, z, 1, w, rw, iw, i3, info )
1230 CALL chkxer(
'CHBEVX', infot, nout, lerr, ok )
1232 CALL chbevx(
'V',
'A',
'U', 2, 1, a, 1, q, 2, 0.0, 0.0, 0, 0,
1233 $ 0.0, m, x, z, 2, w, rw, iw, i3, info )
1234 CALL chkxer(
'CHBEVX', infot, nout, lerr, ok )
1236 CALL chbevx(
'V',
'A',
'U', 2, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1237 $ 0.0, m, x, z, 2, w, rw, iw, i3, info )
1238 CALL chkxer(
'CHBEVX', infot, nout, lerr, ok
1240 CALL chbevx(
'V',
'V',
'U', 1, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1241 $ 0.0, m, x, z, 1, w, rw, iw, i3, info )
1242 CALL chkxer(
'CHBEVX', infot, nout, lerr, ok )
1244 CALL chbevx(
'V',
'I',
'U', 1, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1245 $ 0.0, m, x, z, 1, w, rw, iw, i3, info )
1246 CALL chkxer(
'CHBEVX', infot, nout, lerr, ok )
1248 CALL chbevx(
'V',
'I',
'U', 1, 0, a, 1, q, 1, 0.0, 0.0, 1, 2,
1249 $ 0.0, m, x, z, 1, w, rw, iw, i3, info )
1250 CALL chkxer(
'CHBEVX', infot, nout, lerr, ok )
1252 CALL chbevx(
'V',
'A', 'u', 2, 0, a, 1, q, 2, 0.0, 0.0, 0, 0,
1253 $ 0.0, m, x, z, 1, w, rw, iw, i3, info )
1254 CALL chkxer(
'CHBEVX', infot, nout, lerr, ok )
1261 CALL CHBEVX_2STAGE( '/
', 'a
', 'u
', 0, 0, A, 1, Q, 1,
1262 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
1263 $ M, X, Z, 1, W, 0, RW, IW, I3, INFO )
1265 CALL CHBEVX_2STAGE( 'v
', 'a
', 'u
', 0, 0, A, 1, Q, 1,
1266 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
1267 $ M, X, Z, 1, W, 0, RW, IW, I3, INFO )
1270 CALL CHBEVX_2STAGE( 'n
', '/
', 'u
', 0, 0, A, 1, Q, 1,
1271 $ 0.0D0, 1.0D0, 1, 0, 0.0D0,
1272 $ M, X, Z, 1, W, 0, RW, IW, I3, INFO )
1275 CALL CHBEVX_2STAGE( 'n
', 'a
', '/
', 0, 0, A, 1, Q, 1,
1276 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
1277 $ M, X, Z, 1, W, 0, RW, IW, I3, INFO )
1279 CALL CHBEVX_2STAGE( 'n
', 'a
', 'u
', -1, 0, A, 1, Q, 1,
1280 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
1281 $ M, X, Z, 1, W, 0, RW, IW, I3, INFO )
1284 CALL CHBEVX_2STAGE( 'n
', 'a
', 'u
', 0, -1, A, 1, Q, 1,
1285 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
1286 $ M, X, Z, 1, W, 0, RW, IW, I3, INFO )
1289 CALL CHBEVX_2STAGE( 'n
', 'a
', 'u
', 2, 1, A, 1, Q, 2,
1290 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
1291 $ M, X, Z, 2, W, 0, RW, IW, I3, INFO )
1299 CALL CHBEVX_2STAGE( 'n
', 'v
', 'u
', 1, 0, A, 1, Q, 1,
1300 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
1301 $ M, X, Z, 1, W, 0, RW, IW, I3, INFO )
1304 CALL CHBEVX_2STAGE( 'n
', 'i
', 'u
', 1, 0, A, 1, Q, 1,
1305 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
1306 $ M, X, Z, 1, W, 0, RW, IW, I3, INFO )
1309 CALL CHBEVX_2STAGE( 'n
', 'i
', 'u
', 1, 0, A, 1, Q, 1,
1310 $ 0.0D0, 0.0D0, 1, 2, 0.0D0,
1311 $ M, X, Z, 1, W, 0, RW, IW, I3, INFO )
1314 CALL CHBEVX_2STAGE( 'n
', 'a
', 'u
', 2, 0, A, 1, Q, 2,
1315 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
1316 $ M, X, Z, 0, W, 0, RW, IW, I3, INFO )
1319 CALL CHBEVX_2STAGE( 'n
', 'a
', 'u
', 2, 0, A, 1, Q, 2,
1320 $ 0.0D0, 0.0D0, 0, 0, 0.0D0,
1321 $ M, X, Z, 1, W, 0, RW, IW, I3, INFO )
1329 WRITE( NOUT, FMT = 9999 )PATH, NT
1331 WRITE( NOUT, FMT = 9998 )PATH
1334 9999 FORMAT( 1X, A3, ' routines passed
the tests of
the error exits
',
1335 $ ' (
', I3, ' tests done)
' )
1336 9998 FORMAT( ' ***
', A3, ' routines failed
the tests of
the error
',