SRC/sla_porpvgrw.f(3)
π SRC/sla_porpvgrw.f(3) π
okay so bestie, this page comes straight out of the Library Functions Manual and it's all about one very specific little Fortran function. no cap, she's niche but she's iconic. let's get into it π«Ά
β¨ SYNOPSIS β¨
Functions/Subroutines
real function sla_porpvgrw (uplo, ncols, a, lda, af, ldaf, work)β SLA_PORPVGRW computes the reciprocal pivot growth factor norm(A)/norm(U) for a symmetric or Hermitian positive-definite matrix. that's the whole vibe π
πΈ Function/Subroutine Documentation πΈ
real function sla_porpvgrw (character*1 uplo, integer ncols, real, dimension( lda, * ) a, integer lda, real, dimension( ldaf, * ) af, integer ldaf, real, dimension( * ) work)
k so what does she do? SLA_PORPVGRW computes the reciprocal pivot growth factor norm(A)/norm(U) for a symmetric or Hermitian positive-definite matrix. that's it, that's the function β¨
π Purpose π
!> !> SLA_PORPVGRW computes the reciprocal pivot growth factor !> norm(A)/norm(U). The norm is used. If this is !> much less than 1, the stability of the LU factorization of the !> (equilibrated) matrix A could be poor. This also means that the !> solution X, estimated condition numbers, and error bounds could be > unreliable. !>
(and yes, the original really does just say The norm is used. with a lil gap in there β not me, that's literally what it says, i'm just the messenger π)
okay so translating that from math-speak: if this number comes back much less than 1, the stability of the LU factorization of the (equilibrated) matrix A could be poor, and that means the solution X, the estimated condition numbers, and the error bounds could all be unreliable. so likeβ¦ don't trust the vibes if this number is tiny π
π¦ Parameters π¦
UPLOβUPLOisCHARACTER*1= 'U': Upper triangle of A is stored;= 'L': Lower triangle of A is stored.
NCOLSβNCOLSisINTEGER- The number of columns of the matrix A.
NCOLS >= 0.
- The number of columns of the matrix A.
AβAisREALarray, dimension(LDA,N)- On entry, the N-by-N matrix A.
LDAβLDAisINTEGER- The leading dimension of the array A.
LDA >= max(1,N).
- The leading dimension of the array A.
AFβAFisREALarray, dimension(LDAF,N)- The triangular factor U or L from the Cholesky factorization
A = U**T*UorA = L*L**T, as computed bySPOTRF.
- The triangular factor U or L from the Cholesky factorization
LDAFβLDAFisINTEGER- The leading dimension of the array AF.
LDAF >= max(1,N).
- The leading dimension of the array AF.
WORKβWORKisREALarray, dimension(2*N)- literally just workspace, bestie. give her
(2*N)and she's happy π
- literally just workspace, bestie. give her
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
β Definition at line 103 of file sla_porpvgrw.f. β
π Author π
Generated automatically by Doxygen for LAPACK from the source code. (so nobody hand-wrote this one, the robots did it π€β¨)
and that's the whole page, bestie! this lil guy is SRC/sla_porpvgrw.f(3) from LAPACK Version 3.12.0 β slay, go compute some pivot growth factors responsibly πππ¦
π¦ 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