Previous: chifa Up: ../linpack.html Next: chpco

CHISL

       SUBROUTINE CHISL(A,LDA,N,KPVT,B)
       INTEGER LDA,N,KPVT(1)
       COMPLEX A(LDA,1),B(1)
 C
 C     CHISL SOLVES THE COMPLEX HERMITIAN SYSTEM
 C     A * X = B
 C     USING THE FACTORS COMPUTED BY CHIFA.
 C
 C     ON ENTRY
 C
 C        A       COMPLEX(LDA,N)
 C                THE OUTPUT FROM CHIFA.
 C
 C        LDA     INTEGER
 C                THE LEADING DIMENSION OF THE ARRAY  A .
 C
 C        N       INTEGER
 C                THE ORDER OF THE MATRIX  A .
 C
 C        KPVT    INTEGER(N)
 C                THE PIVOT VECTOR FROM CHIFA.
 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  CHICO  HAS SET RCOND .EQ. 0.0
 C        OR  CHIFA  HAS SET INFO .NE. 0  .
 C
 C     TO COMPUTE  INVERSE(A) * C  WHERE  C  IS A MATRIX
 C     WITH  P  COLUMNS
 C           CALL CHIFA(A,LDA,N,KPVT,INFO)
 C           IF (INFO .NE. 0) GO TO ...
 C           DO 10 J = 1, P
 C              CALL CHISL(A,LDA,N,KPVT,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,CDOTC
 C     FORTRAN CONJG,IABS