TESTING/LIN/zget07.f(3)
πΈ NAME πΈ
TESTING/LIN/zget07.f β okay bestie, this is the page for that file, and it lives in the Library Functions Manual. Itβs giving βlibrary functions manualβ energy from the very first line, periodt. π
β¨ SYNOPSIS β¨
π§ Functions/Subroutines
subroutine zget07 (trans, n, nrhs, a, lda, b, ldb, x, ldx, xact, ldxact,
ferr, chkferr, berr, reslts)
ZGET07 is the routine, and yes, thatβs its whole name. π
π Function/Subroutine Documentation π
subroutine zget07 (character trans, integer n, integer nrhs, complex*16,
dimension( lda, * ) a, integer lda, complex*16, dimension( ldb, * ) b,
integer ldb, complex*16, dimension( ldx, * ) x, integer ldx, complex*16,
dimension( ldxact, * ) xact, integer ldxact, double precision,
dimension( * ) ferr, logical chkferr, double precision, dimension( * )
berr, double precision, dimension( * ) reslts)
ZGET07 β thatβs the main character of this whole page, no cap. π
π Purpose π
Okay so, ZGET07 tests the error bounds from iterative refinement for the computed solution to a system of equations op(A)*X = B, where A is a general n by n matrix and op(A) = A or A**T, depending on TRANS.
The purpose block is literally this:
RESLTS(1) = test of the error bound
= norm(X - XACT) / ( norm(X) * FERR )
A large value is returned if this ratio is not less than one.
RESLTS(2) = residual from the iterative refinement routine
= the maximum of BERR / ( (n+1)*EPS + (*) ), where
(*) = (n+1)*UNFL / (min_i (abs(op(A))*abs(X) +abs(b))_i )
So yeah, itβs checking whether the error-bound ratios are slaying or lowkey sus. π
π¦ Parameters π¦
TRANSβTRANS is CHARACTER*1. It specifies the form of the system of equations:= 'N': A * X = B (No transpose),= 'T': A**T * X = B (Transpose),= 'C': A**H * X = B (Conjugate transpose = Transpose).NβN is INTEGER. This is the number of rows of the matricesXandXACT.N >= 0.NRHSβNRHS is INTEGER. This is the number of columns of the matricesXandXACT.NRHS >= 0.AβA is COMPLEX*16 array, dimension (LDA,N). Itβs the originalnbynmatrixA.LDAβLDA is INTEGER. This is the leading dimension of the arrayA.LDA >= max(1,N).BβB is COMPLEX*16 array, dimension (LDB,NRHS). These are the right hand side vectors for the system of linear equations.LDBβLDB is INTEGER. This is the leading dimension of the arrayB.LDB >= max(1,N).XβX is COMPLEX*16 array, dimension (LDX,NRHS). These are the computed solution vectors. Each vector is stored as a column of the matrixX.LDXβLDX is INTEGER. This is the leading dimension of the arrayX.LDX >= max(1,N).XACTβXACT is COMPLEX*16 array, dimension (LDX,NRHS). These are the exact solution vectors. Each vector is stored as a column of the matrixXACT.LDXACTβLDXACT is INTEGER. This is the leading dimension of the arrayXACT.LDXACT >= max(1,N).FERRβFERR is DOUBLE PRECISION array, dimension (NRHS). These are the estimated forward error bounds for each solution vectorX. IfXTRUEis the true solution,FERRbounds the magnitude of the largest entry in(X - XTRUE)divided by the magnitude of the largest entry inX.CHKFERRβCHKFERR is LOGICAL. Set it to.TRUE.to checkFERR,.FALSE.not to checkFERR. When the test system is ill-conditioned, the solution inXACTmay be incorrect. This is important, bestie β ill-conditioned systems can be messy. πBERRβBERR is DOUBLE PRECISION array, dimension (NRHS). This is the componentwise relative backward error of each solution vector (i.e., the smallest relative change in any entry ofAorBthat makesXan exact solution).RESLTSβRESLTS is DOUBLE PRECISION array, dimension (2). This is the maximum over theNRHSsolution vectors of the ratios:
RESLTS(1) = norm(X - XACT) / ( norm(X) * FERR )
RESLTS(2) = BERR / ( (n+1)*EPS + (*) )
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 164 of file zget07.f. π
π Author (the auto-generated one) π
Generated automatically by Doxygen for LAPACK from the source code.
source/version sign-off, bestie: this iconic little page is LAPACK, Version 3.12.0, and the page is TESTING/LIN/zget07.f(3). Thatβs the footer vibe, and it ate. πβ¨
π¦ 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