TESTING/LIN/dpot06.f(3)
β¨ TESTING/LIN/dpot06.f β¨
okay so bestie, welcome to the group chat of LAPACK test routines π
this is TESTING/LIN/dpot06.f, straight out of the Library Functions Manual, and it's the girl whose whole job is checking that your linear solve actually slayed. no cap.
πΈ NAME πΈ
TESTING/LIN/dpot06.f β yes that's literally the name, she's a file and she's iconic for it π
π¦ SYNOPSIS π¦
Functions/Subroutines
subroutine dpot06 (uplo, n, nrhs, a, lda, x, ldx, b, ldb, rwork, resid)
DPOT06 β that's her. that's the whole vibe. π
π Function/Subroutine Documentation π
subroutine dpot06 (character uplo, integer n, integer nrhs, double
precision, dimension( lda, * ) a, integer lda, double precision,
dimension( ldx, * ) x, integer ldx, double precision, dimension( ldb, *
) b, integer ldb, double precision, dimension( * ) rwork, double
precision resid)
DPOT06
π Purpose:
!>!> DPOT06 computes the residual for a solution of a system of linear!> equations A*x = b :!> RESID = norm(B - A*X,inf) / ( norm(A,inf) * norm(X,inf) * EPS ),!> where EPS is the machine epsilon.!>
so basically she takes your answer, plugs it back in, and measures how far off you are relative to how big everything is β using the infinity norm and machine epsilon. that's the whole plot. π
π Parameters
UPLOβUPLO is CHARACTER*1, and she's the one deciding which half of the symmetric matrixAyou bothered to store:= 'U': Upper triangular= 'L': Lower triangular
NβN is INTEGER, the number of rows and columns of the matrixA.N >= 0. (square matrix behavior, we love to see it β¨)NRHSβNRHS is INTEGER, the number of columns ofB, aka the matrix of right hand sides.NRHS >= 0. that's your lil stack of systems, bestie π«ΆAβA is DOUBLE PRECISION array, dimension (LDA,N). described in the source as "The original M x N matrix A." β yes, it saysM x N, and no,Mis not a parameter, we don't ask questions, we just report the lore πLDAβLDA is INTEGER, the leading dimension of the arrayA.LDA >= max(1,N). don't skimp on it or it's giving undefined behavior πXβX is DOUBLE PRECISION array, dimension (LDX,NRHS), holding the computed solution vectors for the system of linear equations. this is your answer key, the main character π₯LDXβLDX is INTEGER, the leading dimension of the arrayX. and the source says: "IfTRANS = 'N',LDX >= max(1,N); ifTRANS = 'T'or'C',LDX >= max(1,N)." (yep,TRANSβ she's mentioned but she's not invited to the parameter list, ik ik, we're just reading the flyer as printed π)BβB is DOUBLE PRECISION array, dimension (LDB,NRHS). On entry, the right hand side vectors for the system of linear equations. On exit,Bis overwritten with the differenceB - A*X. β οΈ so if you were attached to yourB, too late, she's been transformed πLDBβLDB is INTEGER, the leading dimension of the arrayB. "IFTRANS = 'N',LDB >= max(1,M); ifTRANS = 'T'or'C',LDB >= max(1,N)." (again withMandTRANS, and again we simply report π)RWORKβRWORK is DOUBLE PRECISION array, dimension (N). that's it. that's the whole description. she's a workspace and she's lowkey mysterious π¦RESIDβRESID is DOUBLE PRECISION, the maximum over the number of right hand sides ofnorm(B - A*X) / ( norm(A) * norm(X) * EPS ). the final score, the vibe check, the moment of truth β
π Author
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 125 of file dpot06.f.
π« Author π«
Generated automatically by Doxygen for LAPACK from the source code. (not a human, just a robot doing her little robot thing π€β¨)
and that's the page!! LAPACK, Version 3.12.0, TESTING/LIN/dpot06.f(3) β mwah, go forth and check those residuals bestie π
ππ§
π¦ 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