SUBROUTINE CSPSL(AP,N,KPVT,B) C***BEGIN PROLOGUE CSPSL C***DATE WRITTEN 780814 (YYMMDD) C***REVISION DATE 820801 (YYMMDD) C***CATEGORY NO. D2D1A C***KEYWORDS COMPLEX,LINEAR ALGEBRA,LINPACK,MATRIX,PACKED,SOLVE, C SYMMETRIC C***AUTHOR BUNCH, J., (UCSD) C***PURPOSE Solves the COMPLEX SYMMETRIC system A*X=B using factors C from CSPFA. C***DESCRIPTION C C CSISL solves the complex symmetric system C A * X = B C using the factors computed by CSPFA. C C On Entry C C AP COMPLEX(N*(N+1)/2) C the output from CSPFA. C C N INTEGER C the order of the matrix A . C C KVPT INTEGER(N) C the pivot vector from CSPFA. C C B COMPLEX(N) C the right hand side vector. C C On Return C C B the solution vector X . C C Error Condition C C A division by zero may occur if CSPCO has set RCOND .EQ. 0.0 C or CSPFA has set INFO .NE. 0 . C C To compute INVERSE(A) * C where C is a matrix C with P columns C CALL CSPFA(AP,N,KVPT,INFO) C IF (INFO .NE. 0) GO TO ... C DO 10 J = 1, P C CALL CSPSL(AP,N,KVPT,C(1,J)) C 10 CONTINUE C C LINPACK. This version dated 08/14/78 . C James Bunch, Univ. Calif. San Diego, Argonne Nat. Lab. C C Subroutines and Functions C C BLAS CAXPY,CDOTU C Fortran IABS C***REFERENCES DONGARRA J.J., BUNCH J.R., MOLER C.B., STEWART G.W., C *LINPACK USERS GUIDE*, SIAM, 1979. C***ROUTINES CALLED CAXPY,CDOTU C***END PROLOGUE CSPSL