|
version 1.8, 2010/12/21 13:53:56
|
version 1.9, 2011/07/22 07:38:21
|
|
Line 1
|
Line 1
|
| SUBROUTINE ZTGEX2( WANTQ, WANTZ, N, A, LDA, B, LDB, Q, LDQ, Z, |
SUBROUTINE ZTGEX2( WANTQ, WANTZ, N, A, LDA, B, LDB, Q, LDQ, Z, |
| $ LDZ, J1, INFO ) |
$ LDZ, J1, INFO ) |
| * |
* |
| * -- LAPACK auxiliary routine (version 3.2.2) -- |
* -- LAPACK auxiliary 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..-- |
| * June 2010 |
* -- April 2011 -- |
| * |
* |
| * .. Scalar Arguments .. |
* .. Scalar Arguments .. |
| LOGICAL WANTQ, WANTZ |
LOGICAL WANTQ, WANTZ |
|
Line 28
|
Line 28
|
| * Optionally, the matrices Q and Z of generalized Schur vectors are |
* Optionally, the matrices Q and Z of generalized Schur vectors are |
| * updated. |
* updated. |
| * |
* |
| * Q(in) * A(in) * Z(in)' = Q(out) * A(out) * Z(out)' |
* Q(in) * A(in) * Z(in)**H = Q(out) * A(out) * Z(out)**H |
| * Q(in) * B(in) * Z(in)' = Q(out) * B(out) * Z(out)' |
* Q(in) * B(in) * Z(in)**H = Q(out) * B(out) * Z(out)**H |
| * |
* |
| * |
* |
| * Arguments |
* Arguments |
|
Line 213
|
Line 213
|
| IF( WANDS ) THEN |
IF( WANDS ) THEN |
| * |
* |
| * Strong stability test: |
* Strong stability test: |
| * F-norm((A-QL'*S*QR, B-QL'*T*QR)) <= O(EPS*F-norm((A, B))) |
* F-norm((A-QL**H*S*QR, B-QL**H*T*QR)) <= O(EPS*F-norm((A, B))) |
| * |
* |
| CALL ZLACPY( 'Full', M, M, S, LDST, WORK, M ) |
CALL ZLACPY( 'Full', M, M, S, LDST, WORK, M ) |
| CALL ZLACPY( 'Full', M, M, T, LDST, WORK( M*M+1 ), M ) |
CALL ZLACPY( 'Full', M, M, T, LDST, WORK( M*M+1 ), M ) |
![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to ztgex2.f CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [local] / rpl / lapack / lapack