| @@ -364,7 +364,7 @@ | |||
| *> \verbatim | |||
| *> D1 is REAL array of | |||
| *> dimension( max(NN) ) | |||
| *> The eigenvalues of A, as computed by CSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by CSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> \endverbatim | |||
| @@ -665,8 +665,7 @@ | |||
| EXTERNAL CCOPY, CHET21, CHETRD, CHPT21, CHPTRD, CLACPY, | |||
| $ CLASET, CLATMR, CLATMS, CPTEQR, CSTEDC, CSTEMR, | |||
| $ CSTEIN, CSTEQR, CSTT21, CSTT22, CUNGTR, CUPGTR, | |||
| $ SCOPY, SLABAD, SLASUM, SSTEBZ, SSTECH, SSTERF, | |||
| $ XERBLA | |||
| $ SCOPY, SLASUM, SSTEBZ, SSTECH, SSTERF, XERBLA | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| INTRINSIC ABS, CONJG, INT, LOG, MAX, MIN, REAL, SQRT | |||
| @@ -733,7 +732,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = ONE / UNFL | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| LOG2UI = INT( LOG( ULPINV ) / LOG( TWO ) ) | |||
| @@ -385,7 +385,7 @@ | |||
| *> \verbatim | |||
| *> D1 is REAL array of | |||
| *> dimension( max(NN) ) | |||
| *> The eigenvalues of A, as computed by CSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by CSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> \endverbatim | |||
| @@ -683,10 +683,10 @@ | |||
| EXTERNAL ILAENV, SLAMCH, SLARND, SSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL SCOPY, SLABAD, SLASUM, SSTEBZ, SSTECH, SSTERF, | |||
| $ XERBLA, CCOPY, CHET21, CHETRD, CHPT21, CHPTRD, | |||
| $ CLACPY, CLASET, CLATMR, CLATMS, CPTEQR, CSTEDC, | |||
| $ CSTEMR, CSTEIN, CSTEQR, CSTT21, CSTT22, CUNGTR, | |||
| EXTERNAL SCOPY, SLASUM, SSTEBZ, SSTECH, SSTERF, XERBLA, | |||
| $ CCOPY, CHET21, CHETRD, CHPT21, CHPTRD, CLACPY, | |||
| $ CLASET, CLATMR, CLATMS, CPTEQR, CSTEDC, CSTEMR, | |||
| $ CSTEIN, CSTEQR, CSTT21, CSTT22, CUNGTR, | |||
| $ CUPGTR, CHETRD_2STAGE, SLASET | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| @@ -754,7 +754,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = ONE / UNFL | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| LOG2UI = INT( LOG( ULPINV ) / LOG( TWO ) ) | |||
| @@ -395,7 +395,7 @@ | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, CGESVD, CGET51, CGGESX, CLACPY, CLAKF2, | |||
| $ CLASET, CLATM5, SLABAD, XERBLA | |||
| $ CLASET, CLATM5, XERBLA | |||
| * .. | |||
| * .. Scalars in Common .. | |||
| LOGICAL FS | |||
| @@ -478,7 +478,6 @@ | |||
| ULPINV = ONE / ULP | |||
| SMLNUM = SLAMCH( 'S' ) / ULP | |||
| BIGNUM = ONE / SMLNUM | |||
| CALL SLABAD( SMLNUM, BIGNUM ) | |||
| THRSH2 = TEN*THRESH | |||
| NTESTT = 0 | |||
| NERRS = 0 | |||
| @@ -917,7 +916,7 @@ | |||
| $ / ' 2: A and B are upper triangular matrices, ', | |||
| $ / ' 3: A and B are as type 2, but each second diagonal ', | |||
| $ 'block in A_11 and ', / | |||
| $ ' each third diaongal block in A_22 are 2x2 blocks,', | |||
| $ ' each third diagonal block in A_22 are 2x2 blocks,', | |||
| $ / ' 4: A and B are block diagonal matrices, ', | |||
| $ / ' 5: (A,B) has potentially close or common ', | |||
| $ 'eigenvalues.', / ) | |||
| @@ -236,7 +236,7 @@ | |||
| *> | |||
| *> B COMPLEX array, dimension (LDB , max(NN)) | |||
| *> Used to hold the Hermitian positive definite matrix for | |||
| *> the generailzed problem. | |||
| *> the generalized problem. | |||
| *> On exit, B contains the last matrix actually | |||
| *> used. | |||
| *> Modified. | |||
| @@ -420,7 +420,7 @@ | |||
| * .. External Subroutines .. | |||
| EXTERNAL CHBGV, CHBGVD, CHBGVX, CHEGV, CHEGVD, CHEGVX, | |||
| $ CHPGV, CHPGVD, CHPGVX, CLACPY, CLASET, CLATMR, | |||
| $ CLATMS, CSGT01, SLABAD, SLAFTS, SLASUM, XERBLA | |||
| $ CLATMS, CSGT01, SLAFTS, SLASUM, XERBLA | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| INTRINSIC ABS, MAX, MIN, REAL, SQRT | |||
| @@ -481,7 +481,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = SLAMCH( 'Overflow' ) | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -242,7 +242,7 @@ | |||
| *> | |||
| *> B COMPLEX array, dimension (LDB , max(NN)) | |||
| *> Used to hold the Hermitian positive definite matrix for | |||
| *> the generailzed problem. | |||
| *> the generalized problem. | |||
| *> On exit, B contains the last matrix actually | |||
| *> used. | |||
| *> Modified. | |||
| @@ -426,7 +426,7 @@ | |||
| EXTERNAL LSAME, SLAMCH, SLARND | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL SLABAD, SLAFTS, SLASUM, XERBLA, CHBGV, CHBGVD, | |||
| EXTERNAL SLAFTS, SLASUM, XERBLA, CHBGV, CHBGVD, | |||
| $ CHBGVX, CHEGV, CHEGVD, CHEGVX, CHPGV, CHPGVD, | |||
| $ CHPGVX, CLACPY, CLASET, CLATMR, CLATMS, CSGT01, | |||
| $ CHEGV_2STAGE | |||
| @@ -490,7 +490,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = SLAMCH( 'Overflow' ) | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -204,7 +204,7 @@ | |||
| *> Not modified. | |||
| *> | |||
| *> D1 REAL array, dimension (max(NN)) | |||
| *> The eigenvalues of A, as computed by CSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by CSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> Modified. | |||
| @@ -393,8 +393,8 @@ | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, CHBEV, CHBEVD, CHBEVX, CHEEV, CHEEVD, | |||
| $ CHEEVR, CHEEVX, CHET21, CHET22, CHPEV, CHPEVD, | |||
| $ CHPEVX, CLACPY, CLASET, CLATMR, CLATMS, SLABAD, | |||
| $ SLAFTS, XERBLA | |||
| $ CHPEVX, CLACPY, CLASET, CLATMR, CLATMS, SLAFTS, | |||
| $ XERBLA | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| INTRINSIC ABS, INT, LOG, MAX, MIN, REAL, SQRT | |||
| @@ -451,7 +451,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = SLAMCH( 'Overflow' ) | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -204,7 +204,7 @@ | |||
| *> Not modified. | |||
| *> | |||
| *> D1 REAL array, dimension (max(NN)) | |||
| *> The eigenvalues of A, as computed by CSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by CSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> Modified. | |||
| @@ -391,7 +391,7 @@ | |||
| EXTERNAL SLAMCH, SLARND, SSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, SLABAD, SLAFTS, XERBLA, CHBEV, CHBEVD, | |||
| EXTERNAL ALASVM, SLAFTS, XERBLA, CHBEV, CHBEVD, | |||
| $ CHBEVX, CHEEV, CHEEVD, CHEEVR, CHEEVX, CHET21, | |||
| $ CHET22, CHPEV, CHPEVD, CHPEVX, CLACPY, CLASET, | |||
| $ CHEEVD_2STAGE, CHEEVR_2STAGE, CHEEVX_2STAGE, | |||
| @@ -453,7 +453,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = SLAMCH( 'Overflow' ) | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -628,56 +628,56 @@ | |||
| SRNAMT = 'CHEEVX_2STAGE' | |||
| INFOT = 1 | |||
| CALL CHEEVX_2STAGE( '/', 'A', 'U', 0, A, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 1, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHEEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 1 | |||
| CALL CHEEVX_2STAGE( 'V', 'A', 'U', 0, A, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 1, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHEEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 2 | |||
| CALL CHEEVX_2STAGE( 'N', '/', 'U', 0, A, 1, | |||
| $ 0.0D0, 1.0D0, 1, 0, 0.0D0, | |||
| $ 0.0, 1.0, 1, 0, 0.0, | |||
| $ M, X, Z, 1, W, 1, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHEEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 3 | |||
| CALL CHEEVX_2STAGE( 'N', 'A', '/', 0, A, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 1, RW, IW, I3, INFO ) | |||
| INFOT = 4 | |||
| CALL CHEEVX_2STAGE( 'N', 'A', 'U', -1, A, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 1, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHEEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 6 | |||
| CALL CHEEVX_2STAGE( 'N', 'A', 'U', 2, A, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 2, W, 3, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHEEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 8 | |||
| CALL CHEEVX_2STAGE( 'N', 'V', 'U', 1, A, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 1, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHEEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 9 | |||
| CALL CHEEVX_2STAGE( 'N', 'I', 'U', 1, A, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 1, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHEEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 10 | |||
| CALL CHEEVX_2STAGE( 'N', 'I', 'U', 2, A, 2, | |||
| $ 0.0D0, 0.0D0, 2, 1, 0.0D0, | |||
| $ 0.0, 0.0, 2, 1, 0.0, | |||
| $ M, X, Z, 2, W, 3, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHEEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 15 | |||
| CALL CHEEVX_2STAGE( 'N', 'A', 'U', 2, A, 2, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 0, W, 3, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHEEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 17 | |||
| CALL CHEEVX_2STAGE( 'N', 'A', 'U', 2, A, 2, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 2, W, 0, RW, IW, I1, INFO ) | |||
| CALL CHKXER( 'CHEEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| NT = NT + 11 | |||
| @@ -755,79 +755,79 @@ | |||
| N = 1 | |||
| INFOT = 1 | |||
| CALL CHEEVR_2STAGE( '/', 'A', 'U', 0, A, 1, | |||
| $ 0.0D0, 0.0D0, 1, 1, 0.0D0, | |||
| $ 0.0, 0.0, 1, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), | |||
| $ 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 1 | |||
| CALL CHEEVR_2STAGE( 'V', 'A', 'U', 0, A, 1, | |||
| $ 0.0D0, 0.0D0, 1, 1, 0.0D0, | |||
| $ 0.0, 0.0, 1, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), | |||
| $ 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 2 | |||
| CALL CHEEVR_2STAGE( 'N', '/', 'U', 0, A, 1, | |||
| $ 0.0D0, 0.0D0, 1, 1, 0.0D0, | |||
| $ 0.0, 0.0, 1, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), | |||
| $ 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 3 | |||
| CALL CHEEVR_2STAGE( 'N', 'A', '/', -1, A, 1, | |||
| $ 0.0D0, 0.0D0, 1, 1, 0.0D0, | |||
| $ 0.0, 0.0, 1, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, | |||
| $ IW( 2*N+1 ), 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 4 | |||
| CALL CHEEVR_2STAGE( 'N', 'A', 'U', -1, A, 1, | |||
| $ 0.0D0, 0.0D0, 1, 1, 0.0D0, | |||
| $ 0.0, 0.0, 1, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, | |||
| $ IW( 2*N+1 ), 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 6 | |||
| CALL CHEEVR_2STAGE( 'N', 'A', 'U', 2, A, 1, | |||
| $ 0.0D0, 0.0D0, 1, 1, 0.0D0, | |||
| $ 0.0, 0.0, 1, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), | |||
| $ 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 8 | |||
| CALL CHEEVR_2STAGE( 'N', 'V', 'U', 1, A, 1, | |||
| $ 0.0D0, 0.0D0, 1, 1, 0.0D0, | |||
| $ 0.0, 0.0, 1, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), | |||
| $ 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 9 | |||
| CALL CHEEVR_2STAGE( 'N', 'I', 'U', 1, A, 1, | |||
| $ 0.0D0, 0.0D0, 0, 1, 0.0D0, | |||
| $ 0.0, 0.0, 0, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), | |||
| $ 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 10 | |||
| CALL CHEEVR_2STAGE( 'N', 'I', 'U', 2, A, 2, | |||
| $ 0.0D0, 0.0D0, 2, 1, 0.0D0, | |||
| $ 0.0, 0.0, 2, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), | |||
| $ 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 15 | |||
| CALL CHEEVR_2STAGE( 'N', 'I', 'U', 1, A, 1, | |||
| $ 0.0D0, 0.0D0, 1, 1, 0.0D0, | |||
| $ 0.0, 0.0, 1, 1, 0.0, | |||
| $ M, R, Z, 0, IW, Q, 2*N, RW, 24*N, IW( 2*N+1 ), | |||
| $ 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 18 | |||
| CALL CHEEVR_2STAGE( 'N', 'I', 'U', 1, A, 1, | |||
| $ 0.0D0, 0.0D0, 1, 1, 0.0D0, | |||
| $ 0.0, 0.0, 1, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 2*N-1, RW, 24*N, IW( 2*N+1 ), | |||
| $ 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 20 | |||
| CALL CHEEVR_2STAGE( 'N', 'I', 'U', 1, A, 1, | |||
| $ 0.0D0, 0.0D0, 1, 1, 0.0D0, | |||
| $ 0.0, 0.0, 1, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 26*N, RW, 24*N-1, IW( 2*N-1 ), | |||
| $ 10*N, INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 22 | |||
| CALL CHEEVR_2STAGE( 'N', 'I', 'U', 1, A, 1, | |||
| $ 0.0D0, 0.0D0, 1, 1, 0.0D0, | |||
| $ 0.0, 0.0, 1, 1, 0.0, | |||
| $ M, R, Z, 1, IW, Q, 26*N, RW, 24*N, IW, 10*N-1, | |||
| $ INFO ) | |||
| CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| @@ -1259,65 +1259,65 @@ | |||
| SRNAMT = 'CHBEVX_2STAGE' | |||
| INFOT = 1 | |||
| CALL CHBEVX_2STAGE( '/', 'A', 'U', 0, 0, A, 1, Q, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 0, RW, IW, I3, INFO ) | |||
| INFOT = 1 | |||
| CALL CHBEVX_2STAGE( 'V', 'A', 'U', 0, 0, A, 1, Q, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 0, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHBEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 2 | |||
| CALL CHBEVX_2STAGE( 'N', '/', 'U', 0, 0, A, 1, Q, 1, | |||
| $ 0.0D0, 1.0D0, 1, 0, 0.0D0, | |||
| $ 0.0, 1.0, 1, 0, 0.0, | |||
| $ M, X, Z, 1, W, 0, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHBEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 3 | |||
| CALL CHBEVX_2STAGE( 'N', 'A', '/', 0, 0, A, 1, Q, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 0, RW, IW, I3, INFO ) | |||
| INFOT = 4 | |||
| CALL CHBEVX_2STAGE( 'N', 'A', 'U', -1, 0, A, 1, Q, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 0, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHBEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 5 | |||
| CALL CHBEVX_2STAGE( 'N', 'A', 'U', 0, -1, A, 1, Q, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 0, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHBEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 7 | |||
| CALL CHBEVX_2STAGE( 'N', 'A', 'U', 2, 1, A, 1, Q, 2, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 2, W, 0, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHBEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| * INFOT = 9 | |||
| * CALL CHBEVX_2STAGE( 'V', 'A', 'U', 2, 0, A, 1, Q, 1, | |||
| * $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| * $ 0.0, 0.0, 0, 0, 0.0, | |||
| * $ M, X, Z, 2, W, 0, RW, IW, I3, INFO ) | |||
| * CALL CHKXER( 'CHBEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 11 | |||
| CALL CHBEVX_2STAGE( 'N', 'V', 'U', 1, 0, A, 1, Q, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 0, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHBEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 12 | |||
| CALL CHBEVX_2STAGE( 'N', 'I', 'U', 1, 0, A, 1, Q, 1, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 0, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHBEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 13 | |||
| CALL CHBEVX_2STAGE( 'N', 'I', 'U', 1, 0, A, 1, Q, 1, | |||
| $ 0.0D0, 0.0D0, 1, 2, 0.0D0, | |||
| $ 0.0, 0.0, 1, 2, 0.0, | |||
| $ M, X, Z, 1, W, 0, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHBEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 18 | |||
| CALL CHBEVX_2STAGE( 'N', 'A', 'U', 2, 0, A, 1, Q, 2, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 0, W, 0, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHBEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| INFOT = 20 | |||
| CALL CHBEVX_2STAGE( 'N', 'A', 'U', 2, 0, A, 1, Q, 2, | |||
| $ 0.0D0, 0.0D0, 0, 0, 0.0D0, | |||
| $ 0.0, 0.0, 0, 0, 0.0, | |||
| $ M, X, Z, 1, W, 0, RW, IW, I3, INFO ) | |||
| CALL CHKXER( 'CHBEVX_2STAGE', INFOT, NOUT, LERR, OK ) | |||
| NT = NT + 12 | |||
| @@ -363,7 +363,7 @@ | |||
| *> \verbatim | |||
| *> D1 is DOUBLE PRECISION array of | |||
| *> dimension( max(NN) ) | |||
| *> The eigenvalues of A, as computed by DSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by DSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> \endverbatim | |||
| @@ -645,10 +645,10 @@ | |||
| EXTERNAL ILAENV, DLAMCH, DLARND, DSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL DCOPY, DLABAD, DLACPY, DLASET, DLASUM, DLATMR, | |||
| $ DLATMS, DOPGTR, DORGTR, DPTEQR, DSPT21, DSPTRD, | |||
| $ DSTEBZ, DSTECH, DSTEDC, DSTEMR, DSTEIN, DSTEQR, | |||
| $ DSTERF, DSTT21, DSTT22, DSYT21, DSYTRD, XERBLA | |||
| EXTERNAL DCOPY, DLACPY, DLASET, DLASUM, DLATMR, DLATMS, | |||
| $ DOPGTR, DORGTR, DPTEQR, DSPT21, DSPTRD, DSTEBZ, | |||
| $ DSTECH, DSTEDC, DSTEMR, DSTEIN, DSTEQR, DSTERF, | |||
| $ DSTT21, DSTT22, DSYT21, DSYTRD, XERBLA | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| INTRINSIC ABS, DBLE, INT, LOG, MAX, MIN, SQRT | |||
| @@ -715,7 +715,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = ONE / UNFL | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| LOG2UI = INT( LOG( ULPINV ) / LOG( TWO ) ) | |||
| @@ -384,7 +384,7 @@ | |||
| *> \verbatim | |||
| *> D1 is DOUBLE PRECISION array of | |||
| *> dimension( max(NN) ) | |||
| *> The eigenvalues of A, as computed by DSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by DSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> \endverbatim | |||
| @@ -666,10 +666,10 @@ | |||
| EXTERNAL ILAENV, DLAMCH, DLARND, DSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL DCOPY, DLABAD, DLACPY, DLASET, DLASUM, DLATMR, | |||
| $ DLATMS, DOPGTR, DORGTR, DPTEQR, DSPT21, DSPTRD, | |||
| $ DSTEBZ, DSTECH, DSTEDC, DSTEMR, DSTEIN, DSTEQR, | |||
| $ DSTERF, DSTT21, DSTT22, DSYT21, DSYTRD, XERBLA, | |||
| EXTERNAL DCOPY, DLACPY, DLASET, DLASUM, DLATMR, DLATMS, | |||
| $ DOPGTR, DORGTR, DPTEQR, DSPT21, DSPTRD, DSTEBZ, | |||
| $ DSTECH, DSTEDC, DSTEMR, DSTEIN, DSTEQR, DSTERF, | |||
| $ DSTT21, DSTT22, DSYT21, DSYTRD, XERBLA, | |||
| $ DSYTRD_2STAGE | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| @@ -737,7 +737,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = ONE / UNFL | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| LOG2UI = INT( LOG( ULPINV ) / LOG( TWO ) ) | |||
| @@ -400,7 +400,7 @@ | |||
| EXTERNAL DLCTSX, ILAENV, DLAMCH, DLANGE | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, DGESVD, DGET51, DGET53, DGGESX, DLABAD, | |||
| EXTERNAL ALASVM, DGESVD, DGET51, DGET53, DGGESX, | |||
| $ DLACPY, DLAKF2, DLASET, DLATM5, XERBLA | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| @@ -478,7 +478,6 @@ | |||
| ULPINV = ONE / ULP | |||
| SMLNUM = DLAMCH( 'S' ) / ULP | |||
| BIGNUM = ONE / SMLNUM | |||
| CALL DLABAD( SMLNUM, BIGNUM ) | |||
| THRSH2 = TEN*THRESH | |||
| NTESTT = 0 | |||
| NERRS = 0 | |||
| @@ -984,7 +983,7 @@ | |||
| $ / ' 2: A and B are upper triangular matrices, ', | |||
| $ / ' 3: A and B are as type 2, but each second diagonal ', | |||
| $ 'block in A_11 and ', / | |||
| $ ' each third diaongal block in A_22 are 2x2 blocks,', | |||
| $ ' each third diagonal block in A_22 are 2x2 blocks,', | |||
| $ / ' 4: A and B are block diagonal matrices, ', | |||
| $ / ' 5: (A,B) has potentially close or common ', | |||
| $ 'eigenvalues.', / ) | |||
| @@ -234,7 +234,7 @@ | |||
| *> | |||
| *> B DOUBLE PRECISION array, dimension (LDB , max(NN)) | |||
| *> Used to hold the symmetric positive definite matrix for | |||
| *> the generailzed problem. | |||
| *> the generalized problem. | |||
| *> On exit, B contains the last matrix actually | |||
| *> used. | |||
| *> Modified. | |||
| @@ -399,7 +399,7 @@ | |||
| EXTERNAL LSAME, DLAMCH, DLARND | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL DLABAD, DLACPY, DLAFTS, DLASET, DLASUM, DLATMR, | |||
| EXTERNAL DLACPY, DLAFTS, DLASET, DLASUM, DLATMR, | |||
| $ DLATMS, DSBGV, DSBGVD, DSBGVX, DSGT01, DSPGV, | |||
| $ DSPGVD, DSPGVX, DSYGV, DSYGVD, DSYGVX, XERBLA | |||
| * .. | |||
| @@ -460,7 +460,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = DLAMCH( 'Overflow' ) | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -240,7 +240,7 @@ | |||
| *> | |||
| *> B DOUBLE PRECISION array, dimension (LDB , max(NN)) | |||
| *> Used to hold the symmetric positive definite matrix for | |||
| *> the generailzed problem. | |||
| *> the generalized problem. | |||
| *> On exit, B contains the last matrix actually | |||
| *> used. | |||
| *> Modified. | |||
| @@ -408,7 +408,7 @@ | |||
| EXTERNAL LSAME, DLAMCH, DLARND | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL DLABAD, DLACPY, DLAFTS, DLASET, DLASUM, DLATMR, | |||
| EXTERNAL DLACPY, DLAFTS, DLASET, DLASUM, DLATMR, | |||
| $ DLATMS, DSBGV, DSBGVD, DSBGVX, DSGT01, DSPGV, | |||
| $ DSPGVD, DSPGVX, DSYGV, DSYGVD, DSYGVX, XERBLA, | |||
| $ DSYGV_2STAGE | |||
| @@ -470,7 +470,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = DLAMCH( 'Overflow' ) | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -214,7 +214,7 @@ | |||
| *> Not modified. | |||
| *> | |||
| *> D1 DOUBLE PRECISION array, dimension (max(NN)) | |||
| *> The eigenvalues of A, as computed by DSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by DSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> Modified. | |||
| @@ -502,11 +502,11 @@ | |||
| EXTERNAL DLAMCH, DLARND, DSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, DLABAD, DLACPY, DLAFTS, DLASET, DLATMR, | |||
| $ DLATMS, DSBEV, DSBEVD, DSBEVX, DSPEV, DSPEVD, | |||
| $ DSPEVX, DSTEV, DSTEVD, DSTEVR, DSTEVX, DSTT21, | |||
| $ DSTT22, DSYEV, DSYEVD, DSYEVR, DSYEVX, DSYT21, | |||
| $ DSYT22, XERBLA | |||
| EXTERNAL ALASVM, DLACPY, DLAFTS, DLASET, DLATMR, DLATMS, | |||
| $ DSBEV, DSBEVD, DSBEVX, DSPEV, DSPEVD, DSPEVX, | |||
| $ DSTEV, DSTEVD, DSTEVR, DSTEVX, DSTT21, DSTT22, | |||
| $ DSYEV, DSYEVD, DSYEVR, DSYEVX, DSYT21, DSYT22, | |||
| $ XERBLA | |||
| * .. | |||
| * .. Scalars in Common .. | |||
| CHARACTER*32 SRNAMT | |||
| @@ -574,7 +574,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = DLAMCH( 'Overflow' ) | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -214,7 +214,7 @@ | |||
| *> Not modified. | |||
| *> | |||
| *> D1 DOUBLE PRECISION array, dimension (max(NN)) | |||
| *> The eigenvalues of A, as computed by DSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by DSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> Modified. | |||
| @@ -502,7 +502,7 @@ | |||
| EXTERNAL DLAMCH, DLARND, DSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, DLABAD, DLACPY, DLAFTS, DLASET, DLATMR, | |||
| EXTERNAL ALASVM, DLACPY, DLAFTS, DLASET, DLATMR, | |||
| $ DLATMS, DSBEV, DSBEVD, DSBEVX, DSPEV, DSPEVD, | |||
| $ DSPEVX, DSTEV, DSTEVD, DSTEVR, DSTEVX, DSTT21, | |||
| $ DSTT22, DSYEV, DSYEVD, DSYEVR, DSYEVX, DSYT21, | |||
| @@ -577,7 +577,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = DLAMCH( 'Overflow' ) | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -363,7 +363,7 @@ | |||
| *> \verbatim | |||
| *> D1 is REAL array of | |||
| *> dimension( max(NN) ) | |||
| *> The eigenvalues of A, as computed by SSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by SSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> \endverbatim | |||
| @@ -645,10 +645,10 @@ | |||
| EXTERNAL ILAENV, SLAMCH, SLARND, SSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL SCOPY, SLABAD, SLACPY, SLASET, SLASUM, SLATMR, | |||
| $ SLATMS, SOPGTR, SORGTR, SPTEQR, SSPT21, SSPTRD, | |||
| $ SSTEBZ, SSTECH, SSTEDC, SSTEMR, SSTEIN, SSTEQR, | |||
| $ SSTERF, SSTT21, SSTT22, SSYT21, SSYTRD, XERBLA | |||
| EXTERNAL SCOPY, SLACPY, SLASET, SLASUM, SLATMR, SLATMS, | |||
| $ SOPGTR, SORGTR, SPTEQR, SSPT21, SSPTRD, SSTEBZ, | |||
| $ SSTECH, SSTEDC, SSTEMR, SSTEIN, SSTEQR, SSTERF, | |||
| $ SSTT21, SSTT22, SSYT21, SSYTRD, XERBLA | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| INTRINSIC ABS, INT, LOG, MAX, MIN, REAL, SQRT | |||
| @@ -715,7 +715,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = ONE / UNFL | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| LOG2UI = INT( LOG( ULPINV ) / LOG( TWO ) ) | |||
| @@ -384,7 +384,7 @@ | |||
| *> \verbatim | |||
| *> D1 is REAL array of | |||
| *> dimension( max(NN) ) | |||
| *> The eigenvalues of A, as computed by SSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by SSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> \endverbatim | |||
| @@ -666,10 +666,10 @@ | |||
| EXTERNAL ILAENV, SLAMCH, SLARND, SSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL SCOPY, SLABAD, SLACPY, SLASET, SLASUM, SLATMR, | |||
| $ SLATMS, SOPGTR, SORGTR, SPTEQR, SSPT21, SSPTRD, | |||
| $ SSTEBZ, SSTECH, SSTEDC, SSTEMR, SSTEIN, SSTEQR, | |||
| $ SSTERF, SSTT21, SSTT22, SSYT21, SSYTRD, XERBLA, | |||
| EXTERNAL SCOPY, SLACPY, SLASET, SLASUM, SLATMR, SLATMS, | |||
| $ SOPGTR, SORGTR, SPTEQR, SSPT21, SSPTRD, SSTEBZ, | |||
| $ SSTECH, SSTEDC, SSTEMR, SSTEIN, SSTEQR, SSTERF, | |||
| $ SSTT21, SSTT22, SSYT21, SSYTRD, XERBLA, | |||
| $ SSYTRD_2STAGE | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| @@ -737,7 +737,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = ONE / UNFL | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| LOG2UI = INT( LOG( ULPINV ) / LOG( TWO ) ) | |||
| @@ -400,7 +400,7 @@ | |||
| EXTERNAL SLCTSX, ILAENV, SLAMCH, SLANGE | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, SGESVD, SGET51, SGET53, SGGESX, SLABAD, | |||
| EXTERNAL ALASVM, SGESVD, SGET51, SGET53, SGGESX, | |||
| $ SLACPY, SLAKF2, SLASET, SLATM5, XERBLA | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| @@ -479,7 +479,6 @@ c MINWRK = MAX( 10*( NSIZE+1 ), 5*NSIZE*NSIZE / 2-2 ) | |||
| ULPINV = ONE / ULP | |||
| SMLNUM = SLAMCH( 'S' ) / ULP | |||
| BIGNUM = ONE / SMLNUM | |||
| CALL SLABAD( SMLNUM, BIGNUM ) | |||
| THRSH2 = TEN*THRESH | |||
| NTESTT = 0 | |||
| NERRS = 0 | |||
| @@ -985,7 +984,7 @@ c MINWRK = MAX( 10*( NSIZE+1 ), 5*NSIZE*NSIZE / 2-2 ) | |||
| $ / ' 2: A and B are upper triangular matrices, ', | |||
| $ / ' 3: A and B are as type 2, but each second diagonal ', | |||
| $ 'block in A_11 and ', / | |||
| $ ' each third diaongal block in A_22 are 2x2 blocks,', | |||
| $ ' each third diagonal block in A_22 are 2x2 blocks,', | |||
| $ / ' 4: A and B are block diagonal matrices, ', | |||
| $ / ' 5: (A,B) has potentially close or common ', | |||
| $ 'eigenvalues.', / ) | |||
| @@ -234,7 +234,7 @@ | |||
| *> | |||
| *> B REAL array, dimension (LDB , max(NN)) | |||
| *> Used to hold the symmetric positive definite matrix for | |||
| *> the generailzed problem. | |||
| *> the generalized problem. | |||
| *> On exit, B contains the last matrix actually | |||
| *> used. | |||
| *> Modified. | |||
| @@ -399,7 +399,7 @@ | |||
| EXTERNAL LSAME, SLAMCH, SLARND | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL SLABAD, SLACPY, SLAFTS, SLASET, SLASUM, SLATMR, | |||
| EXTERNAL SLACPY, SLAFTS, SLASET, SLASUM, SLATMR, | |||
| $ SLATMS, SSBGV, SSBGVD, SSBGVX, SSGT01, SSPGV, | |||
| $ SSPGVD, SSPGVX, SSYGV, SSYGVD, SSYGVX, XERBLA | |||
| * .. | |||
| @@ -460,7 +460,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = SLAMCH( 'Overflow' ) | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -240,7 +240,7 @@ | |||
| *> | |||
| *> B REAL array, dimension (LDB , max(NN)) | |||
| *> Used to hold the symmetric positive definite matrix for | |||
| *> the generailzed problem. | |||
| *> the generalized problem. | |||
| *> On exit, B contains the last matrix actually | |||
| *> used. | |||
| *> Modified. | |||
| @@ -408,7 +408,7 @@ | |||
| EXTERNAL LSAME, SLAMCH, SLARND | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL SLABAD, SLACPY, SLAFTS, SLASET, SLASUM, SLATMR, | |||
| EXTERNAL SLACPY, SLAFTS, SLASET, SLASUM, SLATMR, | |||
| $ SLATMS, SSBGV, SSBGVD, SSBGVX, SSGT01, SSPGV, | |||
| $ SSPGVD, SSPGVX, SSYGV, SSYGVD, SSYGVX, XERBLA, | |||
| $ SSYGV_2STAGE | |||
| @@ -470,7 +470,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = SLAMCH( 'Overflow' ) | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -214,7 +214,7 @@ | |||
| *> Not modified. | |||
| *> | |||
| *> D1 REAL array, dimension (max(NN)) | |||
| *> The eigenvalues of A, as computed by SSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by SSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> Modified. | |||
| @@ -502,11 +502,11 @@ | |||
| EXTERNAL SLAMCH, SLARND, SSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, SLABAD, SLACPY, SLAFTS, SLASET, SLATMR, | |||
| $ SLATMS, SSBEV, SSBEVD, SSBEVX, SSPEV, SSPEVD, | |||
| $ SSPEVX, SSTEV, SSTEVD, SSTEVR, SSTEVX, SSTT21, | |||
| $ SSTT22, SSYEV, SSYEVD, SSYEVR, SSYEVX, SSYT21, | |||
| $ SSYT22, XERBLA | |||
| EXTERNAL ALASVM, SLACPY, SLAFTS, SLASET, SLATMR, SLATMS, | |||
| $ SSBEV, SSBEVD, SSBEVX, SSPEV, SSPEVD, SSPEVX, | |||
| $ SSTEV, SSTEVD, SSTEVR, SSTEVX, SSTT21, SSTT22, | |||
| $ SSYEV, SSYEVD, SSYEVR, SSYEVX, SSYT21, SSYT22, | |||
| $ XERBLA | |||
| * .. | |||
| * .. Scalars in Common .. | |||
| CHARACTER*32 SRNAMT | |||
| @@ -574,7 +574,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = SLAMCH( 'Overflow' ) | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -214,7 +214,7 @@ | |||
| *> Not modified. | |||
| *> | |||
| *> D1 REAL array, dimension (max(NN)) | |||
| *> The eigenvalues of A, as computed by SSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by SSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> Modified. | |||
| @@ -502,7 +502,7 @@ | |||
| EXTERNAL SLAMCH, SLARND, SSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, SLABAD, SLACPY, SLAFTS, SLASET, SLATMR, | |||
| EXTERNAL ALASVM, SLACPY, SLAFTS, SLASET, SLATMR, | |||
| $ SLATMS, SSBEV, SSBEVD, SSBEVX, SSPEV, SSPEVD, | |||
| $ SSPEVX, SSTEV, SSTEVD, SSTEVR, SSTEVX, SSTT21, | |||
| $ SSTT22, SSYEV, SSYEVD, SSYEVR, SSYEVX, SSYT21, | |||
| @@ -577,7 +577,6 @@ | |||
| * | |||
| UNFL = SLAMCH( 'Safe minimum' ) | |||
| OVFL = SLAMCH( 'Overflow' ) | |||
| CALL SLABAD( UNFL, OVFL ) | |||
| ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -364,7 +364,7 @@ | |||
| *> \verbatim | |||
| *> D1 is DOUBLE PRECISION array of | |||
| *> dimension( max(NN) ) | |||
| *> The eigenvalues of A, as computed by ZSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by ZSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> \endverbatim | |||
| @@ -662,11 +662,10 @@ | |||
| EXTERNAL ILAENV, DLAMCH, DLARND, DSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL DCOPY, DLABAD, DLASUM, DSTEBZ, DSTECH, DSTERF, | |||
| $ XERBLA, ZCOPY, ZHET21, ZHETRD, ZHPT21, ZHPTRD, | |||
| $ ZLACPY, ZLASET, ZLATMR, ZLATMS, ZPTEQR, ZSTEDC, | |||
| $ ZSTEMR, ZSTEIN, ZSTEQR, ZSTT21, ZSTT22, ZUNGTR, | |||
| $ ZUPGTR | |||
| EXTERNAL DCOPY, DLASUM, DSTEBZ, DSTECH, DSTERF, XERBLA, | |||
| $ ZCOPY, ZHET21, ZHETRD, ZHPT21, ZHPTRD, ZLACPY, | |||
| $ ZLASET, ZLATMR, ZLATMS, ZPTEQR, ZSTEDC, ZSTEMR, | |||
| $ ZSTEIN, ZSTEQR, ZSTT21, ZSTT22, ZUNGTR, ZUPGTR | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| INTRINSIC ABS, DBLE, DCONJG, INT, LOG, MAX, MIN, SQRT | |||
| @@ -733,7 +732,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = ONE / UNFL | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| LOG2UI = INT( LOG( ULPINV ) / LOG( TWO ) ) | |||
| @@ -385,7 +385,7 @@ | |||
| *> \verbatim | |||
| *> D1 is DOUBLE PRECISION array of | |||
| *> dimension( max(NN) ) | |||
| *> The eigenvalues of A, as computed by ZSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by ZSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> \endverbatim | |||
| @@ -683,11 +683,11 @@ | |||
| EXTERNAL ILAENV, DLAMCH, DLARND, DSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL DCOPY, DLABAD, DLASUM, DSTEBZ, DSTECH, DSTERF, | |||
| $ XERBLA, ZCOPY, ZHET21, ZHETRD, ZHPT21, ZHPTRD, | |||
| $ ZLACPY, ZLASET, ZLATMR, ZLATMS, ZPTEQR, ZSTEDC, | |||
| $ ZSTEMR, ZSTEIN, ZSTEQR, ZSTT21, ZSTT22, ZUNGTR, | |||
| $ ZUPGTR, ZHETRD_2STAGE, DLASET | |||
| EXTERNAL DCOPY, DLASUM, DSTEBZ, DSTECH, DSTERF, XERBLA, | |||
| $ ZCOPY, ZHET21, ZHETRD, ZHPT21, ZHPTRD, ZLACPY, | |||
| $ ZLASET, ZLATMR, ZLATMS, ZPTEQR, ZSTEDC, ZSTEMR, | |||
| $ ZSTEIN, ZSTEQR, ZSTT21, ZSTT22, ZUNGTR, ZUPGTR, | |||
| $ ZHETRD_2STAGE, DLASET | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| INTRINSIC ABS, DBLE, DCONJG, INT, LOG, MAX, MIN, SQRT | |||
| @@ -754,7 +754,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = ONE / UNFL | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| LOG2UI = INT( LOG( ULPINV ) / LOG( TWO ) ) | |||
| @@ -394,7 +394,7 @@ | |||
| EXTERNAL ZLCTSX, ILAENV, DLAMCH, ZLANGE | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, DLABAD, XERBLA, ZGESVD, ZGET51, ZGGESX, | |||
| EXTERNAL ALASVM, XERBLA, ZGESVD, ZGET51, ZGGESX, | |||
| $ ZLACPY, ZLAKF2, ZLASET, ZLATM5 | |||
| * .. | |||
| * .. Scalars in Common .. | |||
| @@ -479,7 +479,6 @@ | |||
| ULPINV = ONE / ULP | |||
| SMLNUM = DLAMCH( 'S' ) / ULP | |||
| BIGNUM = ONE / SMLNUM | |||
| CALL DLABAD( SMLNUM, BIGNUM ) | |||
| THRSH2 = TEN*THRESH | |||
| NTESTT = 0 | |||
| NERRS = 0 | |||
| @@ -918,7 +917,7 @@ | |||
| $ / ' 2: A and B are upper triangular matrices, ', | |||
| $ / ' 3: A and B are as type 2, but each second diagonal ', | |||
| $ 'block in A_11 and ', / | |||
| $ ' each third diaongal block in A_22 are 2x2 blocks,', | |||
| $ ' each third diagonal block in A_22 are 2x2 blocks,', | |||
| $ / ' 4: A and B are block diagonal matrices, ', | |||
| $ / ' 5: (A,B) has potentially close or common ', | |||
| $ 'eigenvalues.', / ) | |||
| @@ -236,7 +236,7 @@ | |||
| *> | |||
| *> B COMPLEX*16 array, dimension (LDB , max(NN)) | |||
| *> Used to hold the Hermitian positive definite matrix for | |||
| *> the generailzed problem. | |||
| *> the generalized problem. | |||
| *> On exit, B contains the last matrix actually | |||
| *> used. | |||
| *> Modified. | |||
| @@ -418,7 +418,7 @@ | |||
| EXTERNAL LSAME, DLAMCH, DLARND | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL DLABAD, DLAFTS, DLASUM, XERBLA, ZHBGV, ZHBGVD, | |||
| EXTERNAL DLAFTS, DLASUM, XERBLA, ZHBGV, ZHBGVD, | |||
| $ ZHBGVX, ZHEGV, ZHEGVD, ZHEGVX, ZHPGV, ZHPGVD, | |||
| $ ZHPGVX, ZLACPY, ZLASET, ZLATMR, ZLATMS, ZSGT01 | |||
| * .. | |||
| @@ -481,7 +481,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = DLAMCH( 'Overflow' ) | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -242,7 +242,7 @@ | |||
| *> | |||
| *> B COMPLEX*16 array, dimension (LDB , max(NN)) | |||
| *> Used to hold the Hermitian positive definite matrix for | |||
| *> the generailzed problem. | |||
| *> the generalized problem. | |||
| *> On exit, B contains the last matrix actually | |||
| *> used. | |||
| *> Modified. | |||
| @@ -426,7 +426,7 @@ | |||
| EXTERNAL LSAME, DLAMCH, DLARND | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL DLABAD, DLAFTS, DLASUM, XERBLA, ZHBGV, ZHBGVD, | |||
| EXTERNAL DLAFTS, DLASUM, XERBLA, ZHBGV, ZHBGVD, | |||
| $ ZHBGVX, ZHEGV, ZHEGVD, ZHEGVX, ZHPGV, ZHPGVD, | |||
| $ ZHPGVX, ZLACPY, ZLASET, ZLATMR, ZLATMS, ZSGT01, | |||
| $ ZHEGV_2STAGE | |||
| @@ -490,7 +490,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = DLAMCH( 'Overflow' ) | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -204,7 +204,7 @@ | |||
| *> Not modified. | |||
| *> | |||
| *> D1 DOUBLE PRECISION array, dimension (max(NN)) | |||
| *> The eigenvalues of A, as computed by ZSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by ZSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> Modified. | |||
| @@ -391,10 +391,10 @@ | |||
| EXTERNAL DLAMCH, DLARND, DSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, DLABAD, DLAFTS, XERBLA, ZHBEV, ZHBEVD, | |||
| $ ZHBEVX, ZHEEV, ZHEEVD, ZHEEVR, ZHEEVX, ZHET21, | |||
| $ ZHET22, ZHPEV, ZHPEVD, ZHPEVX, ZLACPY, ZLASET, | |||
| $ ZLATMR, ZLATMS | |||
| EXTERNAL ALASVM, DLAFTS, XERBLA, ZHBEV, ZHBEVD, ZHBEVX, | |||
| $ ZHEEV, ZHEEVD, ZHEEVR, ZHEEVX, ZHET21, ZHET22, | |||
| $ ZHPEV, ZHPEVD, ZHPEVX, ZLACPY, ZLASET, ZLATMR, | |||
| $ ZLATMS | |||
| * .. | |||
| * .. Intrinsic Functions .. | |||
| INTRINSIC ABS, DBLE, INT, LOG, MAX, MIN, SQRT | |||
| @@ -451,7 +451,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = DLAMCH( 'Overflow' ) | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||
| @@ -204,7 +204,7 @@ | |||
| *> Not modified. | |||
| *> | |||
| *> D1 DOUBLE PRECISION array, dimension (max(NN)) | |||
| *> The eigenvalues of A, as computed by ZSTEQR simlutaneously | |||
| *> The eigenvalues of A, as computed by ZSTEQR simultaneously | |||
| *> with Z. On exit, the eigenvalues in D1 correspond with the | |||
| *> matrix in A. | |||
| *> Modified. | |||
| @@ -391,7 +391,7 @@ | |||
| EXTERNAL DLAMCH, DLARND, DSXT1 | |||
| * .. | |||
| * .. External Subroutines .. | |||
| EXTERNAL ALASVM, DLABAD, DLAFTS, XERBLA, ZHBEV, ZHBEVD, | |||
| EXTERNAL ALASVM, DLAFTS, XERBLA, ZHBEV, ZHBEVD, | |||
| $ ZHBEVX, ZHEEV, ZHEEVD, ZHEEVR, ZHEEVX, ZHET21, | |||
| $ ZHET22, ZHPEV, ZHPEVD, ZHPEVX, ZLACPY, ZLASET, | |||
| $ ZHEEVD_2STAGE, ZHEEVR_2STAGE, ZHEEVX_2STAGE, | |||
| @@ -453,7 +453,6 @@ | |||
| * | |||
| UNFL = DLAMCH( 'Safe minimum' ) | |||
| OVFL = DLAMCH( 'Overflow' ) | |||
| CALL DLABAD( UNFL, OVFL ) | |||
| ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' ) | |||
| ULPINV = ONE / ULP | |||
| RTUNFL = SQRT( UNFL ) | |||