arxiv: v2 [math.na] 16 May 2018

Transcription

1 BLOCK DIAGONAL DOMINANCE OF MATRICES REVISITED: BOUNDS FOR THE NORMS OF INVERSES AND EIGENVALUE INCLUSION SETS CARLOS ECHEVERRÍA, JÖRG LIESEN, AND REINHARD NABBEN arxv:72.662v2 [math.na] 6 May 28 Dedcated to Rchard S. Varga Abstract. We generalze the bounds on the nverses of dagonally domnant matrces obtaned n [6] from scalar to block trdagonal matrces. Our dervatons are based on a generalzaton of the classcal condton of block dagonal domnance of matrces gven by Fengold and Varga n []. Based on ths generalzaton, whch was recently presented n [3], we also derve a varant of the Gershgorn Crcle Theorem for general block matrces whch can provde tghter spectral ncluson regons than those obtaned by Fengold and Varga. Key words. block matrces, block dagonal domnance, block trdagonal matrces, decay bounds for the nverse, egenvalue ncluson regons, Gershgorn Crcle Theorem AMS subject classfcatons. A4, 6F. Introducton. Matrces that are characterzed by off-dagonal decay, or more generally localzaton of ther entres, appear n applcatons throughout the mathematcal and computatonal scences. The presence of such localzaton can lead to computatonal savngs, snce t allows to (closely) approxmate a gven matrx by usng ts sgnfcant entres only, and dscardng the neglgble ones accordng to a pre-establshed crteron. In ths context t s then of great practcal nterest to know a pror how many and whch of these entres can be dscarded as nsgnfcant. Many authors have therefore studed decay rates for dfferent matrx classes and functons of matrces; see, e.g., [2, 4,, 6, 7,, 4, 8]. For an excellent survey of the current state-of-the-art we refer to []. An mportant example n ths context s gven by the (nonsymmetrc) dagonally domnant matrces, and n partcular the dagonally domnant trdagonal matrces, whch were studed, e.g., n [, 6]. As shown n these works, the entres of the nverse decay wth an exponental rate along a row or column, dependng on whether the gven matrx s row or column dagonally domnant; see [, Secton 3.2] for a more general treatment of decay bounds for the nverse and further references. Our man goal n ths paper s to generalze results of [6] from scalar to block trdagonal matrces. In order to do so, we use a generalzaton of the classcal defnton of block dagonal domnance of Fengold and Varga [] to derve bounds and decay rates for the block norms of the nverse of block trdagonal matrces. We also show how to mprove these bounds teratvely (Secton 2). Moreover, we obtan a new varant of the Gershgorn Crcle Theorem for general block matrces (Secton 3). Throughout ths paper we assume that s a gven submultplcatve matrx norm. Ths work dated May 6, 28. Insttute of Mathematcs, Techncal Unversty of Berln, Strasse des 7. Jun 36, D-623 Berln, Germany. Carlos Echeverra s work was partally supported by the Ensten Center for Mathematcs, Berln, the Deutsche Akademsche Austausch Denst (DAAD), Germany, and the Consejo Naconal de Cenca y Tecnología (CONACyT), Méxco. Contact e-mals: {echeverra, lesen, nabben}@math.tu-berln.de.

2 2 C. ECHEVERRÍA, J. LIESEN and R. NABBEN 2. Bounds on the nverses of block trdagonal matrces. We start wth a defnton of block dagonally domnant matrces, whch was recently presented n [3] n an applcaton of block dagonal precondtonng. Defnton 2.. Consder a matrx of the form (2.) A = [A j ] wth blocks A j C m m for,j =,...,n. The matrx A s called row block dagonally domnant (wth respect to the matrx norm ) when the dagonal blocks A are nonsngular, and (2.2) j= A A j, for =,...,n. If a strct nequalty holds n (2.2) then A s called row block strctly dagonally domnant (wth respect to the matrx norm ). Obvously, an analogous defnton of column block dagonal domnance s possble. Most of the results n ths paper can be easly rewrtten for that case. Also note that the authors of [3] call a matrx block dagonally domnant when all ts dagonal blocks are nonsngular, and (2.2) or the anologous condtons wth A j A replacng A A j hold (n the -norm). The above defnton of (row) block dagonal domnance generalzes the one of Fengold and Varga gven n [, Defnton ], who consdered a matrx as n (2.) block dagonally domnant when the dagonal blocks A are nonsngular, and (2.3) j= A j ( A ), for =,...,n. It s clear that f a matrx satsfes these condtons, then t also satsfes the condtons gven n Defnton 2.. Accordng to Varga[9, p. 6], the defnton of block dagonal domnance gven n [] s one of the earlest, and t was roughly smultaneously and ndependently consdered also by Ostrowsk [7] and Fedler and Pták [2]. Varga calls ths a Zetgest phenomenon. In the specalcasem =,.e., all the blocksa j areofsze and A j = A j, the nequaltes (2.2) and (2.3) are equvalent, and they can all be wrtten as A j A, for =,...,n, j= whch s the usual defnton of row dagonal domnance. In the rest of ths secton we wll restrct our attenton to block trdagonal matrces of the form (2.4) A = A B C A 2 B C n 2 A n C n B n A n, where A,B,C C m m.

3 BLOCK DIAGONAL DOMINACE AND ITS APPLICATIONS 3 Frst Capovan for the scalar case n [8, 9] and later Ikebe for the block case n [3] (see also [6]), have shown that the nverse of a nonsngular block trdagonal matrx can be descrbed by four sets of matrces. The man result can be stated as follows. Theorem 2.2. Let A be as n (2.4), and suppose that A as well as B and for =,...,n exst. If we wrte A = [Z j ] wth Z j C m m, then there C exst matrces U,V,X,Y C m m wth U V = X Y for =,...,n, and { U V j f j, (2.) Z j = Y X j f j. Moreover, the matrces U,V,X,Y, =,...,n, are recursvely gven by (2.8) (2.6) U = I, U 2 = B U, (2.7) U = B 2U 2 +A U ), for = 3,...,n, V n = (A n U n +C n U n ), V n = V n A n Bn (2.9) V = (V +A + +V +2C +)B, for = n 2,...,. X = I, X 2 = X A C (2.) X = (X 2B 2 +X A )C, for = 3,...,n, (2.2) Y n = (X n A n +X n B n ), Y n = Cn ny n, (2.3) Y = C (A + Y + +B + Y +2 ), for = n 2,...,. The next result s a generalzaton of [, Theorem 3.2]. Lemma 2.3. Let A be a matrx as n Theorem 2.2. Suppose n addton that A s row block dagonally domnant, and that (2.4) A B < and A n C n <. Then the sequence { U } n = s strctly ncreasng, and the sequence { Y } n = s strctly decreasng. Proof. Frst we consder the sequence { U } n =. The defnton of U 2 n (2.6) mples that U = A B U 2. Takng norms and usng the frst nequalty n (2.4) yelds U A B U 2 < U 2. Now suppose that U < U 2 < < U holds for some 3. The equaton for U n (2.7) can be wrtten as A B U = U +A C 2U 2. Rearrangng terms and takng norms we obtan U A C 2 U 2 + A B U < A C 2 U + A B U, wherewehaveused thenducton hypothess,.e., U 2 < U, n ordertoobtan the strct nequalty. Snce A s row block dagonally domnant we have A B + A C 2.

4 4 C. ECHEVERRÍA, J. LIESEN and R. NABBEN Combnng ths wth the prevous nequalty gves U U < A A C 2, B so that ndeed U < U. Next we consder the sequence { Y } n =. The defnton of Y n n (2.2) mples that Y n = An C n Y n. Takng norms and usng the second nequalty n (2.4) yelds Y n A n C n Y n < Y n. Now suppose that Y n < Y n < < Y + holds for some n 2. The equaton for Y n (2.3) can be wrtten as A + C Y = Y + +A + B +Y +2. Rearrangng terms and takng norms we obtan Y + A + C Y + A + B + Y +2 < A + C Y + A + B + Y +, where we have used the nducton hypothess,.e., Y +2 < Y +, n orderto obtan the strct nequalty. Snce A s row block dagonally domnant we have A + C + A + B +. Combnng ths wth the prevous nequalty gves Y + Y < A + C A + B + so that ndeed Y + < Y. Forthe restofths sectonwewll assumethat that AsamatrxasnLemma 2.3. Then the nverse s gven by A = [Z j ] wth Z j = Y X j for j; see Theorem 2.2. Thus, foreachfxed j =,...,n, thestrctdecreaseofthesequence{ Y } n = suggests that the sequence { Z j } n =j decreases as well,.e., that the norms of the blocks of A decay columnwse away from the dagonal. We wll now study ths decay n detal. We set C = B n =, and defne τ = A B A, for =,...,n, C ω = A C A, for =,...,n. B TherowblockdagonaldomnanceofAthenmplesthat τ and ω. Also note that, by assumpton, τ = A B <, ω n = An C n <, and τ n = ω =. In order to obtan bounds on the norms of the block entres A, we wll frst derve alternatve recurrence formulas for the matrces U and Y from Lemma 2.3. To ths end, we ntroduce some ntermedate quanttes and gve bounds on ther norms n the followng result.

5 BLOCK DIAGONAL DOMINACE AND ITS APPLICATIONS Lemma 2.4. The followng assertons hold: (a) The matrces L = T = A B, T 2 = I A 2 C T, and L = T A B, for = 2,...,n, T = I A C L, for = 3,...,n, are all nonsngular, and L τ, for =,...,n. (b) The matrces M n = W n = A n C n, W n = I A n B n W n, and M = W A C, for = n,...,2, W = I A B M +, for = n 2,...,2, are all nonsngular, and M ω, for = 2,...,n. Proof. We only prove (a); the proof of (b) s analogous. The matrces L = T = A B are nonsngular snce both A and B are. Moreover, (2.4) gves L = T = A B = τ <. Now suppose that L τ holds for some 2. Then A C L A C L <, where we have also used that A C A B <. Thus, T = I A C L s nonsngular, and therefore L = T A B s nonsngular. Usng the Neumann seres gves T = (I A C L ) = (A C L ) k A C L k = k= A C L A C L A C, and L = T A B A = τ A C, whch fnshes the proof. Usng Lemma 2.4 we can now derve alternatve recurrences for the matrces U and Y from Lemma 2.3. Lemma 2.. If A s a matrx as n Lemma 2.3, then the correspondng matrces U and Y are gven by B k= (2.) (2.6) U = L U +, for =,...n, Y = M Y, for = n,...2, where the matrces L and M are defned as n Lemma 2.4. Proof. We only prove that (2.) holds; the proof of (2.6) s analogous. From (2.6) and the defnton of T n Lemma 2.3 we obtan We next wrte (2.7) for = 3 as U = A B U 2 = L U 2. A 2 B 2U 3 = A 2 C U +U 2 = A 2 C T U 2 +U 2 = (I A 2 C T )U 2 = T 2 U 2, and hence U 2 = T 2 A 2 B 2U 3 = L 2 U 3.

6 6 C. ECHEVERRÍA, J. LIESEN and R. NABBEN Now suppose that U = L U holds for some 3 n. Then from (2.7) we obtan and hence A B U + = A C U +U = (I A C L )U = T U, U = T A B U + = L U +, whch completes the proof. We are now ready to state and prove our bounds on the norms of the blocks of A, whch generalze [6, Theorems 3. and 3.2] from the scalar to the block case. Theorem 2.6. If A s a matrx as n Lemma 2.3, then A = [Z j ] wth (2.7) (2.8) j Z j Z jj τ k, for all < j, Z j Z jj k= k=j+ ω k, for all > j. Moreover, for =,...,n, (2.9) I A +τ C +ω + B Z I A τ C ω + B, provded that the denomnator of the upper bound s larger than zero, and where we set C = B n =, and τ = ω n+ =. Proof. From Lemma 2. we know that U = L U + holds for =,...,n. Thus, for all < j, ( j ) ( j Z j = U V j = L U + V j = ( ) j L k U j V j = ( ) j L k )Z jj. Takng norms and usng Lemma 2.4 yelds k= j j Z j Z jj L k Z jj τ k. The expresson for > j follows analogously usng the two lemmas. Snce AA = AZ = I we have k= k= C Z, +A Z +B Z +, = I, for =,...,n, where we set C = Z, = B n = Z n+,n =. Usng (2.) and Lemma 2., Z, = U V = L U V = L Z, Z +, = Y + X = M + Y X = M + Z, where we set U = L = Y n+ = M n+ =. Combnng ths wth the prevous equaton yelds (2.2) C L Z B M + Z +A Z = I, for =,...,n, k=

7 BLOCK DIAGONAL DOMINACE AND ITS APPLICATIONS 7 Takng norms and usng agan Lemma 2.4 now gves I = C L Z B M + Z +A Z ( C L + A + B M + ) Z (τ C + A +ω + B ) Z, for =,...,n, where we set τ = ω n+ =, and whch shows the lower bound n (2.9). In order to show the upper bound we wrte (2.2) as Ths yelds I A Z = (C L Z +B M + Z ), for =,...,n. A Z I I A Z = C L Z +B M + Z (τ C +ω + B ) Z. From Z = A A Z A A Z we get A Z Z / A, and combnng ths wth the prevous nequalty yelds (τ C +ω + B ) Z A Z I. When A τ C ω + B > holds,wegettheupperboundn(2.9). Note that the postvty assumpton on the denomnator of the upper bound n (2.9) s ndeed necessary. A smple example for whch the denomnator s equal to zero s gven by the matrx A = trdag(,2, ) R n n wth blocks, whch satsfes all assumptons of Lemma 2.3. Both the off-dagonal bounds(2.7) (2.8) and the dagonal bounds(2.9) depend on the values τ and ω, whch bound L and M, respectvely. We wll now show that by modfyng the proof of Lemma 2.4 the bounds can be mproved n an teratve fashon. Ths s analogous to the teratve mprovement for the case when the blocks of A are scalars, whch was consdered n [6]. We have shown n the nductve proof of Lemma 2.4 that T A C L A C L A C. Ths bound can be mproved by makng use of Lemma 2.4 tself,.e., T and ths yelds A C L A C L A, C τ L = T A B A B A C τ. If we denote the expresson on the rght hand sde by τ,2, then we obtan a modfed verson of Lemma 2.4, where L τ,2 τ 2. Iteratvely we now defne, for all =,...,n and t =,...,n, τ f t =, τ,t = τ,t f t >, A B A C τ,t else.

8 8 C. ECHEVERRÍA, J. LIESEN and R. NABBEN Analogouslywecanproceedforthevalues M,andherewedefne, forall =,...,n and t =,...,n, ω f t =, ω,t = ω,t f n t+ <, A C else. A B ω +,t Usng these defntons we can easly prove the followng modfed verson of Theorem 2.6, whch refnes the bounds (2.7), (2.8) and (2.9) as t ncreases, and whch generalzes [6, Theorems 3.4 and 3.] from the scalar to the block case. Theorem 2.7. If A s a matrx as n Lemma 2.3 wth A = [Z j ], then for each t =,...,n, (2.2) (2.22) j Z j Z jj τ k,t, for all < j, Z j Z jj k= k=j+ ω k,t, for all > j. Moreover, for =,...,n, (2.23) I A +τ,t C +ω +,t B Z I A τ,t C ω +,t B, provded that the denomnator of the upper bound s larger than zero, and where we set C = B n =, and τ,t = ω n+,t =. Note that the statements of Theorem 2.7 wth t = are the same as those n Theorem 2.6. By constructon, the sequences {τ,t } t= n and {ω,t} t= n are decreasng, and hence the bounds (2.2), (2.22) and (2.23) become tghter as t ncreases. However, snce we have used the submultplcatvty property of the matrx norm n the dervaton, t s not guaranteed that the bounds n Theorem 2.7 wth t = n wll gve the exact norms of the blocks of A. Ths s a dfference to the scalar case, where n the last refnement step one obtans the exact nverse; see [6]. Fnally, let us defne ρ,t := max τ,t, ω 2,t := maxω,t, for t =,...,n. Then the off-dagonal bounds (2.2) and (2.22) of Theorem 2.7 mmedately gve the followng result about the decay of the norms Z j ; cf. [6, Corollary 3.7] Corollary 2.8. If A s a matrx as n Theorem 2.7, then and for each t =,...,n. Z j ρ j,t Z jj, for all < j, Z j ρ j 2,t Z jj, for all > j,

9 BLOCK DIAGONAL DOMINACE AND ITS APPLICATIONS 9 In the followng we provde some numercal llustratons of the bounds n Theorem 2.7 for dfferent values of t. We consder dfferent matrces A = [A j ] whch are rowblock dagonallydomnant, and we compute the correspondngmatrcesz = [Z j ] usng the recurrences stated n Theorem 2.2. In all experments we use the matrx 2-norm, 2. For each gven par,j, we denote by u j the value of a computed upper bound (.e., (2.2), (2.22) or (2.23)) on the value Z j 2 and for each we denote by l the value of the computed lower bound for the correspondng dagonal entry (.e., (2.23)). Then the relatve errors n the upper and lower bounds are gven by (2.24) E u j = u j Z j 2 u j and E l = Z 2 l Z 2, respectvely. (Thus, both Ej u and El are between and.) Example 2.. We start wth the symmetrc block Toepltz matrx (2.2) A = T I +I T R 8 8, where T = trdag(,2, ) R 9 9,.e., A s of the form (2.4) wth A = trdag(,4, ), and B = C = dag( ) for all. We have κ 2 (A ) = ,.e., the matrx A s qute well condtoned. For the computed matrx Z = [Z j ] we obtan ZA I 2 = , suggestng that Z s a reasonably accurate approxmaton of the exact nverse A. In the top row of Fgure 2. we show the relatve errors Ej u for the refnement step t = (no refnement) and t = 8 (maxmal refnement). We observe that the upper bounds are qute tght already for t =, and that for t = 8 the maxmal relatve error s on the order 3,.e., the value of the upper bound s almost exact. In the bottom row of Fgure 2. we show the values Z 2 for =,...,9, and the correspondngupper and lowerbounds (2.23) for the refnement steps t = and t = 8. We observe that whle the upper bounds on Z 2 for t = 8 almost exactly match the exact values, the lower bounds do not mprove by the teratve refnement. The maxmal error of the lower bounds for the dagonal block entres of Z n the maxmal refnement step s on the order. The maxmal relatve errors n the upper and lower bounds and all refnement steps are shown n the followng table: t max j E u j max E l

10 C. ECHEVERRÍA, J. LIESEN and R. NABBEN refnement step t= refnement step t=8.8 6 E u j.6.4 E j block ndex j block ndex block ndex j block ndex refnement step t= upper bound (2.23) norm(z ) lower bound (2.23) refnement step t=8 upper bound (2.23) norm(z ) lower bound (2.23) block ndex block ndex Fg. 2.. Relatve errors E u j (top row), upper and lower bounds on Z 2 (bottom row) for the matrx A of Example 2.. Example 2.2. Let A be the nonsymmetrc block Toepltz matrx of the form (2.2) wth T = trdag(,29.999, ) R 9 9,.e., A agan takes the form (2.4) wth A = trdag(,49.999, ), B = dag( ), and C = dag( ). The condton number n ths case s κ 2 (A) = 7.72, and for the computed matrx Z we obtan ZA I 2 =.. The top row of Fgure 2.2 shows the relatve errors for the refnement steps t = and t = 8. We observe that for ths nonsymmetrc example the upper bounds are not as accurate as those gven n the symmetrc case, producng a maxmal relatve error at refnement step t = 8 on the order 3. The bottom row of Fgure 2.2 shows the upper and lower bounds (2.23) as well as the values Z 2 for =,...,9, and refnement steps t = and t = 8. Agan we can observe that whle we obtan a reasonable approxmaton n the upper bounds on Z 2 for t = 8, the lower bounds almost do not mprove by the teratve refnement process. The maxmal relatve errors n the upper and lowerbounds and all refnement steps s shown n the followngtable: t max j Ej u max E l

11 BLOCK DIAGONAL DOMINACE AND ITS APPLICATIONS refnement step t= refnement step t= E u j.4 E u j block ndex j block ndex block ndex j block ndex - refnement step t= upper bound (2.23) norm(z ) lower bound (2.23) - refnement step t=8 upper bound (2.23) norm(z ) lower bound (2.23) block ndex block ndex Fg Relatve errors E u j (top row), and upper and lower bounds on Z 2 (bottom row) for the matrx A of Example 2.2. Example 2.3. We now consder the nonsymmetrc block trdagonal matrx A = (R I)(T I +I T) R 8 8, where T s gven as n Example 2., and R R 9 9 s a random dagonal matrx wth nonzero nteger entres between and and constructed n MATLAB wth the command R = dag(cel(*rand(9,))). Thus, A s of the form (2.4) wth random trdagonal Toepltz matrces A, and random constant dagonal matrces B and C for all. For ths matrx we have κ 2 (A) = , and the computed matrx Z yelds ZA I 2 =.9 9. The relatve errors n the bounds are shown n Fgure 2.3 and n the followng table: t max j E u j max E l

12 2 C. ECHEVERRÍA, J. LIESEN and R. NABBEN refnement step t= -2 2 refnement step t=8.8. E u j E u j. block ndex j block ndex block ndex j block ndex refnement step t= upper bound (2.23) norm(z ) lower bound (2.23) refnement step t=8 upper bound (2.23) norm(z ) lower bound (2.23) block ndex block ndex Fg Relatve errors E u j (top row), and upper and lower bounds on Z 2 (bottom row) for the matrx A of Example 2.3. Example 2.4. Fnally, we consder the nonsymmetrc block trdagonal matrx A = (R I) trdag(trdag(., 2, ), trdag( 2,, 2), trdag(., 2, )), wth A R 8 8, and where R R 9 9 s a random dagonal matrx constructed as n Example 2.3. In ths case A takes the form (2.4) wth A, B and C random trdagonal Toepltz matrces wth nteger entres for all. For ths matrx we have κ 2 (A) = , and ZA I 2 = The relatve errors n the bounds are shown n Fgure 2.4 and the followng table: t max j Ej u max E l

13 BLOCK DIAGONAL DOMINACE AND ITS APPLICATIONS 3 refnement step t= refnement step t= E u j.6.4 E u j.2.2. block ndex j block ndex block ndex j block ndex refnement step t= upper bound (2.23) norm(z ) lower bound (2.23) refnement step t=8 upper bound (2.23) norm(z ) lower bound (2.23) block ndex block ndex Fg Relatve errors E u j (top row), and upper and lower bounds on Z 2 (bottom row) for the matrx A of Example Incluson regons for egenvalues. In ths secton we generalze a result of Fengold and Varga on egenvalue ncluson regons of block matrces. We start wth the followng generalzaton of [, Theorem ]; also cf. [9, Theorem 6.2]. Lemma 3.. If a matrx A as n (2.) s row block strctly dagonally domnant, then A s nonsngular. Proof. The proof closely follows the proof of [, Theorem ]. Suppose that A s row block strctly dagonally domnant but sngular. Then there exsts a nonzero block vector X, parttoned conformally wth respect to the partton of A n (2.), such that Ths s equvalent to A X + A X. X n =. A j X j =, n, j=

14 4 C. ECHEVERRÍA, J. LIESEN and R. NABBEN and, snce the dagonal blocks A are nonsngular, X = j= A A jx j, n, Wthout loss of generalty we can assume that X s normalzed such that X for all n, wth equalty for some = r. For ths ndex we obtan = X r = j= A rr A rj X j n j= j r A rr A rj X j j= j r A rr A rj, whch contradcts the assumpton that A s row block strctly dagonally domnant. Thus, A must be nonsngular. If λ s an egenvalue of A, then A λi s sngular, and hence A λi cannot be block strctly dagonally domnant. Ths mmedately gves the followng result, whch generalzes [, Theorem 2]; also cf. [9, Theorem 6.3]. Corollary 3.2. If a matrx A s as n (2.), and λ s an egenvalue of A, then there exsts at least one {,...,n} wth (3.) (A λi) A j. j= If all the blocks A j of A are of sze, and A j = A j, then ths result reduces to the classcal Gershgorn Crcle Theorem. Corollary 3.2 shows that each egenvalue λ of A must be contaned n the unon of the sets { G new = z C : j= } (A zi) A j, for =,...,n. Due to the submultplcatvty property of the matrx norm, the sets G new are potentally smaller than the ones proposed n [, Defnton 3], { G FV = z C : j= } (A zi) A j,.e., we have G new G FV. We wll llustrate ths fact wth numercal examples. Example 3.. We frst consder the symmetrc matrx A = = [ A A 2 A 2 A 22 whch has the egenvalues.486, 2.382, 4.68, and 7.44 (computed n MATLAB and rounded to fve sgnfcant dgts). The left part of Fgure 3. shows the boundares of the correspondng sets G new j and G FV j for j =,2,.e., the curves for z C where ],

15 BLOCK DIAGONAL DOMINACE AND ITS APPLICATIONS 4 3 egenvalues (A -zi) - A j = (A -zi) - A j = 4 3 egenvalues (A -zi) - A j = (A -zi) - A j = 2 2 Im(z) Im(z) Re(z) Re(z) Fg. 3.. Egenvalue ncluson regons obtaned from the sets G new and G FV n Example 3.. (A zi) A j =, and (A zi) A j =,,j {,2}, j, respectvely. Clearly, the sets G new gve tghter ncluson regons for the egenvalues than the sets G FV as well as the usual Gershgorn crcles for the matrx A, whch are gven by the two crcles centered at z = 4 of radus 3 and 4. We next consder the nonsymmetrc matrx (3.2) A = = [ A A 2 A 2 A 22 whch has the egenvalues.68, 2.99, , and As shown n the rght part of Fgure 3., the sets G new agan gve tghter ncluson regons than the sets as well as the usual Gershgorn crcles. G FV Acknowledgments. Carlos Echeverría s work was partally supported by the Ensten Center for Mathematcs, Berln, the Deutscher Akademscher Austauschdenst (DAAD), Germany, and the Consejo Naconal de Cenca y Tecnología (CONA- CyT), Méxco. We thank Mchele Benz and an anonymous referee for ther helpful comments. ], REFERENCES [] M. Benz, Localzaton n matrx computatons: Theory and applcatons, n Explotng Hdden Structure n Matrx Computatons: Algorthms and Applcatons, vol. 273 of Lecture Notes n Math., Sprnger, Cham, 26, pp [2] M. Benz and P. Boto, Decay propertes for functons of matrces over C -algebras, Lnear Algebra Appl., 46 (24), pp [3] M. Benz, T. M. Evans, S. P. Hamlton, M. Lupo Pasn, and S. R. Slattery, Analyss of Monte Carlo accelerated teratve methods for sparse lnear systems, Numer. Lnear Algebra Appl., 24 (27), pp. e288, 8. [4] M. Benz and G. H. Golub, Bounds for the entres of matrx functons wth applcatons to precondtonng, BIT, 39 (999), pp

16 6 C. ECHEVERRÍA, J. LIESEN and R. NABBEN [] M. Benz and N. Razouk, Decay bounds and O(n) algorthms for approxmatng functons of sparse matrces, Electron. Trans. Numer. Anal., 28 (27/8), pp [6] M. Benz and V. Smoncn, Decay bounds for functons of Hermtan matrces wth banded or Kronecker structure, SIAM J. Matrx Anal. Appl., 36 (2), pp [7] C. Canuto, V. Smoncn, and M. Veran, On the decay of the nverse of matrces that are sum of Kronecker products, Lnear Algebra Appl., 42 (24), pp [8] M. Capovan, Sulla determnazone della nversa delle matrc trdagonal e trdagonal a blocch, Calcolo, 7 (97), pp [9] M. Capovan, Su alcune propretà delle matrc trdagonal e pentadagonal, Calcolo, 8(97), pp [] S. Demko, W. F. Moss, and P. W. Smth, Decay rates for nverses of band matrces, Math. Comp., 43 (984), pp [] D. G. Fengold and R. S. Varga, Block dagonally domnant matrces and generalzatons of the Gerschgorn crcle theorem, Pacfc J. Math., 2 (962), pp [2] M. Fedler and V. Pták, Generalzed norms of matrces and the locaton of the spectrum, Czechoslovak Math. J., 2 (87) (962), pp [3] Y. Ikebe, On nverses of Hessenberg matrces, Lnear Algebra Appl., 24 (979), pp [4] I. Krshtal, T. Strohmer, and T. Wertz, Localzaton of matrx factorzatons, Found. Comput. Math., (2), pp [] R. Nabben, Decay rates of the nverse of nonsymmetrc trdagonal and band matrces, SIAM J. Matrx Anal. Appl., 2 (999), pp [6], Two-sded bounds on the nverses of dagonally domnant trdagonal matrces, Lnear Algebra Appl., 287 (999), pp Specal ssue celebratng the 6th brthday of Ludwg Elsner. [7] A. M. Ostrowsk, On some metrcal propertes of operator matrces and matrces parttoned nto blocks, J. Math. Anal. Appl., 2 (96), pp [8] R. Peluso and T. Polt, Some mprovements for two-sded bounds on the nverse of dagonally domnant trdagonal matrces, Lnear Algebra Appl., 33 (2), pp. 4. [9] R. S. Varga, Geršgorn and hs crcles, vol. 36 of Sprnger Seres n Computatonal Mathematcs, Sprnger-Verlag, Berln, 24.