PDA_DGEDI

Determinant and inverse of a matrix using the factors from PDA_DGEFA.

Origin

SLATEC / CAMSUN
        SUBROUTINE PDA_DGEDI (A, LDA, N, IPVT, DET, WORK, JOB)
  
  
  ***BEGIN PROLOGUE  PDA_DGEDI
  ***PURPOSE  Compute the determinant and inverse of a matrix using the
              factors computed by PDA_DGECO or PDA_DGEFA.
  ***LIBRARY   SLATEC (LINPACK)
  ***CATEGORY  D3A1, D2A1
  ***TYPE      DOUBLE PRECISION (SGEDI-S, PDA_DGEDI-D, CGEDI-C)
  ***KEYWORDS  DETERMINANT, INVERSE, LINEAR ALGEBRA, LINPACK, MATRIX
  ***AUTHOR  Moler, C. B., (U. of New Mexico)
  ***DESCRIPTION
  
       PDA_DGEDI computes the determinant and inverse of a matrix
       using the factors computed by PDA_DGECO or PDA_DGEFA.
  
       On Entry
  
          A       DOUBLE PRECISION(LDA, N)
                  the output from PDA_DGECO or PDA_DGEFA.
  
          LDA     INTEGER
                  the leading dimension of the array  A .
  
          N       INTEGER
                  the order of the matrix  A .
  
          IPVT    INTEGER(N)
                  the pivot vector from PDA_DGECO or PDA_DGEFA.
  
          WORK    DOUBLE PRECISION(N)
                  work vector.  Contents destroyed.
  
          JOB     INTEGER
                  = 11   both determinant and inverse.
                  = 01   inverse only.
                  = 10   determinant only.
  
       On Return
  
          A       inverse of original matrix if requested.
                  Otherwise unchanged.
  
          DET     DOUBLE PRECISION(2)
                  determinant of original matrix if requested.
                  Otherwise not referenced.
                  Determinant = DET(1) * 10.0**DET(2)
                  with  1.0 .LE. ABS(DET(1)) .LT. 10.0
                  or  DET(1) .EQ. 0.0 .
  
       Error Condition
  
          A division by zero will occur if the input factor contains
          a zero on the diagonal and the inverse is requested.
          It will not occur if the subroutines are called correctly
          and if PDA_DGECO has set RCOND .GT. 0.0 or PDA_DGEFA has set
          INFO .EQ. 0 .
  
  ***REFERENCES  J. J. Dongarra, J. R. Bunch, C. B. Moler, and G. W.
                   Stewart, LINPACK Users’ Guide, SIAM, 1979.
  ***ROUTINES CALLED  PDA_DAXPY, PDA_DSCAL, PDA_DSWAP
  ***REVISION HISTORY  (YYMMDD)
     780814  DATE WRITTEN
     890531  Changed all specific intrinsics to generic.  (WRB)
     890831  Modified array declarations.  (WRB)
     890831  REVISION DATE from Version 3.2
     891214  Prologue converted to Version 4.0 format.  (BAB)
     900326  Removed duplicate information from DESCRIPTION section.
             (WRB)
     920501  Reformatted the REFERENCES section.  (WRB)
  ***END PROLOGUE  PDA_DGEDI