dtpsv - solve one of the systems of equations A*x = b, or A'*x = b
SUBROUTINE DTPSV( UPLO, TRANSA, DIAG, N, A, Y, INCY) CHARACTER * 1 UPLO, TRANSA, DIAG INTEGER N, INCY DOUBLE PRECISION A(*), Y(*)
SUBROUTINE DTPSV_64( UPLO, TRANSA, DIAG, N, A, Y, INCY) CHARACTER * 1 UPLO, TRANSA, DIAG INTEGER*8 N, INCY DOUBLE PRECISION A(*), Y(*)
SUBROUTINE TPSV( UPLO, [TRANSA], DIAG, [N], A, Y, [INCY]) CHARACTER(LEN=1) :: UPLO, TRANSA, DIAG INTEGER :: N, INCY REAL(8), DIMENSION(:) :: A, Y
SUBROUTINE TPSV_64( UPLO, [TRANSA], DIAG, [N], A, Y, [INCY]) CHARACTER(LEN=1) :: UPLO, TRANSA, DIAG INTEGER(8) :: N, INCY REAL(8), DIMENSION(:) :: A, Y
#include <sunperf.h>
void dtpsv(char uplo, char transa, char diag, int n, double *a, double *y, int incy);
void dtpsv_64(char uplo, char transa, char diag, long n, double *a, double *y, long incy);
dtpsv solves one of the systems of equations A*x = b, or A'*x = b, where b and x are n element vectors and A is an n by n unit, or non-unit, upper or lower triangular matrix, supplied in packed form.
No test for singularity or near-singularity is included in this routine. Such tests must be performed before calling this routine.
UPLO = 'U' or 'u' A is an upper triangular matrix.
UPLO = 'L' or 'l' A is a lower triangular matrix.
Unchanged on exit.
TRANSA = 'N' or 'n' A*x = b.
TRANSA = 'T' or 't' A'*x = b.
TRANSA = 'C' or 'c' A'*x = b.
Unchanged on exit.
DIAG = 'U' or 'u' A is assumed to be unit triangular.
DIAG = 'N' or 'n' A is not assumed to be unit triangular.
Unchanged on exit.