Subroutine */
| 37207 | } /* dpotrf_ */ |
| 37208 | |
| 37209 | /* Subroutine */ int dpotri_(char *uplo, integer *n, doublereal *a, integer * |
| 37210 | lda, integer *info) |
| 37211 | { |
| 37212 | /* System generated locals */ |
| 37213 | integer a_dim1, a_offset, i__1; |
| 37214 | |
| 37215 | /* Local variables */ |
| 37216 | extern logical lsame_(char *, char *); |
| 37217 | extern /* Subroutine */ int xerbla_(char *, integer *), dlauum_( |
| 37218 | char *, integer *, doublereal *, integer *, integer *), |
| 37219 | dtrtri_(char *, char *, integer *, doublereal *, integer *, |
| 37220 | integer *); |
| 37221 | |
| 37222 | |
| 37223 | /* |
| 37224 | -- LAPACK routine (version 3.2) -- |
| 37225 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 37226 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 37227 | November 2006 |
| 37228 | |
| 37229 | |
| 37230 | Purpose |
| 37231 | ======= |
| 37232 | |
| 37233 | DPOTRI computes the inverse of a real symmetric positive definite |
| 37234 | matrix A using the Cholesky factorization A = U**T*U or A = L*L**T |
| 37235 | computed by DPOTRF. |
| 37236 | |
| 37237 | Arguments |
| 37238 | ========= |
| 37239 | |
| 37240 | UPLO (input) CHARACTER*1 |
| 37241 | = 'U': Upper triangle of A is stored; |
| 37242 | = 'L': Lower triangle of A is stored. |
| 37243 | |
| 37244 | N (input) INTEGER |
| 37245 | The order of the matrix A. N >= 0. |
| 37246 | |
| 37247 | A (input/output) DOUBLE PRECISION array, dimension (LDA,N) |
| 37248 | On entry, the triangular factor U or L from the Cholesky |
| 37249 | factorization A = U**T*U or A = L*L**T, as computed by |
| 37250 | DPOTRF. |
| 37251 | On exit, the upper or lower triangle of the (symmetric) |
| 37252 | inverse of A, overwriting the input factor U or L. |
| 37253 | |
| 37254 | LDA (input) INTEGER |
| 37255 | The leading dimension of the array A. LDA >= max(1,N). |
| 37256 | |
| 37257 | INFO (output) INTEGER |
| 37258 | = 0: successful exit |
| 37259 | < 0: if INFO = -i, the i-th argument had an illegal value |
| 37260 | > 0: if INFO = i, the (i,i) element of the factor U or L is |
| 37261 | zero, and the inverse could not be computed. |
| 37262 | |
| 37263 | ===================================================================== |
| 37264 | |
| 37265 | |
| 37266 | Test the input parameters. |