SRC/ztbrfs.f(3)
okay bestie, welcome to SRC/ztbrfs.f(3) from the Library Functions Manual β this is the LAPACK tea, and itβs lowkey iconic β¨
β¨ NAME β¨
SRC/ztbrfs.f β thatβs the name, and thatβs literally the whole NAME section. no notes π
π SYNOPSIS π
π Functions/Subroutines π
Functions/Subroutines
subroutine ztbrfs (uplo, trans, diag, n, kd, nrhs, ab, ldab, b, ldb, x,
ldx, ferr, berr, work, rwork, info)
ZTBRFS
π¦ Function/Subroutine Documentation π¦
subroutine ztbrfs (character uplo, character trans, character diag, integer
n, integer kd, integer nrhs, complex*16, dimension( ldab, * ) ab,
integer ldab, complex*16, dimension( ldb, * ) b, integer ldb,
complex*16, dimension( ldx, * ) x, integer ldx, double precision,
dimension( * ) ferr, double precision, dimension( * ) berr, complex*16,
dimension( * ) work, double precision, dimension( * ) rwork, integer
info)
ZTBRFS
π Purpose π
ZTBRFS provides error bounds and backward error estimates for the solution to a system of linear equations with a triangular band coefficient matrix.
π¨ bestie, important: the solution matrixXmust be computed byZTBTRSor some other means before entering this routine.ZTBRFSdoes not do iterative refinement because doing so cannot improve the backward error.
πΈ Parameters πΈ
UPLOβUPLO is CHARACTER*1. it tells us which triangle weβre dealing with:= 'U':Ais upper triangular;= 'L':Ais lower triangular.
TRANSβTRANS is 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βDIAG is CHARACTER*1:= 'N':Ais non-unit triangular;= 'U':Ais unit triangular.
NβN is INTEGER. the order of the matrixA.N >= 0.KDβKD is INTEGER. the number of superdiagonals or subdiagonals of the triangular band matrixA.KD >= 0.NRHSβNRHS is INTEGER. the number of right hand sides, i.e., the number of columns of the matricesBandX.NRHS >= 0.ABβAB is COMPLEX*16 array, dimension (LDAB,N). the upper or lower triangular band matrixA, stored in the firstkd+1rows of the array. the j-th column ofAis stored in the j-th column of the arrayABas follows:
if UPLO = 'U', AB(kd+1+i-j,j) = A(i,j) for max(1,j-kd)<=i<=j;
if UPLO = 'L', AB(1+i-j,j) = A(i,j) for j<=i<=min(n,j+kd).
if DIAG = 'U', the diagonal elements of A are not referenced and are assumed to be 1.
LDABβLDAB is INTEGER. the leading dimension of the arrayAB.LDAB >= KD+1.BβB is COMPLEX*16 array, dimension (LDB,NRHS). the right hand side matrixB.LDBβLDB is INTEGER. the leading dimension of the arrayB.LDB >= max(1,N).XβX is COMPLEX*16 array, dimension (LDX,NRHS). the solution matrixX.LDXβLDX is INTEGER. the leading dimension of the arrayX.LDX >= max(1,N).FERRβFERR is DOUBLE PRECISION array, dimension (NRHS). the estimated forward error bound for each solution vectorX(j)(the j-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βBERR is DOUBLE PRECISION 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βWORK is COMPLEX*16 array, dimension (2*N).RWORKβRWORK is DOUBLE PRECISION array, dimension (N).INFOβINFO is INTEGER:= 0: successful exit< 0: ifINFO = -i, the i-th argument had an illegal value
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 186 of file ztbrfs.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
okay thatβs the whole page, bestie β source/version: LAPACK Version 3.12.0, and the page is SRC/ztbrfs.f(3). slay and stay triangular π
β¨
π¦ 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