Subroutine */
| 32852 | } /* dlatrd_ */ |
| 32853 | |
| 32854 | /* Subroutine */ int dlauu2_(char *uplo, integer *n, doublereal *a, integer * |
| 32855 | lda, integer *info) |
| 32856 | { |
| 32857 | /* System generated locals */ |
| 32858 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 32859 | |
| 32860 | /* Local variables */ |
| 32861 | static integer i__; |
| 32862 | static doublereal aii; |
| 32863 | extern doublereal ddot_(integer *, doublereal *, integer *, doublereal *, |
| 32864 | integer *); |
| 32865 | extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, |
| 32866 | integer *); |
| 32867 | extern logical lsame_(char *, char *); |
| 32868 | extern /* Subroutine */ int dgemv_(char *, integer *, integer *, |
| 32869 | doublereal *, doublereal *, integer *, doublereal *, integer *, |
| 32870 | doublereal *, doublereal *, integer *); |
| 32871 | static logical upper; |
| 32872 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 32873 | |
| 32874 | |
| 32875 | /* |
| 32876 | -- LAPACK auxiliary routine (version 3.2) -- |
| 32877 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 32878 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 32879 | November 2006 |
| 32880 | |
| 32881 | |
| 32882 | Purpose |
| 32883 | ======= |
| 32884 | |
| 32885 | DLAUU2 computes the product U * U' or L' * L, where the triangular |
| 32886 | factor U or L is stored in the upper or lower triangular part of |
| 32887 | the array A. |
| 32888 | |
| 32889 | If UPLO = 'U' or 'u' then the upper triangle of the result is stored, |
| 32890 | overwriting the factor U in A. |
| 32891 | If UPLO = 'L' or 'l' then the lower triangle of the result is stored, |
| 32892 | overwriting the factor L in A. |
| 32893 | |
| 32894 | This is the unblocked form of the algorithm, calling Level 2 BLAS. |
| 32895 | |
| 32896 | Arguments |
| 32897 | ========= |
| 32898 | |
| 32899 | UPLO (input) CHARACTER*1 |
| 32900 | Specifies whether the triangular factor stored in the array A |
| 32901 | is upper or lower triangular: |
| 32902 | = 'U': Upper triangular |
| 32903 | = 'L': Lower triangular |
| 32904 | |
| 32905 | N (input) INTEGER |
| 32906 | The order of the triangular factor U or L. N >= 0. |
| 32907 | |
| 32908 | A (input/output) DOUBLE PRECISION array, dimension (LDA,N) |
| 32909 | On entry, the triangular factor U or L. |
| 32910 | On exit, if UPLO = 'U', the upper triangle of A is |
| 32911 | overwritten with the upper triangle of the product U * U'; |