Preconditioning strategies for 2D Finite Difference matrix sequences

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).