CPPSL(3F)
CPPSL - CPPSL solves the complex Hermitian positive definite system A * X = B using the factors computed by CPPCO or CPPFA.
As shipped in IRIX 6.5.30. Unchanged since IRIX 6.5.
NAME CPPSL - CPPSL solves the complex Hermitian positive definite system A * X = B using the factors computed by CPPCO or CPPFA. SYNOPSYS SUBROUTINE CPPSL(AP,N,B) DESCRIPTION On Entry AP COMPLEX (N*(N+1)/2) the output from CPPCO or CPPFA. 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 CPPCO(AP,N,RCOND,Z,INFO) IF (RCOND is too small .OR. INFO .NE. 0) GO TO ... DO 10 J = 1, P CALL CPPSL(AP,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