CSPSL(3)

SYNOPSIS
SUBROUTINE CSPSL (AP, N, KPVT, B)

NAME
cspsl - Solve a complex symmetric system using the factors obtained from CSPFA.

LIBRARY
SLATEC (LINPACK)

CATEGORY
D2C1

TYPE
COMPLEX (SSPSL-S, DSPSL-D, CHPSL-C, CSPSL-C)

KEYWORDS
LINEAR ALGEBRA, LINPACK, MATRIX, PACKED, SOLVE, SYMMETRIC

AUTHOR
Bunch, J., (UCSD)

DESCRIPTION
CSISL solves the complex symmetric system A * X = B using the factors computed by CSPFA. On Entry AP COMPLEX(N*(N+1)/2) the output from CSPFA. N INTEGER the order of the matrix A . KVPT INTEGER(N) the pivot vector from CSPFA. B COMPLEX(N) the right hand side vector. On Return B the solution vector X . Error Condition A division by zero may occur if CSPCO has set RCOND .EQ. 0.0 or CSPFA has set INFO .NE. 0 . To compute INVERSE(A) * C where C is a matrix with P columns CALL CSPFA(AP,N,KVPT,INFO) IF (INFO .NE. 0) GO TO ... DO 10 J = 1, P CALL CSPSL(AP,N,KVPT,C(1,J)) 10 CONTINUE

REFERENCES
J. J. Dongarra, J. R. Bunch, C. B. Moler, and G. W. Stewart, LINPACK Users' Guide, SIAM, 1979.

ROUTINES CALLED
CAXPY, CDOTU

REVISION HISTORY (YYMMDD)
780814 DATE WRITTEN
890531 Changed all specific intrinsics to generic. (WRB)
890831 Modified array declarations. (WRB)
891107 Corrected category and modified routine equivalence list. (WRB)
891107 REVISION DATE from Version 3.2
891214 Prologue converted to Version 4.0 format. (BAB)
900326 Removed duplicate information from DESCRIPTION section. (WRB)
920501 Reformatted the REFERENCES section. (WRB)