dptsl

NAME

dptsl - (obsolete) solve the linear system Ax = b for a symmetric positive definite tridiagonal matrix A and vectors b and x.

SYNOPSIS

  SUBROUTINE DPTSL( N, DIAG, OFFD, B)
  INTEGER N
  DOUBLE PRECISION DIAG(*), OFFD(*), B(*)
 
  SUBROUTINE DPTSL_64( N, DIAG, OFFD, B)
  INTEGER*8 N
  DOUBLE PRECISION DIAG(*), OFFD(*), B(*)

C INTERFACE

#include <sunperf.h>

void dptsl(int n, double *diag, double *offd, double *b);

void dptsl_64(long n, double *diag, double *offd, double *b);

ARGUMENTS

* N (input): Order of the matrix A. N >= 0.
* DIAG (input): Main diagonal elements of A.
* OFFD (input): Off diagonal elements of A. OFFD(1) through OFFD(N*1) contains the off diagonal elements; OFFD(N) is not referenced.
* B (input/output): On entry, the right-hand side vector b. On exit, the solution vector x.