DPPRFS(3)
okay bestie, welcome to DPPRFS(3) from the LAPACK pages in the Library Functions Manual β this is the whole page, and it's giving "precision linear algebra but make it cute" energy. β¨
π NAME π
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/dpprfs.f
β¨ SYNOPSIS β¨
π§ Functions/Subroutines π§
okay so the function lineup is short and iconic:
subroutine DPPRFS (uplo, n, nrhs, ap, afp, b, ldb, x, ldx, ferr, berr,
work, iwork, info)
DPPRFS
πΈ Function/Subroutine Documentation πΈ
here's the full signature, copied exactly because we do not freestyle code in this group chat:
subroutine DPPRFS (character uplo, integer n, integer nrhs, double
precision, dimension( * ) ap, double precision, dimension( * ) afp,
double precision, dimension( ldb, * ) b, integer ldb, double precision,
dimension( ldx, * ) x, integer ldx, double precision, dimension( * )
ferr, double precision, dimension( * ) berr, double precision,
dimension( * ) work, integer, dimension( * ) iwork, integer info)
DPPRFS
π Purpose π
okay so DPPRFS improves the computed solution to a system of linear equations when the coefficient matrix is symmetric positive definite and packed, and provides error bounds and backward error estimates for the solution. that's the assignment, and it ate. π
π¦ Parameters π¦
UPLOβUPLO is CHARACTER*1;= 'U': upper triangle ofAis stored;= 'L': lower triangle ofAis stored. pick your triangle, bestie. πNβN is INTEGER; the order of the matrixA.N >= 0.NRHSβNRHS is INTEGER; the number of right hand sides, i.e., the number of columns of the matricesBandX.NRHS >= 0.APβAP is DOUBLE PRECISION array, dimension (N*(N+1)/2); the upper or lower triangle of the symmetric matrixA, packed columnwise in a linear array. Thej-th column ofAis stored in the arrayAPas follows: ifUPLO = 'U',AP(i + (j-1)*j/2) = A(i,j)for1<=i<=j; ifUPLO = 'L',AP(i + (j-1)*(2n-j)/2) = A(i,j)forj<=i<=n. packing, no cap. πAFPβAFP is DOUBLE PRECISION array, dimension (N*(N+1)/2); the triangular factorUorLfrom the Cholesky factorizationA = U**T*UorA = L*L**T, as computed byDPPTRF/ZPPTRF, packed columnwise in a linear array in the same format asA(seeAP).BβB is DOUBLE PRECISION array, dimension (LDB,NRHS); the right hand side matrixB.LDBβLDB is INTEGER; the leading dimension of the arrayB.LDB >= max(1,N).XβX is DOUBLE PRECISION array, dimension (LDX,NRHS); on entry, the solution matrixX, as computed byDPPTRS. On exit, the improved solution matrixX.LDXβLDX is INTEGER; the leading dimension of the arrayX.LDX >= max(1,N).FERRβFERR is DOUBLE PRECISION array, dimension (NRHS); the estimated forward error bound for each solution vectorX(j)(thej-th column of the solution matrixX). IfXTRUEis the true solution corresponding toX(j),FERR(j)is an estimated upper bound for the magnitude of the largest element in(X(j) - XTRUE)divided by the magnitude of the largest element inX(j). The estimate is as reliable as the estimate forRCOND, and is almost always a slight overestimate of the true error. she's honest about being a tiny overestimate. πBERRβBERR is DOUBLE PRECISION array, dimension (NRHS); the componentwise relative backward error of each solution vectorX(j)(i.e., the smallest relative change in any element ofAorBthat makesX(j)an exact solution).WORKβWORK is DOUBLE PRECISION array, dimension (3*N).IWORKβIWORK is INTEGER array, dimension (N).INFOβINFO is INTEGER;= 0: successful exit;< 0: ifINFO = -i, thei-th argument had an illegal value. illegal-value vibe check, bestie. π
π Internal Parameters π
ITMAXis the maximum number of steps of iterative refinement. that's the internal tea, and it's lowkey important. β
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 169 of file dpprfs.f.
π¦ Author π¦
Generated automatically by Doxygen for LAPACK from the source code. not me writing this by hand, Doxygen did that. π
bye bestie! this page is from LAPACK, Version 3.12.0, with footer path /home/abui...RC/dpprfs.f(3) β that's the source/version info, and there's no separate date line in the footer, periodt. π«Άβ¨
π¦ where this page comes from
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms