CSPDI(3F)CSPDI(3F)NAMECSPDI - CSPDI computes the determinant and inverse of a complex
symmetric matrix using the factors from CSPFA, where the matrix is stored
in packed form.
SYNOPSYS
SUBROUTINE CSPDI(AP,N,KPVT,DET,WORK,JOB)
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.
WORK COMPLEX(N)
work vector. Contents ignored.
JOB INTEGER
JOB has the decimal expansion AB where
if B .NE. 0, the inverse is computed,
if A .NE. 0, the determinant is computed.
For example, JOB = 11 gives both. On Return Variables not requested
by JOB are not used.
AP contains the upper triangle of the inverse of
the original matrix, stored in packed form.
The columns of the upper triangle are stored
sequentially in a one-dimensional array.
DET COMPLEX(2)
determinant of original matrix.
Determinant = DET(1) * 10.0**DET(2)
with 1.0 .LE. ABS(DET(1)) .LT. 10.0
or DET(1) = 0.0. Error Condition
A division by zero will occur if the inverse is requested and CSPCO has
set RCOND .EQ. 0.0 or CSPFA has set INFO .NE. 0 . LINPACK. This
version dated 08/14/78 . James Bunch, Univ. Calif. San Diego, Argonne
Nat. Lab. Subroutines and Functions BLAS CAXPY,CCOPY,CDOTU,CSWAP Fortran
ABS,CMPLX,IABS,MOD,REAL
Page 1