SUBROUTINE DTBTRS( UPLO, TRANSA, DIAG, N, NDIAG, NRHS, A, LDA, B, * LDB, INFO) CHARACTER * 1 UPLO, TRANSA, DIAG INTEGER N, NDIAG, NRHS, LDA, LDB, INFO DOUBLE PRECISION A(LDA,*), B(LDB,*) SUBROUTINE DTBTRS_64( UPLO, TRANSA, DIAG, N, NDIAG, NRHS, A, LDA, B, * LDB, INFO) CHARACTER * 1 UPLO, TRANSA, DIAG INTEGER*8 N, NDIAG, NRHS, LDA, LDB, INFO DOUBLE PRECISION A(LDA,*), B(LDB,*)
SUBROUTINE TBTRS( UPLO, TRANSA, DIAG, N, NDIAG, NRHS, A, [LDA], B, * [LDB], [INFO]) CHARACTER(LEN=1) :: UPLO, TRANSA, DIAG INTEGER :: N, NDIAG, NRHS, LDA, LDB, INFO REAL(8), DIMENSION(:,:) :: A, B SUBROUTINE TBTRS_64( UPLO, TRANSA, DIAG, N, NDIAG, NRHS, A, [LDA], * B, [LDB], [INFO]) CHARACTER(LEN=1) :: UPLO, TRANSA, DIAG INTEGER(8) :: N, NDIAG, NRHS, LDA, LDB, INFO REAL(8), DIMENSION(:,:) :: A, B
void dtbtrs(char uplo, char transa, char diag, int n, int ndiag, int nrhs, double *a, int lda, double *b, int ldb, int *info);
void dtbtrs_64(char uplo, char transa, char diag, long n, long ndiag, long nrhs, double *a, long lda, double *b, long ldb, long *info);
where A is a triangular band matrix of order N, and B is an N-by NRHS matrix. A check is made to verify that A is nonsingular.