CTRT05(3)
β¨ NAME β¨
okay bestie, this is CTRT05(3) from the Library Functions Manual β header path is /home/abui...IN/ctrt05.f(3), and the full file path is:
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/LIN/ctrt05.f
that path is literally the NAME, no cap π
π¦ SYNOPSIS π¦
β¨ Functions/Subroutines β¨
subroutine CTRT05 (uplo, trans, diag, n, nrhs, a, lda, b, ldb, x, ldx,
xact, ldxact, ferr, berr, reslts)
CTRT05
π Function/Subroutine Documentation π
subroutine CTRT05 (character uplo, character trans, character diag, integer
n, integer nrhs, complex, dimension( lda, * ) a, integer lda, complex,
dimension( ldb, * ) b, integer ldb, complex, dimension( ldx, * ) x,
integer ldx, complex, dimension( ldxact, * ) xact, integer ldxact, real,
dimension( * ) ferr, real, dimension( * ) berr, real, dimension( * )
reslts)
CTRT05
π¦ Purpose π¦
okay so the whole vibe of CTRT05 is: it tests the error bounds from iterative refinement for the computed solution to a system of equations A*X = B, where A is a triangular n by n matrix. iconic honestly.
RESLTS(1) = test of the error bound
= norm(X - XACT) / ( norm(X) * FERR )
if that ratio is not less than one, a large value is returned β thatβs the vibe check failing, bestie π
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 )
π Parameters π
UPLOβUPLO is CHARACTER*1. tells you whether the matrixAis upper or lower triangular.= 'U':Upper triangular.= 'L':Lower triangular. πTRANSβTRANS is CHARACTER*1. tells you the form of the system of equations.= 'N':A * X = B(No transpose).= 'T':A'* X = B(Transpose).= 'C':A'* X = B(Conjugate transpose = Transpose). β¨DIAGβDIAG is CHARACTER*1. specifies whether or not the matrixAis unit triangular.= 'N':Non-unit triangular.= 'U':Unit triangular. πNβN is INTEGER. the number of rows of the matricesX,B, andXACT, and the order of the matrixA.N >= 0. πΈNRHSβNRHS is INTEGER. the number of columns of the matricesX,B, andXACT.NRHS >= 0. π«ΆAβA is COMPLEX array, dimension (LDA,N). the triangular matrixA. ifUPLO = 'U', the leadingnbynupper triangular part of the arrayAcontains the upper triangular matrix, and the strictly lower triangular part ofAis not referenced. ifUPLO = 'L', the leadingnbynlower triangular part of the arrayAcontains the lower triangular matrix, and the strictly upper triangular part ofAis not referenced. ifDIAG = 'U', the diagonal elements ofAare also not referenced and are assumed to be1. πLDAβLDA is INTEGER. the leading dimension of the arrayA.LDA >= max(1,N). πBβB is COMPLEX array, dimension (LDB,NRHS). the right hand side vectors for the system of linear equations. πLDBβLDB is INTEGER. the leading dimension of the arrayB.LDB >= max(1,N). πXβX is COMPLEX array, dimension (LDX,NRHS). the computed solution vectors. each vector is stored as a column of the matrixX. π¦LDXβLDX is INTEGER. the leading dimension of the arrayX.LDX >= max(1,N). πXACTβXACT is COMPLEX array, dimension (LDX,NRHS). the exact solution vectors. each vector is stored as a column of the matrixXACT. π§LDXACTβLDXACT is INTEGER. the leading dimension of the arrayXACT.LDXACT >= max(1,N). π₯FERRβFERR is REAL 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βBERR is REAL 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βRESLTS is REAL array, dimension (2). 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 180 of file ctrt05.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
cute sign-off time, bestie: this footer says LAPACK Version 3.12.0 and /home/abui...IN/ctrt05.f(3) β no separate date is shown, just those header/footer goodies. 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