ZPSTRF(3)

πŸ“¦ version 3.12.0, in opensuse-leap 16.0

πŸ’… ZPSTRF(3) β€” Library Functions Manual πŸ’…

okay bestie, this is ZPSTRF(3) from the Library Functions Manual. The header path is /home/abui...RC/zpstrf.f(3) β€” and yes, it shows up twice in the raw header, iconic behavior πŸ’…βœ¨


πŸŽ€ NAME πŸŽ€

/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/zpstrf.f

That whole path is the name section. No cap, that's the source file.


πŸ¦‹ SYNOPSIS πŸ¦‹

🧚 Functions/Subroutines 🧚

subroutine ZPSTRF (uplo, n, a, lda, piv, rank, tol, work, info)

ZPSTRF computes the Cholesky factorization with complete pivoting of a complex Hermitian positive semidefinite matrix.


🌸 Function/Subroutine Documentation 🌸

subroutine ZPSTRF (character uplo, integer n, complex*16, dimension( lda, * ) a, integer lda, integer, dimension( n ) piv, integer rank, double precision tol, double precision, dimension( 2*n ) work, integer info)

ZPSTRF computes the Cholesky factorization with complete pivoting of a complex Hermitian positive semidefinite matrix.

✨ Purpose ✨

ZPSTRF computes the Cholesky factorization with complete pivoting of a complex Hermitian positive semidefinite matrix A.

The factorization has the form

P**T * A * P = U**H * U ,  if UPLO = 'U',
P**T * A * P = L  * L**H,  if UPLO = 'L',

where U is an upper triangular matrix and L is lower triangular, and P is stored as vector PIV.

Heads up, bestie: this algorithm does not attempt to check that A is positive semidefinite. This version of the algorithm calls level 3 BLAS.

🌈 Parameters 🌈

  • UPLO β€” UPLO is CHARACTER*1. Specifies whether the upper or lower triangular part of the symmetric matrix A is stored.
    • = 'U': Upper triangular
    • = 'L': Lower triangular
  • N β€” N is INTEGER. The order of the matrix A. N >= 0.
  • A β€” A is COMPLEX*16 array, dimension (LDA,N). On entry, the symmetric matrix A. If UPLO = 'U', the leading n by n upper triangular part of A contains the upper triangular part of the matrix A, and the strictly lower triangular part of A is not referenced. If UPLO = 'L', the leading n by n lower triangular part of A contains the lower triangular part of the matrix A, and the strictly upper triangular part of A is not referenced. On exit, if INFO = 0, the factor U or L from the Cholesky factorization as above.
  • LDA β€” LDA is INTEGER. The leading dimension of the array A. LDA >= max(1,N).
  • PIV β€” PIV is INTEGER array, dimension (N). PIV is such that the nonzero entries are P( PIV(K), K ) = 1.
  • RANK β€” RANK is INTEGER. The rank of A given by the number of steps the algorithm completed.
  • TOL β€” TOL is DOUBLE PRECISION. User defined tolerance. If TOL < 0, then N*U*MAX( A(K,K) ) will be used. The algorithm terminates at the (K-1)st step if the pivot <= TOL.
  • WORK β€” WORK is DOUBLE PRECISION array, dimension (2*N). Work space.
  • INFO β€” INFO is INTEGER.
    • < 0: If INFO = -K, the K-th argument had an illegal value,
    • = 0: algorithm completed successfully, and
    • > 0: the matrix A is either rank deficient with computed rank as returned in RANK, or is not positive semidefinite. See Section 7 of LAPACK Working Note #161 for further information.

🫢 Author 🫢

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Definition at line 141 of file zpstrf.f.


πŸ‘‘ Author πŸ‘‘

Generated automatically by Doxygen for LAPACK from the source code.


mwah πŸ’‹ β€” source/version is LAPACK Version 3.12.0, and the footer path is /home/abui...RC/zpstrf.f(3). That's the whole vibe, periodt 🫢✨

πŸ“¦ where this page comes from