F08YFF (DTGEXC) reorders the generalized Schur factorization of a matrix pair in real generalized Schur form.
SUBROUTINE F08YFF ( |
WANTQ, WANTZ, N, A, LDA, B, LDB, Q, LDQ, Z, LDZ, IFST, ILST, WORK, LWORK, INFO) |
INTEGER |
N, LDA, LDB, LDQ, LDZ, IFST, ILST, LWORK, INFO |
REAL (KIND=nag_wp) |
A(LDA,*), B(LDB,*), Q(LDQ,*), Z(LDZ,*), WORK(max(1,LWORK)) |
LOGICAL |
WANTQ, WANTZ |
|
F08YFF (DTGEXC) reorders the generalized real
by
matrix pair
in real generalized Schur form, so that the diagonal element or block of
with row index
is moved to row
, using an orthogonal equivalence transformation. That is,
and
are factorized as
where
are also in real generalized Schur form.
The pair
are in real generalized Schur form if
is block upper triangular with
by
and
by
diagonal blocks and
is upper triangular as returned, for example, by
F08XAF (DGGES), or
F08XEF (DHGEQZ) with
.
If
and
are the result of a generalized Schur factorization of a matrix pair
then, optionally, the matrices
and
can be updated as
and
.
Anderson E, Bai Z, Bischof C, Blackford S, Demmel J, Dongarra J J, Du Croz J J, Greenbaum A, Hammarling S, McKenney A and Sorensen D (1999)
LAPACK Users' Guide (3rd Edition) SIAM, Philadelphia
http://www.netlib.org/lapack/lugThe computed generalized Schur form is nearly the exact generalized Schur form for nearby matrices
and
, where
and
is the
machine precision. See Section 4.11 of
Anderson et al. (1999) for further details of error bounds for the generalized nonsymmetric eigenproblem.
The complex analogue of this routine is
F08YTF (ZTGEXC).
This example exchanges blocks
and
of the matrix pair
, where