version 1.1, 2010/01/26 15:22:46
|
version 1.8, 2011/07/22 07:38:21
|
Line 1
|
Line 1
|
SUBROUTINE ZTRSEN( JOB, COMPQ, SELECT, N, T, LDT, Q, LDQ, W, M, S, |
SUBROUTINE ZTRSEN( JOB, COMPQ, SELECT, N, T, LDT, Q, LDQ, W, M, S, |
$ SEP, WORK, LWORK, INFO ) |
$ SEP, WORK, LWORK, INFO ) |
* |
* |
* -- LAPACK routine (version 3.2) -- |
* -- LAPACK routine (version 3.3.1) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* November 2006 |
* -- April 2011 -- |
* |
* |
* Modified to call ZLACN2 in place of ZLACON, 10 Feb 03, SJH. |
* Modified to call ZLACN2 in place of ZLACON, 10 Feb 03, SJH. |
* |
* |
Line 119
|
Line 119
|
* transformation Z to move them to the top left corner of T. In other |
* transformation Z to move them to the top left corner of T. In other |
* words, the selected eigenvalues are the eigenvalues of T11 in: |
* words, the selected eigenvalues are the eigenvalues of T11 in: |
* |
* |
* Z'*T*Z = ( T11 T12 ) n1 |
* Z**H * T * Z = ( T11 T12 ) n1 |
* ( 0 T22 ) n2 |
* ( 0 T22 ) n2 |
* n1 n2 |
* n1 n2 |
* |
* |
* where N = n1+n2 and Z' means the conjugate transpose of Z. The first |
* where N = n1+n2. The first |
* n1 columns of Z span the specified invariant subspace of T. |
* n1 columns of Z span the specified invariant subspace of T. |
* |
* |
* If T has been obtained from the Schur factorization of a matrix |
* If T has been obtained from the Schur factorization of a matrix |
* A = Q*T*Q', then the reordered Schur factorization of A is given by |
* A = Q*T*Q**H, then the reordered Schur factorization of A is given by |
* A = (Q*Z)*(Z'*T*Z)*(Q*Z)', and the first n1 columns of Q*Z span the |
* A = (Q*Z)*(Z**H*T*Z)*(Q*Z)**H, and the first n1 columns of Q*Z span the |
* corresponding invariant subspace of A. |
* corresponding invariant subspace of A. |
* |
* |
* The reciprocal condition number of the average of the eigenvalues of |
* The reciprocal condition number of the average of the eigenvalues of |
Line 331
|
Line 331
|
$ IERR ) |
$ IERR ) |
ELSE |
ELSE |
* |
* |
* Solve T11'*R - R*T22' = scale*X. |
* Solve T11**H*R - R*T22**H = scale*X. |
* |
* |
CALL ZTRSYL( 'C', 'C', -1, N1, N2, T, LDT, |
CALL ZTRSYL( 'C', 'C', -1, N1, N2, T, LDT, |
$ T( N1+1, N1+1 ), LDT, WORK, N1, SCALE, |
$ T( N1+1, N1+1 ), LDT, WORK, N1, SCALE, |