gtcon(3)
okay bestie, welcome to gtcon(3) from the Library Functions Manual π
This page is all about condition number estimates, and it comes in four flavors β c, d, s, and z β because of course it does. Letβs get into it β¨
πΈ NAME πΈ
gtcon β gtcon: condition number estimate π
β¨ SYNOPSIS β¨
π¦ Functions π¦
subroutine cgtcon (norm, n, dl, d, du, du2, ipiv, anorm, rcond, work,
info)
CGTCON
subroutine dgtcon (norm, n, dl, d, du, du2, ipiv, anorm, rcond, work,
iwork, info)
DGTCON
subroutine sgtcon (norm, n, dl, d, du, du2, ipiv, anorm, rcond, work,
iwork, info)
SGTCON
subroutine zgtcon (norm, n, dl, d, du, du2, ipiv, anorm, rcond, work,
info)
ZGTCON
π Detailed Description π
This part is literally just the doorway into the function docs, bestie β no extra side quest here π«
π Function Documentation π
π¦ cgtcon / CGTCON π¦
subroutine cgtcon (character norm, integer n, complex, dimension( * ) dl,
complex, dimension( * ) d, complex, dimension( * ) du, complex,
dimension( * ) du2, integer, dimension( * ) ipiv, real anorm, real
rcond, complex, dimension( * ) work, integer info)
π― Purpose
CGTCON estimates the reciprocal of the condition number of a complex tridiagonal matrix A using the LU factorization as computed by CGTTRF.
An estimate is obtained for norm(inv(A)), and the reciprocal of the condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))). π₯
π§Έ Parameters
NORMβ isCHARACTER*1. It specifies whether the 1-norm condition number or the infinity-norm condition number is required:= '1' or 'O': 1-norm;= 'I': Infinity-norm.
Nβ isINTEGER. The order of the matrixA.N >= 0.DLβ isCOMPLEX array, dimension(N-1). The(n-1)multipliers that define the matrixLfrom the LU factorization ofAas computed byCGTTRF.Dβ isCOMPLEX array, dimension(N). Thendiagonal elements of the upper triangular matrixUfrom the LU factorization ofA.DUβ isCOMPLEX array, dimension(N-1). The(n-1)elements of the first superdiagonal ofU.DU2β isCOMPLEX array, dimension(N-2). The(n-2)elements of the second superdiagonal ofU.IPIVβ isINTEGER array, dimension(N). The pivot indices; for1 <= i <= n, rowiof the matrix was interchanged with rowIPIV(i).IPIV(i)will always be eitheriori+1;IPIV(i) = iindicates a row interchange was not required.ANORMβ isREAL. IfNORM = '1' or 'O', the 1-norm of the original matrixA. IfNORM = 'I', the infinity-norm of the original matrixA.RCONDβ isREAL. The reciprocal of the condition number of the matrixA, computed asRCOND = 1/(ANORM * AINVNM), whereAINVNMis an estimate of the 1-norm ofinv(A)computed in this routine.WORKβ isCOMPLEX array, dimension(2*N).INFOβ isINTEGER:= 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 139 of file cgtcon.f. π
π dgtcon / DGTCON π
subroutine dgtcon (character norm, integer n, double precision, dimension( *
) dl, double precision, dimension( * ) d, double precision, dimension( *
) du, double precision, dimension( * ) du2, integer, dimension( * )
ipiv, double precision anorm, double precision rcond, double precision,
dimension( * ) work, integer, dimension( * ) iwork, integer info)
π― Purpose
DGTCON estimates the reciprocal of the condition number of a real tridiagonal matrix A using the LU factorization as computed by DGTTRF.
An estimate is obtained for norm(inv(A)), and the reciprocal of the condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))). β
π§Έ Parameters
NORMβ isCHARACTER*1. It specifies whether the 1-norm condition number or the infinity-norm condition number is required:= '1' or 'O': 1-norm;= 'I': Infinity-norm.
Nβ isINTEGER. The order of the matrixA.N >= 0.DLβ isDOUBLE PRECISION array, dimension(N-1). The(n-1)multipliers that define the matrixLfrom the LU factorization ofAas computed byDGTTRF.Dβ isDOUBLE PRECISION array, dimension(N). Thendiagonal elements of the upper triangular matrixUfrom the LU factorization ofA.DUβ isDOUBLE PRECISION array, dimension(N-1). The(n-1)elements of the first superdiagonal ofU.DU2β isDOUBLE PRECISION array, dimension(N-2). The(n-2)elements of the second superdiagonal ofU.IPIVβ isINTEGER array, dimension(N). The pivot indices; for1 <= i <= n, rowiof the matrix was interchanged with rowIPIV(i).IPIV(i)will always be eitheriori+1;IPIV(i) = iindicates a row interchange was not required.ANORMβ isDOUBLE PRECISION. IfNORM = '1' or 'O', the 1-norm of the original matrixA. IfNORM = 'I', the infinity-norm of the original matrixA.RCONDβ isDOUBLE PRECISION. The reciprocal of the condition number of the matrixA, computed asRCOND = 1/(ANORM * AINVNM), whereAINVNMis an estimate of the 1-norm ofinv(A)computed in this routine.WORKβ isDOUBLE PRECISION array, dimension(2*N).IWORKβ isINTEGER array, dimension(N).INFOβ isINTEGER:= 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 144 of file dgtcon.f. π
π sgtcon / SGTCON π
subroutine sgtcon (character norm, integer n, real, dimension( * ) dl, real,
dimension( * ) d, real, dimension( * ) du, real, dimension( * ) du2,
integer, dimension( * ) ipiv, real anorm, real rcond, real, dimension( *
) work, integer, dimension( * ) iwork, integer info)
π― Purpose
SGTCON estimates the reciprocal of the condition number of a real tridiagonal matrix A using the LU factorization as computed by SGTTRF.
An estimate is obtained for norm(inv(A)), and the reciprocal of the condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))). π
π§Έ Parameters
NORMβ isCHARACTER*1. It specifies whether the 1-norm condition number or the infinity-norm condition number is required:= '1' or 'O': 1-norm;= 'I': Infinity-norm.
Nβ isINTEGER. The order of the matrixA.N >= 0.DLβ isREAL array, dimension(N-1). The(n-1)multipliers that define the matrixLfrom the LU factorization ofAas computed bySGTTRF.Dβ isREAL array, dimension(N). Thendiagonal elements of the upper triangular matrixUfrom the LU factorization ofA.DUβ isREAL array, dimension(N-1). The(n-1)elements of the first superdiagonal ofU.DU2β isREAL array, dimension(N-2). The(n-2)elements of the second superdiagonal ofU.IPIVβ isINTEGER array, dimension(N). The pivot indices; for1 <= i <= n, rowiof the matrix was interchanged with rowIPIV(i).IPIV(i)will always be eitheriori+1;IPIV(i) = iindicates a row interchange was not required.ANORMβ isREAL. IfNORM = '1' or 'O', the 1-norm of the original matrixA. IfNORM = 'I', the infinity-norm of the original matrixA.RCONDβ isREAL. The reciprocal of the condition number of the matrixA, computed asRCOND = 1/(ANORM * AINVNM), whereAINVNMis an estimate of the 1-norm ofinv(A)computed in this routine.WORKβ isREAL array, dimension(2*N).IWORKβ isINTEGER array, dimension(N).INFOβ isINTEGER:= 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 144 of file sgtcon.f. π
π¦ zgtcon / ZGTCON π¦
subroutine zgtcon (character norm, integer n, complex*16, dimension( * ) dl,
complex*16, dimension( * ) d, complex*16, dimension( * ) du, complex*16,
dimension( * ) du2, integer, dimension( * ) ipiv, double precision
anorm, double precision rcond, complex*16, dimension( * ) work, integer
info)
π― Purpose
ZGTCON estimates the reciprocal of the condition number of a complex tridiagonal matrix A using the LU factorization as computed by ZGTTRF.
An estimate is obtained for norm(inv(A)), and the reciprocal of the condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))). π
π§Έ Parameters
NORMβ part 2 picks up right here, bestie; this is where part 1 ends π«ΆNORMisCHARACTER*1β tells us which condition-number flavor we're serving: the 1-norm condition number or the infinity-norm condition number. The recipe:'1'or'O': 1-norm;'I': infinity-norm.
NisINTEGERβ the order of matrixA.N >= 0, bestie, no negative-order nonsense.DLisCOMPLEX*16 array, dimension (N-1)β the(n-1)multipliers that define the matrixLfrom the LU factorization ofAas computed byZGTTRF. Slay, that's the lower-triangular tea.DisCOMPLEX*16 array, dimension (N)β thendiagonal elements of the upper triangular matrixUfrom the LU factorization ofA.DUisCOMPLEX*16 array, dimension (N-1)β the(n-1)elements of the first superdiagonal ofU.DU2isCOMPLEX*16 array, dimension (N-2)β the(n-2)elements of the second superdiagonal ofU.IPIVisINTEGER array, dimension (N)β the pivot indices; for1 <= i <= n, rowiof the matrix was interchanged with rowIPIV(i).IPIV(i)will always be eitheriori+1;IPIV(i) = iindicates a row interchange was not required. No cap, that's the whole pivot story.ANORMisDOUBLE PRECISIONβ ifNORM = '1'or'O', this is the 1-norm of the original matrixA. IfNORM = 'I', it's the infinity-norm of the original matrixA.RCONDisDOUBLE PRECISIONβ the reciprocal of the condition number of matrixA, computed asRCOND = 1/(ANORM * AINVNM), whereAINVNMis an estimate of the 1-norm ofinv(A)computed in this routine.WORKisCOMPLEX*16 array, dimension (2*N)β workspace, doing its lil behind-the-scenes thing. β¨INFOisINTEGERβ the vibe check:= 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 139 of file zgtcon.f.
π Author
Generated automatically by Doxygen for LAPACK from the source code.
π Okay bestie, that's the whole page: sourceLAPACK, version3.12.0, pagegtcon(3), and the footer didn't hand us a date so we're simply timeless. Bye bestie! π
π¦ 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