Subroutine */
| 36801 | } /* dormtr_ */ |
| 36802 | |
| 36803 | /* Subroutine */ int dpotf2_(char *uplo, integer *n, doublereal *a, integer * |
| 36804 | lda, integer *info) |
| 36805 | { |
| 36806 | /* System generated locals */ |
| 36807 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 36808 | doublereal d__1; |
| 36809 | |
| 36810 | /* Local variables */ |
| 36811 | static integer j; |
| 36812 | static doublereal ajj; |
| 36813 | extern doublereal ddot_(integer *, doublereal *, integer *, doublereal *, |
| 36814 | integer *); |
| 36815 | extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, |
| 36816 | integer *); |
| 36817 | extern logical lsame_(char *, char *); |
| 36818 | extern /* Subroutine */ int dgemv_(char *, integer *, integer *, |
| 36819 | doublereal *, doublereal *, integer *, doublereal *, integer *, |
| 36820 | doublereal *, doublereal *, integer *); |
| 36821 | static logical upper; |
| 36822 | extern logical disnan_(doublereal *); |
| 36823 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 36824 | |
| 36825 | |
| 36826 | /* |
| 36827 | -- LAPACK routine (version 3.2) -- |
| 36828 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 36829 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 36830 | November 2006 |
| 36831 | |
| 36832 | |
| 36833 | Purpose |
| 36834 | ======= |
| 36835 | |
| 36836 | DPOTF2 computes the Cholesky factorization of a real symmetric |
| 36837 | positive definite matrix A. |
| 36838 | |
| 36839 | The factorization has the form |
| 36840 | A = U' * U , if UPLO = 'U', or |
| 36841 | A = L * L', if UPLO = 'L', |
| 36842 | where U is an upper triangular matrix and L is lower triangular. |
| 36843 | |
| 36844 | This is the unblocked version of the algorithm, calling Level 2 BLAS. |
| 36845 | |
| 36846 | Arguments |
| 36847 | ========= |
| 36848 | |
| 36849 | UPLO (input) CHARACTER*1 |
| 36850 | Specifies whether the upper or lower triangular part of the |
| 36851 | symmetric matrix A is stored. |
| 36852 | = 'U': Upper triangular |
| 36853 | = 'L': Lower triangular |
| 36854 | |
| 36855 | N (input) INTEGER |
| 36856 | The order of the matrix A. N >= 0. |
| 36857 | |
| 36858 | A (input/output) DOUBLE PRECISION array, dimension (LDA,N) |
| 36859 | On entry, the symmetric matrix A. If UPLO = 'U', the leading |
| 36860 | n by n upper triangular part of A contains the upper |