BZPEER
More

Discover

About BZPEERSearchGraph

Research tools

HomeCollectionsResearchCollaborationPublish

Network

ExpertsCommunities

Opportunities

Listings

Account

Log inCreate account
Log inCreate account

Source author record

Marco Caliari

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

ResearcherUnclaimed source record
math.NAmath.APNumerical Analysis

Catalog footprint

Research overview

7works
3topics
4close collaborators

Actions

Is this your work?

Log in to claim
Open graphBrowse papers

Published work

7 publications

preprint2022arXiv

A μμ-mode BLAS approach for multidimensional tensor-structured problems

In this manuscript, we present a common tensor framework which can be used to generalize one-dimensional numerical tasks to arbitrary dimension dd by means of tensor product formulas. This is useful, for example, in the context of multivariate interpolation, multidimensional function approximation using pseudospectral expansions and solution of stiff differential equations on tensor product domains. The key point to obtain an efficient-to-implement BLAS formulation consists in the suitable usage of the μμ-mode product (also known as tensor-matrix product or mode- nn product) and related operations, such as the Tucker operator. Their MathWorks MATLAB/GNU Octave implementations are discussed in the paper, and collected in the package KronPACK. We present numerical results on experiments up to dimension six from different fields of numerical analysis, which show the effectiveness of the approach.

Open access
preprint2021arXiv

A μμ-mode integrator for solving evolution equations in Kronecker form

In this paper, we propose a μμ-mode integrator for computing the solution of stiff evolution equations. The integrator is based on a dd-dimensional splitting approach and uses exact (usually precomputed) one-dimensional matrix exponentials. We show that the action of the exponentials, i.e. the corresponding batched matrix-vector products, can be implemented efficiently on modern computer systems. We further explain how μμ-mode products can be used to compute spectral transforms efficiently even if no fast transform is available. We illustrate the performance of the new integrator by solving, among the others, three-dimensional linear and nonlinear Schrödinger equations, and we show that the μμ-mode integrator can significantly outperform numerical methods well established in the field. We also discuss how to efficiently implement this integrator on both multi-core CPUs and GPUs. Finally, the numerical experiments show that using GPUs results in performance improvements between a factor of 1010 and 2020, depending on the problem.

Open access
preprint2016arXiv

A splitting approach for the magnetic Schrödinger equation

The Schrödinger equation in the presence of an external electromagnetic field is an important problem in computational quantum mechanics. It also provides a nice example of a differential equation whose flow can be split with benefit into three parts. After presenting a splitting approach for three operators with two of them being unbounded, we exemplarily prove first-order convergence of Lie splitting in this framework. The result is then applied to the magnetic Schrödinger equation, which is split into its potential, kinetic and advective parts. The latter requires special treatment in order not to lose the conservation properties of the scheme. We discuss several options. Numerical examples in one, two and three space dimensions show that the method of characteristics coupled with a nonequispaced fast Fourier transform (NFFT) provides a fast and reliable technique for achieving mass conservation at the discrete level.

Open access
preprint2016arXiv

Reliability of the time splitting Fourier method for singular solutions in quantum fluids

We extensively study the numerical accuracy of the well-known time splitting Fourier spectral method for the approximation of singular solutions of the Gross-Pitaevskii equation. In particular, we explore its capability of preserving a steady-state vortex solution, whose density profile is approximated by a very accurate diagonal Padé expansion of order 8, here explicitly derived for the first time. Although the Fourier spectral method turns out to be only slightly more accurate than a time splitting finite difference scheme, the former is reliable and efficient. Moreover, at a post-processing stage, it allows an accurate evaluation of the solution outside grid points, thus becoming particularly appealing when high resolution is needed, such as in the study of quantum vortex interactions.

Open access
preprint2016arXiv

The Leja method revisited: backward error analysis for the matrix exponential

The Leja method is a polynomial interpolation procedure that can be used to compute matrix functions. In particular, computing the action of the matrix exponential on a given vector is a typical application. This quantity is required, e.g., in exponential integrators. The Leja method essentially depends on three parameters: the scaling parameter, the location of the interpolation points, and the degree of interpolation. We present here a backward error analysis that allows us to determine these three parameters as a function of the prescribed accuracy. Additional aspects that are required for an efficient and reliable implementation are discussed. Numerical examples that illustrate the performance of our Matlab code are included.

Open access

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.