LAPACK  3.4.2
LAPACK: Linear Algebra PACKage
 All Files Functions Groups
zgbt02.f
Go to the documentation of this file.
1 *> \brief \b ZGBT02
2 *
3 * =========== DOCUMENTATION ===========
4 *
5 * Online html documentation available at
6 * http://www.netlib.org/lapack/explore-html/
7 *
8 * Definition:
9 * ===========
10 *
11 * SUBROUTINE ZGBT02( TRANS, M, N, KL, KU, NRHS, A, LDA, X, LDX, B,
12 * LDB, RESID )
13 *
14 * .. Scalar Arguments ..
15 * CHARACTER TRANS
16 * INTEGER KL, KU, LDA, LDB, LDX, M, N, NRHS
17 * DOUBLE PRECISION RESID
18 * ..
19 * .. Array Arguments ..
20 * COMPLEX*16 A( LDA, * ), B( LDB, * ), X( LDX, * )
21 * ..
22 *
23 *
24 *> \par Purpose:
25 * =============
26 *>
27 *> \verbatim
28 *>
29 *> ZGBT02 computes the residual for a solution of a banded system of
30 *> equations A*x = b or A'*x = b:
31 *> RESID = norm( B - A*X ) / ( norm(A) * norm(X) * EPS).
32 *> where EPS is the machine precision.
33 *> \endverbatim
34 *
35 * Arguments:
36 * ==========
37 *
38 *> \param[in] TRANS
39 *> \verbatim
40 *> TRANS is CHARACTER*1
41 *> Specifies the form of the system of equations:
42 *> = 'N': A *x = b
43 *> = 'T': A'*x = b, where A' is the transpose of A
44 *> = 'C': A'*x = b, where A' is the transpose of A
45 *> \endverbatim
46 *>
47 *> \param[in] M
48 *> \verbatim
49 *> M is INTEGER
50 *> The number of rows of the matrix A. M >= 0.
51 *> \endverbatim
52 *>
53 *> \param[in] N
54 *> \verbatim
55 *> N is INTEGER
56 *> The number of columns of the matrix A. N >= 0.
57 *> \endverbatim
58 *>
59 *> \param[in] KL
60 *> \verbatim
61 *> KL is INTEGER
62 *> The number of subdiagonals within the band of A. KL >= 0.
63 *> \endverbatim
64 *>
65 *> \param[in] KU
66 *> \verbatim
67 *> KU is INTEGER
68 *> The number of superdiagonals within the band of A. KU >= 0.
69 *> \endverbatim
70 *>
71 *> \param[in] NRHS
72 *> \verbatim
73 *> NRHS is INTEGER
74 *> The number of columns of B. NRHS >= 0.
75 *> \endverbatim
76 *>
77 *> \param[in] A
78 *> \verbatim
79 *> A is COMPLEX*16 array, dimension (LDA,N)
80 *> The original matrix A in band storage, stored in rows 1 to
81 *> KL+KU+1.
82 *> \endverbatim
83 *>
84 *> \param[in] LDA
85 *> \verbatim
86 *> LDA is INTEGER
87 *> The leading dimension of the array A. LDA >= max(1,KL+KU+1).
88 *> \endverbatim
89 *>
90 *> \param[in] X
91 *> \verbatim
92 *> X is COMPLEX*16 array, dimension (LDX,NRHS)
93 *> The computed solution vectors for the system of linear
94 *> equations.
95 *> \endverbatim
96 *>
97 *> \param[in] LDX
98 *> \verbatim
99 *> LDX is INTEGER
100 *> The leading dimension of the array X. If TRANS = 'N',
101 *> LDX >= max(1,N); if TRANS = 'T' or 'C', LDX >= max(1,M).
102 *> \endverbatim
103 *>
104 *> \param[in,out] B
105 *> \verbatim
106 *> B is COMPLEX*16 array, dimension (LDB,NRHS)
107 *> On entry, the right hand side vectors for the system of
108 *> linear equations.
109 *> On exit, B is overwritten with the difference B - A*X.
110 *> \endverbatim
111 *>
112 *> \param[in] LDB
113 *> \verbatim
114 *> LDB is INTEGER
115 *> The leading dimension of the array B. IF TRANS = 'N',
116 *> LDB >= max(1,M); if TRANS = 'T' or 'C', LDB >= max(1,N).
117 *> \endverbatim
118 *>
119 *> \param[out] RESID
120 *> \verbatim
121 *> RESID is DOUBLE PRECISION
122 *> The maximum over the number of right hand sides of
123 *> norm(B - A*X) / ( norm(A) * norm(X) * EPS ).
124 *> \endverbatim
125 *
126 * Authors:
127 * ========
128 *
129 *> \author Univ. of Tennessee
130 *> \author Univ. of California Berkeley
131 *> \author Univ. of Colorado Denver
132 *> \author NAG Ltd.
133 *
134 *> \date November 2011
135 *
136 *> \ingroup complex16_lin
137 *
138 * =====================================================================
139  SUBROUTINE zgbt02( TRANS, M, N, KL, KU, NRHS, A, LDA, X, LDX, B,
140  $ ldb, resid )
141 *
142 * -- LAPACK test routine (version 3.4.0) --
143 * -- LAPACK is a software package provided by Univ. of Tennessee, --
144 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
145 * November 2011
146 *
147 * .. Scalar Arguments ..
148  CHARACTER trans
149  INTEGER kl, ku, lda, ldb, ldx, m, n, nrhs
150  DOUBLE PRECISION resid
151 * ..
152 * .. Array Arguments ..
153  COMPLEX*16 a( lda, * ), b( ldb, * ), x( ldx, * )
154 * ..
155 *
156 * =====================================================================
157 *
158 * .. Parameters ..
159  DOUBLE PRECISION zero, one
160  parameter( zero = 0.0d+0, one = 1.0d+0 )
161  COMPLEX*16 cone
162  parameter( cone = ( 1.0d+0, 0.0d+0 ) )
163 * ..
164 * .. Local Scalars ..
165  INTEGER i1, i2, j, kd, n1
166  DOUBLE PRECISION anorm, bnorm, eps, xnorm
167 * ..
168 * .. External Functions ..
169  LOGICAL lsame
170  DOUBLE PRECISION dlamch, dzasum
171  EXTERNAL lsame, dlamch, dzasum
172 * ..
173 * .. External Subroutines ..
174  EXTERNAL zgbmv
175 * ..
176 * .. Intrinsic Functions ..
177  INTRINSIC max, min
178 * ..
179 * .. Executable Statements ..
180 *
181 * Quick return if N = 0 pr NRHS = 0
182 *
183  IF( m.LE.0 .OR. n.LE.0 .OR. nrhs.LE.0 ) THEN
184  resid = zero
185  return
186  END IF
187 *
188 * Exit with RESID = 1/EPS if ANORM = 0.
189 *
190  eps = dlamch( 'Epsilon' )
191  kd = ku + 1
192  anorm = zero
193  DO 10 j = 1, n
194  i1 = max( kd+1-j, 1 )
195  i2 = min( kd+m-j, kl+kd )
196  anorm = max( anorm, dzasum( i2-i1+1, a( i1, j ), 1 ) )
197  10 continue
198  IF( anorm.LE.zero ) THEN
199  resid = one / eps
200  return
201  END IF
202 *
203  IF( lsame( trans, 'T' ) .OR. lsame( trans, 'C' ) ) THEN
204  n1 = n
205  ELSE
206  n1 = m
207  END IF
208 *
209 * Compute B - A*X (or B - A'*X )
210 *
211  DO 20 j = 1, nrhs
212  CALL zgbmv( trans, m, n, kl, ku, -cone, a, lda, x( 1, j ), 1,
213  $ cone, b( 1, j ), 1 )
214  20 continue
215 *
216 * Compute the maximum over the number of right hand sides of
217 * norm(B - A*X) / ( norm(A) * norm(X) * EPS ).
218 *
219  resid = zero
220  DO 30 j = 1, nrhs
221  bnorm = dzasum( n1, b( 1, j ), 1 )
222  xnorm = dzasum( n1, x( 1, j ), 1 )
223  IF( xnorm.LE.zero ) THEN
224  resid = one / eps
225  ELSE
226  resid = max( resid, ( ( bnorm / anorm ) / xnorm ) / eps )
227  END IF
228  30 continue
229 *
230  return
231 *
232 * End of ZGBT02
233 *
234  END