dpocon dpotrf dpotri dpotrs dppcon dpptrf dpptri dpptrs dsposv

dppcon() void dppcon ( char  uplo, int  n, double  ap[], double  anorm, double *  rcond, double  work[], int  iwork[], int *  info  ) Condition number of a symmetric positive definite matrix in packed form Purpose This routine estimates the reciprocal of the condition number (in the 1-norm) of a real symmetric positive definite packed matrix using the Cholesky factorization A = U^T*U or A = L*L^T computed by dpptrf. 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] uplo Specifies whether the factor U or L is stored. = 'U': Upper triangular factor U from the Cholesky factorization A = U^T*U. = 'L': Lower triangular factor L from the Cholesky factorization A = L*L^T. [in] n Order of the matrix A. (n >= 0) (If n = 0, returns rcond = 1) [in] ap[] Array ap[lap] (lap >= n(n + 1)/2) The triangular factor U or L in packed form from the Cholesky factorization A = U^T*U or A = L*L^T, as computed by dpptrf. [in] anorm The 1-norm (or infinity-norm) of the symmetric matrix A. (anorm >= 0) [out] rcond The 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 >= 3*n) Work array. [out] iwork[] Array iwork[liwork] (liwork >= 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) = -4: The argument anorm had an illegal value (anorm < 0) Reference LAPACK