CPBSL(3F)CPBSL(3F)NAMECPBSL - CPBSL solves the complex Hermitian positive definite band
system A*X = B using the factors computed by CPBCO or CPBFA.
SYNOPSYS
SUBROUTINE CPBSL(ABD,LDA,N,M,B)
DESCRIPTION
On Entry
ABD COMPLEX(LDA, N)
the output from CPBCO or CPBFA.
LDA INTEGER
the leading dimension of the array ABD .
N INTEGER
the order of the matrix A .
M INTEGER
the number of diagonals above the main diagonal.
B COMPLEX(N)
the right hand side vector. On Return
B the solution vector X . Error Condition
A division by zero will occur if the input factor contains a zero on the
diagonal. Technically this indicates singularity but it is usually
caused by improper subroutine arguments. It will not occur if the
subroutines are called correctly and INFO .EQ. 0 . To compute
INVERSE(A) * C where C is a matrix with P columns
CALL CPBCO(ABD,LDA,N,RCOND,Z,INFO)
IF (RCOND is too small .OR. INFO .NE. 0) GO TO ...
DO 10 J = 1, P
CALL CPBSL(ABD,LDA,N,C(1,J)) 10 CONTINUE LINPACK. This version dated
08/14/78 . Cleve Moler, University of New Mexico, Argonne National Lab.
Subroutines and Functions BLAS CAXPY,CDOTC Fortran MIN0
Page 1