 |
LAPACK 3.11.0
LAPACK: Linear Algebra PACKage
|
|
LAPACK 3.11.0
LAPACK: Linear Algebra PACKage
|
| recursive subroutine dlaqz3 | ( | logical, intent(in) | ILSCHUR, |
| logical, intent(in) | ILQ, | ||
| logical, intent(in) | ILZ, | ||
| integer, intent(in) | N, | ||
| integer, intent(in) | ILO, | ||
| integer, intent(in) | IHI, | ||
| integer, intent(in) | NW, | ||
| double precision, dimension( lda, * ), intent(inout) | A, | ||
| integer, intent(in) | LDA, | ||
| double precision, dimension( ldb, * ), intent(inout) | B, | ||
| integer, intent(in) | LDB, | ||
| double precision, dimension( ldq, * ), intent(inout) | Q, | ||
| integer, intent(in) | LDQ, | ||
| double precision, dimension( ldz, * ), intent(inout) | Z, | ||
| integer, intent(in) | LDZ, | ||
| integer, intent(out) | NS, | ||
| integer, intent(out) | ND, | ||
| double precision, dimension( * ), intent(inout) | ALPHAR, | ||
| double precision, dimension( * ), intent(inout) | ALPHAI, | ||
| double precision, dimension( * ), intent(inout) | BETA, | ||
| double precision, dimension( ldqc, * ) | QC, | ||
| integer, intent(in) | LDQC, | ||
| double precision, dimension( ldzc, * ) | ZC, | ||
| integer, intent(in) | LDZC, | ||
| double precision, dimension( * ) | WORK, | ||
| integer, intent(in) | LWORK, | ||
| integer, intent(in) | REC, | ||
| integer, intent(out) | INFO | ||
| ) |
DLAQZ3
Download DLAQZ3 + dependencies [TGZ] [ZIP] [TXT]
DLAQZ3 performs AED
| [in] | ILSCHUR | ILSCHUR is LOGICAL
Determines whether or not to update the full Schur form |
| [in] | ILQ | ILQ is LOGICAL
Determines whether or not to update the matrix Q |
| [in] | ILZ | ILZ is LOGICAL
Determines whether or not to update the matrix Z |
| [in] | N | N is INTEGER
The order of the matrices A, B, Q, and Z. N >= 0. |
| [in] | ILO | ILO is INTEGER |
| [in] | IHI | IHI is INTEGER
ILO and IHI mark the rows and columns of (A,B) which
are to be normalized |
| [in] | NW | NW is INTEGER
The desired size of the deflation window. |
| [in,out] | A | A is DOUBLE PRECISION array, dimension (LDA, N) |
| [in] | LDA | LDA is INTEGER
The leading dimension of the array A. LDA >= max( 1, N ). |
| [in,out] | B | B is DOUBLE PRECISION array, dimension (LDB, N) |
| [in] | LDB | LDB is INTEGER
The leading dimension of the array B. LDB >= max( 1, N ). |
| [in,out] | Q | Q is DOUBLE PRECISION array, dimension (LDQ, N) |
| [in] | LDQ | LDQ is INTEGER |
| [in,out] | Z | Z is DOUBLE PRECISION array, dimension (LDZ, N) |
| [in] | LDZ | LDZ is INTEGER |
| [out] | NS | NS is INTEGER
The number of unconverged eigenvalues available to
use as shifts. |
| [out] | ND | ND is INTEGER
The number of converged eigenvalues found. |
| [out] | ALPHAR | ALPHAR is DOUBLE PRECISION array, dimension (N)
The real parts of each scalar alpha defining an eigenvalue
of GNEP. |
| [out] | ALPHAI | ALPHAI is DOUBLE PRECISION array, dimension (N)
The imaginary parts of each scalar alpha defining an
eigenvalue of GNEP.
If ALPHAI(j) is zero, then the j-th eigenvalue is real; if
positive, then the j-th and (j+1)-st eigenvalues are a
complex conjugate pair, with ALPHAI(j+1) = -ALPHAI(j). |
| [out] | BETA | BETA is DOUBLE PRECISION array, dimension (N)
The scalars beta that define the eigenvalues of GNEP.
Together, the quantities alpha = (ALPHAR(j),ALPHAI(j)) and
beta = BETA(j) represent the j-th eigenvalue of the matrix
pair (A,B), in one of the forms lambda = alpha/beta or
mu = beta/alpha. Since either lambda or mu may overflow,
they should not, in general, be computed. |
| [in,out] | QC | QC is DOUBLE PRECISION array, dimension (LDQC, NW) |
| [in] | LDQC | LDQC is INTEGER |
| [in,out] | ZC | ZC is DOUBLE PRECISION array, dimension (LDZC, NW) |
| [in] | LDZC | LDZ is INTEGER |
| [out] | WORK | WORK is DOUBLE PRECISION array, dimension (MAX(1,LWORK))
On exit, if INFO >= 0, WORK(1) returns the optimal LWORK. |
| [in] | LWORK | LWORK is INTEGER
The dimension of the array WORK. LWORK >= max(1,N).
If LWORK = -1, then a workspace query is assumed; the routine
only calculates the optimal size of the WORK array, returns
this value as the first entry of the WORK array, and no error
message related to LWORK is issued by XERBLA. |
| [in] | REC | REC is INTEGER
REC indicates the current recursion level. Should be set
to 0 on first call. |
| [out] | INFO | INFO is INTEGER
= 0: successful exit
< 0: if INFO = -i, the i-th argument had an illegal value |
1.9.4