1DTPSV(1)                         BLAS routine                         DTPSV(1)
2
3
4

NAME

6       DTPSV - solves one of the systems of equations   A*x = b, or A'*x = b,
7

SYNOPSIS

9       SUBROUTINE DTPSV(UPLO,TRANS,DIAG,N,AP,X,INCX)
10
11           INTEGER                                   INCX,N
12
13           CHARACTER                                 DIAG,TRANS,UPLO
14
15           DOUBLE                                    PRECISION AP(*),X(*)
16

PURPOSE

18       DTPSV  solves one of the systems of equations
19
20       where  b  and  x are n element vectors and A is an n by n unit, or non-
21       unit, upper or lower triangular matrix, supplied in packed form.
22
23       No test for singularity or near-singularity is included  in  this  rou‐
24       tine. Such tests must be performed before calling this routine.
25
26

ARGUMENTS

28       UPLO   - CHARACTER*1.
29              On entry, UPLO specifies whether the matrix is an upper or lower
30              triangular matrix as follows:
31
32              UPLO = 'U' or 'u'   A is an upper triangular matrix.
33
34              UPLO = 'L' or 'l'   A is a lower triangular matrix.
35
36              Unchanged on exit.
37
38       TRANS  - CHARACTER*1.
39              On entry, TRANS specifies the equations to be solved as follows:
40
41              TRANS = 'N' or 'n'   A*x = b.
42
43              TRANS = 'T' or 't'   A'*x = b.
44
45              TRANS = 'C' or 'c'   A'*x = b.
46
47              Unchanged on exit.
48
49       DIAG   - CHARACTER*1.
50              On entry, DIAG specifies whether or not A is unit triangular  as
51              follows:
52
53              DIAG = 'U' or 'u'   A is assumed to be unit triangular.
54
55              DIAG = 'N' or 'n'   A is not assumed to be unit triangular.
56
57              Unchanged on exit.
58
59       N      - INTEGER.
60              On  entry,  N specifies the order of the matrix A.  N must be at
61              least zero.  Unchanged on exit.
62
63       AP     - DOUBLE PRECISION array of DIMENSION at least
64              ( ( n*( n + 1 ) )/2 ).  Before entry with  UPLO =  'U'  or  'u',
65              the  array  AP  must  contain the upper triangular matrix packed
66              sequentially, column by column, so that AP( 1 ) contains a( 1, 1
67              ),  AP(  2 ) and AP( 3 ) contain a( 1, 2 ) and a( 2, 2 ) respec‐
68              tively, and so on.  Before entry with UPLO =  'L'  or  'l',  the
69              array AP must contain the lower triangular matrix packed sequen‐
70              tially, column by column, so that AP( 1 ) contains a(  1,  1  ),
71              AP(  2  )  and  AP(  3 ) contain a( 2, 1 ) and a( 3, 1 ) respec‐
72              tively, and so on.  Note that when  DIAG = 'U' or 'u', the diag‐
73              onal  elements  of  A  are not referenced, but are assumed to be
74              unity.  Unchanged on exit.
75
76       X      - DOUBLE PRECISION array of dimension at least
77              ( 1 + ( n - 1 )*abs( INCX ) ).  Before  entry,  the  incremented
78              array  X must contain the n element right-hand side vector b. On
79              exit, X is overwritten with the solution vector x.
80
81       INCX   - INTEGER.
82              On entry, INCX specifies the increment for the  elements  of  X.
83              INCX must not be zero.  Unchanged on exit.
84

FURTHER DETAILS

86       Level 2 Blas routine.
87
88       -- Written on 22-October-1986.
89          Jack Dongarra, Argonne National Lab.
90          Jeremy Du Croz, Nag Central Office.
91          Sven Hammarling, Nag Central Office.
92          Richard Hanson, Sandia National Labs.
93
94
95
96
97BLAS routine                     November 2008                        DTPSV(1)
Impressum