1DLAQGB(1) LAPACK auxiliary routine (version 3.1) DLAQGB(1)
2
3
4
6 DLAQGB - a general M by N band matrix A with KL subdiagonals and KU
7 superdiagonals using the row and scaling factors in the vectors R and C
8
10 SUBROUTINE DLAQGB( M, N, KL, KU, AB, LDAB, R, C, ROWCND, COLCND, AMAX,
11 EQUED )
12
13 CHARACTER EQUED
14
15 INTEGER KL, KU, LDAB, M, N
16
17 DOUBLE PRECISION AMAX, COLCND, ROWCND
18
19 DOUBLE PRECISION AB( LDAB, * ), C( * ), R( * )
20
22 DLAQGB equilibrates a general M by N band matrix A with KL subdiagonals
23 and KU superdiagonals using the row and scaling factors in the vectors
24 R and C.
25
26
28 M (input) INTEGER
29 The number of rows of the matrix A. M >= 0.
30
31 N (input) INTEGER
32 The number of columns of the matrix A. N >= 0.
33
34 KL (input) INTEGER
35 The number of subdiagonals within the band of A. KL >= 0.
36
37 KU (input) INTEGER
38 The number of superdiagonals within the band of A. KU >= 0.
39
40 AB (input/output) DOUBLE PRECISION array, dimension (LDAB,N)
41 On entry, the matrix A in band storage, in rows 1 to KL+KU+1.
42 The j-th column of A is stored in the j-th column of the array
43 AB as follows: AB(ku+1+i-j,j) = A(i,j) for max(1,j-
44 ku)<=i<=min(m,j+kl)
45
46 On exit, the equilibrated matrix, in the same storage format as
47 A. See EQUED for the form of the equilibrated matrix.
48
49 LDAB (input) INTEGER
50 The leading dimension of the array AB. LDA >= KL+KU+1.
51
52 R (input) DOUBLE PRECISION array, dimension (M)
53 The row scale factors for A.
54
55 C (input) DOUBLE PRECISION array, dimension (N)
56 The column scale factors for A.
57
58 ROWCND (input) DOUBLE PRECISION
59 Ratio of the smallest R(i) to the largest R(i).
60
61 COLCND (input) DOUBLE PRECISION
62 Ratio of the smallest C(i) to the largest C(i).
63
64 AMAX (input) DOUBLE PRECISION
65 Absolute value of largest matrix entry.
66
67 EQUED (output) CHARACTER*1
68 Specifies the form of equilibration that was done. = 'N': No
69 equilibration
70 = 'R': Row equilibration, i.e., A has been premultiplied by
71 diag(R). = 'C': Column equilibration, i.e., A has been post‐
72 multiplied by diag(C). = 'B': Both row and column equilibra‐
73 tion, i.e., A has been replaced by diag(R) * A * diag(C).
74
76 THRESH is a threshold value used to decide if row or column scaling
77 should be done based on the ratio of the row or column scaling factors.
78 If ROWCND < THRESH, row scaling is done, and if COLCND < THRESH, column
79 scaling is done.
80
81 LARGE and SMALL are threshold values used to decide if row scaling
82 should be done based on the absolute size of the largest matrix ele‐
83 ment. If AMAX > LARGE or AMAX < SMALL, row scaling is done.
84
85
86
87 LAPACK auxiliary routine (versionNo3v.e1m)ber 2006 DLAQGB(1)