CHPSL(3F) CHPSL(3F)
CHPSL - CHISL solves the complex Hermitian system A * X = B using the
factors computed by CHPFA.
SUBROUTINE CHPSL(AP,N,KPVT,B)
On Entry
AP COMPLEX(N*(N+1)/2)
the output from CHPFA.
N INTEGER
the order of the matrix A . KVPT INTEGER(N)
the pivot vector from CHPFA.
B COMPLEX(N)
the right hand side vector. On Return
B the solution vector X . Error Condition
A division by zero may occur if CHPCO has set RCOND .EQ. 0.0 or CHPFA
has set INFO .NE. 0 . To compute INVERSE(A) * C where C is a matrix
with P columns
CALL CHPFA(AP,N,KVPT,INFO)
IF (INFO .NE. 0) GO TO ...
DO 10 J = 1, P
CALL CHPSL(AP,N,KVPT,C(1,J)) 10 CONTINUE LINPACK. This version dated
08/14/78 . James Bunch, Univ. Calif. San Diego, Argonne Nat. Lab.
Subroutines and Functions BLAS CAXPY,CDOTC Fortran CONJG,IABS
PPPPaaaaggggeeee 1111 [ Back ]
|