BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Stefano Serra-Capizzano

Stefano 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.NANumerical Analysismath.FA

Catalog footprint

What is connected

15works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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.

Published work

15 published item(s)

preprint2022arXiv

Fast algebraic multigrid for block-structured dense and Toeplitz-like-plus-Cross systems arising from nonlocal diffusion problems

Algebraic multigrid (AMG) is one of the most efficient iterative methods for solving large sparse system of equations. However, how to build/check restriction and prolongation operators in practical of AMG methods for nonsymmetric {\em sparse} systems is still an interesting open question [Brezina, Manteuffel, McCormick, Runge, and Sanders, SIAM J. Sci. Comput. (2010); Manteuffel and Southworth, SIAM J. Sci. Comput. (2019)]. This paper deals with the block-structured dense and Toeplitz-like-plus-Cross systems, including {\em nonsymmetric} indefinite, symmetric positive definite (SPD), arising from nonlocal diffusion problem and peridynamic problem. The simple (traditional) restriction operator and prolongation operator are employed in order to handle such block-structured dense and Toeplitz-like-plus-Cross systems, which is convenient and efficient when employing a fast AMG. We focus our efforts on providing the detailed proof of the convergence of the two-grid method for such SPD situations. The numerical experiments are performed in order to verify the convergence with a computational cost of only $\mathcal{O}(N \mbox{log} N)$ arithmetic operations, by using few fast Fourier transforms, where $N$ is the number of the grid points. To the best of our knowledge, this is the first contribution regarding Toeplitz-like-plus-Cross linear systems solved by means of a fast AMG.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rectangular GLT Sequences

The theory of generalized locally Toeplitz (GLT) sequences is a powerful apparatus for computing the asymptotic spectral distribution of square matrices $A_n$ arising from the discretization of differential problems. Indeed, as the mesh fineness parameter $n$ increases to $\infty$, the sequence $\{A_n\}_n$ often turns out to be a GLT sequence. In this paper, motivated by recent applications, we further enhance the GLT apparatus by developing a full theory of rectangular GLT sequences as an extension of the theory of classical square GLT sequences. We also detail an example of application as an illustration of the potential impact of the theory presented herein.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Theoretical results for eigenvalues, singular values, and eigenvectors of (flipped) Toeplitz matrices and related computational proposals

In a series of recent papers the spectral behavior of the matrix sequence $\{Y_nT_n(f)\}$ is studied in the sense of the spectral distribution, where $Y_n$ is the main antidiagonal (or flip matrix) and $T_n(f)$ is the Toeplitz matrix generated by the function $f$, with $f$ being Lebesgue integrable and with real Fourier coefficients. This kind of study is also motivated by computational purposes for the solution of the related large linear systems using the (preconditioned) MINRES algorithm. Here we complement the spectral study with more results holding both asymptotically and for a fixed dimension $n$, and with regard to eigenvalues, singular values, and eigenvectors of $T_n(f), Y_nT_n(f)$ and to several relationships among them: beside fast linear solvers, a further target is the design of ad hoc procedures for the computation of the related spectra via matrix-less algorithms, with a cost being linear in the number of computed eigenvalues. We emphasize that the challenge of the case of non-monotone generating functions is considered in the current work, for which the previous matrix-less algorithms fail. Numerical experiments are reported and commented, with the aim of showing in a visual way the theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Toeplitz Momentary Symbols: definition, results, and limitations in the spectral analysis of Structured Matrices

A powerful tool for analyzing and approximating the singular values and eigenvalues of structured matrices is the theory of GLT sequences. By the GLT theory one can derive a function, which describes the singular value or the eigenvalue distribution of the sequence, the latter under precise assumptions. However, for small values of the matrix size of the considered sequence, the approximations may not be as good as it is desirable, since in the construction of the GLT symbol one disregards small norm and low-rank perturbations. On the other hand, LFA can be used to construct polynomial symbols in a similar manner for discretizations, where the geometric information is present, but the small norm perturbations are retained. The main focus of this paper is the introduction of the concept of sequence of "Toeplitz momentary symbols", associated with a given sequence of truncated Toeplitz-like matrices. We construct the symbol in the same way as in the GLT theory, but we keep the information of the small norm contributions. The low-rank contributions are still disregarded, and we give an idea on the reason why this is negligible in certain cases and why it is not in other cases, being aware that in presence of high nonnormality the same low-rank perturbation can produce a dramatic change in the eigenvalue distribution. Moreover, a difference with respect to the LFA symbols is that GLT symbols and Toeplitz momentary symbols are more general and are applicable to a larger class of matrices. We show the applicability of the approach which leads to higher accuracy in some cases when compared with the GLT symbol. Finally, since for many applications and their analysis it is often necessary to consider non-square Toeplitz matrices, we formalize and provide some useful definitions, applicable for non-square Toeplitz momentary symbols.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asymptotic spectra of large (grid) graphs with a uniform local structure

We are concerned with sequences of graphs having a grid geometry, with a uniform local structure in a bounded domain $Ω\subset {\mathbb R}^d$, $d\ge 1$. We assume $Ω$ to be Lebesgue measurable with regular boundary and contained, for convenience, in the cube $[0,1]^d$. When $Ω=[0,1]$, such graphs include the standard Toeplitz graphs and, for $Ω=[0,1]^d$, the considered class includes $d$-level Toeplitz graphs. In the general case, the underlying sequence of adjacency matrices has a canonical eigenvalue distribution, in the Weyl sense, and we show that we can associate to it a symbol $f$. The knowledge of the symbol and of its basic analytical features provide many informations on the eigenvalue structure, of localization, spectral gap, clustering, and distribution type. Few generalizations are also considered in connection with the notion of generalized locally Toeplitz sequences and applications are discussed, stemming e.g. from the approximation of differential operators via numerical schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Iterated fractional Tikhonov regularization

Fractional Tikhonov regularization methods have been recently proposed to reduce the oversmoothing property of the Tikhonov regularization in standard form, in order to preserve the details of the approximated solution. Their regularization and convergence properties have been previously investigated showing that they are of optimal order. This paper provides saturation and converse results on their convergence rates. Using the same iterative refinement strategy of iterated Tikhonov regularization, new iterated fractional Tikhonov regularization methods are introduced. We show that these iterated methods are of optimal order and overcome the previous saturation results. Furthermore, nonstationary iterated fractional Tikhonov regularization methods are investigated, establishing their convergence rate under general conditions on the iteration parameters. Numerical results confirm the effectiveness of the proposed regularization iterations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimal order conditions for filter based regularization methods, with applications to variants of the Tikhonov method

We study filter based regularization methods for linear ill-posed problems between Hilbert spaces. We derive optimal order conditions under a-priori choice rules for the regularization parameter. Such analysis is applied to the fractional Tikhonov method and the weighted Tikhonov method recently investigated in [9] and [8], respectively. These variants of the classical Tikhonov method can be viewed as classes of methods depending on a further parameter, controlling the smoothness of the computed solution. From our previous analysis, we derive sufficient conditions on such parameter in order to have filter based regularization methods that agrees with the previous results in the literature. On the other hand, our analysis shows sufficient and necessary conditions to have an optimal order fractional Tikhonov method, which are stronger than those in [9]. The same analysis is performed also for the weighted Tikhonov method.

0 citations0 reviewsNo code linkedOpen access
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 $I_k:=(-π,π)^k$, let $\mathcal M_s$ be the linear space of complex $s\times s$ matrices, and let $f,g:I_k\to\mathcal M_s$ be functions whose components $f_{ij},\,g_{ij}:I_k\to\mathbb C,\ i,j=1,\ldots,s,$ belong to $L^\infty$. Consider the matrices $T_n^{-1}(g)T_n(f)$, where $n:=(n_1,\ldots,n_k)$ varies in $\mathbb N^k$ and $T_n(f),T_n(g)$ are the multilevel block Toeplitz matrices of size $n_1\cdots n_ks$ generated by $f,g$. Then $\{T_n^{-1}(g)T_n(f)\}_{n\in\mathbb N^k}\sim_λg^{-1}f$, i.e. the family of matrices $\{T_n^{-1}(g)T_n(f)\}_{n\in\mathbb N^k}$ has a global (asymptotic) spectral distribution described by the function $g^{-1}f$, provided $g$ possesses certain properties (which ensure in particular the invertibility of $T_n^{-1}(g)$ for all $n$) and the following topological conditions are met: the essential range of $g^{-1}f$, defined as the union of the essential ranges of the eigenvalue functions $λ_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=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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Fast Alternating Minimization Algorithm for Total Variation Deblurring Without Boundary Artifacts

Recently, a fast alternating minimization algorithm for total variation image deblurring (FTVd) has been presented by Wang, Yang, Yin, and Zhang [{\em SIAM J. Imaging Sci.}, 1 (2008), pp. 248--272]. The method in a nutshell consists of a discrete Fourier transform-based alternating minimization algorithm with periodic boundary conditions and in which two fast Fourier transforms (FFTs) are required per iteration. In this paper, we propose an alternating minimization algorithm for the continuous version of the total variation image deblurring problem. We establish convergence of the proposed continuous alternating minimization algorithm. The continuous setting is very useful to have a unifying representation of the algorithm, independently of the discrete approximation of the deconvolution problem, in particular concerning the strategies for dealing with boundary artifacts. Indeed, an accurate restoration of blurred and noisy images requires a proper treatment of the boundary. A discrete version of our continuous alternating minimization algorithm is obtained following two different strategies: the imposition of appropriate boundary conditions and the enlargement of the domain. The first one is computationally useful in the case of a symmetric blur, while the second one can be efficiently applied for a nonsymmetric blur. Numerical tests show that our algorithm generates higher quality images in comparable running times with respect to the Fast Total Variation deconvolution algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Proposal of Multigrid Methods for Hermitian Positive Definite Linear Systems enjoying an order relation

Given a multigrid procedure for linear systems with coefficient matrices $A_n$, we discuss the optimality of a related multigrid procedure with the same smoother and the same projector, when applied to properly related algebraic problems with coefficient matrices $B_n$: we assume that both $A_n$ and $B_n$ are positive definite with $A_n\le \vartheta B_n$, for some positive $\vartheta$ independent of $n$. In this context we prove the Two-Grid method optimality. We apply this elementary strategy for designing a multigrid solution for modifications of multilevel structured (Toeplitz, circulants, Hartley, sine ($τ$ class) and cosine algebras) linear systems, in which the coefficient matrix is banded in a multilevel sense and Hermitian positive definite. In such a way, several linear systems arising from the approximation of integro-differential equations with various boundary conditions can be efficiently solved in linear time (with respect to the size of the algebraic problem). Some numerical experiments are presented and discussed, both with respect to Two-Grid and multigrid procedures.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Korovkin results and Frobenius optimal approximants for infinite dimensional bounded linear operators

The classical as well as non commutative Korovkin-type theorems deal with convergence of positive linear maps with respect to modes of convergences such as norm convergence and weak operator convergence. In this article, Korovkin-type theorems are proved for convergence of completely positive maps with respect to weak, strong and uniform clustering of sequences of matrices of growing order. Such modes of convergence were originally considered for Toeplitz matrices (see [23],[26]). As an application, we translate the Korovkin-type approach used in the finite dimensional case, in the setting of preconditioning large linear systems with Toeplitz structure, into the infinite dimensional context of operators acting on separable Hilbert spaces. The asymptotic of these pre-conditioners are obtained and analyzed using the concept of completely positive maps. It is observed that any two limit points of the same sequence of pre-conditioners are the same modulo compact operators. Finally, we prove the generalized versions of the Korovkin type theorems in [23]. Keywords: Completely positive maps, Frobenius norm, Pre-conditioners.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal preconditioning for image deblurring with Anti-Reflective boundary conditions

Inspired by the theoretical results on optimal preconditioning stated by Ng, R.Chan, and Tang in the framework of Reflective boundary conditions (BCs), in this paper we present analogous results for Anti-Reflective BCs, where an additional technical difficulty is represented by the non orthogonal character of the Anti-Reflective transform and indeed the technique of Ng, R.Chan, and Tang can not be used. Nevertheless, in both cases, the optimal preconditioner is the blurring matrix associated to the symmetrized Point Spread Function (PSF). The geometrical idea on which our proof is based is very simple and general, so it may be useful in the future to prove theoretical results for new proposed boundary conditions. Computational results show that the preconditioning strategy is effective and it is able to give rise to a meaningful acceleration both for slightly and highly non-symmetric PSFs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fast Preconditioners for Total Variation Deblurring with Anti-Reflective Boundary Conditions

In recent works several authors have proposed the use of precise boundary conditions (BCs) for blurring models and they proved that the resulting choice (Neumann or reflective, anti-reflective) leads to fast algorithms both for deblurring and for detecting the regularization parameters in presence of noise. When considering a symmetric point spread function, the crucial fact is that such BCs are related to fast trigonometric transforms. In this paper we combine the use of precise BCs with the Total Variation (TV) approach in order to preserve the jumps of the given signal (edges of the given image) as much as possible. We consider a classic fixed point method with a preconditioned Krylov method (usually the conjugate gradient method) for the inner iteration. Based on fast trigonometric transforms, we propose some preconditioning strategies which are suitable for reflective and anti-reflective BCs. A theoretical analysis motivates the choice of our preconditioners and an extensive numerical experimentation is reported and critically discussed. The latter shows that the TV regularization with anti-reflective BCs implies not only a reduced analytical error, but also a lower computational cost of the whole restoration procedure over the other BCs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multigrid and preconditioning strategies for implicit PDE solvers for degenerate parabolic equations

The novel contribution of this paper relies in the proposal of a fully implicit numerical method designed for nonlinear degenerate parabolic equations, in its convergence/stability analysis, and in the study of the related computational cost. In fact, due to the nonlinear nature of the underlying mathematical model, the use of a fixed point scheme is required and every step implies the solution of large, locally structured, linear systems. A special effort is devoted to the spectral analysis of the relevant matrices and to the design of appropriate iterative or multi-iterative solvers, with special attention to preconditioned Krylov methods and to multigrid procedures: in particular we investigate the mutual benefit of combining in various ways suitable preconditioners with V-cycle algorithms. Numerical experiments in one and two spatial dimensions for the validation of our multi-facet analysis complement this contribution.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access