Researcher profile

Marco Donatelli

Marco Donatelli contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
3topics
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

6 published item(s)

preprint2022arXiv

Symbol based convergence analysis in multigrid methods for saddle point problems

Saddle point problems arise in a variety of applications, e.g., when solving the Stokes equations. They can be formulated such that the system matrix is symmetric, but indefinite, so the variational convergence theory that is usually used to prove multigrid convergence cannot be applied. In a 2016 paper in Numerische Mathematik Notay has presented a different algebraic approach that analyzes properly preconditioned saddle point problems, proving convergence of the Two-Grid method. In the present paper we analyze saddle point problems where the blocks are circulant within this framework. We are able to derive sufficient conditions for convergence and provide optimal parameters for the preconditioning of the saddle point problem and for the point smoother that is used. The analysis is based on the generating symbols of the circulant blocks. Further, we show that the structure can be kept on the coarse level, allowing for a recursive application of the approach in a W- or V-cycle and proving the "level independency" property. Numerical results demonstrate the efficiency of the proposed method in the circulant and the Toeplitz case.

preprint2021arXiv

Compatibility, embedding and regularization of non-local random walks on graphs

Several variants of the graph Laplacian have been introduced to model non-local diffusion processes, which allow a random walker to {\textquotedblleft jump\textquotedblright} to non-neighborhood nodes, most notably the transformed path graph Laplacians and the fractional graph Laplacian. From a rigorous point of view, this new dynamics is made possible by having replaced the original graph $G$ with a weighted complete graph $G'$ on the same node-set, that depends on $G$ and wherein the presence of new edges allows a direct passage between nodes that were not neighbors in $G$. We show that, in general, the graph $G'$ is not compatible with the dynamics characterizing the original model graph $G$: the random walks on $G'$ subjected to move on the edges of $G$ are not stochastically equivalent, in the wide sense, to the random walks on $G$. From a purely analytical point of view, the incompatibility of $G'$ with $G$ means that the normalized graph $\hat{G}$ can not be embedded into the normalized graph $\hat{G}'$. Eventually, we provide a regularization method to guarantee such compatibility and preserving at the same time all the nice properties granted by $G'$.

preprint2021arXiv

Graph Laplacian for image deblurring

Image deblurring is relevant in many fields of science and engineering. To solve this problem, many different approaches have been proposed and among the various methods, variational ones are extremely popular. These approaches are characterized by substituting the original problem with a minimization one where the functional is composed of two terms, a data fidelity term and a regularization term. In this paper we propose, in the classical $\ell^2-\ell^1$ minimization with the non-negativity constraint of the solution, the use of the graph Laplacian as regularization operator. Firstly, we describe how to construct the graph Laplacian from the observed noisy and blurred image. Once the graph Laplacian has been built, we solve efficiently the proposed minimization problem splitting the convolution operator and the graph Laplacian by the alternating direction method of multipliers (ADMM). Some selected numerical examples show the good performances of the proposed algorithm.

preprint2020arXiv

An inexact non stationary Tikhonov procedure for large-scale nonlinear ill-posed problems

In this work we consider the stable numerical solution of large-scale ill-posed nonlinear least squares problems with nonzero residual. We propose a non-stationary Tikhonov method with inexact step computation, specially designed for large-scale problems. At each iteration the method requires the solution of an elliptical trust-region subproblem to compute the step. This task is carried out employing a Lanczos approach, by which an approximated solution is computed. The trust region radius is chosen to ensure the resulting Tikhonov regularization parameter to satisfy a prescribed condition on the model, which is proved to ensure regularizing properties to the method. The proposed approach is tested on a parameter identification problem and on an image registration problem, and it is shown to provide important computational savings with respect to its exact counterpart.

preprint2015arXiv

Regularization preconditioners for frame-based image deblurring with reduced boundary artifacts

Thresholding iterative methods are recently successfully applied to image deblurring problems. In this paper, we investigate the modified linearized Bregman algorithm (MLBA) used in image deblurring problems, with a proper treatment of the boundary artifacts. We consider two standard approaches: the imposition of boundary conditions and the use of the rectangular blurring matrix. The fast convergence of the MLBA depends on a regularizing preconditioner that could be computationally expensive and hence it is usually chosen as a block circulant circulant block (BCCB) matrix, diagonalized by discrete Fourier transform. We show that the standard approach based on the BCCB preconditioner may provide low quality restored images and we propose different preconditioning strategies, that improve the quality of the restoration and save some computational cost at the same time. Motivated by a recent nonstationary preconditioned iteration, we propose a new algorithm that combines such method with the MLBA.We prove that it is a regularizing and convergent method. A variant with a stationary preconditioner is also considered. Finally, a large number of numerical experiments shows that our methods provide accurate and fast restorations, when compared with the state of the art.

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.