DSTEIN(3)
π
DSTEIN(3) β Library Functions Manual energy π
Okay so this is DSTEIN(3) from the Library Functions Manual (manual name folded in, bestie π«Ά), and the page path in the header/footer is /home/abui...RC/dstein.f(3). This is the group chat where we compute eigenvectors with inverse iteration. Slay. β¨
πΈ NAME πΈ
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/dstein.f
That's the full source path. Iconic. π
β¨ SYNOPSIS β¨
π§ Functions/Subroutines
subroutine DSTEIN (n, d, e, m, w, iblock, isplit, z, ldz, work, iwork,
ifail, info)
DSTEIN
π Function/Subroutine Documentation π
subroutine DSTEIN (integer n, double precision, dimension( * ) d, double
precision, dimension( * ) e, integer m, double precision, dimension( * )
w, integer, dimension( * ) iblock, integer, dimension( * ) isplit,
double precision, dimension( ldz, * ) z, integer ldz, double precision,
dimension( * ) work, integer, dimension( * ) iwork, integer, dimension(
* ) ifail, integer info)
DSTEIN
π Purpose π
Bestie, listen: DSTEIN computes the eigenvectors of a real symmetric tridiagonal matrix T corresponding to specified eigenvalues, using inverse iteration. The maximum number of iterations allowed for each eigenvector is specified by an internal parameter MAXITS (currently set to 5). That's the whole purpose, no cap. π₯
π«Ά Parameters π«Ά
NβNisINTEGER. The order of the matrix.N >= 0.DβDisDOUBLE PRECISION array, dimension (N). Thendiagonal elements of the tridiagonal matrixT.EβEisDOUBLE PRECISION array, dimension (N-1). The(n-1)subdiagonal elements of the tridiagonal matrixT, in elements1toN-1.MβMisINTEGER. The number of eigenvectors to be found.0 <= M <= N.WβWisDOUBLE PRECISION array, dimension (N). The firstMelements ofWcontain the eigenvalues for which eigenvectors are to be computed. The eigenvalues should be grouped by split-off block and ordered from smallest to largest within the block. ( The output arrayWfromDSTEBZwithORDER = 'B'is expected here. )IBLOCKβIBLOCKisINTEGER array, dimension (N). The submatrix indices associated with the corresponding eigenvalues inW;IBLOCK(i)=1if eigenvalueW(i)belongs to the first submatrix from the top,=2ifW(i)belongs to the second submatrix, etc. ( The output arrayIBLOCKfromDSTEBZis expected here. )ISPLITβISPLITisINTEGER array, dimension (N). The splitting points, at whichTbreaks up into submatrices. The first submatrix consists of rows/columns1toISPLIT( 1 ), the second of rows/columnsISPLIT( 1 )+1throughISPLIT( 2 ), etc. ( The output arrayISPLITfromDSTEBZis expected here. )ZβZisDOUBLE PRECISION array, dimension (LDZ, M). The computed eigenvectors. The eigenvector associated with the eigenvalueW(i)is stored in thei-th column ofZ. Any vector which fails to converge is set to its current iterate afterMAXITSiterations.LDZβLDZisINTEGER. The leading dimension of the arrayZ.LDZ >= max(1,N).WORKβWORKisDOUBLE PRECISION array, dimension (5*N).IWORKβIWORKisINTEGER array, dimension (N).IFAILβIFAILisINTEGER array, dimension (M). On normal exit, all elements ofIFAILare zero. If one or more eigenvectors fail to converge afterMAXITSiterations, then their indices are stored in arrayIFAIL.INFOβINFOisINTEGER.= 0: successful exit.< 0: ifINFO = -i, thei-th argument had an illegal value.> 0: ifINFO = i, thenieigenvectors failed to converge inMAXITSiterations. Their indices are stored in arrayIFAIL.
π heads up, bestie: ifINFO > 0, that means some eigenvectors did not converge inMAXITSiterations, and their indices are chilling inIFAIL.
π§ Internal Parameters π§
MAXITSβINTEGER, default =5. The maximum number of iterations performed.EXTRAβINTEGER, default =2. The number of iterations performed after norm growth criterion is satisfied, should be at least1.
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 172 of file dstein.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
π footer tea, bestie:LAPACK,Version 3.12.0, and/home/abui...RC/dstein.f(3). That's the source/version + path sign-off, periodt. Bye! π β¨
π¦ 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