HP MLIB User's Guide Vol. 1 7th Ed.

400 HP MLIB User’s Guide
F_STBSV/F_DTBSV/F_CTBSV/F_ZTBSV Triangular banded solve
Name F_STBSV/F_DTBSV/F_CTBSV/F_ZTBSV
Triangular banded solve
Purpose F_xTBSV solves one of the following equations:
where x is a vector and the matrix T is a unit, non-unit, upper, or lower
triangular banded matrix, T
−T
is the inverse of the transpose of T, and is
the inverse of the conjugate transpose of T.
Refer to “STBSV/DTBSV/CTBSV/ZTBSV” on page 301 for a description of the
equivalent HP MLIB legacy BLAS subprograms.
Matrix
Storage
For these subprograms, you supply A in a two-dimensional array large enough
to hold a square matrix. The other triangle of the array is not referenced. If A
has an unstored unit diagonal (see input argument DIAG), then the diagonal
elements of the array also are not referenced.
Usage VECLIB
INTEGER*4 DIAG, INCX, K, N, TRANS, UPLO
REAL*4 ALPHA, A( LDA, * ), X( * )
SUBROUTINE F_STBSV (UPLO, TRANS, DIAG, N, K, ALPHA, A, LDA, X,
INCX)
INTEGER*4 DIAG, INCX, K, N, TRANS, UPLO
REAL*8 ALPHA, A( LDA, * ), X( * )
SUBROUTINE F_DTBSV (UPLO, TRANS, DIAG, N, K, ALPHA, A, LDA, X,
INCX)
INTEGER*4 DIAG, INCX, K, N, TRANS, UPLO
COMPLEX*8 ALPHA, A( LDA, * ), X( * )
SUBROUTINE F_CTBSV (UPLO, TRANS, DIAG, N, K, ALPHA, A, LDA, X,
INCX)
INTEGER*4 DIAG, INCX, K, N, TRANS, UPLO
COMPLEX*16 ALPHA, A( LDA, * ), X( * )
SUBROUTINE F_ZTBSV (UPLO, TRANS, DIAG, N, K, ALPHA, A, LDA, X,
INCX)
VECLIB8
x αT
1–
x←
x αT
T–
x←
x αT
*–
x←
T
*–