ZLAQR3(3)
π ZLAQR3(3) π
okay so bestie, welcome to the Library Functions Manual entry for /home/abui...RC/zlaqr3.f(3) β yes that's the real header, we're not copying it, we're honoring it. This is the page for ZLAQR3, an absolute icon of a subroutine that does aggressive early deflation. no cap. π
π¦ NAME π¦
The file in question, in all its glorious path-ness:
/home/abuild/rpmbuild/BUILD/lapack-
static-3.12.0-build/lapack-3.12.0/SRC/zlaqr3.f
β¨ SYNOPSIS β¨
πΈ Functions/Subroutines πΈ
subroutine ZLAQR3 (wantt, wantz, n, ktop, kbot, nw, h, ldh, iloz, ihiz,
z, ldz, ns, nd, sh, v, ldv, nh, t, ldt, nv, wv, ldwv, work, lwork)
ZLAQR3 performs the unitary similarity transformation of a Hessenberg matrix to detect and deflate fully converged eigenvalues from a trailing principal submatrix (aggressive early deflation). it's giving main character energy, honestly π
π Function/Subroutine Documentation π
subroutine ZLAQR3 (logical wantt, logical wantz, integer n, integer ktop,
integer kbot, integer nw, complex*16, dimension( ldh, * ) h, integer
ldh, integer iloz, integer ihiz, complex*16, dimension( ldz, * ) z,
integer ldz, integer ns, integer nd, complex*16, dimension( * ) sh,
complex*16, dimension( ldv, * ) v, integer ldv, integer nh, complex*16,
dimension( ldt, * ) t, integer ldt, integer nv, complex*16, dimension(
ldwv, * ) wv, integer ldwv, complex*16, dimension( * ) work, integer
lwork)
ZLAQR3 performs the unitary similarity transformation of a Hessenberg matrix to detect and deflate fully converged eigenvalues from a trailing principal submatrix (aggressive early deflation). same energy, second mention, because she deserves it π«Ά
π₯ Purpose π₯
okay so, aggressive early deflation, bestie:
ZLAQR3 accepts as input an upper Hessenberg matrix H and performs an unitary similarity transformation designed to detect and deflate fully converged eigenvalues from a trailing principal submatrix. On output H has been over-written by a new Hessenberg matrix that is a perturbation of an unitary similarity transformation of H. It is to be hoped that the final version of H has many zero subdiagonal entries. π€ like, manifesting that for you fr fr.
π Parameters π
WANTTβWANTTisLOGICAL. If.TRUE., then the Hessenberg matrix H is fully updated so that the triangular Schur factor may be computed (in cooperation with the calling subroutine). If.FALSE., then only enough of H is updated to preserve the eigenvalues. π«WANTZβWANTZisLOGICAL. If.TRUE., then the unitary matrix Z is updated so that the unitary Schur factor may be computed (in cooperation with the calling subroutine). If.FALSE., then Z is not referenced. she's simply not invited to the function πNβNisINTEGER. The order of the matrix H and (ifWANTZis.TRUE.) the order of the unitary matrix Z. πKTOPβKTOPisINTEGER. It is assumed that eitherKTOP = 1orH(KTOP,KTOP-1)=0.KBOTandKTOPtogether determine an isolated block along the diagonal of the Hessenberg matrix. the duo we didn't know we needed π«ΆKBOTβKBOTisINTEGER. It is assumed without a check that eitherKBOT = NorH(KBOT+1,KBOT)=0.KBOTandKTOPtogether determine an isolated block along the diagonal of the Hessenberg matrix. ngl "without a check" is so real of her πNWβNWisINTEGER. Deflation window size.1 <= NW <= (KBOT-KTOP+1). π keep it in range bestie.HβHisCOMPLEX*16array, dimension(LDH,N). On input the initial N-by-N section of H stores the Hessenberg matrix undergoing aggressive early deflation. On output H has been transformed by a unitary similarity transformation, perturbed, and the returned to Hessenberg form that (it is to be hoped) has some zero subdiagonal entries. glow-up arc complete β¨LDHβLDHisINTEGER. Leading dimension of H just as declared in the calling subroutine.N <= LDH. π§ILOZβILOZisINTEGER. that's it, that's the whole description. she said minimal effort, maximum vibe πIHIZβIHIZisINTEGER. Specify the rows of Z to which transformations must be applied ifWANTZis.TRUE..1 <= ILOZ <= IHIZ <= N. π¦ZβZisCOMPLEX*16array, dimension(LDZ,N). IFWANTZis.TRUE., then on output, the unitary similarity transformation mentioned above has been accumulated intoZ(ILOZ:IHIZ,ILOZ:IHIZ)from the right. IfWANTZis.FALSE., then Z is unreferenced. not read, not written, not perceived πLDZβLDZisINTEGER. The leading dimension of Z just as declared in the calling subroutine.1 <= LDZ. βNSβNSisINTEGER. The number of unconverged (ie approximate) eigenvalues returned inSRandSIthat may be used as shifts by the calling subroutine. (yes,SRandSI, we're quoting her exactly, don't come for me π)NDβNDisINTEGER. The number of converged eigenvalues uncovered by this subroutine. found them, they were converged the whole time πSHβSHisCOMPLEX*16array, dimension(KBOT). On output, approximate eigenvalues that may be used for shifts are stored inSH(KBOT-ND-NS+1)throughSR(KBOT-ND). Converged eigenvalues are stored inSH(KBOT-ND+1)throughSH(KBOT). (again with theSR, that's literally what the original says, keeping it verbatim π«Ά)VβVisCOMPLEX*16array, dimension(LDV,NW). An NW-by-NW work array. just vibing in the background πLDVβLDVisINTEGER. The leading dimension of V just as declared in the calling subroutine.NW <= LDV. πNHβNHisINTEGER. The number of columns of T.NH >= NW. πTβTisCOMPLEX*16array, dimension(LDT,NW). no further notes, she's mysterious like that π¦LDTβLDTisINTEGER. The leading dimension of T just as declared in the calling subroutine.NW <= LDT. π§NVβNVisINTEGER. The number of rows of work arrayWVavailable for workspace.NV >= NW. π«WVβWVisCOMPLEX*16array, dimension(LDWV,NW). (that's it, that's all she wrote π )LDWVβLDWVisINTEGER. The leading dimension of W just as declared in the calling subroutine.NW <= LDV. β οΈ note: it really does say "W" and it really does sayLDVβ I'm not touching that, that's canon.WORKβWORKisCOMPLEX*16array, dimension(LWORK). On exit,WORK(1)is set to an estimate of the optimal value ofLWORKfor the given values ofN,NW,KTOPandKBOT. πLWORKβLWORKisINTEGER. The dimension of the work arrayWORK.LWORK = 2*NWsuffices, but greater efficiency may result from larger values ofLWORK.
IfLWORK = -1, then a workspace query is assumed;ZLAQR3only estimates the optimal workspace size for the given values ofN,NW,KTOPandKBOT. The estimate is returned inWORK(1). No error message related toLWORKis issued byXERBLA. Neither H nor Z are accessed. π΅οΈ like, they're not even looked at. iconique.
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Contributors: π
Karen Braman and Ralph Byers, Department of Mathematics, University of Kansas, USA. absolute legends, we stan π«Ά
β Definition at line 264 of file zlaqr3.f. β
π Author π
Generated automatically by Doxygen for LAPACK from the source code. (a bot wrote this, and honestly? it ate π€β¨)
and that's the whole page babes β LAPACK Version 3.12.0, /home/abui...RC/zlaqr3.f(3), no facts dropped, no flags merged, periodt π
π¦ where this page comes from
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms