dpbsl
dpbsl - (obsolete) section solve the linear system Ax = b for a symmetric positive definite
matrix A in banded storage, which has been Cholesky-factored by SPBCO or
SPBFA, and vectors b and x.
SUBROUTINE DPBSL( A, LDA, N, NDIAG, B)
INTEGER LDA, N, NDIAG
DOUBLE PRECISION A(LDA,*), B(*)
SUBROUTINE DPBSL_64( A, LDA, N, NDIAG, B)
INTEGER*8 LDA, N, NDIAG
DOUBLE PRECISION A(LDA,*), B(*)
#include <sunperf.h>
void dpbsl(double *a, int lda, int n, int ndiag, double *b);
void dpbsl_64(double *a, long lda, long n, long ndiag, double *b);
-
* A (input)
-
Cholesky factorization of the matrix A, as computed by SPBCO or SPBFA.
-
* LDA (input)
-
Leading dimension of the array A as specified in a dimension or type
statement. LDA >= NDIAG + 1.
-
* N (input)
-
Order of the matrix A. N >= 0.
-
* NDIAG (input)
-
Number of diagonals of matrix A. N-1 >= NDIAG >= 0 but if
N = 0 then NDIAG = 0.
-
* B (input/output)
-
On entry, the right-hand side vector b.
On exit, the solution vector x.