1SLAQSP(1)           LAPACK auxiliary routine (version 3.1)           SLAQSP(1)
2
3
4

NAME

6       SLAQSP - a symmetric matrix A using the scaling factors in the vector S
7

SYNOPSIS

9       SUBROUTINE SLAQSP( UPLO, N, AP, S, SCOND, AMAX, EQUED )
10
11           CHARACTER      EQUED, UPLO
12
13           INTEGER        N
14
15           REAL           AMAX, SCOND
16
17           REAL           AP( * ), S( * )
18

PURPOSE

20       SLAQSP  equilibrates  a symmetric matrix A using the scaling factors in
21       the vector S.
22
23

ARGUMENTS

25       UPLO    (input) CHARACTER*1
26               Specifies whether the upper or lower  triangular  part  of  the
27               symmetric matrix A is stored.  = 'U':  Upper triangular
28               = 'L':  Lower triangular
29
30       N       (input) INTEGER
31               The order of the matrix A.  N >= 0.
32
33       AP      (input/output) REAL array, dimension (N*(N+1)/2)
34               On  entry,  the upper or lower triangle of the symmetric matrix
35               A, packed columnwise in a linear array.  The j-th column  of  A
36               is  stored  in  the  array AP as follows: if UPLO = 'U', AP(i +
37               (j-1)*j/2) =  A(i,j)  for  1<=i<=j;  if  UPLO  =  'L',  AP(i  +
38               (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.
39
40               On  exit,  the  equilibrated matrix:  diag(S) * A * diag(S), in
41               the same storage format as A.
42
43       S       (input) REAL array, dimension (N)
44               The scale factors for A.
45
46       SCOND   (input) REAL
47               Ratio of the smallest S(i) to the largest S(i).
48
49       AMAX    (input) REAL
50               Absolute value of largest matrix entry.
51
52       EQUED   (output) CHARACTER*1
53               Specifies whether or not equilibration was done.   =  'N':   No
54               equilibration.
55               =  'Y':   Equilibration  was done, i.e., A has been replaced by
56               diag(S) * A * diag(S).
57

PARAMETERS

59       THRESH is a threshold value used to decide if scaling  should  be  done
60       based  on the ratio of the scaling factors.  If SCOND < THRESH, scaling
61       is done.
62
63       LARGE and SMALL are threshold values used to decide if  scaling  should
64       be  done  based on the absolute size of the largest matrix element.  If
65       AMAX > LARGE or AMAX < SMALL, scaling is done.
66
67
68
69 LAPACK auxiliary routine (versionNo3v.e1m)ber 2006                       SLAQSP(1)
Impressum