Subroutine */
| 23762 | } /* zpotrf_ */ |
| 23763 | |
| 23764 | /* Subroutine */ int zpotri_(char *uplo, integer *n, doublecomplex *a, |
| 23765 | integer *lda, integer *info) |
| 23766 | { |
| 23767 | /* System generated locals */ |
| 23768 | integer a_dim1, a_offset, i__1; |
| 23769 | |
| 23770 | /* Local variables */ |
| 23771 | extern logical lsame_(char *, char *); |
| 23772 | extern /* Subroutine */ int xerbla_(char *, integer *), zlauum_( |
| 23773 | char *, integer *, doublecomplex *, integer *, integer *), |
| 23774 | ztrtri_(char *, char *, integer *, doublecomplex *, integer *, |
| 23775 | integer *); |
| 23776 | |
| 23777 | |
| 23778 | /* |
| 23779 | -- LAPACK routine (version 3.2) -- |
| 23780 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 23781 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 23782 | November 2006 |
| 23783 | |
| 23784 | |
| 23785 | Purpose |
| 23786 | ======= |
| 23787 | |
| 23788 | ZPOTRI computes the inverse of a complex Hermitian positive definite |
| 23789 | matrix A using the Cholesky factorization A = U**H*U or A = L*L**H |
| 23790 | computed by ZPOTRF. |
| 23791 | |
| 23792 | Arguments |
| 23793 | ========= |
| 23794 | |
| 23795 | UPLO (input) CHARACTER*1 |
| 23796 | = 'U': Upper triangle of A is stored; |
| 23797 | = 'L': Lower triangle of A is stored. |
| 23798 | |
| 23799 | N (input) INTEGER |
| 23800 | The order of the matrix A. N >= 0. |
| 23801 | |
| 23802 | A (input/output) COMPLEX*16 array, dimension (LDA,N) |
| 23803 | On entry, the triangular factor U or L from the Cholesky |
| 23804 | factorization A = U**H*U or A = L*L**H, as computed by |
| 23805 | ZPOTRF. |
| 23806 | On exit, the upper or lower triangle of the (Hermitian) |
| 23807 | inverse of A, overwriting the input factor U or L. |
| 23808 | |
| 23809 | LDA (input) INTEGER |
| 23810 | The leading dimension of the array A. LDA >= max(1,N). |
| 23811 | |
| 23812 | INFO (output) INTEGER |
| 23813 | = 0: successful exit |
| 23814 | < 0: if INFO = -i, the i-th argument had an illegal value |
| 23815 | > 0: if INFO = i, the (i,i) element of the factor U or L is |
| 23816 | zero, and the inverse could not be computed. |
| 23817 | |
| 23818 | ===================================================================== |
| 23819 | |
| 23820 | |
| 23821 | Test the input parameters. |