Subroutine */
| 36666 | } /* sormtr_ */ |
| 36667 | |
| 36668 | /* Subroutine */ int spotf2_(char *uplo, integer *n, real *a, integer *lda, |
| 36669 | integer *info) |
| 36670 | { |
| 36671 | /* System generated locals */ |
| 36672 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 36673 | real r__1; |
| 36674 | |
| 36675 | /* Local variables */ |
| 36676 | static integer j; |
| 36677 | static real ajj; |
| 36678 | extern doublereal sdot_(integer *, real *, integer *, real *, integer *); |
| 36679 | extern logical lsame_(char *, char *); |
| 36680 | extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *), |
| 36681 | sgemv_(char *, integer *, integer *, real *, real *, integer *, |
| 36682 | real *, integer *, real *, real *, integer *); |
| 36683 | static logical upper; |
| 36684 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 36685 | extern logical sisnan_(real *); |
| 36686 | |
| 36687 | |
| 36688 | /* |
| 36689 | -- LAPACK routine (version 3.2) -- |
| 36690 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 36691 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 36692 | November 2006 |
| 36693 | |
| 36694 | |
| 36695 | Purpose |
| 36696 | ======= |
| 36697 | |
| 36698 | SPOTF2 computes the Cholesky factorization of a real symmetric |
| 36699 | positive definite matrix A. |
| 36700 | |
| 36701 | The factorization has the form |
| 36702 | A = U' * U , if UPLO = 'U', or |
| 36703 | A = L * L', if UPLO = 'L', |
| 36704 | where U is an upper triangular matrix and L is lower triangular. |
| 36705 | |
| 36706 | This is the unblocked version of the algorithm, calling Level 2 BLAS. |
| 36707 | |
| 36708 | Arguments |
| 36709 | ========= |
| 36710 | |
| 36711 | UPLO (input) CHARACTER*1 |
| 36712 | Specifies whether the upper or lower triangular part of the |
| 36713 | symmetric matrix A is stored. |
| 36714 | = 'U': Upper triangular |
| 36715 | = 'L': Lower triangular |
| 36716 | |
| 36717 | N (input) INTEGER |
| 36718 | The order of the matrix A. N >= 0. |
| 36719 | |
| 36720 | A (input/output) REAL array, dimension (LDA,N) |
| 36721 | On entry, the symmetric matrix A. If UPLO = 'U', the leading |
| 36722 | n by n upper triangular part of A contains the upper |
| 36723 | triangular part of the matrix A, and the strictly lower |
| 36724 | triangular part of A is not referenced. If UPLO = 'L', the |
| 36725 | leading n by n lower triangular part of A contains the lower |