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

NAME

6       ZPBEQU  -  computes  row  and column scalings intended to equilibrate a
7       Hermitian positive definite band matrix A and reduce its condition num‐
8       ber (with respect to the two-norm)
9

SYNOPSIS

11       SUBROUTINE ZPBEQU( UPLO, N, KD, AB, LDAB, S, SCOND, AMAX, INFO )
12
13           CHARACTER      UPLO
14
15           INTEGER        INFO, KD, LDAB, N
16
17           DOUBLE         PRECISION AMAX, SCOND
18
19           DOUBLE         PRECISION S( * )
20
21           COMPLEX*16     AB( LDAB, * )
22

PURPOSE

24       ZPBEQU  computes row and column scalings intended to equilibrate a Her‐
25       mitian positive definite band matrix A and reduce its condition  number
26       (with  respect  to the two-norm).  S contains the scale factors, S(i) =
27       1/sqrt(A(i,i)), chosen so that the scaled matrix B with elements B(i,j)
28       = S(i)*A(i,j)*S(j) has ones on the diagonal.  This choice of S puts the
29       condition number of B within a factor N of the smallest possible condi‐
30       tion number over all possible diagonal scalings.
31

ARGUMENTS

33       UPLO    (input) CHARACTER*1
34               = 'U':  Upper triangular of A is stored;
35               = 'L':  Lower triangular of A is stored.
36
37       N       (input) INTEGER
38               The order of the matrix A.  N >= 0.
39
40       KD      (input) INTEGER
41               The  number of superdiagonals of the matrix A if UPLO = 'U', or
42               the number of subdiagonals if UPLO = 'L'.  KD >= 0.
43
44       AB      (input) COMPLEX*16 array, dimension (LDAB,N)
45               The upper or lower triangle of the  Hermitian  band  matrix  A,
46               stored in the first KD+1 rows of the array.  The j-th column of
47               A is stored in the j-th column of the array AB as  follows:  if
48               UPLO  =  'U', AB(kd+1+i-j,j) = A(i,j) for max(1,j-kd)<=i<=j; if
49               UPLO = 'L', AB(1+i-j,j)    = A(i,j) for j<=i<=min(n,j+kd).
50
51       LDAB     (input) INTEGER
52                The leading dimension of the array A.  LDAB >= KD+1.
53
54       S       (output) DOUBLE PRECISION array, dimension (N)
55               If INFO = 0, S contains the scale factors for A.
56
57       SCOND   (output) DOUBLE PRECISION
58               If INFO = 0, S contains the ratio of the smallest S(i)  to  the
59               largest  S(i).   If  SCOND >= 0.1 and AMAX is neither too large
60               nor too small, it is not worth scaling by S.
61
62       AMAX    (output) DOUBLE PRECISION
63               Absolute value of largest matrix  element.   If  AMAX  is  very
64               close to overflow or very close to underflow, the matrix should
65               be scaled.
66
67       INFO    (output) INTEGER
68               = 0:  successful exit
69               < 0:  if INFO = -i, the i-th argument had an illegal value.
70               > 0:  if INFO = i, the i-th diagonal element is nonpositive.
71
72
73
74 LAPACK routine (version 3.2)    November 2008                       ZPBEQU(1)
Impressum