Subroutine */
| 41330 | } /* strexc_ */ |
| 41331 | |
| 41332 | /* Subroutine */ int strti2_(char *uplo, char *diag, integer *n, real *a, |
| 41333 | integer *lda, integer *info) |
| 41334 | { |
| 41335 | /* System generated locals */ |
| 41336 | integer a_dim1, a_offset, i__1, i__2; |
| 41337 | |
| 41338 | /* Local variables */ |
| 41339 | static integer j; |
| 41340 | static real ajj; |
| 41341 | extern logical lsame_(char *, char *); |
| 41342 | extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *); |
| 41343 | static logical upper; |
| 41344 | extern /* Subroutine */ int strmv_(char *, char *, char *, integer *, |
| 41345 | real *, integer *, real *, integer *), |
| 41346 | xerbla_(char *, integer *); |
| 41347 | static logical nounit; |
| 41348 | |
| 41349 | |
| 41350 | /* |
| 41351 | -- LAPACK routine (version 3.2) -- |
| 41352 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 41353 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 41354 | November 2006 |
| 41355 | |
| 41356 | |
| 41357 | Purpose |
| 41358 | ======= |
| 41359 | |
| 41360 | STRTI2 computes the inverse of a real upper or lower triangular |
| 41361 | matrix. |
| 41362 | |
| 41363 | This is the Level 2 BLAS version of the algorithm. |
| 41364 | |
| 41365 | Arguments |
| 41366 | ========= |
| 41367 | |
| 41368 | UPLO (input) CHARACTER*1 |
| 41369 | Specifies whether the matrix A is upper or lower triangular. |
| 41370 | = 'U': Upper triangular |
| 41371 | = 'L': Lower triangular |
| 41372 | |
| 41373 | DIAG (input) CHARACTER*1 |
| 41374 | Specifies whether or not the matrix A is unit triangular. |
| 41375 | = 'N': Non-unit triangular |
| 41376 | = 'U': Unit triangular |
| 41377 | |
| 41378 | N (input) INTEGER |
| 41379 | The order of the matrix A. N >= 0. |
| 41380 | |
| 41381 | A (input/output) REAL array, dimension (LDA,N) |
| 41382 | On entry, the triangular matrix A. If UPLO = 'U', the |
| 41383 | leading n by n upper triangular part of the array A contains |
| 41384 | the upper triangular matrix, and the strictly lower |
| 41385 | triangular part of A is not referenced. If UPLO = 'L', the |
| 41386 | leading n by n lower triangular part of the array A contains |
| 41387 | the lower triangular matrix, and the strictly upper |
| 41388 | triangular part of A is not referenced. If DIAG = 'U', the |
| 41389 | diagonal elements of A are also not referenced and are |