75 parameter( nmax = 3, liw = 12*nmax, lw = 20*nmax )
79 INTEGER I, INFO, J, 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 )
161 srnamt =
'CHETRD_2STAGE'
165 CALL chkxer(
'CHETRD_2STAGE', infot, nout, lerr, ok )
169 CALL chkxer(
'CHETRD_2STAGE', infot, nout, lerr, ok )
173 CALL chkxer(
'CHETRD_2STAGE', infot, nout, lerr, ok )
177 CALL chkxer(
'CHETRD_2STAGE', infot, nout, lerr, ok )
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 )
219 CALL CHETRD_HB2ST( '/
', 'n
', 'u
', 0, 0, A, 1, D, E,
223 CALL CHETRD_HB2ST( 'y
', '/
', 'u
', 0, 0, A, 1, D, E,
227 CALL CHETRD_HB2ST( 'y
', 'h
', 'u
', 0, 0, A, 1, D, E,
231 CALL CHETRD_HB2ST( 'y
', 'n
', '/', 0, 0, a, 1, d, e,
233 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
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
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 )
352 CALL chkxer(
'CUPMTR', infot, nout, lerr, ok )
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 )
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 )
377 CALL CHKXER( 'cstein', INFOT, NOUT, LERR, OK )
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 )
435 CALL CHKXER( 'cstedc', INFOT, NOUT, LERR, OK )
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, lerr, ok )
507 $ rw, 1, iw, 1, info )
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 )
575 CALL CHEEV_2STAGE( 'n
', 'u
', 2, A, 1, X, W, 3, RW, INFO )
578 CALL CHEEV_2STAGE( 'n
', 'u
', 2, A, 2, X, W, 2, RW, INFO )
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 )
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 )
603 CALL CHKXER( 'cheevx', INFOT, NOUT, LERR, OK )
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 )
628 srnamt =
'CHEEVX_2STAGE'
632 $ m, x, z, 1, w, 1, rw, iw, i3, info )
633 CALL chkxer(
'CHEEVX_2STAGE', infot, nout, lerr, ok )
636 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
637 $ m, x, z, 1, w, 1, rw, iw, i3, info )
638 CALL chkxer(
'CHEEVX_2STAGE', infot, nout, lerr, ok )
641 $ 0.0d0, 1.0d0, 1, 0, 0.0d0,
642 $ m, x, z, 1, w, 1, rw, iw, i3, info )
643 CALL chkxer(
'CHEEVX_2STAGE', infot, nout, lerr, ok )
646 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
647 $ m, x, z, 1, w, 1, rw, iw, i3, info )
650 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
651 $ m, x, z, 1, w, 1, rw, iw, i3, info )
652 CALL chkxer(
'CHEEVX_2STAGE', infot, nout, lerr, ok )
655 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
656 $ m, x, z, 2, w, 3, rw, iw, i3, info )
657 CALL chkxer(
'CHEEVX_2STAGE', infot, nout, lerr, ok )
660 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
661 $ m, x, z, 1, w, 1, rw, iw, i3, info )
662 CALL chkxer(
'CHEEVX_2STAGE', infot, nout, lerr, ok )
665 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
666 $ m, x, z, 1, w, 1, rw, iw, i3, info )
667 CALL chkxer(
'CHEEVX_2STAGE', infot, nout, lerr, ok )
670 $ 0.0d0, 0.0d0, 2, 1, 0.0d0,
671 $ m, x, z, 2, w, 3, rw, iw, i3, info )
672 CALL chkxer(
'CHEEVX_2STAGE', infot, nout, lerr, ok )
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 ),
761 CALL chkxer(
'CHEEVR_2STAGE', infot, nout, lerr, ok )
764 $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
765 $ m, r, z, 1, iw, q, 2*n, rw, 24*n, iw( 2*n+
767 CALL chkxer(
'CHEEVR_2STAGE', infot, nout, lerr, ok )
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 ),
773 CALL chkxer(
'CHEEVR_2STAGE', infot, nout, lerr, ok )
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,
866 CALL CHKXER( 'chpevd', INFOT, NOUT, LERR, OK )
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'
980 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
984 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
988 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
992 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
996 CALL chkxer(
'CHETRD_HB2ST', infot, nout, lerr, ok )
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 )
1104 CALL CHBEVD_2STAGE( 'n
', 'u
', 2, 1, A, 1, X, Z, 1,
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
',