CSISL(3F)
CSISL - CSISL solves the complex symmetric system A * X = B using the factors computed by CSIFA.
As shipped in IRIX 6.5.5. Unchanged since IRIX 6.5.
NAME CSISL - CSISL solves the complex symmetric system A * X = B using the factors computed by CSIFA. SYNOPSYS SUBROUTINE CSISL(A,LDA,N,KPVT,B) DESCRIPTION On Entry A COMPLEX(LDA,N) the output from CSIFA. LDA INTEGER the leading dimension of the array A . N INTEGER the order of the matrix A . KVPT INTEGER(N) the pivot vector from CSIFA. B COMPLEX(N) the right hand side vector. On Return B the solution vector X . Error Condition A division by zero may occur if CSICO has set RCOND .EQ. 0.0 or CSIFA has set INFO .NE. 0 . To compute INVERSE(A) * C where C is a matrix with P columns CALL CSIFA(A,LDA,N,KVPT,INFO) If (INFO .NE. 0) GO TO ... DO 10 J = 1, P CALL CSISL(A,LDA,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