SRC/zpbcon.f(3)
π SRC/zpbcon.f(3) π
okay so this page is from the Library Functions Manual, and itβs giving βI estimate reciprocal condition numbers for complex Hermitian positive definite band matricesβ energy. The name on the cover is SRC/zpbcon.f β periodt.
π NAME π
SRC/zpbcon.f β thatβs the name, bestie. This is the manual page for the ZPBCON subroutine, and it is lowkey all about condition-number estimates. π¦
β¨ SYNOPSIS β¨
π Functions/Subroutines π
subroutine zpbcon (uplo, n, kd, ab, ldab, anorm, rcond, work, rwork,
info)
ZPBCON β yes, thatβs the main character. π
π Function/Subroutine Documentation π
subroutine zpbcon (character uplo, integer n, integer kd, complex*16,
dimension( ldab, * ) ab, integer ldab, double precision anorm, double
precision rcond, complex*16, dimension( * ) work, double precision,
dimension( * ) rwork, integer info)
ZPBCON β say it with your chest. π
πΈ Purpose πΈ
bestie, ZPBCON estimates the reciprocal of the condition number (in the 1-norm) of a complex Hermitian positive definite band matrix using the Cholesky factorization A = U**H*U or A = L*L**H computed by ZPBTRF.
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 glow-up. β¨
π¦ Parameters π¦
UPLOβ isCHARACTER*1.= 'U': Upper triangular factor stored inAB;= 'L': Lower triangular factor stored inAB.
Nβ isINTEGER. The order of the matrixA.N >= 0.KDβ isINTEGER. The number of superdiagonals of the matrixAifUPLO = 'U', or the number of sub-diagonals ifUPLO = 'L'.KD >= 0.ABβ isCOMPLEX*16 array, dimension (LDAB,N). The triangular factorUorLfrom the Cholesky factorizationA = U**H*UorA = L*L**Hof the band matrixA, stored in the firstKD+1rows of the array. Thej-th column ofUorLis stored in thej-th column of the arrayABas follows:- if
UPLO ='U',AB(kd+1+i-j,j) = U(i,j)formax(1,j-kd)<=i<=j; - if
UPLO ='L',AB(1+i-j,j) = L(i,j)forj<=i<=min(n,j+kd).
- if
LDABβ isINTEGER. The leading dimension of the arrayAB.LDAB >= KD+1.ANORMβ isDOUBLE PRECISION. The 1-norm (or infinity-norm) of the Hermitian band 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β isCOMPLEX*16 array, dimension (2*N).RWORKβ isDOUBLE PRECISION 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 131 of file zpbcon.f. π
π§ Author π§
Generated automatically by Doxygen for LAPACK from the source code. no cap, the robots handled the paperwork. π€β¨
π xoxo, your LAPACK bestie β source/version/date vibes from the footer: LAPACK, Version 3.12.0, page SRC/zpbcon.f(3). The footer doesnβt list a date, so Iβm not inventing one, 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