ZGTT05(3)
π ZGTT05(3) β the LAPACK test routine bestie π
okay so this page is from the Library Functions Manual, and we're talking about ZGTT05. It's a subroutine from LAPACK, and it tests the error bounds from iterative refinement for the computed solution to a system of equations A*X = B, where A is a general tridiagonal matrix of order n and op(A) = A or A**T, depending on TRANS. (Yes, the tridiagonal girlies are eating.) β¨
π NAME π
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/LIN/zgtt05.f
That's the file path, no cap. π«Ά
π«Ά SYNOPSIS π«Ά
Functions/Subroutines
ZGTT05β the function/subroutine you came for, in all its icon behavior.
subroutine ZGTT05 (trans, n, nrhs, dl, d, du, b, ldb, x, ldx, xact,
ldxact, ferr, berr, reslts)
π§ Function/Subroutine Documentation π§
subroutine ZGTT05
subroutine ZGTT05 (character trans, integer n, integer nrhs, complex*16,
dimension( * ) dl, complex*16, dimension( * ) d, complex*16, dimension(
* ) du, 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, double
precision, dimension( * ) berr, double precision, dimension( * ) reslts)
ZGTT05 β yes bestie, that's the whole vibe.
π― Purpose: π―
okay so the purpose is literally: ZGTT05 tests the error bounds from iterative refinement for the computed solution to a system of equations A*X = B, where A is a general tridiagonal matrix of order n and op(A) = A or A**T, depending on TRANS.
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 / ( NZ*EPS + (*) ), where
(*) = NZ*UNFL / (min_i (abs(op(A))*abs(X) +abs(b))_i )
and NZ = max. number of nonzeros in any row of A, plus 1
That's the formula glow-up. π
πΈ Parameters πΈ
TRANSβCHARACTER*1. 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βINTEGER. The number of rows of the matricesXandXACT.N >= 0.NRHSβINTEGER. The number of columns of the matricesXandXACT.NRHS >= 0.DLβCOMPLEX*16 array, dimension (N-1). The(n-1)sub-diagonal elements ofA.DβCOMPLEX*16 array, dimension (N). The diagonal elements ofA.DUβCOMPLEX*16 array, dimension (N-1). The(n-1)super-diagonal elements ofA.BβCOMPLEX*16 array, dimension (LDB,NRHS). The right hand side vectors for the system of linear equations.LDBβINTEGER. The leading dimension of the arrayB.LDB >= max(1,N).XβCOMPLEX*16 array, dimension (LDX,NRHS). The computed solution vectors. Each vector is stored as a column of the matrixX.LDXβINTEGER. The leading dimension of the arrayX.LDX >= max(1,N).XACTβCOMPLEX*16 array, dimension (LDX,NRHS). The exact solution vectors. Each vector is stored as a column of the matrixXACT.LDXACTβINTEGER. The leading dimension of the arrayXACT.LDXACT >= max(1,N).FERRβDOUBLE PRECISION array, dimension (NRHS). 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.BERRβDOUBLE PRECISION array, dimension (NRHS). The componentwise relative backward error of each solution vector (i.e., the smallest relative change in any entry ofAorBthat makesXan exact solution).RESLTSβDOUBLE PRECISION array, dimension (2). The maximum over theNRHSsolution vectors of the ratios:
RESLTS(1) = norm(X - XACT) / ( norm(X) * FERR )
RESLTS(2) = BERR / ( NZ*EPS + (*) )
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Definition at line 163 of file zgtt05.f.
π Author (auto-generated edition) π
Generated automatically by Doxygen for LAPACK from the source code.
π okay bye bestie π β this footer energy is LAPACK Version 3.12.0, from the path /home/abui...IN/zgtt05.f(3). There's literally no date in the footer, so no date to give you; don't come for me π. Mwah! π«Ά
π¦ 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