CPOCO(3F)CPOCO(3F)NAMECPOCO - 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.
Page 1
CPOCO(3F)CPOCO(3F)
= 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
Page 2