CPOSL(3F)
CPOSL - CPOSL solves the COMPLEX Hermitian positive definite system A * X = B using the factors computed by CPOCO or CPOFA.
As shipped in IRIX 6.5.19. Unchanged since IRIX 6.5.
NAME CPOSL - CPOSL solves the COMPLEX Hermitian positive definite system A * X = B using the factors computed by CPOCO or CPOFA. SYNOPSYS SUBROUTINE CPOSL(A,LDA,N,B) DESCRIPTION On Entry A COMPLEX(LDA, N) the output from CPOCO or CPOFA. LDA INTEGER the leading dimension of the array A . N INTEGER the order of the matrix A . B COMPLEX(N) the right hand side vector. On Return B the solution vector X . Error Condition A division by zero will occur if the input factor contains a zero on the diagonal. Technically this indicates singularity but it is usually caused by improper subroutine arguments. It will not occur if the subroutines are called correctly and INFO .EQ. 0 . To compute INVERSE(A) * C where C is a matrix with P columns CALL CPOCO(A,LDA,N,RCOND,Z,INFO) IF (RCOND is too small .OR. INFO .NE. 0) GO TO ... DO 10 J = 1, P CALL CPOSL(A,LDA,N,C(1,J)) 10 CONTINUE LINPACK. This version dated 08/14/78 . Cleve Moler, University of New Mexico, Argonne National Lab. Subroutines and Functions BLAS CAXPY,CDOTC