SGEDI(3F)

NAME
     SGEDI   - SGEDI computes the determinant and inverse of a matrix using
     the factors computed by SGECO or SGEFA.


SYNOPSYS
      SUBROUTINE SGEDI(A,LDA,N,IPVT,DET,WORK,JOB)

DESCRIPTION
     On Entry

     A REAL(LDA, N)
        the output from SGECO or SGEFA.

     LDA INTEGER
        the leading dimension of the array  A .

     N INTEGER
        the order of the matrix  A .

     IPVT INTEGER(N)
        the pivot vector from SGECO or SGEFA.

     WORK REAL(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 REAL(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 SGECO has set RCOND .GT. 0.0 or
     SGEFA has set INFO .EQ. 0 .  LINPACK.  This version dated 08/14/78 .
     Cleve Moler, University of New Mexico, Argonne National Lab.  Subroutines
     and Functions BLAS SAXPY,SSCAL,SSWAP Fortran ABS,MOD