dlaqz3.f - Man Page

recursive subroutine dlaqz3 (ilschur, ilq, ilz, n, ilo, ihi, nw, a, lda, b, ldb, q, ldq, z, ldz, ns, nd, alphar, alphai, beta, qc, ldqc, zc, ldzc, work, lwork, rec, info)
DLAQZ3

DLAQZ3  

Purpose:

 DLAQZ3 performs AED

Parameters

ILSCHUR

          ILSCHUR is LOGICAL
              Determines whether or not to update the full Schur form

ILQ

          ILQ is LOGICAL
              Determines whether or not to update the matrix Q

ILZ

          ILZ is LOGICAL
              Determines whether or not to update the matrix Z

N

          N is INTEGER
          The order of the matrices A, B, Q, and Z.  N >= 0.

ILO

          ILO is INTEGER

IHI

          IHI is INTEGER
          ILO and IHI mark the rows and columns of (A,B) which
          are to be normalized

NW

          NW is INTEGER
          The desired size of the deflation window.

A

          A is DOUBLE PRECISION array, dimension (LDA, N)

LDA

          LDA is INTEGER
          The leading dimension of the array A.  LDA >= max( 1, N ).

B

          B is DOUBLE PRECISION array, dimension (LDB, N)

LDB

          LDB is INTEGER
          The leading dimension of the array B.  LDB >= max( 1, N ).

Q

          Q is DOUBLE PRECISION array, dimension (LDQ, N)

LDQ

          LDQ is INTEGER

Z

          Z is DOUBLE PRECISION array, dimension (LDZ, N)

LDZ

          LDZ is INTEGER

NS

          NS is INTEGER
          The number of unconverged eigenvalues available to
          use as shifts.

ND

          ND is INTEGER
          The number of converged eigenvalues found.

ALPHAR

          ALPHAR is DOUBLE PRECISION array, dimension (N)
          The real parts of each scalar alpha defining an eigenvalue
          of GNEP.

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).

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.

QC

          QC is DOUBLE PRECISION array, dimension (LDQC, NW)

LDQC

          LDQC is INTEGER

ZC

          ZC is DOUBLE PRECISION array, dimension (LDZC, NW)

LDZC

          LDZ is INTEGER

WORK

          WORK is DOUBLE PRECISION array, dimension (MAX(1,LWORK))
          On exit, if INFO >= 0, WORK(1) returns the optimal LWORK.

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.

REC

          REC is INTEGER
             REC indicates the current recursion level. Should be set
             to 0 on first call.

INFO

          INFO is INTEGER
          = 0: successful exit
          < 0: if INFO = -i, the i-th argument had an illegal value

Author

Thijs Steel, KU Leuven

Date

May 2020

Definition at line 235 of file dlaqz3.f.