gbcon(3)
π NAME π
okay so, bestie: this is gbcon(3) from the Library Functions Manual, and it is literally gbcon: condition number estimate. That's the name, that's the vibe, periodt. π
β¨ SYNOPSIS β¨
π Functions π
Here's the whole squad of subroutine besties β four of them, all lined up:
subroutine cgbcon (norm, n, kl, ku, ab, ldab, ipiv, anorm, rcond, work,
rwork, info)
CGBCON
subroutine dgbcon (norm, n, kl, ku, ab, ldab, ipiv, anorm, rcond, work,
iwork, info)
DGBCON
subroutine sgbcon (norm, n, kl, ku, ab, ldab, ipiv, anorm, rcond, work,
iwork, info)
SGBCON
subroutine zgbcon (norm, n, kl, ku, ab, ldab, ipiv, anorm, rcond, work,
rwork, info)
ZGBCON
π¦ Detailed Description π¦
This is the detailed-description zone, and it's basically the runway before the function docs walk out. No extra lore here β the real tea is in Function Documentation below. π
π Function Documentation π
π subroutine cgbcon β CGBCON π
subroutine cgbcon (character norm, integer n, integer kl, integer ku,
complex, dimension( ldab, * ) ab, integer ldab, integer, dimension( * )
ipiv, real anorm, real rcond, complex, dimension( * ) work, real,
dimension( * ) rwork, integer info)
CGBCON
π― Purpose π―
CGBCON estimates the reciprocal of the condition number of a complex general band matrix A, in either the 1-norm or the infinity-norm, using the LU factorization computed by CGBTRF.
An estimate is obtained for norm(inv(A)), and the reciprocal of the condition number is computed as:
RCOND = 1 / ( norm(A) * norm(inv(A)) ).
πΈ Parameters πΈ
NORMβNORMisCHARACTER*1. It says whether you want the 1-norm condition number or the infinity-norm condition number:= '1' or 'O': 1-norm;= 'I': Infinity-norm.
NβNisINTEGER. The order of the matrix A.N >= 0.KLβKLisINTEGER. The number of subdiagonals within the band of A.KL >= 0.KUβKUisINTEGER. The number of superdiagonals within the band of A.KU >= 0.ABβABisCOMPLEX array, dimension (LDAB,N). Details of the LU factorization of the band matrix A, as computed byCGBTRF. U is stored as an upper triangular band matrix withKL+KUsuperdiagonals in rows1toKL+KU+1, and the multipliers used during the factorization are stored in rowsKL+KU+2to2*KL+KU+1.LDABβLDABisINTEGER. The leading dimension of the array AB.LDAB >= 2*KL+KU+1.IPIVβIPIVisINTEGER array, dimension (N). The pivot indices; for1 <= i <= N, rowiof the matrix was interchanged with rowIPIV(i).ANORMβANORMisREAL. IfNORM = '1' or 'O', the 1-norm of the original matrix A. IfNORM = 'I', the infinity-norm of the original matrix A.RCONDβRCONDisREAL. The reciprocal of the condition number of the matrix A, computed asRCOND = 1/(norm(A) * norm(inv(A))).WORKβWORKisCOMPLEX array, dimension (2*N).RWORKβRWORKisREAL array, dimension (N).INFOβINFOisINTEGER:= 0: successful exit< 0: if INFO = -i, the i-th argument had an illegal value
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Definition π
Definition at line 145 of file cgbcon.f.
π subroutine dgbcon β DGBCON π
subroutine dgbcon (character norm, integer n, integer kl, integer ku, double
precision, dimension( ldab, * ) ab, integer ldab, integer, dimension( *
) ipiv, double precision anorm, double precision rcond, double
precision, dimension( * ) work, integer, dimension( * ) iwork, integer
info)
DGBCON
π― Purpose π―
DGBCON estimates the reciprocal of the condition number of a real general band matrix A, in either the 1-norm or the infinity-norm, using the LU factorization computed by DGBTRF.
An estimate is obtained for norm(inv(A)), and the reciprocal of the condition number is computed as:
RCOND = 1 / ( norm(A) * norm(inv(A)) ).
πΈ Parameters πΈ
NORMβNORMisCHARACTER*1. It says whether you want the 1-norm condition number or the infinity-norm condition number:= '1' or 'O': 1-norm;= 'I': Infinity-norm.
NβNisINTEGER. The order of the matrix A.N >= 0.KLβKLisINTEGER. The number of subdiagonals within the band of A.KL >= 0.KUβKUisINTEGER. The number of superdiagonals within the band of A.KU >= 0.ABβABisDOUBLE PRECISION array, dimension (LDAB,N). Details of the LU factorization of the band matrix A, as computed byDGBTRF. U is stored as an upper triangular band matrix withKL+KUsuperdiagonals in rows1toKL+KU+1, and the multipliers used during the factorization are stored in rowsKL+KU+2to2*KL+KU+1.LDABβLDABisINTEGER. The leading dimension of the array AB.LDAB >= 2*KL+KU+1.IPIVβIPIVisINTEGER array, dimension (N). The pivot indices; for1 <= i <= N, rowiof the matrix was interchanged with rowIPIV(i).ANORMβANORMisDOUBLE PRECISION. IfNORM = '1' or 'O', the 1-norm of the original matrix A. IfNORM = 'I', the infinity-norm of the original matrix A.RCONDβRCONDisDOUBLE PRECISION. The reciprocal of the condition number of the matrix A, computed asRCOND = 1/(norm(A) * norm(inv(A))).WORKβWORKisDOUBLE PRECISION array, dimension (3*N).IWORKβIWORKisINTEGER array, dimension (N).INFOβINFOisINTEGER:= 0: successful exit< 0: if INFO = -i, the i-th argument had an illegal value
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Definition π
Definition at line 144 of file dgbcon.f.
π subroutine sgbcon β SGBCON π
subroutine sgbcon (character norm, integer n, integer kl, integer ku, real,
dimension( ldab, * ) ab, integer ldab, integer, dimension( * ) ipiv,
real anorm, real rcond, real, dimension( * ) work, integer, dimension( *
) iwork, integer info)
SGBCON
π― Purpose π―
SGBCON estimates the reciprocal of the condition number of a real general band matrix A, in either the 1-norm or the infinity-norm, using the LU factorization computed by SGBTRF.
An estimate is obtained for norm(inv(A)), and the reciprocal of the condition number is computed as:
RCOND = 1 / ( norm(A) * norm(inv(A)) ).
πΈ Parameters πΈ
NORMβNORMisCHARACTER*1. It says whether you want the 1-norm condition number or the infinity-norm condition number:= '1' or 'O': 1-norm;= 'I': Infinity-norm.
NβNisINTEGER. The order of the matrix A.N >= 0.KLβKLisINTEGER. The number of subdiagonals within the band of A.KL >= 0.KUβKUisINTEGER. The number of superdiagonals within the band of A.KU >= 0.ABβABisREAL array, dimension (LDAB,N). Details of the LU factorization of the band matrix A, as computed bySGBTRF. U is stored as an upper triangular band matrix withKL+KUsuperdiagonals in rows1toKL+KU+1, and the multipliers used during the factorization are stored in rowsKL+KU+2to2*KL+KU+1.LDABβLDABisINTEGER. The leading dimension of the array AB.LDAB >= 2*KL+KU+1.IPIVβIPIVisINTEGER array, dimension (N). The pivot indices; for1 <= i <= N, rowiof the matrix was interchanged with rowIPIV(i).ANORMβANORMisREAL. IfNORM = '1' or 'O', the 1-norm of the original matrix A. IfNORM = 'I', the infinity-norm of the original matrix A.RCONDβRCONDisREAL. The reciprocal of the condition number of the matrix A, computed asRCOND = 1/(norm(A) * norm(inv(A))).WORKβWORKisREAL array, dimension (3*N).IWORKβIWORKisINTEGER array, dimension (N).INFOβINFOisINTEGER:= 0: successful exit< 0: if INFO = -i, the i-th argument had an illegal value
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Definition π
Definition at line 144 of file sgbcon.f.
π subroutine zgbcon β ZGBCON π
subroutine zgbcon (character norm, integer n, integer kl, integer ku,
complex*16, dimension( ldab, * ) ab, integer ldab, integer, dimension( *
) ipiv, double precision anorm, double precision rcond, complex*16,
dimension( * ) work, double precision, dimension( * ) rwork, integer
info)
ZGBCON
π― Purpose π―
ZGBCON estimates the reciprocal of the condition number of a complex general band matrix A, in either the 1-norm or the infinity-norm, using the LU factorization computed by ZGBTRF.
An estimate is obtained for norm(inv(A)), and the reciprocal of the condition number is computed as:
RCOND = 1 / ( norm(A) * norm(inv(A)) ).
πΈ Parameters πΈ
NORMβ okay bestie, this is where part 1 literally cuts off. NORM's full description is about to walk in with part 2 πNORMβCHARACTER*1, bestie, listen up: this one picks which condition number we're computing, okay?
= '1' or 'O': 1-norm;
= 'I': Infinity-norm.
NβINTEGER. The order of the matrixA. Gotta beN >= 0β no negative girlies allowed πKLβINTEGER. The number of subdiagonals within the band ofA.KL >= 0, always. That's the rule πKUβINTEGER. The number of superdiagonals within the band ofA. AlsoKU >= 0. Symmetry, we love to see it β¨ABβCOMPLEX*16array, dimension(LDAB,N). This is the whole tea β β the details of the LU factorization of the band matrixA, exactly as computed byZGBTRF. So like:Uis stored as an upper triangular band matrix withKL+KUsuperdiagonals sitting in rows1toKL+KU+1, and then the multipliers used during the factorization are chilling in rowsKL+KU+2to2*KL+KU+1. Bookkeeping but make it iconic πLDABβINTEGER. The leading dimension of the arrayAB. Constraint:LDAB >= 2*KL+KU+1. Respect the boundary pls π§IPIVβINTEGERarray, dimension(N). The pivot indices, and no cap this part actually matters: for1 <= i <= N, rowiof the matrix was interchanged with rowIPIV(i). That's the receipts π§ΎANORMβDOUBLE PRECISION. The norm of the original matrixA, and it depends on your vibe check up top: ifNORM = '1' or 'O', it's the 1-norm of the original matrixA; ifNORM = 'I', it's the infinity-norm of the original matrixA. Match it up correctly or it's giving chaos πRCONDβDOUBLE PRECISION. The reciprocal of the condition number of the matrixA, computed as:
RCOND = 1/(norm(A) * norm(inv(A)))
Literally the main character output of this whole routine, iykyk π
WORKβCOMPLEX*16array, dimension(2*N). Scratch space, doing the background work π«ΆRWORKβDOUBLE PRECISIONarray, dimension(N). More scratch space, real-valued this time πINFOβINTEGER. The final verdict, periodt:
= 0: successful exit
< 0: if INFO = -i, the i-th argument had an illegal value
So if you get a negative, some arg showed up to the function looking a mess π
πΈ Author πΈ
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 145 of file zgbcon.f. Not me knowing exactly where it lives π
π¦ Author π¦
Generated automatically by Doxygen for LAPACK from the source code. The robots wrote the docs, we just made them cute π§
okay that's the whole page bestie, go forth and compute condition numbers responsibly π
Shoutout to LAPACK, Version 3.12.0, and this gbcon(3) page β that's the source and version, and the date lives on the footer, but we're saying bye in style instead. 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