ZPTCON(3)
πΈ ZPTCON(3) β the LAPACK girlies' guide to condition numbers πΈ
okay so bestie, welcome to the zptcon.f(3) page π
this one lives in the Library Functions Manual and it's about the very specific, very iconic act of asking "how badly-conditioned is my Hermitian positive definite tridiagonal matrix, actually?" π¦
it lives at /home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/zptcon.f, no cap.
β¨ SYNOPSIS β¨
Functions/Subroutines π
subroutine ZPTCON (n, d, e, anorm, rcond, rwork, info)
ZPTCON
π Function/Subroutine Documentation π
subroutine ZPTCON (integer n, double precision, dimension( * ) d,
complex*16, dimension( * ) e, double precision anorm, double precision
rcond, double precision, dimension( * ) rwork, integer info)
ZPTCON
π§ Purpose π§
ok so, literally the whole point:
!>
!> ZPTCON computes the reciprocal of the condition number (in the
!> 1-norm) of a complex Hermitian positive definite tridiagonal matrix
!> using the factorization A = L*D*L**H or A = U**H*D*U computed by
!> ZPTTRF.
!>
!> Norm(inv(A)) is computed by a direct method, and the reciprocal of
!> the condition number is computed as
!> RCOND = 1 / (ANORM * norm(inv(A))).
!>
like, it takes the factorization your girl ZPTTRF already produced and does the β¨norm of the inverseβ¨ thing β direct method, no games β and then flips it into a reciprocal. it's giving measurement, periodt. π
π Parameters π
NβN is INTEGER. the order of the matrix A. and listen:N >= 0. that's the vibe check. π
!> N is INTEGER
!> The order of the matrix A. N >= 0.
!>
DβD is DOUBLE PRECISION array, dimension (N). thendiagonal elements of the diagonal matrix D from the factorization of A, as computed byZPTTRF. it already did the work, bestie, we're just reading it back. π
!> D is DOUBLE PRECISION array, dimension (N)
!> The n diagonal elements of the diagonal matrix D from the
!> factorization of A, as computed by ZPTTRF.
!>
EβE is COMPLEX*16 array, dimension (N-1). the(n-1)off-diagonal elements of the unit bidiagonal factor U or L from the factorization of A, as computed byZPTTRF. note theN-1β it's one shorter, don't come for me, that's just how tridiagonal off-diagonals are. π
!> E is COMPLEX*16 array, dimension (N-1)
!> The (n-1) off-diagonal elements of the unit bidiagonal factor
!> U or L from the factorization of A, as computed by ZPTTRF.
!>
ANORMβANORM is DOUBLE PRECISION. the 1-norm of the original matrix A. you gotta bring this in yourself, it's not computed here, iykyk. π
!> ANORM is DOUBLE PRECISION
!> The 1-norm of the original matrix A.
!>
RCONDβRCOND is DOUBLE PRECISION. the reciprocal of the condition number of the matrix A, computed asRCOND = 1/(ANORM * AINVNM), whereAINVNMis the 1-norm ofinv(A)computed in this routine. this is the output. this is the moment. main character energy fr fr. β
!> RCOND is DOUBLE PRECISION
!> The reciprocal of the condition number of the matrix A,
!> computed as RCOND = 1/(ANORM * AINVNM), where AINVNM is the
!> 1-norm of inv(A) computed in this routine.
!>
RWORKβRWORK is DOUBLE PRECISION array, dimension (N). workspace! give itNdoubles and let it cook. π«Ά
!> RWORK is DOUBLE PRECISION array, dimension (N)
!>
INFOβINFO is INTEGER. the receipts:
!> INFO is INTEGER
!> = 0: successful exit
!> < 0: if INFO = -i, the i-th argument had an illegal value
!>
so 0 means slay, and anything negative means argument number i was a mess. π
π Author π
these legends:
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π₯ Further Details: π₯
!>
!> The method used is described in Nicholas J. Higham, , SIAM J. Sci. Stat. Comput., Vol. 7, No. 1, January 1986.
!>
so if you wanna read the lore, that's the citation. Higham really said "let me explain" and we respect it. π
Definition at line 118 of file zptcon.f. π¦
πΈ Author πΈ
Generated automatically by Doxygen for LAPACK from the source code.
okay and finally, the little sign-off: this whole thing is LAPACK, Version 3.12.0, and the page is /home/abui...RC/zptcon.f(3) β that's from the source/version and date line at the bottom, mwah π go forth and compute your condition numbers, bestie, you're gonna do amazing πβ¨
π¦ 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