STPSV/DTPSV/CTPSV/ZTPSV

Solve triangular system

Lower triangular matrix

If A is

11

0

0

0

21

22

0

0

31

32

33

0

41

42

43

44

then A is packed column by column into an array ap as follows:

k

1

2

3

4

5

6

7

8

9

10

ap(k)

11

21

31

41

22

32

42

33

43

44

 

 

 

 

 

 

 

 

 

 

 

Lower-triangular matrix element aij is stored in array element ap(i+((j−1)⋅(2n−j))/2).

Usage

VECLIB:

 

 

CHARACTER*1

uplo, trans, diag

 

INTEGER*4

n, incx

 

REAL*4

ap(lenap), x(lenx)

CALL STPSV(uplo, trans, diag, n, ap, x, incx)

CHARACTER*1

uplo, trans, diag

INTEGER*4

n, incx

REAL*8

ap(lenap), x(lenx)

CALL DTPSV(uplo, trans, diag, n, ap, x, incx)

CHARACTER*1

uplo, trans, diag

INTEGER*4

n, incx

COMPLEX*8

ap(lenap), x(lenx)

CALL CTPSV(uplo, trans, diag, n, ap, x, incx)

CHARACTER*1

uplo, trans, diag

INTEGER*4

n, incx

COMPLEX*16

ap(lenap), x(lenx)

CALL ZTPSV(uplo, trans, diag, n, ap, x, incx)

VECLIB8:

CHARACTER*1

uplo, trans, diag

INTEGER*8

n, incx

REAL*4

ap(lenap), x(lenx)

CALL STPSV(uplo, trans, diag, n, ap, x, incx)

314HP MLIB User’s Guide