home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / CSPDI.z / CSPDI
Encoding:
Text File  |  1998-10-30  |  2.0 KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCSSSSPPPPDDDDIIII((((3333FFFF))))                                                            CCCCSSSSPPPPDDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CSPDI   - CSPDI computes the determinant and inverse of a complex
  10.      symmetric matrix using the factors from CSPFA, where the matrix is stored
  11.      in packed form.
  12.  
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  15.       SUBROUTINE CSPDI(AP,N,KPVT,DET,WORK,JOB)
  16.  
  17. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  18.      On Entry
  19.  
  20.      AAAAPPPP COMPLEX (N*(N+1)/2)
  21.         the output from CSPFA.
  22.  
  23.      NNNN INTEGER
  24.         the order of the matrix A .  KVPT    INTEGER(N)
  25.         the pivot vector from CSPFA.
  26.  
  27.      WWWWOOOORRRRKKKK COMPLEX(N)
  28.         work vector.  Contents ignored.
  29.  
  30.      JJJJOOOOBBBB INTEGER
  31.         JOB has the decimal expansion  AB  where
  32.         if  B .NE. 0, the inverse is computed,
  33.         if  A .NE. 0, the determinant is computed.
  34.         For example, JOB = 11  gives both.  On Return Variables not requested
  35.      by JOB are not used.
  36.  
  37.      AAAAPPPP contains the upper triangle of the inverse of
  38.         the original matrix, stored in packed form.
  39.         The columns of the upper triangle are stored
  40.         sequentially in a one-dimensional array.
  41.  
  42.      DDDDEEEETTTT COMPLEX(2)
  43.         determinant of original matrix.
  44.         Determinant = DET(1) * 10.0**DET(2)
  45.         with 1.0 .LE. ABS(DET(1)) .LT. 10.0
  46.         or DET(1) = 0.0.  Error Condition
  47.  
  48.      AAAA division by zero will occur if the inverse is requested and  CSPCO  has
  49.      set RCOND .EQ. 0.0 or  CSPFA  has set  INFO .NE. 0 .  LINPACK.  This
  50.      version dated 08/14/78 .  James Bunch, Univ. Calif. San Diego, Argonne
  51.      Nat. Lab.  Subroutines and Functions BLAS CAXPY,CCOPY,CDOTU,CSWAP Fortran
  52.      ABS,CMPLX,IABS,MOD,REAL
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.