QZ (version 0.1-7)

QZ Decomposition Reordering: Reordering QZ Decomposition

Description

This function performs QZ decomposition on input paired matrices (A,B) or a single matrix A with reordering.

Usage

ordqz(A, B = NULL, cluster = NULL,
        keyword = c("lhp", "rhp", "udi", "udo", "ref", "cef"), ...)

Arguments

A

a 'complex/real' matrix, dim = c(N, N).

B

a 'complex/real' matrix, dim = c(N, N).

cluster

specifies the eigenvalues in the selected cluster.

keyword

as similarly used in MATLAB.

...

options to qz.* functions.

Value

Returns a list from the call.

Details

Either cluster or keyword should be specified.

cluster actually is the same as select in all qz.* functions.

keywork actually is similar as MATLAB.

keyword	Selected Region
'lhp'	Left-half plane (real(E) < 0)
'rhp'	Right-half plane (real(E) > 0)
'udi'	Interior of unit disk (abs(E) < 1)
'udo'	Exterior of unit disk (abs(E) < 1)
'ref'	Real eigenvalues first (top-left conner)
'cef'	Complex eigenvalues first (top-left conner)

References

Anderson, E., et al. (1999) LAPACK User's Guide, 3rd edition, SIAM, Philadelphia.

http://en.wikipedia.org/wiki/Schur_decomposition

See Also

qz, geigen.

Examples

# NOT RUN {
<!-- % \dontrun{ -->
# }
# NOT RUN {
library(QZ, quiet = TRUE)

# Reordering eigenvalues
(ret <- ordqz(exAB1$A, exAB1$B, keyword = "lhp"))
(ret <- ordqz(exAB1$A, exAB1$B, keyword = "rhp"))
(ret <- ordqz(exAB1$A, exAB1$B, keyword = "udi"))
(ret <- ordqz(exAB1$A, exAB1$B, keyword = "udo"))
(ret <- ordqz(exAB1$A, exAB1$B, keyword = "ref"))
(ret <- ordqz(exAB1$A, exAB1$B, keyword = "cef"))

# }
# NOT RUN {
<!-- % } -->
# }