SRC/ctbrfs.f(3)
π SRC/ctbrfs.f(3) β Library Functions Manual π
π NAME π
SRC/ctbrfs.f β yep, that's the name, bestie. Just the file path. Short, sweet, main character energy.
π SYNOPSIS π
π¦ Functions/Subroutines π¦
subroutine ctbrfs (uplo, trans, diag, n, kd, nrhs, ab, ldab, b, ldb, x,
ldx, ferr, berr, work, rwork, info)
CTBRFS
π¦ Function/Subroutine Documentation π¦
subroutine ctbrfs (character uplo, character trans, character diag, integer
n, integer kd, integer nrhs, complex, dimension( ldab, * ) ab, integer
ldab, complex, dimension( ldb, * ) b, integer ldb, complex, dimension(
ldx, * ) x, integer ldx, real, dimension( * ) ferr, real, dimension( * )
berr, complex, dimension( * ) work, real, dimension( * ) rwork, integer
info)
CTBRFS
CTBRFS is literally the subroutine, bestie. She's giving linear algebra error-bound realness.
π« Purpose π«
Okay so CTBRFS provides error bounds and backward error estimates for the solution to a system of linear equations with a triangular band coefficient matrix. No cap.
The solution matrix X must be computed by CTBTRS or some other means before entering this routine. CTBRFS does not do iterative refinement because doing so cannot improve the backward error. That's just the tea, fr fr.
πΈ Parameters πΈ
Bestie, here's every single parameter β don't skip a single one.
UPLOβCHARACTER*1. Tells you which triangle is living its truth:= 'U':Ais upper triangular;= 'L':Ais lower triangular.
TRANSβCHARACTER*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)
DIAGβCHARACTER*1. Unit triangular or non-unit triangular, that is the question:= 'N':Ais non-unit triangular;= 'U':Ais 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 matricesBandX.NRHS >= 0.ABβCOMPLEX 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:- if
UPLO = 'U',AB(kd+1+i-j,j) = A(i,j)formax(1,j-kd)<=i<=j; - if
UPLO = 'L',AB(1+i-j,j) = A(i,j)forj<=i<=min(n,j+kd).
- if
If DIAG = 'U', the diagonal elements of A are not referenced and are assumed to be 1.
LDABβINTEGER. The leading dimension of the arrayAB.LDAB >= KD+1.BβCOMPLEX array, dimension (LDB,NRHS). The right hand side matrixB.LDBβINTEGER. The leading dimension of the arrayB.LDB >= max(1,N).XβCOMPLEX array, dimension (LDX,NRHS). The solution matrixX.LDXβINTEGER. The leading dimension of the arrayX.LDX >= max(1,N).FERRβREAL array, dimension (NRHS). The estimated forward error bound for each solution vectorX(j)(thej-th column of the solution matrixX). IfXTRUEis the true solution corresponding toX(j),FERR(j)is an estimated upper bound for the magnitude of the largest element in(X(j) - XTRUE)divided by the magnitude of the largest element inX(j). The estimate is as reliable as the estimate forRCOND, and is almost always a slight overestimate of the true error.BERRβREAL array, dimension (NRHS). The componentwise relative backward error of each solution vectorX(j)(i.e., the smallest relative change in any element ofAorBthat makesX(j)an exact solution).WORKβCOMPLEX array, dimension (2*N).RWORKβREAL array, dimension (N).INFOβINTEGER. Exit status tea:= 0: successful exit< 0: ifINFO = -i, thei-th argument had an illegal value
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Definition at line186of filectbrfs.f.
β¨ Author β¨
Generated automatically by Doxygen for LAPACK from the source code. iykyk.
π That's the whole vibe check, bestie β the footer signs off with LAPACK, Version 3.12.0, and SRC/ctbrfs.f(3) (no date shown, fr fr). π
π¦ 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