NAG Library Routine Document
F08QTF (ZTREXC)
1 Purpose
F08QTF (ZTREXC) reorders the Schur factorization of a complex general matrix.
2 Specification
INTEGER |
N, LDT, LDQ, IFST, ILST, INFO |
COMPLEX (KIND=nag_wp) |
T(LDT,*), Q(LDQ,*) |
CHARACTER(1) |
COMPQ |
|
The routine may be called by its
LAPACK
name ztrexc.
3 Description
F08QTF (ZTREXC) reorders the Schur factorization of a complex general matrix
, so that the diagonal element of
with row index
IFST is moved to row
ILST.
The reordered Schur form is computed by a unitary similarity transformation: . Optionally the updated matrix of Schur vectors is computed as , giving .
4 References
Golub G H and Van Loan C F (1996) Matrix Computations (3rd Edition) Johns Hopkins University Press, Baltimore
5 Parameters
- 1: – CHARACTER(1)Input
-
On entry: indicates whether the matrix
of Schur vectors is to be updated.
- The matrix of Schur vectors is updated.
- No Schur vectors are updated.
Constraint:
or .
- 2: – INTEGERInput
-
On entry: , the order of the matrix .
Constraint:
.
- 3: – COMPLEX (KIND=nag_wp) arrayInput/Output
-
Note: the second dimension of the array
T
must be at least
.
On entry: the
by
upper triangular matrix
, as returned by
F08PSF (ZHSEQR).
On exit:
T is overwritten by the updated matrix
.
- 4: – INTEGERInput
-
On entry: the first dimension of the array
T as declared in the (sub)program from which F08QTF (ZTREXC) is called.
Constraint:
.
- 5: – COMPLEX (KIND=nag_wp) arrayInput/Output
-
Note: the second dimension of the array
Q
must be at least
if
and at least
if
.
On entry: if
,
Q must contain the
by
unitary matrix
of Schur vectors.
On exit: if
,
Q contains the updated matrix of Schur vectors.
If
,
Q is not referenced.
- 6: – INTEGERInput
-
On entry: the first dimension of the array
Q as declared in the (sub)program from which F08QTF (ZTREXC) is called.
Constraints:
- if , ;
- if , .
- 7: – INTEGERInput
- 8: – INTEGERInput
-
On entry:
IFST and
ILST must specify the reordering of the diagonal elements of
. The element with row index
IFST is moved to row
ILST by a sequence of exchanges between adjacent elements.
Constraint:
and .
- 9: – INTEGEROutput
On exit:
unless the routine detects an error (see
Section 6).
6 Error Indicators and Warnings
-
If , argument had an illegal value. An explanatory message is output, and execution of the program is terminated.
7 Accuracy
The computed matrix
is exactly similar to a matrix
, where
and
is the
machine precision.
The values of the eigenvalues are never changed by the reordering.
8 Parallelism and Performance
F08QTF (ZTREXC) is not threaded by NAG in any implementation.
F08QTF (ZTREXC) makes calls to BLAS and/or LAPACK routines, which may be threaded within the vendor library used by this implementation. Consult the documentation for the vendor library for further information.
Please consult the
X06 Chapter Introduction for information on how to control and interrogate the OpenMP environment used within this routine. Please also consult the
Users' Note for your implementation for any additional implementation-specific information.
The total number of real floating-point operations is approximately if , and if , where .
The real analogue of this routine is
F08QFF (DTREXC).
10 Example
This example reorders the Schur factorization of the matrix
so that element
is moved to
, where
10.1 Program Text
Program Text (f08qtfe.f90)
10.2 Program Data
Program Data (f08qtfe.d)
10.3 Program Results
Program Results (f08qtfe.r)