Source author record

S. Serra-Capizzano

S. Serra-Capizzano appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA math.FA math.SP Numerical Analysis

Catalog footprint

What is connected

5works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Fast Toeplitz eigenvalue computations, joining interpolation-extrapolation matrix-less algorithms and simple-loop theory

Under appropriate technical assumptions, the simple-loop theory allows to deduce various types of asymptotic expansions for the eigenvalues of Toeplitz matrices generated by a function $f$. Independently and under the milder hypothesis that $f$ is even and monotonic over $[0,π]$, matrix-less algorithms have been developed for the fast eigenvalue computation of large Toeplitz matrices, within a linear complexity in the matrix order: behind the high efficiency of such algorithms there are the expansions predicted by the simple-loop theory, combined with the extrapolation idea. Here we focus our attention on a change of variable, followed by the asymptotic expansion of the new variable, and we adapt the matrix-less algorithm to the considered new setting. Numerical experiments show a higher precision (till machine precision) and the same linear computation cost, when compared with the matrix-less procedures already presented in the relevant literature. Among the advantages, we concisely mention the following: a) when the coefficients of the simple-loop function are analytically known, the algorithm computes them perfectly; b) while the proposed algorithm is better or at worst comparable to the previous ones for computing the inner eigenvalues, it is extremely better for the computation of the extreme eigenvalues.

preprint2016arXiv

Staggered discontinuous Galerkin methods for the incompressible Navier-Stokes equations: spectral analysis and computational results

The goal of this paper is to create a fruitful bridge between the numerical methods for approximating partial differential equations (PDEs) in fluid dynamics and the (iterative) numerical methods for dealing with the resulting large linear systems. Among the main objectives are the design of new efficient iterative solvers and a rigorous analysis of their convergence speed. The link we have in mind is either the structure or the hidden structure that the involved coefficient matrices inherit, both from the continuous PDE and from the approximation scheme: in turn, the resulting structure is used for deducing spectral information, crucial for the conditioning and convergence analysis, and for the design of more efficient solvers. As specific problem we consider the incompressible Navier-Stokes equations, as numerical technique we consider a novel family of high order accurate Discontinuous Galerkin methods on staggered meshes, and as tools we use the theory of Toeplitz matrices generated by a function (in the most general block, multi-level form) and the more recent theory of Generalized Locally Toeplitz matrix-sequences. We arrive at a quite complete picture of the spectral features of the underlying matrices and this information is employed for giving a forecast of the convergence history of the conjugate gradient method, together with a discussion on new more advanced techniques (involving preconditioning, multigrid, multi-iterative solvers). Several numerical tests are provided and critically illustrated in order to show the validity and the potential of our analysis.

preprint2014arXiv

Essential spectral equivalence via multiple step preconditioning and applications to ill conditioned Toeplitz matrices

In this note, we study the fast solution of Toeplitz linear systems with coefficient matrix $T_n(f)$, where the generating function $f$ is nonnegative and has a unique zero at zero of any real positive order $θ$. As preconditioner we choose a matrix $τ_n(f)$ belonging to the so-called $τ$ algebra, which is diagonalized by the sine transform associated to the discrete Laplacian. In previous works, the spectral equivalence of the matrix sequences $\{τ_n(f)\}_n $ and $\{T_n(f) \}_n$ was proven under the assumption that the order of the zero is equal to $2$: in other words the preconditioned matrix sequence $\{τ^{-1}_n(f)T_n(f) \}_n $ has eigenvalues, which are uniformly away from zero and from infinity. Here we prove a generalization of the above result when $θ<2$. Furthermore, by making use of multiple step preconditioning, we show that the matrix sequences $\{τ_n(f)\}_n $ and $\{T_n(f) \}_n$ are essentially spectrally equivalent for every $θ>2$, i.e., for every $θ>2$, there exist $m_θ$ and a positive interval $[α_θ,β_θ]$ such that all the eigenvalues of $\{τ^{-1}_n(f)T_n(f) \}_n $ belong to this interval, except at most $m_θ$ outliers larger than $β_θ$. Such a nice property, already known only when $θ$ is an even positive integer greater than 2, is coupled with the fact that the preconditioned sequence has an eigenvalue cluster at one, so that the convergence rate of the associated preconditioned conjugate gradient method is optimal. As a conclusion we discuss possible generalizations and we present selected numerical experiments.

preprint2012arXiv

Perturbation of operators and approximation of spectrum

Let A(x) be a holomorphic family of bounded self-adjoint operators on a separable Hilbert space H and let A(x)_n be the orthogonal compressions of A(x) to the span of first n elements of an orthonormal basis of H. The problem considered here is to approximate the spectrum of A(x) using the sequence of eigenvalues of A(x)_n. We show that the bounds of the essential spectrum and the discrete spectral values outside the bounds of essential spectrum of A(x) can be approximated uniformly on all compact subsets by the sequence of eigenvalue functions of A(x)_n. The known results for a bounded selfadjoint operator, are translated into the case of a holomorphic family of operators. Also an attempt is made to predict the existence of spectral gaps that may occur between the bounds of essential spectrum of A(0) = A and study the effect of holomorphic perturbation of operators in the prediction of spectral gaps. As an example, gap issues of some block Toeplitz-Laurent operators are discussed. The pure linear algebraic approach is the main advantage of the results here.

preprint2008arXiv

A Note on Algebraic Multigrid Methods for the Discrete Weighted Laplacian

In recent contributions, algebraic multigrid methods have been designed and studied from the viewpoint of the spectral complementarity. In this note we focus our efforts on specific applications and, more precisely, on large linear systems arising from the approximation of weighted Laplacian with various boundary conditions. We adapt the multigrid idea to this specific setting and we present and critically discuss a wide numerical experimentation showing the potentiality of the considered approach.