CSPSL(3F)
CSPSL - CSISL solves the complex symmetric system A * X = B using the factors computed by CSPFA.
As shipped in IRIX 6.5.7. Unchanged since IRIX 6.5.
NAME CSPSL - CSISL solves the complex symmetric system A * X = B using the factors computed by CSPFA. SYNOPSYS SUBROUTINE CSPSL(AP,N,KPVT,B) DESCRIPTION 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 LINPACK. This version dated 08/14/78 . James Bunch, Univ. Calif. San Diego, Argonne Nat. Lab. Subroutines and Functions BLAS CAXPY,CDOTU Fortran IABS