Subroutine */
| 23330 | } /* zlauum_ */ |
| 23331 | |
| 23332 | /* Subroutine */ int zpotf2_(char *uplo, integer *n, doublecomplex *a, |
| 23333 | integer *lda, integer *info) |
| 23334 | { |
| 23335 | /* System generated locals */ |
| 23336 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 23337 | doublereal d__1; |
| 23338 | doublecomplex z__1, z__2; |
| 23339 | |
| 23340 | /* Local variables */ |
| 23341 | static integer j; |
| 23342 | static doublereal ajj; |
| 23343 | extern logical lsame_(char *, char *); |
| 23344 | extern /* Double Complex */ VOID zdotc_(doublecomplex *, integer *, |
| 23345 | doublecomplex *, integer *, doublecomplex *, integer *); |
| 23346 | extern /* Subroutine */ int zgemv_(char *, integer *, integer *, |
| 23347 | doublecomplex *, doublecomplex *, integer *, doublecomplex *, |
| 23348 | integer *, doublecomplex *, doublecomplex *, integer *); |
| 23349 | static logical upper; |
| 23350 | extern logical disnan_(doublereal *); |
| 23351 | extern /* Subroutine */ int xerbla_(char *, integer *), zdscal_( |
| 23352 | integer *, doublereal *, doublecomplex *, integer *), zlacgv_( |
| 23353 | integer *, doublecomplex *, integer *); |
| 23354 | |
| 23355 | |
| 23356 | /* |
| 23357 | -- LAPACK routine (version 3.2) -- |
| 23358 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 23359 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 23360 | November 2006 |
| 23361 | |
| 23362 | |
| 23363 | Purpose |
| 23364 | ======= |
| 23365 | |
| 23366 | ZPOTF2 computes the Cholesky factorization of a complex Hermitian |
| 23367 | positive definite matrix A. |
| 23368 | |
| 23369 | The factorization has the form |
| 23370 | A = U' * U , if UPLO = 'U', or |
| 23371 | A = L * L', if UPLO = 'L', |
| 23372 | where U is an upper triangular matrix and L is lower triangular. |
| 23373 | |
| 23374 | This is the unblocked version of the algorithm, calling Level 2 BLAS. |
| 23375 | |
| 23376 | Arguments |
| 23377 | ========= |
| 23378 | |
| 23379 | UPLO (input) CHARACTER*1 |
| 23380 | Specifies whether the upper or lower triangular part of the |
| 23381 | Hermitian matrix A is stored. |
| 23382 | = 'U': Upper triangular |
| 23383 | = 'L': Lower triangular |
| 23384 | |
| 23385 | N (input) INTEGER |
| 23386 | The order of the matrix A. N >= 0. |
| 23387 | |
| 23388 | A (input/output) COMPLEX*16 array, dimension (LDA,N) |
| 23389 | On entry, the Hermitian matrix A. If UPLO = 'U', the leading |