naginterfaces.library.lapacklin.zpftrf

naginterfaces.library.lapacklin.zpftrf(transr, uplo, n, ar)

zpftrf computes the Cholesky factorization of a complex Hermitian positive definite matrix stored in Rectangular Full Packed (RFP) format.

For full information please refer to the NAG Library document for f07wr

https://support.nag.com/numeric/nl/nagdoc_30.2/flhtml/f07/f07wrf.html

Parameters

transrstr, length 1

Specifies whether the normal RFP representation of $A$ or its conjugate transpose is stored.

$transr=\text{'N'}$

The matrix $A$ is stored in normal RFP format.

$transr=\text{'C'}$

The conjugate transpose of the RFP representation of the matrix $A$ is stored.

uplostr, length 1

Specifies whether the upper or lower triangular part of $A$ is stored.

$uplo=\text{'U'}$

The upper triangular part of $A$ is stored, and $A$ is factorized as $U^{H}U$, where $U$ is upper triangular.

$uplo=\text{'L'}$

The lower triangular part of $A$ is stored, and $A$ is factorized as $L{L}^H$, where $L$ is lower triangular.

nint

$n$, the order of the matrix $A$.

arcomplex, array-like, shape $(n\times (n+1)/2)$

The upper or lower triangular part (as specified by $uplo$) of the $n\times n$ Hermitian matrix $A$, in either normal or transposed RFP format (as specified by $transr$). The storage format is described in detail in the F07 Introduction.

Returns

arcomplex, ndarray, shape $(n\times (n+1)/2)$

If no exception or warning is raised, the factor $U$ or $L$ from the Cholesky factorization $A=U^{H}U$ or $A=L{L}^H$, in the same storage format as $A$.

Raises

NagValueError

(errno $-1$)

On entry, error in parameter $transr$.

Constraint: $transr=\text{'N'}$ or $\text{'C'}$.

(errno $-2$)

On entry, error in parameter $uplo$.

Constraint: $uplo=\text{'U'}$ or $\text{'L'}$.

(errno $-3$)

On entry, error in parameter $n$.

Constraint: $n\ge 0$.

(errno $i>0$)

The leading minor of order $⟨value⟩$ is not positive definite and the factorization could not be completed. Hence $A$ itself is not positive definite. This may indicate an error in forming the matrix $A$. There is no function specifically designed to factorize a Hermitian matrix stored in RFP format which is not positive definite; the matrix must be treated as a full Hermitian matrix, by calling zhetrf().

Notes

zpftrf forms the Cholesky factorization of a complex Hermitian positive definite matrix $A$ either as $A=U^{H}U$ if $uplo=\text{'U'}$ or $A=L{L}^H$ if $uplo=\text{'L'}$, where $U$ is an upper triangular matrix and $L$ is a lower triangular, stored in RFP format. The RFP storage format is described in the F07 Introduction.

References

Demmel, J W, 1989, On floating-point errors in Cholesky, LAPACK Working Note No. 14, University of Tennessee, Knoxville, https://www.netlib.org/lapack/lawnspdf/lawn14.pdf

Gustavson, F G, Waśniewski, J, Dongarra, J J and Langou, J, 2010, Rectangular full packed format for Cholesky’s algorithm: factorization, solution, and inversion, ACM Trans. Math. Software (37, 2)