CSPDI(3F)							     CSPDI(3F)

NAME    [Toc]    [Back]

     CSPDI   - 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    [Toc]    [Back]

      SUBROUTINE CSPDI(AP,N,KPVT,DET,WORK,JOB)

DESCRIPTION    [Toc]    [Back]

     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


									PPPPaaaaggggeeee 1111