CLAQR5(3)
π CLAQR5(3) β Library Functions Manual π
Manual: Library Functions Manual Β· page path: /home/abui...RC/claqr5.f(3)
πΈ NAME πΈ
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/claqr5.f
Yes bestie, that's the source file path this page is serving. Keep it cute, keep it exact. π
β¨ SYNOPSIS β¨
π Functions/Subroutines π
subroutine CLAQR5 (wantt, wantz, kacc22, n, ktop, kbot, nshfts, s, h, ldh, iloz, ihiz, z, ldz, v, ldv, u, ldu, nv, wv, ldwv, nh, wh, ldwh)β CLAQR5 performs a single small-bulge multi-shift QR sweep.
π¦ Function/Subroutine Documentation π¦
π subroutine CLAQR5 π
subroutine CLAQR5 (logical wantt, logical wantz, integer kacc22, integer n, integer ktop, integer kbot, integer nshfts, complex, dimension( * ) s, complex, dimension( ldh, * ) h, integer ldh, integer iloz, integer ihiz, complex, dimension( ldz, * ) z, integer ldz, complex, dimension( ldv, * ) v, integer ldv, complex, dimension( ldu, * ) u, integer ldu, integer nv, complex, dimension( ldwv, * ) wv, integer ldwv, integer nh, complex, dimension( ldwh, * ) wh, integer ldwh)
CLAQR5 performs a single small-bulge multi-shift QR sweep. That's the whole vibe of this subroutine, and honestly? Very main character energy. β¨
π― Purpose π―
CLAQR5 called by CLAQR0 performs a single small-bulge multi-shift QR sweep. Periodt. π
π Parameters π
WANTTβWANTTis LOGICAL.WANTT = .true.if the triangular Schur factor is being computed.WANTTis set to.false.otherwise.WANTZβWANTZis LOGICAL.WANTZ = .true.if the unitary Schur factor is being computed.WANTZis set to.false.otherwise.KACC22βKACC22is INTEGER with value0,1, or2. Specifies the computation mode of far-from-diagonal orthogonal updates.= 0:CLAQR5does not accumulate reflections and does not use matrix-matrix multiply to update far-from-diagonal matrix entries.= 1:CLAQR5accumulates reflections and uses matrix-matrix multiply to update the far-from-diagonal matrix entries.= 2: Same asKACC22 = 1. This option used to enable exploiting the 2-by-2 structure during matrix multiplications, but this is no longer supported.
NβNis INTEGER.Nis the order of the Hessenberg matrixHupon which this subroutine operates.KTOPβKTOPis INTEGER.KBOTβKBOTis INTEGER. These are the first and last rows and columns of an isolated diagonal block upon which the QR sweep is to be applied.
It is assumed without a check that eitherKTOP = 1orH(KTOP,KTOP-1) = 0, and eitherKBOT = NorH(KBOT+1,KBOT) = 0.
NSHFTSβNSHFTSis INTEGER.NSHFTSgives the number of simultaneous shifts.NSHFTSmust be positive and even.SβSis COMPLEX array, dimension(NSHFTS).Scontains the shifts of origin that define the multi-shift QR sweep. On outputSmay be reordered.HβHis COMPLEX array, dimension(LDH,N). On inputHcontains a Hessenberg matrix. On output a multi-shift QR sweep with shiftsSR(J)+i*SI(J)is applied to the isolated diagonal block in rows and columnsKTOPthroughKBOT.LDHβLDHis INTEGER.LDHis the leading dimension ofHjust as declared in the calling procedure.LDH >= MAX(1,N).ILOZβILOZis INTEGER.IHIZβIHIZis INTEGER. Specify the rows ofZto which transformations must be applied ifWANTZis.TRUE..1 <= ILOZ <= IHIZ <= N.ZβZis COMPLEX array, dimension(LDZ,IHIZ). IfWANTZ = .TRUE., then the QR Sweep unitary similarity transformation is accumulated intoZ(ILOZ:IHIZ,ILOZ:IHIZ)from the right. IfWANTZ = .FALSE., thenZis unreferenced.LDZβLDZis INTEGER. The original literally saysLDAis the leading dimension ofZjust as declared in the calling procedure β keeping that exact wording, bestie.LDZ >= N.VβVis COMPLEX array, dimension(LDV,NSHFTS/2).LDVβLDVis INTEGER.LDVis the leading dimension ofVas declared in the calling procedure.LDV >= 3.UβUis COMPLEX array, dimension(LDU,2*NSHFTS).LDUβLDUis INTEGER.LDUis the leading dimension ofUjust as declared in the calling subroutine.LDU >= 2*NSHFTS.NVβNVis INTEGER.NVis the number of rows inWVagailable for workspace. (Yes, the page says βagailableβ β we're keeping it π)NV >= 1.WVβWVis COMPLEX array, dimension(LDWV,2*NSHFTS).LDWVβLDWVis INTEGER.LDWVis the leading dimension ofWVas declared in the calling subroutine.LDWV >= NV.NHβNHis INTEGER.NHis the number of columns in arrayWHavailable for workspace.NH >= 1.WHβWHis COMPLEX array, dimension(LDWH,NH).LDWHβLDWHis INTEGER. Leading dimension ofWHjust as declared in the calling procedure.LDWH >= 2*NSHFTS.
π 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
- Lars Karlsson, Daniel Kressner, and Bruno Lang
- Thijs Steel, Department of Computer science, KU Leuven, Belgium
π References π
- K. Braman, R. Byers and R. Mathias, The Multi-Shift QR Algorithm Part I: Maintaining Well Focused Shifts, and Level 3 Performance, SIAM Journal of Matrix Analysis, volume 23, pages 929--947, 2002.
- Lars Karlsson, Daniel Kressner, and Bruno Lang, Optimally packed chains of bulges in multishift QR algorithms. ACM Trans. Math. Softw. 40, 2, Article 12 (February 2014).
Definition at line 254 of file claqr5.f.
ποΈ Author (the auto-gen one) ποΈ
Generated automatically by Doxygen for LAPACK from the source code.
Sign-off: LAPACK Version 3.12.0 Β· page path /home/abui...RC/claqr5.f(3). No date is shown in this footer, so I'm not inventing one, bestie. Mwah ππ
π¦ 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