TESTING/LIN/cgbt02.f(3)
π TESTING/LIN/cgbt02.f π
okay bestie, welcome to the <pink>TESTING/LIN/cgbt02.f(3)</pink> page β yes, straight from the <purple>Library Functions Manual</purple>, and she is here to explain CGBT02 in full, no skips, no cap. π
β¨ NAME β¨
TESTING/LIN/cgbt02.f β the LAPACK testing routine CGBT02, and yes, sheβs iconic. π
π¦ SYNOPSIS π¦
Functions/Subroutines
subroutine cgbt02 (trans, m, n, kl, ku, nrhs, a, lda, x, ldx, b, ldb,
rwork, resid)
The routine name is CGBT02. π
π§ Function/Subroutine Documentation π§
Hereβs the full signature, character for character:
subroutine cgbt02 (character trans, integer m, integer n, integer kl,
integer ku, integer nrhs, complex, dimension( lda, * ) a, integer lda,
complex, dimension( ldx, * ) x, integer ldx, complex, dimension( ldb, *
) b, integer ldb, real, dimension( * ) rwork, real resid)
CGBT02
π Purpose π
okay so the whole point is this, bestie:
CGBT02 computes the residual for a solution of a banded system of
equations op(A)*X = B:
RESID = norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ),
where op(A) = A, A**T, or A**H, depending on TRANS, and EPS is the
machine epsilon.
So sheβs computing that residual, exactly with the formula above. β¨
π Parameters π
TRANSβ isCHARACTER*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)
Mβ isINTEGER. The number of rows of the matrix A.M >= 0.Nβ isINTEGER. The number of columns of the matrix A.N >= 0.KLβ isINTEGER. The number of subdiagonals within the band of A.KL >= 0.KUβ isINTEGER. The number of superdiagonals within the band of A.KU >= 0.NRHSβ isINTEGER. The number of columns of B.NRHS >= 0.Aβ isCOMPLEX array, dimension (LDA,N). The original matrix A in band storage, stored in rows 1 toKL+KU+1.LDAβ isINTEGER. The leading dimension of the array A.LDA >= max(1,KL+KU+1).Xβ isCOMPLEX array, dimension (LDX,NRHS). The computed solution vectors for the system of linear equations.LDXβ isINTEGER. The leading dimension of the array X. IfTRANS = 'N',LDX >= max(1,N); ifTRANS = 'T'or'C',LDX >= max(1,M).Bβ isCOMPLEX array, dimension (LDB,NRHS). On entry, the right hand side vectors for the system of linear equations. On exit, B is overwritten with the differenceB - A*X.LDBβ isINTEGER. The leading dimension of the array B. IFTRANS = 'N',LDB >= max(1,M); ifTRANS = 'T'or'C',LDB >= max(1,N).RWORKβ isREAL array, dimension (MAX(1,LRWORK)), whereLRWORK >= MwhenTRANS = 'T'or'C'; otherwise, RWORK is not referenced.RESIDβ isREAL. The maximum over the number of right hand sides ofnorm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ).
lowkey important:RWORKonly gets referenced whenTRANS = 'T'or'C'β otherwise sheβs literally not invited. π«Ά
πΈ Author πΈ
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 146 of file cgbt02.f. π
π Author π
Generated automatically by Doxygen for LAPACK from the source code. π
okay thatβs the whole page, bestie! final footer vibe: <rainbow>LAPACK Version 3.12.0</rainbow>, source TESTING/LIN/cgbt02.f(3), no date printed, generated automatically by Doxygen for LAPACK from the source code. bye! πβ¨π
π¦ 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