Researcher profile

S. Serra-Capizzano

S. Serra-Capizzano contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

4 published item(s)

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.

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.