Subroutine */
| 33004 | } /* dlauu2_ */ |
| 33005 | |
| 33006 | /* Subroutine */ int dlauum_(char *uplo, integer *n, doublereal *a, integer * |
| 33007 | lda, integer *info) |
| 33008 | { |
| 33009 | /* System generated locals */ |
| 33010 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4; |
| 33011 | |
| 33012 | /* Local variables */ |
| 33013 | static integer i__, ib, nb; |
| 33014 | extern /* Subroutine */ int dgemm_(char *, char *, integer *, integer *, |
| 33015 | integer *, doublereal *, doublereal *, integer *, doublereal *, |
| 33016 | integer *, doublereal *, doublereal *, integer *); |
| 33017 | extern logical lsame_(char *, char *); |
| 33018 | extern /* Subroutine */ int dtrmm_(char *, char *, char *, char *, |
| 33019 | integer *, integer *, doublereal *, doublereal *, integer *, |
| 33020 | doublereal *, integer *); |
| 33021 | static logical upper; |
| 33022 | extern /* Subroutine */ int dsyrk_(char *, char *, integer *, integer *, |
| 33023 | doublereal *, doublereal *, integer *, doublereal *, doublereal *, |
| 33024 | integer *), dlauu2_(char *, integer *, |
| 33025 | doublereal *, integer *, integer *), xerbla_(char *, |
| 33026 | integer *); |
| 33027 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 33028 | integer *, integer *, ftnlen, ftnlen); |
| 33029 | |
| 33030 | |
| 33031 | /* |
| 33032 | -- LAPACK auxiliary routine (version 3.2) -- |
| 33033 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 33034 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 33035 | November 2006 |
| 33036 | |
| 33037 | |
| 33038 | Purpose |
| 33039 | ======= |
| 33040 | |
| 33041 | DLAUUM computes the product U * U' or L' * L, where the triangular |
| 33042 | factor U or L is stored in the upper or lower triangular part of |
| 33043 | the array A. |
| 33044 | |
| 33045 | If UPLO = 'U' or 'u' then the upper triangle of the result is stored, |
| 33046 | overwriting the factor U in A. |
| 33047 | If UPLO = 'L' or 'l' then the lower triangle of the result is stored, |
| 33048 | overwriting the factor L in A. |
| 33049 | |
| 33050 | This is the blocked form of the algorithm, calling Level 3 BLAS. |
| 33051 | |
| 33052 | Arguments |
| 33053 | ========= |
| 33054 | |
| 33055 | UPLO (input) CHARACTER*1 |
| 33056 | Specifies whether the triangular factor stored in the array A |
| 33057 | is upper or lower triangular: |
| 33058 | = 'U': Upper triangular |
| 33059 | = 'L': Lower triangular |
| 33060 | |
| 33061 | N (input) INTEGER |
| 33062 | The order of the triangular factor U or L. N >= 0. |
| 33063 |