laqr2#
Functions
-
void slaqr2(const INT wantt, const INT wantz, const INT n, const INT ktop, const INT kbot, const INT nw, f32 *H, const INT ldh, const INT iloz, const INT ihiz, f32 *Z, const INT ldz, INT *ns, INT *nd, f32 *sr, f32 *si, f32 *V, const INT ldv, const INT nh, f32 *T, const INT ldt, const INT nv, f32 *WV, const INT ldwv, f32 *work, const INT lwork)#
SLAQR2 performs the orthogonal similarity transformation of a Hessenberg matrix to detect and deflate fully converged eigenvalues from a trailing principal submatrix (aggressive early deflation).
SLAQR2 is identical to SLAQR3 except that it avoids recursion by calling SLAHQR instead of SLAQR4.
Parameters
inwanttIf nonzero, the Hessenberg matrix H is fully updated.
inwantzIf nonzero, the orthogonal matrix Z is updated.
innThe order of the matrix H. n >= 0.
inktopFirst row/column of isolated block (0-based).
inkbotLast row/column of isolated block (0-based).
innwDeflation window size. 1 <= nw <= (kbot - ktop + 1).
inoutHDouble precision array, dimension (ldh, n).
inldhLeading dimension of H. ldh >= n.
inilozFirst row of Z to update (0-based).
inihizLast row of Z to update (0-based).
inoutZDouble precision array, dimension (ldz, n).
inldzLeading dimension of Z. ldz >= 1.
outnsNumber of unconverged eigenvalues (shifts).
outndNumber of converged (deflated) eigenvalues.
outsrDouble precision array, dimension (kbot+1). Real parts.
outsiDouble precision array, dimension (kbot+1). Imaginary parts.
outVDouble precision array, dimension (ldv, nw).
inldvLeading dimension of V. ldv >= nw.
innhNumber of columns of T. nh >= nw.
outTDouble precision array, dimension (ldt, nw).
inldtLeading dimension of T. ldt >= nw.
innvNumber of rows of WV. nv >= nw.
outWVDouble precision array, dimension (ldwv, nw).
inldwvLeading dimension of WV. ldwv >= nv.
outworkDouble precision array, dimension (lwork).
inlworkDimension of work array. lwork >= 2*nw. If lwork = -1, workspace query is assumed.
void slaqr2(
const INT wantt,
const INT wantz,
const INT n,
const INT ktop,
const INT kbot,
const INT nw,
f32* H,
const INT ldh,
const INT iloz,
const INT ihiz,
f32* Z,
const INT ldz,
INT* ns,
INT* nd,
f32* sr,
f32* si,
f32* V,
const INT ldv,
const INT nh,
f32* T,
const INT ldt,
const INT nv,
f32* WV,
const INT ldwv,
f32* work,
const INT lwork
);
Functions
-
void dlaqr2(const INT wantt, const INT wantz, const INT n, const INT ktop, const INT kbot, const INT nw, f64 *H, const INT ldh, const INT iloz, const INT ihiz, f64 *Z, const INT ldz, INT *ns, INT *nd, f64 *sr, f64 *si, f64 *V, const INT ldv, const INT nh, f64 *T, const INT ldt, const INT nv, f64 *WV, const INT ldwv, f64 *work, const INT lwork)#
DLAQR2 performs the orthogonal similarity transformation of a Hessenberg matrix to detect and deflate fully converged eigenvalues from a trailing principal submatrix (aggressive early deflation).
DLAQR2 is identical to DLAQR3 except that it avoids recursion by calling DLAHQR instead of DLAQR4.
Parameters
inwanttIf nonzero, the Hessenberg matrix H is fully updated.
inwantzIf nonzero, the orthogonal matrix Z is updated.
innThe order of the matrix H. n >= 0.
inktopFirst row/column of isolated block (0-based).
inkbotLast row/column of isolated block (0-based).
innwDeflation window size. 1 <= nw <= (kbot - ktop + 1).
inoutHDouble precision array, dimension (ldh, n).
inldhLeading dimension of H. ldh >= n.
inilozFirst row of Z to update (0-based).
inihizLast row of Z to update (0-based).
inoutZDouble precision array, dimension (ldz, n).
inldzLeading dimension of Z. ldz >= 1.
outnsNumber of unconverged eigenvalues (shifts).
outndNumber of converged (deflated) eigenvalues.
outsrDouble precision array, dimension (kbot+1). Real parts.
outsiDouble precision array, dimension (kbot+1). Imaginary parts.
outVDouble precision array, dimension (ldv, nw).
inldvLeading dimension of V. ldv >= nw.
innhNumber of columns of T. nh >= nw.
outTDouble precision array, dimension (ldt, nw).
inldtLeading dimension of T. ldt >= nw.
innvNumber of rows of WV. nv >= nw.
outWVDouble precision array, dimension (ldwv, nw).
inldwvLeading dimension of WV. ldwv >= nv.
outworkDouble precision array, dimension (lwork).
inlworkDimension of work array. lwork >= 2*nw. If lwork = -1, workspace query is assumed.
void dlaqr2(
const INT wantt,
const INT wantz,
const INT n,
const INT ktop,
const INT kbot,
const INT nw,
f64* H,
const INT ldh,
const INT iloz,
const INT ihiz,
f64* Z,
const INT ldz,
INT* ns,
INT* nd,
f64* sr,
f64* si,
f64* V,
const INT ldv,
const INT nh,
f64* T,
const INT ldt,
const INT nv,
f64* WV,
const INT ldwv,
f64* work,
const INT lwork
);
Functions
-
void claqr2(const INT wantt, const INT wantz, const INT n, const INT ktop, const INT kbot, const INT nw, c64 *H, const INT ldh, const INT iloz, const INT ihiz, c64 *Z, const INT ldz, INT *ns, INT *nd, c64 *SH, c64 *V, const INT ldv, const INT nh, c64 *T, const INT ldt, const INT nv, c64 *WV, const INT ldwv, c64 *work, const INT lwork)#
CLAQR2 performs the unitary similarity transformation of a Hessenberg matrix to detect and deflate fully converged eigenvalues from a trailing principal submatrix (aggressive early deflation).
CLAQR2 is identical to CLAQR3 except that it avoids recursion by calling CLAHQR instead of CLAQR4.
Parameters
inwanttIf nonzero, the Hessenberg matrix H is fully updated.
inwantzIf nonzero, the unitary matrix Z is updated.
innThe order of the matrix H. n >= 0.
inktopFirst row/column of isolated block (0-based).
inkbotLast row/column of isolated block (0-based).
innwDeflation window size. 1 <= nw <= (kbot - ktop + 1).
inoutHComplex array, dimension (ldh, n).
inldhLeading dimension of H. ldh >= n.
inilozFirst row of Z to update (0-based).
inihizLast row of Z to update (0-based).
inoutZComplex array, dimension (ldz, n).
inldzLeading dimension of Z. ldz >= 1.
outnsNumber of unconverged eigenvalues (shifts).
outndNumber of converged (deflated) eigenvalues.
outSHComplex array, dimension (kbot+1). Eigenvalues/shifts.
outVComplex array, dimension (ldv, nw).
inldvLeading dimension of V. ldv >= nw.
innhNumber of columns of T. nh >= nw.
outTComplex array, dimension (ldt, nw).
inldtLeading dimension of T. ldt >= nw.
innvNumber of rows of WV. nv >= nw.
outWVComplex array, dimension (ldwv, nw).
inldwvLeading dimension of WV. ldwv >= nv.
outworkComplex array, dimension (lwork).
inlworkDimension of work array. lwork >= 2*nw. If lwork = -1, workspace query is assumed.
void claqr2(
const INT wantt,
const INT wantz,
const INT n,
const INT ktop,
const INT kbot,
const INT nw,
c64* H,
const INT ldh,
const INT iloz,
const INT ihiz,
c64* Z,
const INT ldz,
INT* ns,
INT* nd,
c64* SH,
c64* V,
const INT ldv,
const INT nh,
c64* T,
const INT ldt,
const INT nv,
c64* WV,
const INT ldwv,
c64* work,
const INT lwork
);
Functions
-
void zlaqr2(const INT wantt, const INT wantz, const INT n, const INT ktop, const INT kbot, const INT nw, c128 *H, const INT ldh, const INT iloz, const INT ihiz, c128 *Z, const INT ldz, INT *ns, INT *nd, c128 *SH, c128 *V, const INT ldv, const INT nh, c128 *T, const INT ldt, const INT nv, c128 *WV, const INT ldwv, c128 *work, const INT lwork)#
ZLAQR2 performs the unitary similarity transformation of a Hessenberg matrix to detect and deflate fully converged eigenvalues from a trailing principal submatrix (aggressive early deflation).
ZLAQR2 is identical to ZLAQR3 except that it avoids recursion by calling ZLAHQR instead of ZLAQR4.
Parameters
inwanttIf nonzero, the Hessenberg matrix H is fully updated.
inwantzIf nonzero, the unitary matrix Z is updated.
innThe order of the matrix H. n >= 0.
inktopFirst row/column of isolated block (0-based).
inkbotLast row/column of isolated block (0-based).
innwDeflation window size. 1 <= nw <= (kbot - ktop + 1).
inoutHComplex array, dimension (ldh, n).
inldhLeading dimension of H. ldh >= n.
inilozFirst row of Z to update (0-based).
inihizLast row of Z to update (0-based).
inoutZComplex array, dimension (ldz, n).
inldzLeading dimension of Z. ldz >= 1.
outnsNumber of unconverged eigenvalues (shifts).
outndNumber of converged (deflated) eigenvalues.
outSHComplex array, dimension (kbot+1). Eigenvalues/shifts.
outVComplex array, dimension (ldv, nw).
inldvLeading dimension of V. ldv >= nw.
innhNumber of columns of T. nh >= nw.
outTComplex array, dimension (ldt, nw).
inldtLeading dimension of T. ldt >= nw.
innvNumber of rows of WV. nv >= nw.
outWVComplex array, dimension (ldwv, nw).
inldwvLeading dimension of WV. ldwv >= nv.
outworkComplex array, dimension (lwork).
inlworkDimension of work array. lwork >= 2*nw. If lwork = -1, workspace query is assumed.
void zlaqr2(
const INT wantt,
const INT wantz,
const INT n,
const INT ktop,
const INT kbot,
const INT nw,
c128* H,
const INT ldh,
const INT iloz,
const INT ihiz,
c128* Z,
const INT ldz,
INT* ns,
INT* nd,
c128* SH,
c128* V,
const INT ldv,
const INT nh,
c128* T,
const INT ldt,
const INT nv,
c128* WV,
const INT ldwv,
c128* work,
const INT lwork
);