CPOCO(3F)
CPOCO - CPOCO factors a complex Hermitian positive definite matrix and estimates the condition of the matrix.
Showing IRIX 6.5.30 (default release). Unchanged since IRIX 6.5.
NAME CPOCO - CPOCO factors a complex Hermitian positive definite matrix and estimates the condition of the matrix. If RCOND is not needed, CPOFA is slightly faster. To solve A*X = B , follow CPOCO by CPOSL. To compute INVERSE(A)*C , follow CPOCO by CPOSL. To compute DETERMINANT(A) , follow CPOCO by CPODI. To compute INVERSE(A) , follow CPOCO by CPODI. SYNOPSYS SUBROUTINE CPOCO(A,LDA,N,RCOND,Z,INFO) DESCRIPTION On Entry A COMPLEX(LDA, N) the Hermitian matrix to be factored. Only the diagonal and upper triangle are used. LDA INTEGER the leading dimension of the array A . N INTEGER the order of the matrix A . On Return A an upper triangular matrix R so that A = CTRANS(R)*R where CTRANS(R) is the conjugate transpose. The strict lower triangle is unaltered. If INFO .NE. 0 , the factorization is not complete. RCOND REAL an estimate of the reciprocal condition of A . For the system A*X = B , relative perturbations in A and B of size EPSILON may cause relative perturbations in X of size EPSILON/RCOND . If RCOND is so small that the logical expression 1.0 + RCOND .EQ. 1.0 is true, then A may be singular to working precision. In particular, RCOND is zero if exact singularity is detected or the estimate underflows. If INFO .NE. 0 , RCOND is unchanged. Z COMPLEX(N) a work vector whose contents are usually unimportant. If A is close to a singular matrix, then Z is an approximate null vector in the sense that NORM(A*Z) = RCOND*NORM(A)*NORM(Z) . If INFO .NE. 0 , Z is unchanged. INFO INTEGER = 0 for normal return. = K signals an error condition. The leading minor of order K is not positive definite. LINPACK. This version dated 08/14/78 . Cleve Moler, University of New Mexico, Argonne National Lab. Subroutines and Functions LINPACK CPOFA BLAS CAXPY,CDOTC,CSSCAL,SCASUM Fortran ABS,AIMAG,AMAX1,CMPLX,CONJG,REAL