TESTING/LIN/dtbt02.f(3)
okay bestie, welcome to TESTING/LIN/dtbt02.f(3) from the Library Functions Manual β the LAPACK page for a subroutine that is literally the vibe check for triangular band solves. πβ¨
πΈ NAME πΈ
TESTING/LIN/dtbt02.f β that's the manual page, no cap. π
β¨ SYNOPSIS β¨
π¦ Functions/Subroutines π¦
Here's the signature, bestie:
subroutine dtbt02 (uplo, trans, diag, n, kd, nrhs, ab, ldab, x, ldx, b,
ldb, work, resid)
and the name is DTBT02 β iconic. π
π FUNCTION/SUBROUTINE DOCUMENTATION π
subroutine dtbt02 (character uplo, character trans, character diag, integer
n, integer kd, integer nrhs, double precision, dimension( ldab, * ) ab,
integer ldab, double precision, dimension( ldx, * ) x, integer ldx,
double precision, dimension( ldb, * ) b, integer ldb, double precision,
dimension( * ) work, double precision resid)
DTBT02 β the full documentation moment. π«Ά
π Purpose π
okay so DTBT02 computes the residual for the computed solution to a triangular system of linear equations op(A)*X = B, when A is a triangular band matrix. The test ratio is the maximum over
norm(b - op(A)*x) / ( ||op(A)||_1 * norm(x) * EPS ),
where op(A) = A or A**T, b is the column of B, x is the solution vector, and EPS is the machine epsilon. The norm used is the 1-norm. slay. π
π§ Parameters π§
Every parameter, every constraint, every little detail β here we go, bestie:
UPLOβCHARACTER*1. Specifies whether the matrixAis upper or lower triangular.= 'U': Upper triangular= 'L': Lower triangular
TRANSβCHARACTER*1. Specifies the operation applied toA.= 'N': A * X = B (No transpose)= 'T': A**T * X = B (Transpose)= 'C': A**H * X = B (Conjugate transpose = Transpose)
DIAGβCHARACTER*1. Specifies whether or not the matrixAis unit triangular.= 'N': Non-unit triangular= 'U': Unit triangular
NβINTEGER. The order of the matrixA.N >= 0.KDβINTEGER. The number of superdiagonals or subdiagonals of the triangular band matrixA.KD >= 0.NRHSβINTEGER. The number of right hand sides, i.e., the number of columns of the matricesXandB.NRHS >= 0.ABβDOUBLE PRECISION array, dimension (LDAB,N). The upper or lower triangular band matrixA, stored in the firstkd+1rows of the array. Thej-th column ofAis stored in thej-th column of the arrayABas follows: ifUPLO = 'U',AB(kd+1+i-j,j) = A(i,j)formax(1,j-kd)<=i<=j; ifUPLO = 'L',AB(1+i-j,j) = A(i,j)forj<=i<=min(n,j+kd).LDABβINTEGER. The leading dimension of the arrayAB.LDAB >= KD+1.XβDOUBLE PRECISION array, dimension (LDX,NRHS). The computed solution vectors for the system of linear equations.LDXβINTEGER. The leading dimension of the arrayX.LDX >= max(1,N).BβDOUBLE PRECISION 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).WORKβDOUBLE PRECISION array, dimension (N).RESIDβDOUBLE PRECISION. The maximum over the number of right hand sides ofnorm(op(A)*x - b) / ( norm(op(A)) * norm(x) * EPS ).
π Author π
the academic credits, periodt:
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 152 of file dtbt02.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
and that's the whole page, bestie. footer vibe: LAPACK Version 3.12.0, source TESTING/LIN/dtbt02.f(3). date? the footer doesn't give one, so I'm not making one up β no cap. π
β¨
π¦ 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