dptsv - compute the solution to a real system of linear equations A*X = B, where A is an N-by-N symmetric positive definite tridiagonal matrix, and X and B are N-by-NRHS matrices.
SUBROUTINE DPTSV( N, NRHS, DIAG, SUB, B, LDB, INFO) INTEGER N, NRHS, LDB, INFO DOUBLE PRECISION DIAG(*), SUB(*), B(LDB,*)
SUBROUTINE DPTSV_64( N, NRHS, DIAG, SUB, B, LDB, INFO) INTEGER*8 N, NRHS, LDB, INFO DOUBLE PRECISION DIAG(*), SUB(*), B(LDB,*)
SUBROUTINE PTSV( [N], [NRHS], DIAG, SUB, B, [LDB], [INFO]) INTEGER :: N, NRHS, LDB, INFO REAL(8), DIMENSION(:) :: DIAG, SUB REAL(8), DIMENSION(:,:) :: B
SUBROUTINE PTSV_64( [N], [NRHS], DIAG, SUB, B, [LDB], [INFO]) INTEGER(8) :: N, NRHS, LDB, INFO REAL(8), DIMENSION(:) :: DIAG, SUB REAL(8), DIMENSION(:,:) :: B
#include <sunperf.h>
void dptsv(int n, int nrhs, double *diag, double *sub, double *b, int ldb, int *info);
void dptsv_64(long n, long nrhs, double *diag, double *sub, double *b, long ldb, long *info);
dptsv computes the solution to a real system of linear equations A*X = B, where A is an N-by-N symmetric positive definite tridiagonal matrix, and X and B are N-by-NRHS matrices.
A is factored as A = L*D*L**T, and the factored form of A is then used to solve the system of equations.
= 0: successful exit
< 0: if INFO = -i, the i-th argument had an illegal value
> 0: if INFO = i, the leading minor of order i is not positive definite, and the solution has not been computed. The factorization has not been completed unless i = N.