SRC/ztbcon.f(3)
π SRC/ztbcon.f
okay so bestie, this is a Fortran subroutine page from LAPACK β it lives in the Library Functions Manual, and it's literally the whole glow-up of SRC/ztbcon.f. no cap, it's giving "reference documentation that ate". let's get into it π
β¨ SYNOPSIS β¨
Functions/Subroutines
subroutine ztbcon (norm, uplo, diag, n, kd, ab, ldab, rcond, work, rwork, info)β that's ZTBCON, the main character π
subroutine ztbcon (norm, uplo, diag, n, kd, ab, ldab, rcond, work, rwork, info)
π¦ Function/Subroutine Documentation π¦
subroutine ztbcon
here's the full signature, character for character, no edits:
subroutine ztbcon (character norm, character uplo, character diag, integer n,
integer kd, complex*16, dimension( ldab, * ) ab, integer ldab,
double precision rcond, complex*16, dimension( * ) work,
double precision, dimension( * ) rwork, integer info)
ZTBCON π
π§ Purpose
so what does this icon actually do? ZTBCON estimates the reciprocal of the condition number of a triangular band matrix A, in either the 1-norm or the infinity-norm. periodt.
the norm of A gets computed and an estimate is obtained for norm(inv(A)), and then the reciprocal of the condition number is computed as:
RCOND = 1 / ( norm(A) * norm(inv(A)) )
that's the whole vibe π₯
πΈ Parameters
NORMβ it'sCHARACTER*1, and it specifies whether you want the 1-norm condition number or the infinity-norm condition number:= '1' or 'O': 1-norm;= 'I': Infinity-norm.
UPLOβ it'sCHARACTER*1:= 'U':Ais upper triangular;= 'L':Ais lower triangular.
DIAGβ it'sCHARACTER*1:= 'N':Ais non-unit triangular;= 'U':Ais unit triangular.
Nβ it'sINTEGER, the order of the matrixA.N >= 0.KDβ it'sINTEGER, the number of superdiagonals or subdiagonals of the triangular band matrixA.KD >= 0.ABβ it'sCOMPLEX*16 array, dimension (LDAB,N). this is the upper or lower triangular band matrixA, stored in the firstkd+1rows of the array. thej-th column ofAis stored in thej-th column of the arrayABas follows:
if UPLO = 'U', AB(kd+1+i-j,j) = A(i,j) for max(1,j-kd)<=i<=j;
if UPLO = 'L', AB(1+i-j,j) = A(i,j) for j<=i<=min(n,j+kd).
and hey β if DIAG = 'U', the diagonal elements of A are not referenced and are assumed to be 1. noted π
LDABβ it'sINTEGER, the leading dimension of the arrayAB.LDAB >= KD+1.RCONDβ it'sDOUBLE PRECISION, the reciprocal of the condition number of the matrixA, computed asRCOND = 1/(norm(A) * norm(inv(A))).WORKβ it'sCOMPLEX*16 array, dimension (2*N).RWORKβ it'sDOUBLE PRECISION array, dimension (N).INFOβ it'sINTEGER:= 0: successful exit β¨< 0: ifINFO = -i, thei-th argument had an illegal value π
π Author
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 141 of file ztbcon.f. π«Ά
π Author (again, but the meta one) π
Generated automatically by Doxygen for LAPACK from the source code. it's giving "the docs wrote themselves", iconic honestly π
that's the whole page bestie, every last flag and formula intact π this one's LAPACK Version 3.12.0, SRC/ztbcon.f(3) β go forth and condition-number responsibly, mwah π
π¦ where this page comes from
- arch rolling: lapack-doc 3.12.1-2 Β· custom
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms