Source author record

Debora Sesana

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

3works
2topics
4close collaborators

Published work

3 publications

preprint2014arXiv

Spectral behavior of preconditioned non-Hermitian multilevel block Toeplitz matrices with matrix-valued symbol

This note is devoted to preconditioning strategies for non-Hermitian multilevel block Toeplitz linear systems associated with a multivariate Lebesgue integrable matrix-valued symbol. In particular, we consider special preconditioned matrices, where the preconditioner has a band multilevel block Toeplitz structure, and we complement known results on the localization of the spectrum with global distribution results for the eigenvalues of the preconditioned matrices. In this respect, our main result is as follows. Let Ik:=(−π,π)kI_k:=(-π,π)^k, let Ms\mathcal M_s be the linear space of complex s×ss\times s matrices, and let f,g:Ik→Msf,g:I_k\to\mathcal M_s be functions whose components fij, gij:Ik→C, i,j=1,…,s,f_{ij},\,g_{ij}:I_k\to\mathbb C,\ i,j=1,\ldots,s, belong to L∞L^\infty. Consider the matrices Tn−1(g)Tn(f)T_n^{-1}(g)T_n(f), where n:=(n1,…,nk)n:=(n_1,\ldots,n_k) varies in Nk\mathbb N^k and Tn(f),Tn(g)T_n(f),T_n(g) are the multilevel block Toeplitz matrices of size n1⋯nksn_1\cdots n_ks generated by f,gf,g. Then {Tn−1(g)Tn(f)}n∈Nk∼λg−1f\{T_n^{-1}(g)T_n(f)\}_{n\in\mathbb N^k}\sim_λg^{-1}f, i.e. the family of matrices {Tn−1(g)Tn(f)}n∈Nk\{T_n^{-1}(g)T_n(f)\}_{n\in\mathbb N^k} has a global (asymptotic) spectral distribution described by the function g−1fg^{-1}f, provided gg possesses certain properties (which ensure in particular the invertibility of Tn−1(g)T_n^{-1}(g) for all nn) and the following topological conditions are met: the essential range of g−1fg^{-1}f, defined as the union of the essential ranges of the eigenvalue functions λj(g−1f), j=1,…,sλ_j(g^{-1}f),\ j=1,\ldots,s, does not disconnect the complex plane and has empty interior. This result generalizes the one obtained by Donatelli, Neytcheva, Serra-Capizzano in a previous work, concerning the non-preconditioned case g=1g=1. The last part of this note is devoted to numerical experiments, which confirm the theoretical analysis and suggest the choice of optimal GMRES preconditioning techniques to be used for the considered linear systems.

preprint2010arXiv

Multigrid methods for Toeplitz linear systems with different size reduction

Starting from the spectral analysis of g-circulant matrices, we consider a new multigrid method for circulant and Toeplitz matrices with given generating function. We assume that the size n of the coefficient matrix is divisible by g \geq 2 such that at the lower level the system is reduced to one of size n/g by employing g-circulant based projectors. We perform a rigorous two-grid convergence analysis in the circulant case and we extend experimentally the results to the Toeplitz setting, by employing structure preserving projectors. The optimality of the proposed two-grid method and of the multigrid method is proved, when the number theta \in N of recursive calls is such that 1 < theta < g. The previous analysis is used also to overcome some pathological cases, in which the generating function has zeros located at "mirror points" and the standard two-grid method with g = 2 is not optimal. The numerical experiments show the correctness and applicability of the proposed ideas both for circulant and Toeplitz matrices.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.