NAG Library Routine Document

f08qyf  (ztrsna)

1

Purpose

2

Specification

3

Description

4

References

5

Arguments

1:     $\mathbf{job}$ – Character(1)Input

On entry: indicates whether condition numbers are required for eigenvalues and/or eigenvectors.

$\mathbf{job}=\text{'E'}$

Condition numbers for eigenvalues only are computed.

$\mathbf{job}=\text{'V'}$

Condition numbers for eigenvectors only are computed.

$\mathbf{job}=\text{'B'}$

Condition numbers for both eigenvalues and eigenvectors are computed.

Constraint: $\mathbf{job}=\text{'E'}$, $\text{'V'}$ or $\text{'B'}$.

2:     $\mathbf{howmny}$ – Character(1)Input

On entry: indicates how many condition numbers are to be computed.

$\mathbf{howmny}=\text{'A'}$

Condition numbers for all eigenpairs are computed.

$\mathbf{howmny}=\text{'S'}$

Condition numbers for selected eigenpairs (as specified by select) are computed.

Constraint: $\mathbf{howmny}=\text{'A'}$ or $\text{'S'}$.

3:     $\mathbf{select}\left(*\right)$ – Logical arrayInput

Note: the dimension of the array select must be at least $\max \left(1,\mathbf{n}\right)$ if $\mathbf{howmny}=\text{'S'}$, and at least $1$ otherwise.

On entry: specifies the eigenpairs for which condition numbers are to be computed if $\mathbf{howmny}=\text{'S'}$. To select condition numbers for the eigenpair corresponding to the eigenvalue ${\lambda }_j$, $\mathbf{select}\left(j\right)$ must be set to .TRUE..

If $\mathbf{howmny}=\text{'A'}$, select is not referenced.

4:     $\mathbf{n}$ – IntegerInput

On entry: $n$, the order of the matrix $T$.

Constraint: $\mathbf{n}\ge 0$.

5:     $\mathbf{t}\left(\mathbf{ldt},*\right)$ – Complex (Kind=nag_wp) arrayInput

Note: the second dimension of the array t must be at least $\max \left(1,\mathbf{n}\right)$.

On entry: the $n$ by $n$ upper triangular matrix $T$, as returned by f08psf (zhseqr).

6:     $\mathbf{ldt}$ – IntegerInput

On entry: the first dimension of the array t as declared in the (sub)program from which f08qyf (ztrsna) is called.

Constraint: $\mathbf{ldt}\ge \max \left(1,\mathbf{n}\right)$.

7:     $\mathbf{vl}\left(\mathbf{ldvl},*\right)$ – Complex (Kind=nag_wp) arrayInput

Note: the second dimension of the array vl must be at least $\max \left(1,\mathbf{mm}\right)$ if $\mathbf{job}=\text{'E'}$ or $\text{'B'}$ and at least $1$ if $\mathbf{job}=\text{'V'}$.

On entry: if $\mathbf{job}=\text{'E'}$ or $\text{'B'}$, vl must contain the left eigenvectors of $T$ (or of any matrix $QT{Q}^{\mathrm{H}}$ with $Q$ unitary) corresponding to the eigenpairs specified by howmny and select. The eigenvectors must be stored in consecutive columns of vl, as returned by f08pxf (zhsein) or f08qxf (ztrevc).

If $\mathbf{job}=\text{'V'}$, vl is not referenced.

8:     $\mathbf{ldvl}$ – IntegerInput

On entry: the first dimension of the array vl as declared in the (sub)program from which f08qyf (ztrsna) is called.

Constraints:

• if $\mathbf{job}=\text{'E'}$ or $\text{'B'}$, $\mathbf{ldvl}\ge \max \left(1,\mathbf{n}\right)$;
• if $\mathbf{job}=\text{'V'}$, $\mathbf{ldvl}\ge 1$.

9:     $\mathbf{vr}\left(\mathbf{ldvr},*\right)$ – Complex (Kind=nag_wp) arrayInput

Note: the second dimension of the array vr must be at least $\max \left(1,\mathbf{mm}\right)$ if $\mathbf{job}=\text{'E'}$ or $\text{'B'}$ and at least $1$ if $\mathbf{job}=\text{'V'}$.

On entry: if $\mathbf{job}=\text{'E'}$ or $\text{'B'}$, vr must contain the right eigenvectors of $T$ (or of any matrix $QT{Q}^{\mathrm{H}}$ with $Q$ unitary) corresponding to the eigenpairs specified by howmny and select. The eigenvectors must be stored in consecutive columns of vr, as returned by f08pxf (zhsein) or f08qxf (ztrevc).

If $\mathbf{job}=\text{'V'}$, vr is not referenced.

10:   $\mathbf{ldvr}$ – IntegerInput

On entry: the first dimension of the array vr as declared in the (sub)program from which f08qyf (ztrsna) is called.

Constraints:

• if $\mathbf{job}=\text{'E'}$ or $\text{'B'}$, $\mathbf{ldvr}\ge \max \left(1,\mathbf{n}\right)$;
• if $\mathbf{job}=\text{'V'}$, $\mathbf{ldvr}\ge 1$.

11:   $\mathbf{s}\left(*\right)$ – Real (Kind=nag_wp) arrayOutput

Note: the dimension of the array s must be at least $\max \left(1,\mathbf{mm}\right)$ if $\mathbf{job}=\text{'E'}$ or $\text{'B'}$, and at least $1$ otherwise.

On exit: the reciprocal condition numbers of the selected eigenvalues if $\mathbf{job}=\text{'E'}$ or $\text{'B'}$, stored in consecutive elements of the array. Thus $\mathbf{s}\left(j\right)$, $\mathbf{sep}\left(j\right)$ and the $j$th columns of vl and vr all correspond to the same eigenpair (but not in general the $j$th eigenpair unless all eigenpairs have been selected).

If $\mathbf{job}=\text{'V'}$, s is not referenced.

12:   $\mathbf{sep}\left(*\right)$ – Real (Kind=nag_wp) arrayOutput

Note: the dimension of the array sep must be at least $\max \left(1,\mathbf{mm}\right)$ if $\mathbf{job}=\text{'V'}$ or $\text{'B'}$, and at least $1$ otherwise.

On exit: the estimated reciprocal condition numbers of the selected right eigenvectors if $\mathbf{job}=\text{'V'}$ or $\text{'B'}$, stored in consecutive elements of the array.

If $\mathbf{job}=\text{'E'}$, sep is not referenced.

13:   $\mathbf{mm}$ – IntegerInput

On entry: the number of elements in the arrays s and sep, and the number of columns in the arrays vl and vr (if used). The precise number required, $\mathit{m}$, is $n$ if $\mathbf{howmny}=\text{'A'}$; if $\mathbf{howmny}=\text{'S'}$, $\mathit{m}$ is the number of selected eigenpairs (see select), in which case $0\le \mathit{m}\le n$.

Constraints:

• if $\mathbf{howmny}=\text{'A'}$, $\mathbf{mm}\ge \mathbf{n}$;
• otherwise $\mathbf{mm}\ge \mathit{m}$.

14:   $\mathbf{m}$ – IntegerOutput

On exit: $\mathit{m}$, the number of selected eigenpairs. If $\mathbf{howmny}=\text{'A'}$, m is set to $n$.

15:   $\mathbf{work}\left(\mathbf{ldwork},*\right)$ – Complex (Kind=nag_wp) arrayWorkspace

Note: the second dimension of the array work must be at least $\max \left(1,\mathbf{n}+1\right)$ if $\mathbf{job}=\text{'V'}$ or $\text{'B'}$ and at least $1$ if $\mathbf{job}=\text{'E'}$.

If $\mathbf{job}=\text{'E'}$, work is not referenced.

16:   $\mathbf{ldwork}$ – IntegerInput

On entry: the first dimension of the array work as declared in the (sub)program from which f08qyf (ztrsna) is called.

Constraints:

• if $\mathbf{job}=\text{'V'}$ or $\text{'B'}$, $\mathbf{ldwork}\ge \max \left(1,\mathbf{n}\right)$;
• if $\mathbf{job}=\text{'E'}$, $\mathbf{ldwork}\ge 1$.

17:   $\mathbf{rwork}\left(*\right)$ – Real (Kind=nag_wp) arrayWorkspace

Note: the dimension of the array rwork must be at least $\max \left(1,\mathbf{n}\right)$.

18:   $\mathbf{info}$ – IntegerOutput

On exit: $\mathbf{info}=0$ unless the routine detects an error (see Section 6).

6

Error Indicators and Warnings

$\mathbf{info}<0$

If $\mathbf{info}=-i$, argument $i$ had an illegal value. An explanatory message is output, and execution of the program is terminated.

7

Accuracy

8

Parallelism and Performance

9

Further Comments

10

Example

10.1

Program Text

Program Text (f08qyfe.f90)

10.2

Program Data

Program Data (f08qyfe.d)

10.3

Program Results

Program Results (f08qyfe.r)