SRC/cgtcon.f(3)
Bestie, welcome to SRC/cgtcon.f(3) β yes, thatβs the header title β and it lives in the Library Functions Manual. Buckle up, because this page is lowkey a whole vibe: itβs the tea on CGTCON, and itβs giving numerical linear algebra excellence. π
β¨
π NAME π
SRC/cgtcon.f
π SYNOPSIS π
π§© Functions/Subroutines π§©
subroutine cgtcon (norm, n, dl, d, du, du2, ipiv, anorm, rcond, work, info)βCGTCON
π οΈ Function/Subroutine Documentation π οΈ
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)
CGTCON
π― Purpose π―
Okay so 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))). Thatβs the whole main-character formula, no cap. π
ποΈ Parameters ποΈ
NORMβNORM is CHARACTER*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βN is INTEGER. This is the order of the matrixA.N >= 0.DLβDL is COMPLEX array, dimension (N-1). These are the(n-1)multipliers that define the matrixLfrom the LU factorization ofAas computed byCGTTRF.DβD is COMPLEX array, dimension (N). These are thendiagonal elements of the upper triangular matrixUfrom the LU factorization ofA.DUβDU is COMPLEX array, dimension (N-1). These are the(n-1)elements of the first superdiagonal ofU.DU2βDU2 is COMPLEX array, dimension (N-2). These are the(n-2)elements of the second superdiagonal ofU.IPIVβIPIV is INTEGER array, dimension (N). These are 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βANORM is REAL. IfNORM = '1' or 'O', this is the 1-norm of the original matrixA. IfNORM = 'I', this is the infinity-norm of the original matrixA.RCONDβRCOND is REAL. This is 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βWORK is COMPLEX array, dimension (2*N).INFOβINFO is INTEGER:= 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.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
Stay reciprocal and radiant, bestie π«Άπ β the footer is giving LAPACK Version 3.12.0 and SRC/cgtcon.f(3), with the date not shown on this page. Periodt. π
β¨
π¦ 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