pptri(3)
β¨ pptri(3) β¨ β the triangular inverse girlies
okay so bestie, this is the manual page for pptri, living in section 3 of the Library Functions Manual π. it's the LAPACK crew that inverts your positive definite matrices after you've already done the Cholesky glow-up. four flavors, one vibe. let's go π
πΈ NAME πΈ
pptri β pptri: triangular inverse
π SYNOPSIS π
π Functions
subroutine cpptri (uplo, n, ap, info)β CPPTRIsubroutine dpptri (uplo, n, ap, info)β DPPTRIsubroutine spptri (uplo, n, ap, info)β SPPTRIsubroutine zpptri (uplo, n, ap, info)β ZPPTRI
π¦ Detailed Description π¦
(yes babe, the heading is literally "Detailed Description" and there's nothing under it in the source β it's giving mysterious main character energy. the real tea is in the function docs below π)
π Function Documentation π
π§ subroutine cpptri (character uplo, integer n, complex, dimension( * ) ap, integer info) π§
CPPTRI
β Purpose
CPPTRI computes the inverse of a complex Hermitian positive definite matrix A using the Cholesky factorization A = U**H*U or A = L*L**H computed by CPPTRF. literally the inverse, handed to you on a platter β¨
π₯ Parameters
UPLOβUPLOisCHARACTER*1= 'U': Upper triangular factor is stored inAP;= 'L': Lower triangular factor is stored inAP.
NβNisINTEGER. The order of the matrixA.N >= 0.APβAPisCOMPLEXarray, dimension(N*(N+1)/2).
On entry, the triangular factor U or L from the Cholesky factorization A = U**H*U or A = L*L**H, packed columnwise as a linear array. The j-th column of U or L is stored in the array AP as follows:
if UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for 1<=i<=j;
if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for j<=i<=n.
On exit, the upper or lower triangle of the (Hermitian) inverse of A, overwriting the input factor U or L. she ate and left no crumbs π
INFOβINFOisINTEGER
= 0: successful exit
< 0: if INFO = -i, the i-th argument had an illegal value
> 0: if INFO = i, the (i,i) element of the factor U or L is
zero, and the inverse could not be computed.
π Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
π Definition at line 92 of file cpptri.f.
π subroutine dpptri (character uplo, integer n, double precision, dimension( * ) ap, integer info) π
DPPTRI
β Purpose
DPPTRI computes the inverse of a real symmetric positive definite matrix A using the Cholesky factorization A = U**T*U or A = L*L**T computed by DPPTRF.
π₯ Parameters
UPLOβUPLOisCHARACTER*1= 'U': Upper triangular factor is stored inAP;= 'L': Lower triangular factor is stored inAP.
NβNisINTEGER. The order of the matrixA.N >= 0.APβAPisDOUBLE PRECISIONarray, dimension(N*(N+1)/2).
On entry, the triangular factor U or L from the Cholesky factorization A = U**T*U or A = L*L**T, packed columnwise as a linear array. The j-th column of U or L is stored in the array AP as follows:
if UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for 1<=i<=j;
if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for j<=i<=n.
On exit, the upper or lower triangle of the (symmetric) inverse of A, overwriting the input factor U or L.
INFOβINFOisINTEGER
= 0: successful exit
< 0: if INFO = -i, the i-th argument had an illegal value
> 0: if INFO = i, the (i,i) element of the factor U or L is
zero, and the inverse could not be computed.
π Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
π Definition at line 92 of file dpptri.f.
π subroutine spptri (character uplo, integer n, real, dimension( * ) ap, integer info) π
SPPTRI
β Purpose
SPPTRI computes the inverse of a real symmetric positive definite matrix A using the Cholesky factorization A = U**T*U or A = L*L**T computed by SPPTRF.
π₯ Parameters
UPLOβUPLOisCHARACTER*1= 'U': Upper triangular factor is stored inAP;= 'L': Lower triangular factor is stored inAP.
NβNisINTEGER. The order of the matrixA.N >= 0.APβAPisREALarray, dimension(N*(N+1)/2).
On entry, the triangular factor U or L from the Cholesky factorization A = U**T*U or A = L*L**T, packed columnwise as a linear array. The j-th column of U or L is stored in the array AP as follows:
if UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for 1<=i<=j;
if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for j<=i<=n.
On exit, the upper or lower triangle of the (symmetric) inverse of A, overwriting the input factor U or L.
INFOβINFOisINTEGER
= 0: successful exit
< 0: if INFO = -i, the i-th argument had an illegal value
> 0: if INFO = i, the (i,i) element of the factor U or L is
zero, and the inverse could not be computed.
π Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
π Definition at line 92 of file spptri.f.
π subroutine zpptri (character uplo, integer n, complex*16, dimension( * ) ap, integer info) π
ZPPTRI
β Purpose
ZPPTRI computes the inverse of a complex Hermitian positive definite matrix A using the Cholesky factorization A = U**H*U or A = L*L**H computed by ZPPTRF.
π₯ Parameters
UPLOβUPLOisCHARACTER*1= 'U': Upper triangular factor is stored inAP;= 'L': Lower triangular factor is stored inAP.
NβNisINTEGER. The order of the matrixA.N >= 0.APβAPisCOMPLEX*16array, dimension(N*(N+1)/2).
On entry, the triangular factor U or L from the Cholesky factorization A = U**H*U or A = L*L**H, packed columnwise as a linear array. The j-th column of U or L is stored in the array AP as follows:
if UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for 1<=i<=j;
if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for j<=i<=n.
On exit, the upper or lower triangle of the (Hermitian) inverse of A, overwriting the input factor U or L. no cap, the input gets fully replaced π«Ά
INFOβINFOisINTEGER
= 0: successful exit
< 0: if INFO = -i, the i-th argument had an illegal value
> 0: if INFO = i, the (i,i) element of the factor U or L is
zero, and the inverse could not be computed.
π Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
π Definition at line 92 of file zpptri.f.
π real talk: ifINFO > 0, that(i,i)element of your factorUorLis literally zero and the inverse could not be computed. that's not a vibe, that's a singular situation. and ifINFO < 0, you passed a bad argument βINFO = -itells you exactly which one. check yourself bestie π
π«Ά Author π«Ά
Generated automatically by Doxygen for LAPACK from the source code.
and that's the whole page, bestie πΈ the source is LAPACK, Version 3.12.0, and this cutie is pptri(3) β go invert something iconic πβ¨
π¦ where this page comes from
- arch rolling: lapack-doc 3.12.1-2 Β· custom
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms