1CLA_SYAMV(1)LAPACK routine (version 3.2)                           CLA_SYAMV(1)
2
3
4

NAME

6       CLA_SYAMV    -    performs   the   matrix-vector   operation     y   :=
7       alpha*abs(A)*abs(x) + beta*abs(y),
8

SYNOPSIS

10       SUBROUTINE CLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY )
11
12           IMPLICIT          NONE
13
14           REAL              ALPHA, BETA
15
16           INTEGER           INCX, INCY, LDA, N
17
18           INTEGER           UPLO
19
20           COMPLEX           A( LDA, * ), X( * )
21
22           REAL              Y( * )
23

PURPOSE

25       CLA_SYAMV  performs the matrix-vector operation where  alpha  and  beta
26       are scalars, x and y are vectors and A is an n by n symmetric matrix.
27       This  function  is primarily used in calculating error bounds.  To pro‐
28       tect against underflow during evaluation, components in  the  resulting
29       vector  are  perturbed  away  from  zero  by  (N+1) times the underflow
30       threshold.  To prevent unnecessarily large errors  for  block-structure
31       embedded in general matrices,
32       "symbolically" zero components are not perturbed.  A zero entry is con‐
33       sidered "symbolic" if all multiplications involved  in  computing  that
34       entry have at least one zero multiplicand.
35

ARGUMENTS

37       UPLO   - INTEGER
38              On  entry,  UPLO specifies whether the upper or lower triangular
39              part of the array A is to  be  referenced  as  follows:  UPLO  =
40              BLAS_UPPER   Only the upper triangular part of A is to be refer‐
41              enced.  UPLO = BLAS_LOWER   Only the lower triangular part of  A
42              is to be referenced.  Unchanged on exit.
43
44       N      - INTEGER.
45              On  entry, N specifies the number of columns of the matrix A.  N
46              must be at least zero.  Unchanged on exit.
47
48       ALPHA  - REAL            .
49              On entry, ALPHA specifies the scalar alpha.  Unchanged on exit.
50
51       A      - COMPLEX             array of DIMENSION ( LDA, n ).
52              Before entry, the leading m by n part of the array A  must  con‐
53              tain the matrix of coefficients.  Unchanged on exit.
54
55       LDA    - INTEGER.
56              On  entry, LDA specifies the first dimension of A as declared in
57              the calling (sub) program. LDA must be at least  max(  1,  n  ).
58              Unchanged on exit.
59
60       X      - COMPLEX             array of DIMENSION at least
61              (  1  +  (  n  - 1 )*abs( INCX ) ) Before entry, the incremented
62              array X must contain the vector x.  Unchanged on exit.
63
64       INCX   - INTEGER.
65              On entry, INCX specifies the increment for the  elements  of  X.
66              INCX must not be zero.  Unchanged on exit.
67
68       BETA   - REAL            .
69              On  entry, BETA specifies the scalar beta. When BETA is supplied
70              as zero then Y need not be set on input.  Unchanged on exit.
71
72       Y      - REAL             array of DIMENSION at least
73              ( 1 + ( n - 1 )*abs( INCY ) ) Before entry with  BETA  non-zero,
74              the incremented array Y must contain the vector y. On exit, Y is
75              overwritten by the updated vector y.
76
77       INCY   - INTEGER.
78              On entry, INCY specifies the increment for the  elements  of  Y.
79              INCY must not be zero.  Unchanged on exit.
80

FURTHER DETAILS

82       Level 2 Blas routine.
83       -- Written on 22-October-1986.
84          Jack Dongarra, Argonne National Lab.
85          Jeremy Du Croz, Nag Central Office.
86          Sven Hammarling, Nag Central Office.
87          Richard Hanson, Sandia National Labs.
88       -- Modified for the absolute-value product, April 2006
89          Jason Riedy, UC Berkeley
90
91
92
93    LAPACK routine (version 3.2) November 2008                    CLA_SYAMV(1)
Impressum