Subroutine */
| 23218 | } /* clauum_ */ |
| 23219 | |
| 23220 | /* Subroutine */ int cpotf2_(char *uplo, integer *n, complex *a, integer *lda, |
| 23221 | integer *info) |
| 23222 | { |
| 23223 | /* System generated locals */ |
| 23224 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 23225 | real r__1; |
| 23226 | complex q__1, q__2; |
| 23227 | |
| 23228 | /* Local variables */ |
| 23229 | static integer j; |
| 23230 | static real ajj; |
| 23231 | extern /* Complex */ VOID cdotc_(complex *, integer *, complex *, integer |
| 23232 | *, complex *, integer *); |
| 23233 | extern logical lsame_(char *, char *); |
| 23234 | extern /* Subroutine */ int cgemv_(char *, integer *, integer *, complex * |
| 23235 | , complex *, integer *, complex *, integer *, complex *, complex * |
| 23236 | , integer *); |
| 23237 | static logical upper; |
| 23238 | extern /* Subroutine */ int clacgv_(integer *, complex *, integer *), |
| 23239 | csscal_(integer *, real *, complex *, integer *), xerbla_(char *, |
| 23240 | integer *); |
| 23241 | extern logical sisnan_(real *); |
| 23242 | |
| 23243 | |
| 23244 | /* |
| 23245 | -- LAPACK routine (version 3.2) -- |
| 23246 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 23247 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 23248 | November 2006 |
| 23249 | |
| 23250 | |
| 23251 | Purpose |
| 23252 | ======= |
| 23253 | |
| 23254 | CPOTF2 computes the Cholesky factorization of a complex Hermitian |
| 23255 | positive definite matrix A. |
| 23256 | |
| 23257 | The factorization has the form |
| 23258 | A = U' * U , if UPLO = 'U', or |
| 23259 | A = L * L', if UPLO = 'L', |
| 23260 | where U is an upper triangular matrix and L is lower triangular. |
| 23261 | |
| 23262 | This is the unblocked version of the algorithm, calling Level 2 BLAS. |
| 23263 | |
| 23264 | Arguments |
| 23265 | ========= |
| 23266 | |
| 23267 | UPLO (input) CHARACTER*1 |
| 23268 | Specifies whether the upper or lower triangular part of the |
| 23269 | Hermitian matrix A is stored. |
| 23270 | = 'U': Upper triangular |
| 23271 | = 'L': Lower triangular |
| 23272 | |
| 23273 | N (input) INTEGER |
| 23274 | The order of the matrix A. N >= 0. |
| 23275 | |
| 23276 | A (input/output) COMPLEX array, dimension (LDA,N) |
| 23277 | On entry, the Hermitian matrix A. If UPLO = 'U', the leading |