next up previous contents index
Next: Examples Up: Standard Nonsymmetric Eigenvalue Problems Previous: Purpose   Contents   Index

Arguments

A
(input/output) REAL or COMPLEX square array, shape $(:,:)$.
On entry, the matrix $A$.
On exit, the contents of A are destroyed.

w
(output) REAL or COMPLEX array, shape $(:)$ with $size$(w) $= size$(A,1).
The computed eigenvalues.
${\it w(:)}$ ::= WR(:), WI(:) $\mid$ W(:),
where
WR(:), WI(:) are of REAL type (for the real and imaginary parts) and W(:) is of COMPLEX type.
Note: If $A$ is real, then a complex-conjugate pair appear consecutively, with the eigenvalue having the positive imaginary part appearing first.

VL
Optional (output) REAL or COMPLEX square array, shape $(:,:)$ with $size$(VL,1) $= size$(A,1).
The left eigenvectors $u_j$ are stored in the columns of VL in the order of their eigenvalues. Each eigenvector is scaled so that the Euclidean norm is 1 and the largest component is real.
Note: If $A$ is real then complex eigenvectors, like their eigenvalues, occur in complex conjugate pairs. The real and imaginary parts of the first eigenvector of the pair are stored in VL$_{:,j}$ and VL$_{:,j+1}$, respectively. Thus a complex conjugate pair is given by

\begin{displaymath}u_j = {\bf VL}_{:,j}\,+\,i\, {\bf VL}_{:,j+1}\,,\;\;\;
u_{j+1} = {\bf VL}_{:,j}\,-\,i\, {\bf VL}_{:,j+1} \end{displaymath}

VR
Optional (output) REAL or COMPLEX square array, shape $(:,:)$ with $size$(VR,1) $= size$(A,1).
The right eigenvectors $v_j$ are stored in the columns of VR in the order of their eigenvalues. Each eigenvector is scaled so that the Euclidean norm is 1 and the largest component is real.
Note: If $A$ is real then complex eigenvectors, like their eigenvalues, occur in complex conjugate pairs. The real and imaginary parts of the first eigenvector of the pair are stored in VR$_{:,j}$ and VR$_{:,j+1}$, respectively. Thus a complex conjugate pair is given by

\begin{displaymath}v_j = {\bf VR}_{:,j}\,+\,i\, {\bf VR}_{:,j+1}\,,\;\;\;
v_{j+1} = {\bf VR}_{:,j}\,-\,i\, {\bf VR}_{:,j+1} \end{displaymath}

INFO
Optional (output) INTEGER.

\begin{infoarg}
\item[{$=$\ 0:}] successful exit.
\item[{$<$\ 0:}] if {\bf INF...
...ain eigenvalues which have converged. \\
$n$\ is the order of $A$ \end{infoarg}
If INFO is not present and an error occurs, then the program is terminated with an error message.
References: [1] and [17,9,20].
next up previous contents index
Next: Examples Up: Standard Nonsymmetric Eigenvalue Problems Previous: Purpose   Contents   Index
Susan Blackford 2001-08-19