23_qz: QZ Decomposition

Description Usage Arguments Details Value Author(s) References See Also Examples

Description

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

Usage

1
  qz(A, B = NULL, select = NULL, only.values = FALSE, ...)

Arguments

A

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

B

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

select

specifies the eigenvalues in the selected cluster.

only.values

if 'TRUE', only the eigenvalues are computed and returned, otherwise both eigenvalues and eigenvectors are returned.

...

options to qz.* functions.

Details

If select is NULL, then call one of qz.zgges, qz.dgges, qz.zgees, or qz.dgees depending on the input arguments and types.

If select is not NULL, then call one of qz.zgges + qz.ztgsen, qz.dgges + qz.dtgsen,
qz.zgees + qz.ztrsen, or qz.dgees + qz.dtrsen depending on the input arguments and types.

Value

Returns a list from the call.

Author(s)

Wei-Chen Chen [email protected]

References

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

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

See Also

ordqz, geigen.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
library(QZ, quiet = TRUE)

### http://www.nag.com/lapack-ex/node124.html
(ret <- qz(exAB1$A, exAB1$B))

### http://www.nag.com/lapack-ex/node119.html
(ret <- qz(exAB2$A, exAB2$B))

### http://www.nag.com/lapack-ex/node94.html
(ret <- qz(exA1$A))

### http://www.nag.com/lapack-ex/node89.html
(ret <- qz(exA2$A))

# Reordering eigenvalues
select1 <- c(TRUE, FALSE, FALSE, TRUE)
select2 <- c(FALSE, TRUE, TRUE, FALSE)
(ret <- qz(exAB1$A, exAB1$B, select = select1))
(ret <- qz(exAB2$A, exAB2$B, select = select2))
(ret <- qz(exA1$A, select = select1))
(ret <- qz(exA2$A, select = select1))

snoweye/QZ documentation built on Dec. 8, 2019, 6:15 p.m.