TESTING/LIN/ctbt02.f(3)
πΈ NAME πΈ
TESTING/LIN/ctbt02.f β that's the manual-page main character, bestie π and yes, this is the TESTING/LIN/ctbt02.f(3) page from the Library Functions Manual.
β¨ SYNOPSIS β¨
okay so the synopsis is giving "Functions/Subroutines" energy, periodt:
Functions/Subroutines
subroutine ctbt02 (uplo, trans, diag, n, kd, nrhs, ab, ldab, x, ldx, b,
ldb, work, rwork, resid)
CTBT02
π Function/Subroutine Documentation π
the declaration is right here, no cap:
subroutine ctbt02 (character uplo, character trans, character diag, integer
n, integer kd, integer nrhs, complex, dimension( ldab, * ) ab, integer
ldab, complex, dimension( ldx, * ) x, integer ldx, complex, dimension(
ldb, * ) b, integer ldb, complex, dimension( * ) work, real, dimension(
* ) rwork, real resid)
CTBT02
π― Purpose π―
CTBT02 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, A**T, or A**H, b is the column of B, x is the solution vector, and EPS is the machine epsilon. slay, math bestie π
ποΈ Parameters ποΈ
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)
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βCOMPLEX array, dimension (LDA,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) for max(1,j-kd)<=i<=j; ifUPLO = 'L',AB(1+i-j,j) = A(i,j) for j<=i<=min(n,j+kd).Keeping thatLDA,Nexactly as the source says, bestie πLDABβINTEGER. The leading dimension of the arrayAB.LDAB >= max(1,KD+1).XβCOMPLEX 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βCOMPLEX 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βCOMPLEX array, dimension (N).RWORKβREAL array, dimension (N).RESIDβREAL. The maximum over the number of right hand sides ofnorm(op(A)*x - b) / ( norm(op(A)) * norm(x) * EPS ).
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 157 of file ctbt02.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
and that's the page, bestie π«Άπ β source/version sign-off: LAPACK Version 3.12.0, page TESTING/LIN/ctbt02.f(3), from the Library Functions Manual. The footer doesn't show a date on this one, so no date to report, periodt πβ¨
π¦ 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