In this paper we are concerned with the spectral analysis of the sequence of preconditioned matrices {P-1 n An(a, m1, m2, k)}n, where n = (n1, n2), N(n) = n1n2 and where An (a, m1, m2, k) ∈ ℝN(n) × N(n) is the symmetric two-level matrix coming from a high-order Finite Difference (FD) discretization of the problem (equation presented) with ν denoting the unit outward normal direction and where m1 and m2 are parameters identifying the precision order of the used FD schemes. We assume that the coefficient a(x, y) is nonnegative and that the set of the possible zeros can be represented by a finite collection of curves. The proposed preconditioning matrix sequences correspond to two different choices: the Toeplitz sequence {An(1, m1, m2, k)}n and a Toeplitz based sequence that adds to the Toeplitz structure the informative content given by the suitable scaled diagonal part of An(a, m1, m2 k). The former case gives rise to optimal preconditioning sequences under the assumption of positivity and boundedness of a. With respect to the latter, the main result is the proof of the asymptotic clustering at unity of the eigenvalues of the preconditioned matrices, where the "strength" of the cluster depends on the order k, on the regularity features of a(x, y) and on the presence of zeros of a(x, y).

Preconditioning strategies for 2D Finite Difference matrix sequences

SERRA CAPIZZANO, STEFANO;
2003

Abstract

In this paper we are concerned with the spectral analysis of the sequence of preconditioned matrices {P-1 n An(a, m1, m2, k)}n, where n = (n1, n2), N(n) = n1n2 and where An (a, m1, m2, k) ∈ ℝN(n) × N(n) is the symmetric two-level matrix coming from a high-order Finite Difference (FD) discretization of the problem (equation presented) with ν denoting the unit outward normal direction and where m1 and m2 are parameters identifying the precision order of the used FD schemes. We assume that the coefficient a(x, y) is nonnegative and that the set of the possible zeros can be represented by a finite collection of curves. The proposed preconditioning matrix sequences correspond to two different choices: the Toeplitz sequence {An(1, m1, m2, k)}n and a Toeplitz based sequence that adds to the Toeplitz structure the informative content given by the suitable scaled diagonal part of An(a, m1, m2 k). The former case gives rise to optimal preconditioning sequences under the assumption of positivity and boundedness of a. With respect to the latter, the main result is the proof of the asymptotic clustering at unity of the eigenvalues of the preconditioned matrices, where the "strength" of the cluster depends on the order k, on the regularity features of a(x, y) and on the presence of zeros of a(x, y).
Clustering and preconditioning; Finite differences; Spectral distribution; Toeplitz and Vandermonde matrices
SERRA CAPIZZANO, Stefano; Tablino Possio, C.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: http://hdl.handle.net/11383/4781
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 7
  • ???jsp.display-item.citation.isi??? 8
social impact