Previous: sppco Up: ../linpack.html Next: sppfa
SUBROUTINE SPPDI(AP,N,DET,JOB)
INTEGER N,JOB
REAL AP(1)
REAL DET(2)
C
C SPPDI COMPUTES THE DETERMINANT AND INVERSE
C OF A REAL SYMMETRIC POSITIVE DEFINITE MATRIX
C USING THE FACTORS COMPUTED BY SPPCO OR SPPFA .
C
C ON ENTRY
C
C AP REAL (N*(N+1)/2)
C THE OUTPUT FROM SPPCO OR SPPFA.
C
C N INTEGER
C THE ORDER OF THE MATRIX A .
C
C JOB INTEGER
C = 11 BOTH DETERMINANT AND INVERSE.
C = 01 INVERSE ONLY.
C = 10 DETERMINANT ONLY.
C
C ON RETURN
C
C AP THE UPPER TRIANGULAR HALF OF THE INVERSE .
C THE STRICT LOWER TRIANGLE IS UNALTERED.
C
C DET REAL(2)
C DETERMINANT OF ORIGINAL MATRIX IF REQUESTED.
C OTHERWISE NOT REFERENCED.
C DETERMINANT = DET(1) * 10.0**DET(2)
C WITH 1.0 .LE. DET(1) .LT. 10.0
C OR DET(1) .EQ. 0.0 .
C
C ERROR CONDITION
C
C A DIVISION BY ZERO WILL OCCUR IF THE INPUT FACTOR CONTAINS
C A ZERO ON THE DIAGONAL AND THE INVERSE IS REQUESTED.
C IT WILL NOT OCCUR IF THE SUBROUTINES ARE CALLED CORRECTLY
C AND IF SPOCO OR SPOFA HAS SET INFO .EQ. 0 .
C
C LINPACK. THIS VERSION DATED 08/14/78 .
C CLEVE MOLER, UNIVERSITY OF NEW MEXICO, ARGONNE NATIONAL LAB.
C
C SUBROUTINES AND FUNCTIONS
C
C BLAS SAXPY,SSCAL
C FORTRAN MOD
C