ZTPT05(3)
πΈ NAME πΈ
okay so bestie, this page is ztpt05.f(3) from the Library Functions Manual, and the NAME is literally this whole path:
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/LIN/ztpt05.f
β¨ SYNOPSIS β¨
π Functions/Subroutines π
subroutine ZTPT05 (uplo, trans, diag, n, nrhs, ap, b, ldb, x, ldx, xact, ldxact, ferr, berr, reslts)βZTPT05
subroutine ZTPT05 (uplo, trans, diag, n, nrhs, ap, b, ldb, x, ldx, xact,
ldxact, ferr, berr, reslts)
ZTPT05
π«Ά Function/Subroutine Documentation π«Ά
subroutine ZTPT05 (character uplo, character trans, character diag, integer
n, integer nrhs, complex*16, dimension( * ) ap, 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)
ZTPT05
ZTPT05
π Purpose π
bestie, ZTPT05 tests the error bounds from iterative refinement for the computed solution to a system of equations A*X = B, where A is a triangular matrix in packed storage format. no cap, that's the whole job π
RESLTS(1) = test of the error bound
= norm(X - XACT) / ( norm(X) * FERR )
RESLTS(2) = residual from the iterative refinement routine
= the maximum of BERR / ( (n+1)*EPS + (*) ), where
(*) = (n+1)*UNFL / (min_i (abs(A)*abs(X) +abs(b))_i )
A large value is returned if this ratio is not less than one. so yeah, that's the vibe check π₯
π Parameters π
UPLOβUPLOisCHARACTER*1and it says whether matrixAis upper or lower triangular.= 'U':Upper triangular= 'L':Lower triangular
TRANSβTRANSisCHARACTER*1and it says what form the system of equations is in.= 'N':A * X = B(No transpose)= 'T':A'* X = B(Transpose)= 'C':A'* X = B(Conjugate transpose = Transpose)
DIAGβDIAGisCHARACTER*1and it says whether or not matrixAis unit triangular.= 'N':Non-unit triangular= 'U':Unit triangular
NβNisINTEGER. It's the number of rows of the matricesX,B, andXACT, and the order of matrixA.N >= 0.NRHSβNRHSisINTEGER. It's the number of columns of the matricesX,B, andXACT.NRHS >= 0.APβAPisCOMPLEX*16 array, dimension (N*(N+1)/2). It's the upper or lower triangular matrixA, packed columnwise in a linear array. Thej-th column ofAis stored in the arrayAPas follows:
if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.
If DIAG = 'U', the diagonal elements of A are not referenced and are assumed to be 1.
BβBisCOMPLEX*16 array, dimension (LDB,NRHS). These are the right hand side vectors for the system of linear equations.LDBβLDBisINTEGER. It's the leading dimension of the arrayB.LDB >= max(1,N).XβXisCOMPLEX*16 array, dimension (LDX,NRHS). These are the computed solution vectors. Each vector is stored as a column of the matrixX.LDXβLDXisINTEGER. It's the leading dimension of the arrayX.LDX >= max(1,N).XACTβXACTisCOMPLEX*16 array, dimension (LDX,NRHS). These are the exact solution vectors. Each vector is stored as a column of the matrixXACT.LDXACTβLDXACTisINTEGER. It's the leading dimension of the arrayXACT.LDXACT >= max(1,N).FERRβFERRisDOUBLE 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.BERRβBERRisDOUBLE 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βRESLTSisDOUBLE PRECISION array, dimension (2). It's 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 173 of file ztpt05.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
okay that's the whole page, bestie π«Ά it's from LAPACK Version 3.12.0, source path /home/abui...IN/ztpt05.f(3), and the footer didn't give us a date so we're not inventing one β slay responsibly π
β¨
π¦ 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