CSPSL(3F)

CSPSL - CSISL solves the complex symmetric system A * X = B using the factors computed by CSPFA.

As shipped in IRIX 6.5.22. 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