Subroutine */
| 41655 | } /* dtrti2_ */ |
| 41656 | |
| 41657 | /* Subroutine */ int dtrtri_(char *uplo, char *diag, integer *n, doublereal * |
| 41658 | a, integer *lda, integer *info) |
| 41659 | { |
| 41660 | /* System generated locals */ |
| 41661 | address a__1[2]; |
| 41662 | integer a_dim1, a_offset, i__1, i__2[2], i__3, i__4, i__5; |
| 41663 | char ch__1[2]; |
| 41664 | |
| 41665 | /* Local variables */ |
| 41666 | static integer j, jb, nb, nn; |
| 41667 | extern logical lsame_(char *, char *); |
| 41668 | extern /* Subroutine */ int dtrmm_(char *, char *, char *, char *, |
| 41669 | integer *, integer *, doublereal *, doublereal *, integer *, |
| 41670 | doublereal *, integer *), dtrsm_( |
| 41671 | char *, char *, char *, char *, integer *, integer *, doublereal * |
| 41672 | , doublereal *, integer *, doublereal *, integer *); |
| 41673 | static logical upper; |
| 41674 | extern /* Subroutine */ int dtrti2_(char *, char *, integer *, doublereal |
| 41675 | *, integer *, integer *), xerbla_(char *, integer |
| 41676 | *); |
| 41677 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 41678 | integer *, integer *, ftnlen, ftnlen); |
| 41679 | static logical nounit; |
| 41680 | |
| 41681 | |
| 41682 | /* |
| 41683 | -- LAPACK routine (version 3.2) -- |
| 41684 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 41685 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 41686 | November 2006 |
| 41687 | |
| 41688 | |
| 41689 | Purpose |
| 41690 | ======= |
| 41691 | |
| 41692 | DTRTRI computes the inverse of a real upper or lower triangular |
| 41693 | matrix A. |
| 41694 | |
| 41695 | This is the Level 3 BLAS version of the algorithm. |
| 41696 | |
| 41697 | Arguments |
| 41698 | ========= |
| 41699 | |
| 41700 | UPLO (input) CHARACTER*1 |
| 41701 | = 'U': A is upper triangular; |
| 41702 | = 'L': A is lower triangular. |
| 41703 | |
| 41704 | DIAG (input) CHARACTER*1 |
| 41705 | = 'N': A is non-unit triangular; |
| 41706 | = 'U': A is unit triangular. |
| 41707 | |
| 41708 | N (input) INTEGER |
| 41709 | The order of the matrix A. N >= 0. |
| 41710 | |
| 41711 | A (input/output) DOUBLE PRECISION array, dimension (LDA,N) |
| 41712 | On entry, the triangular matrix A. If UPLO = 'U', the |
| 41713 | leading N-by-N upper triangular part of the array A contains |
| 41714 | the upper triangular matrix, and the strictly lower |