SLAHQR(3)
π SLAHQR(3) β the LAPACK Hessenberg bestie page π
okay so we're in the Library Functions Manual, bestie, and this is the slahqr.f(3) page β yes, the one whose source path is /home/abui...RC/slahqr.f(3). It's serving all the SLAHQR details, and we are keeping every single flag, dimension, formula, warning, and INFO emergency, no cap π.
πΈ NAME πΈ
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/slahqr.f
That's the whole NAME section: the source-file path. She's literally just the path, bestie ππ.
β¨ SYNOPSIS β¨
π Functions/Subroutines π
subroutine SLAHQR (wantt, wantz, n, ilo, ihi, h, ldh, wr, wi, iloz,
ihiz, z, ldz, info)
SLAHQR computes the eigenvalues and Schur factorization of an upper Hessenberg matrix, using the double-shift/single-shift QR algorithm.
π Function/Subroutine Documentation π
π¦ subroutine SLAHQR π¦
subroutine SLAHQR (logical wantt, logical wantz, integer n, integer ilo,
integer ihi, real, dimension( ldh, * ) h, integer ldh, real, dimension(
* ) wr, real, dimension( * ) wi, integer iloz, integer ihiz, real,
dimension( ldz, * ) z, integer ldz, integer info)
SLAHQR computes the eigenvalues and Schur factorization of an upper Hessenberg matrix, using the double-shift/single-shift QR algorithm.
π Purpose π
SLAHQR is an auxiliary routine called by SHSEQR to update the eigenvalues and Schur decomposition already computed by SHSEQR, by dealing with the Hessenberg submatrix in rows and columns ILO to IHI.
π«Ά Parameters π«Ά
WANTTβWANTT is LOGICAL= .TRUE.: the full Schur formTis required;= .FALSE.: only eigenvalues are required.
WANTZβWANTZ is LOGICAL= .TRUE.: the matrix of Schur vectorsZis required;= .FALSE.: Schur vectors are not required.
NβN is INTEGER- The order of the matrix
H.N >= 0.
- The order of the matrix
ILOβILO is INTEGER- That's the whole description, bestie β no extra lore here π.
IHIβIHI is INTEGER- It is assumed that
His already upper quasi-triangular in rows and columnsIHI+1:N, and thatH(ILO,ILO-1) = 0(unlessILO = 1).SLAHQRworks primarily with the Hessenberg submatrix in rows and columnsILOtoIHI, but applies transformations to all ofHifWANTTis.TRUE.. 1 <= ILO <= max(1,IHI); IHI <= N.
- It is assumed that
HβH is REAL array, dimension (LDH,N)- On entry, the upper Hessenberg matrix
H. - On exit, if
INFOis zero and ifWANTTis.TRUE.,His upper quasi-triangular in rows and columnsILO:IHI, with any 2-by-2 diagonal blocks in standard form. IfINFOis zero andWANTTis.FALSE., the contents ofHare unspecified on exit. The output state ofHifINFOis nonzero is given below under the description ofINFO.
- On entry, the upper Hessenberg matrix
LDHβLDH is INTEGER- The leading dimension of the array
H.LDH >= max(1,N).
- The leading dimension of the array
WRβWR is REAL array, dimension (N)- That's the whole line for
WRon her own, bestie π.
- That's the whole line for
WIβWI is REAL array, dimension (N)- The real and imaginary parts, respectively, of the computed eigenvalues
ILOtoIHIare stored in the corresponding elements ofWRandWI. If two eigenvalues are computed as a complex conjugate pair, they are stored in consecutive elements ofWRandWI, say the i-th and (i+1)th, withWI(i) > 0andWI(i+1) < 0. IfWANTTis.TRUE., the eigenvalues are stored in the same order as on the diagonal of the Schur form returned inH, withWR(i) = H(i,i), and, ifH(i:i+1,i:i+1)is a 2-by-2 diagonal block,WI(i) = sqrt(H(i+1,i)*H(i,i+1))andWI(i+1) = -WI(i).
- The real and imaginary parts, respectively, of the computed eigenvalues
ILOZβILOZ is INTEGER- Again: that's all the page says, bestie π.
IHIZβIHIZ is INTEGER- Specify the rows of
Zto which transformations must be applied ifWANTZis.TRUE.. 1 <= ILOZ <= ILO; IHI <= IHIZ <= N.
- Specify the rows of
ZβZ is REAL array, dimension (LDZ,N)- If
WANTZis.TRUE., on entryZmust contain the current matrixZof transformations accumulated bySHSEQR, and on exitZhas been updated; transformations are applied only to the submatrixZ(ILOZ:IHIZ,ILO:IHI). - If
WANTZis.FALSE.,Zis not referenced.
- If
LDZβLDZ is INTEGER- The leading dimension of the array
Z.LDZ >= max(1,N).
- The leading dimension of the array
INFOβINFO is INTEGER= 0: successful exit> 0: IfINFO = i,SLAHQRfailed to compute all the eigenvaluesILOtoIHIin a total of 30 iterations per eigenvalue; elementsi+1:ihiofWRandWIcontain those eigenvalues which have been successfully computed.- If
INFO > 0andWANTTis.FALSE., then on exit, the remaining unconverged eigenvalues are the eigenvalues of the upper Hessenberg matrix rows and columnsILOthroughINFOof the final, output value ofH. - If
INFO > 0andWANTTis.TRUE., then on exit(*) (initial value of H)*U = U*(final value of H)whereUis an orthogonal matrix. The final value ofHis upper Hessenberg and triangular in rows and columnsINFO+1throughIHI. - If
INFO > 0andWANTZis.TRUE., then on exit(final value of Z) = (initial value of Z)*UwhereUis the orthogonal matrix in(*)(regardless of the value ofWANTT).
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π₯ Further Details: π₯
02-96 Based on modifications by David Day, Sandia National Laboratory, USA
12-04 Further modifications by Ralph Byers, University of Kansas, USA
This is a modified version of SLAHQR from LAPACK version 3.0. It is (1) more robust against overflow and underflow and (2) adopts the more conservative Ahues & Tisseur stopping criterion (LAWN 122, 1997).
π Definition at line 205 of file slahqr.f. π
π€ Author π€
Generated automatically by Doxygen for LAPACK from the source code.
okay that's the whole page, bestie π«Ά β LAPACK Version 3.12.0, source/footer path /home/abui...RC/slahqr.f(3), and the footer doesn't print a separate date, so no invented date energy here β just version + path realness ππ.
π¦ 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