XLPack 7.0
XLPack Numerical Library (C API) Reference Manual
Loading...
Searching...
No Matches

◆ zpocon()

void zpocon ( char  uplo,
int  n,
int  lda,
doublecomplex  a[],
double  anorm,
double *  rcond,
doublecomplex  work[],
double  rwork[],
int *  info 
)

Condition number of a Hermitian positive definite matrix

Purpose
This routine estimates the reciprocal of the condition number (in the 1-norm) of a Hermitian positive definite matrix using the Cholesky factorization A = U^H*U or A = L*L^H computed by zpotrf.

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
[in]uploSpecifies whether the factor U or L is stored.
= 'U': Upper triangular factor U from the Cholesky factorization A = U^H*U.
= 'L': Lower triangular factor L from the Cholesky factorization A = L*L^H.
[in]nOrder of the matrix A. (n >= 0) (If n = 0, returns rcond = 1)
[in]ldaLeading dimension of the two dimensional array a[][]. (lda >= max(1, n))
[in]a[][]Array a[la][lda] (la >= n)
The triangular factor U or L from the Cholesky factorization A = U^H*U or A = L*L^H, as computed by zpotrf.
[in]anormThe 1-norm (or infinity-norm) of the Hermitian matrix A. (anorm >= 0)
[out]rcondThe reciprocal of the condition number of the matrix A, computed as rcond = 1/(anorm * ainvnm), where ainvnm is an estimate of the 1-norm of inv(A) computed in this routine.
[out]work[]Array work[lwork] (lwork >= 2*n)
Work array.
[out]rwork[]Array rwork[lrwork] (lrwork >= n)
Work array.
[out]info= 0: Successful exit
= -1: The argument uplo had an illegal value (uplo != 'U' nor 'L')
= -2: The argument n had an illegal value (n < 0)
= -3: The argument lda had an illegal value (lda < max(1, n))
= -5: The argument anorm had an illegal value (anorm < 0)
Reference
LAPACK