Source author record

Mariarosa Mazza

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

ResearcherUnclaimed source record

math.NA Numerical Analysis

Catalog footprint

Research overview

4works

2topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

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.

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'$ .

preprint2014arXiv

On the constrained mock-Chebyshev least-squares

The algebraic polynomial interpolation on uniformly distributed nodes is affected by the Runge phenomenon, also when the function to be interpolated is analytic. Among all techniques that have been proposed to defeat this phenomenon, there is the mock-Chebyshev interpolation which is an interpolation made on a subset of the given nodes whose elements mimic as well as possible the Chebyshev-Lobatto points. In this work we use the simultaneous approximation theory to combine the previous technique with a polynomial regression in order to increase the accuracy of the approximation of a given analytic function. We give indications on how to select the degree of the simultaneous regression in order to obtain polynomial approximant good in the uniform norm and provide a sufficient condition to improve, in that norm, the accuracy of the mock-Chebyshev interpolation with a simultaneous regression. Numerical results are provided.

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.

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.