Source author record

Franco Zivcovich

Franco Zivcovich appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis

Catalog footprint

What is connected

3works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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 $d$ 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- $n$ 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.

preprint2022arXiv

A second-order low-regularity correction of Lie splitting for the semilinear Klein--Gordon equation

The numerical approximation of the semilinear Klein--Gordon equation in the $d$-dimensional space, with $d=1,2,3$, is studied by analyzing the consistency errors in approximating the solution. By discovering and utilizing a new cancellation structure in the semilinear Klein--Gordon equation, a low-regularity correction of the Lie splitting method is constructed, which can have second-order convergence in the energy space under the regularity condition $(u,\partial_tu)\in L^\infty(0,T;H^{1+\frac{d}{4}}\times H^{\frac{d}{4}})$, where $d=1,2,3$ denotes the dimension of space. In one dimension, the proposed method is shown to have a convergence order arbitrarily close to $\frac53$ in the energy space for solutions in the same space, i.e. no additional regularity in the solution is required. Rigorous error estimates are presented for a fully discrete spectral method with the proposed low-regularity time-stepping scheme. Numerical examples are provided to support the theoretical analysis and to illustrate the performance of the proposed method in approximating both nonsmooth and smooth solutions of the semilinear Klein--Gordon equation.

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 $d$-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 $10$ and $20$, depending on the problem.