BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Paris Vassalos

Paris Vassalos contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysis

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author
3works
0followers
2topics
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

3 published item(s)

preprint2022arXiv

Fast Toeplitz eigenvalue computations, joining interpolation-extrapolation matrix-less algorithms and simple-loop conjectures: the preconditioned setting

Under appropriate technical assumptions, the simple-loop theory allows to deduce various types of asymptotic expansions for the eigenvalues of Toeplitz matrices $T_{n}(f)$ generated by a function $f$, unfortunately, such a theory is not available in the preconditioning setting, that is for matrices of the form $T_{n}^{-1}(g)T_{n}(l)$ with $l,g$ real-valued, $g$ nonnnegative and not identically zero almost everywhere. Independently and under the milder hypothesis that $f=\frac{l}{g}$ is even and monotonic over $[0,π]$, matrix-less algorithms have been developed for the fast eigenvalue computation of large preconditioned matrices of the type above, within a linear complexity in the matrix order: behind the high efficiency of such algorithms there are the expansions as in the case $g\equiv 1$, combined with the extrapolation idea, and hence we conjecture that the simple-loop theory has to be extended in such a new setting, as the numerics strongly suggest.Here we focus our attention on a change of variable, followed by the asymptotic expansion of the new variable, and we consider new matrix-less algorithms ad hoc for the current case. Numerical experiments show a much higher precision till machine precision and the same linear computation cost, when compared with the matrix-less procedures already proposed in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Preconditioners for Fractional Diffusion Equations Based on the Spectral Symbol

It is well known that the discretization of fractional diffusion equations (FDEs) with fractional derivatives $α\in(1,2)$, using the so-called weighted and shifted Grünwald formula, leads to linear systems whose coefficient matrices show a Toeplitz-like structure. More precisely, in the case of variable coefficients, the related matrix sequences belong to the so-called Generalized Locally Toeplitz (GLT) class. Conversely, when the given FDE have constant coefficients, using a suitable discretization, we encounter a Toeplitz structure associated to a nonnegative function $\mathcal{F}_α$, called the spectral symbol, having a unique zero at zero of real positive order between one and two. For the fast solution of such systems by preconditioned Krylov methods, several preconditioning techniques have been proposed in both the one and two dimensional cases. In this note we propose a new preconditioner denoted by $\mathcal{P}_{\mathcal{F}_α}$ which belongs to the $τ$ algebra and it is based on the spectral symbol $\mathcal{F}_α$. Comparing with some of the previously proposed preconditioners, we show that although the low band structure preserving preconditioners are more effective in the one-dimensional case, the new preconditioner performs better in the more challenging multi-dimensional setting.

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